模糊自适应PID控制器

模糊自适应PID控制器的设计模糊自适应PID控制器的设计一、模糊自适应原理模糊控制是以模糊集合论、模糊语言变量和模糊逻辑推理为基础的一种计算机控制方法，作为智能控制的一个重要分支，在控制领域获得了广泛应用，模糊控制与传统控制方式相比具有以下突出优点:·不需要精确的被控对象的数学模型；·使用自然语言方法，控制方法易于掌握；·鲁棒性好，能够较大范围的适应参数变化；·与常规PID控制相比，动态响应品质优良。常规模糊控制器的原理如图1所示:图1模糊控制系统框图PID控制规律:101()[()()()]pDIdutketetdtTetTdt式中：pk---比例系数；IT---积分时间常数；DT---微分时间常数。在工业生产中过程中，许多被控对象随着负荷变化或干扰因素影响，其对象特性参数或结构发生改变。自适应控制运用现代控制理论在线辨识对象特征参数，实时改变其控制策略，使控制系统品质指标保持在最佳范围内，但其控制效果的好坏取决于辨识模型的精确度，这对于复杂系统是非常困难的。因此，在工业生产中过程中，大量采用的仍然是PID算法，PID参数的整定方法很多，但大多数都以对象特性为基础。随着计算机技术的以展，人们利用人工智能的方法将操作人员的调整经验作为知识存入计算机中，根据现场实际情况，计算机能自动调整PID参数，这样就出现了智能PID.这种控制器把古典的PID控制与先进的专家系统相结合，实现系统的最佳控制。这种控制必须精确地确定模型，首先将操作人员长期实践积累Fuzzy化Fuzzy控制算法Fuzzy判决被控过程d/dt+r-eecuyec的经验知识用控制然后运用推理便可对PID参数实现最佳调整。由于操作者经验不易精确描述，控制过程中各种信号量及评价指标不易定量表示，模糊理念是解决这一问题的有效途径，所以人们运用学的基本理论和方法，把规则的条件、操作用模糊集表示，并把这些模糊有关信息作为知识存入计算机知识库中，然后计算机根据控制系统的实际响应情况，运用模糊推理，即可自动实现对PID参数的最佳调整，这就是模糊自适应PID控制，目前模糊自适应PID控制器有多种结构，但其工作原理基本一致。自适应模糊PID以误差e和误差变化率ec作为输入，可以满足不同时刻的e和ec对PID参数自整定的要求，利用模糊控制规则在线对PID参数进行修改，便构成了自适应模糊PID控制器，其结构如图2所示.图2自适应模糊控制器结构PID参数模糊自整定是找出PID三个参数与e和ec之间的模糊关系，在运行中通过不断检测e和ec，根据模糊控制原理对三个参数进行在线修改，以满足不同e和ec对控制参数的不同要求，而使被控对象有良好的动、静态性能。从系统的稳定性、响应速度、超调量和精度等各方面来考虑，kp.ki,kd的作用如下：（1）比例系数kp的作用是加快系统的响应速度，提高系统的调节精度。kp越大，系统的响应速度越快，系统的调节精度越高，但易产生超调，甚至会导致系统不稳定。kp取值过小，则会降低调节精度，使响应速度缓慢，从而延长调节时间，使系统的动态静态特性变坏。（2）积分作用系数ki的作用是消除系统的稳态误差。ki越大，的静态误差消除越快，但ki过大，在响应过程的初期会产生积分饱和现象，从而引起响应过程的较大超调。若ki过小，将使系统静态误差难以消除，影响系统的调节精度。（3）微分作用系数kd的作用是改善系统的动态特性，其作用主要是在响应过程中抑制偏差向任何方向的变化，对偏差变化进行提前预报。但kd过大，会使响应过程提前制动，从而延长调节时间，而且会降低系统的抗干扰性能二、自适应模糊PID控制的设计思路de/dtPID调节器模糊推理对象rinerrorkpkikdkdyoutPID参数的整定必须考虑到在不同时刻三个参数的作用及相互之间的关系。在线实时模糊自整定PID算法的基础上，通过计算当前系统e和误差变化率ec，利用模糊规则进行模糊推理，查询模糊矩阵表进行参数调整。核心是总结工程设计人员的技术知识和实际操作经验，建立合适的模糊规则表，得到针对kp,ki,kd三个参数分别整定的模糊控制表（表1，表2，和表3）。表1kp模糊控制规则表kpeceNBNMNSZOPSPMPBNBPBPBPMPMPSZOZONMPBPBPMPSPSZONSNSPMPMPMPSZONSNSZOPMPMPSZONSNMNMPSPSPSZONSNSNMNMPMPSZONSNMNMNMNBPBZOZONMNMNMNBNB表2ki模糊控制规则表kieceNBNMNSZOPSPMPBNBNBNBNMNMNSZOZONMNBNBNMNSNSZOZONSNBNMPSNSZOPSPSZONMNMNSZOPSPMPMPSNMNSZOPSPSPMPBPMZOZOPSPSPMPBPBPBZOZOPSPMPMPBPB表3kd模糊控制规则表kdeceNBNMNSZOPSPMPBNBPSNSNBNBNBNMPSNMPSNSNBNMNMNSZONSZONSNMNMNSNSZOZOZONSNSNSNSNSZOPSZOZOZOZOZOZOZOPMPMNSPSPSPSPSPBPBPBPMPMPMPSPSPBkp,ki,kd的控制规则表建立好后，可根据如下方法kp,ki,kd的自适应校正。将系统e和误差变化率ec变化范围定义为模糊集上的论域。e,ec={-5,-4,-3,-2,-1,0,1,2,3,4,5}其模糊子集为e,ec={NB,NM,NS,ZO,PS,PM,PB}子集中元素分别代表负大，负中，负小，零，正小，正中，正大。根据各模糊子集的隶属度赋值表和各参数模糊控制模型，应用模糊合成推理设计PID参数的模糊矩阵表，查出修正参数，代入下式计算：()(1){,}iipkpkkpkeec()(1){,}iiikikkikeec()(1){,}iidkdkkdkeec在线运行过程中，控制系统通过模糊逻辑规则的结果处理、查表和运算，完成对PID参数的在线自校正。其工作流程如图3所示：图3在线自校正工作流程图三、仿真结果被控对象为：32523500()87.3510470Gpssss采样时间为1ms，采用模糊PID控制进行跃响应，在第300个采样时间控制器输出加1.0的干扰，相应的响应结果如图所示：取当前采样值入口e(k)=r(k)-y(k)ec(k)=e(k)-e(k-1)返回PID控制器输出计算当前kp,ki,kd模糊整定△kp,△ki,△kde(k),ec(k)的模糊化e(k-1)=e(k)00.050.10.150.20.250.30.350.40.450.500.20.40.60.811.21.4time(s)rin,yout模糊PID控制阶跃响应00.050.10.150.20.250.30.350.40.450.5-0.200.20.40.60.811.2time(s)error模糊PID控制误差响应a)模糊PID控制阶跃响应b)模糊PID控制误差响应00.050.10.150.20.250.30.350.40.450.5-0.200.20.40.60.811.21.4time(s)u控制器输出u00.050.10.150.20.250.30.350.40.450.50.330.340.350.360.370.380.390.40.410.420.43time(s)kpkp的自适应调整c)控制器输出ud)kp的自适应调整00.050.10.150.20.250.30.350.40.450.5-20246810121416x10-3time(s)kiki的自适应调整00.050.10.150.20.250.30.350.40.450.50.20.30.40.50.60.70.80.91time(s)kdkd的自适应调整e)ki的自适应控制调整e)kd的自适应控制调整图4仿真结果其中，e,de,kp,ki,kd的隶属度函数：-3-2-1012300.20.40.60.81eDegreeofmembershipNBNMNSZPSPMPBe的隶属度函数-3-2-1012300.20.40.60.81ecDegreeofmembershipNBNMNSZPSPMPBec的隶属度函数a)误差的隶属度函数b)误差变化率的隶属度函数-0.25-0.2-0.15-0.1-0.0500.050.10.150.20.2500.20.40.60.81kpDegreeofmembershipNBNMNSZPSPMPBkp的隶属度函数-0.06-0.04-0.0200.020.040.0600.20.40.60.81kiDegreeofmembershipNBNMNSZPSPMPBki的隶属度函数c)kp的隶属度函数d)ki的隶属度函数-3-2-1012300.20.40.60.81kdDegreeofmembershipNBNMNSZPSPMPBkd的隶属度函数e)kd的隶属度函数图5各变量的隶属度函数Systemfuzzpid:2inputs,3outputs,49rulese(7)ec(7)kp(7)ki(7)kd(7)fuzzpid(mamdani)49rules图6模糊PID控制系统构成得到的模糊规则有以下49条：MATLAB下运行showrule(a)，可得到以下49条模糊规则：1.If(eisNB)and(ecisNB)then(kpisPB)(kiisNB)(kdisPS)(1)2.If(eisNB)and(ecisNM)then(kpisPB)(kiisNB)(kdisNS)(1)3.If(eisNB)and(ecisNS)then(kpisPM)(kiisNM)(kdisNB)(1)4.If(eisNB)and(ecisZ)then(kpisPM)(kiisNM)(kdisNB)(1)5.If(eisNB)and(ecisPS)then(kpisPS)(kiisNS)(kdisNB)(1)6.If(eisNB)and(ecisPM)then(kpisZ)(kiisZ)(kdisNM)(1)7.If(eisNB)and(ecisPB)then(kpisZ)(kiisZ)(kdisPS)(1)8.If(eisNM)and(ecisNB)then(kpisPB)(kiisNB)(kdisPS)(1)9.If(eisNM)and(ecisNM)then(kpisPB)(kiisNB)(kdisNS)(1)10.If(eisNM)and(ecisNS)then(kpisPM)(kiisNM)(kdisNB)(1)11.If(eisNM)and(ecisZ)then(kpisPS)(kiisNS)(kdisNM)(1)12.If(eisNM)and(ecisPS)then(kpisPS)(kiisNS)(kdisNM)(1)13.If(eisNM)and(ecisPM)then(kpisZ)(kiisZ)(kdisNS)(1)14.If(eisNM)and(ecisPB)then(kpisNS)(kiisZ)(kdisZ)(1)15.If(eisNS)and(ecisNB)then(kpisPM)(kiisNB)(kdisZ)(1)16.If(eisNS)and(ecisNM)then(kpisPM)(kiisNM)(kdisNS)(1)17.If(eisNS)and(ecisNS)then(kpisPM)(kiisNS)(kdisNM)(1)18.If(eisNS)and(ecisZ)then(kpisPS)(kiisNS)(kdisNM)(1)19.If(eisNS)and(ecisPS)then(kpisZ)(kiisZ)(kdisNS)(1)20.If(eisNS)and(ecisPM)then(kpisNS)(kiisPS)(kdisNS)(1)21.If(eisNS)and(ecisPB)then(kpisNS)(kiisPS)(kdisZ)(1)22.If(eisZ)and(ecisNB)the