（参考）模糊.docVIP

编程思路：根据PID控制器的三个参数与偏差e和偏差的变化ec之间的模糊关系，在运行时不断检测e及ec，通过事先确定的关系，利用模糊推理的方法，在线修改PID控制器的三个参数，让PID参数可自整定。 模糊控制器的设计： 本例模糊控制器的设计直接采用matlab中自带的FIS工具箱，控制器的输入为e和(e ,输出为PID控制器的控制三个参数Kp、Ki和Kd。隶属度函数分别采用Z型，三角形和S型隶属度函数。清晰化方法采用重心法。因为e和ec各有7个模糊分割，故模糊规则最多有49条规则,模糊规则采用AND方法。 图一：模糊自适应整定PID控制原理图 PID算法： PID参数的作用 kp：增加kp可加快响应速度，减少稳态误差，但易产生超调，影响稳定性； ki：增加ki有利于减少稳态误差，但易使系统响应滞后，动态品质变差； kd：增加kd可改善系统动态品质，但会对躁声敏感，影响系统稳定性。 被控对象： 输入信号：采用采样周期为0.1s，幅值为1的单位阶跃信号。 模糊控制设计的核心是总结工程设计人员的技术知识和实际操作经验，建立合适的模糊规则表，下表是Kp、Ki和Kd这3个参数分别整定的模糊控制表。 KP ec NB NM NS ZE PS PM PB e NB PB PB PM PM PS ZO ZO NM PB PB PM PS PS ZO NS NS PM PM PM PS ZO NS NS ZE PM PM PS ZO NS NM NM PS PS PS ZO NS NS NM NM PM PS ZO NS NS NM NM NB PB ZO ZO NM NM NM NB NB 表1 Kp模糊控制规则表 Ki ec NB NM NS ZE PS PM PB e NB NB NB NM NM NS ZO ZO NM NB NB NM NS NS ZO ZO NS NB NM NS NS ZO PS PS ZE NM NM NS ZO PS PM PM PS NM NS ZO PS PS PM PB PM ZO ZO PS PS PM PB PB PB ZO ZO PS PM PM PB PB 表2 Ki模糊控制规则表 Kd ec NB NM NS ZE PS PM PB e NB PS NS NB NB NB NM PS NM PS NS NB NM NM NS ZO NS ZO NS NM NM NS NS ZO ZE ZO NS NS NS NS NS ZO PS PB ZO ZO ZO ZO ZO ZO PM PB PS PS PS PS PS PB PB PB PM PM PS PS PS PB 表3 Kd模糊控制规则表 Kp、Ki、Kd的模糊规则表建立好后，可根据如下方法进行Kp、Ki、Kd的自适应校正。将系统误差e和误差变化率ec变化范围定义为模糊集上的论域，既 e，ec={-3，-2，-1，0，1，2，3 } 其模糊子集为： e，ec={NB，NM，NS，ZO ,PS ,PM,PB} 子集中元素分别代表负大、负中、负小、零、正小、正中、正大。设e、ec和Kp、Ki、Kd均服从正态分布因此可得出各模糊子集的隶属度，根据各模糊子集的隶属度赋值表和各参数模糊控制模型，应用模糊合成推理设计PID参数的模糊规则表，查出修正参数带入下式：Kp=Kp’+△Kp Ki=Ki’+△Ki Kd=Kd’+△Kd 在线运行过程中,控制系统通过对模糊规则的结果处理\查表和运算,完成对PID参数的在线自校正。 仿真结果如下： Matlab程序如下： clear a=newfis(fuzzypid);%创建新的FIS %*******设定e的范围、隶属度******** a=addvar(a,input,e,[-3 3]);%向模糊推理系统（FIS）添加语言变量e %********向模糊推理系统（FIS）的语言变量e添加隶属度函数**** a=addmf(a,input,1,NB,zmf,[-3 -1]); %z形隶属度函数 a=addmf(a,input,1,NM,trimf,[-3 -2 0]);%三角形隶属度函数 a=addmf(a,input,1,Ns,trimf,[-3 -1 1]);%三角形隶属度函数 a=addmf(a,input,1,Z,trimf,[-2 0 2]);%三角形隶属度函数 a=addmf(a,input,1,PS,trimf,[-1 1 3]);%三角形隶属度函数 a=addmf(a,input,1,PM,trimf,[0 2 3]);%三角形隶属度函数 a=addmf(a,input,1,PB,smf,[1 3]