模糊PID控制器设计

1、模糊PID控制器o引言：PID控制作为一种典型的传统反馈控制器，以其结构简单，易于实现和鲁棒性好等特点在工业过程控制中广泛应用。但是传统PID控制器的参数需要被控对象的数学模型来进行调整，而控制过程 中的滞后性、控制参数的非线性和高阶陛增加了对Kp、Ki、Kd三个参数的调整难度。所以对确定的控制系统通过复杂的计算后，其三个参数的值在控制运行中一般是固定的，不易进行在线的调整。而在实际的工业生产过程中，许多被控对象受到负荷变化和干扰因素的作用，其对象参数的特征和结构易发生改变，这就需要对参数进行动态的调整。同样因为被控系 统的复杂性和不确定性，其精确的数学模型难以建立，甚至无法建立模型，所以需要

2、利用模糊控制技术等方 法来解决。模糊PID无需考虑被控系统的模型，而只根据其误差e和误差变化ec等检测数据来自适应调整Kp、Ki、Kd的值，最终使被控系统处于稳定工作态。 1、传统PID控制器：PID参数模糊自整定是找出PID中3个参数与e和ec之间的模糊关系，在运行中通过不断检测e和ec,根据模糊控制原理来对3个参数进行在线修改，以满足不同e和ec时对控制参数的不同要求，而使被控对象有良好的动稳态性能。从系统的稳定性、响应速度、超调量和稳态精度 等方面来考虑Kp,Ki,Kd的作用如下：(1)比例系数Kp的作用是：加快系统的响应速度，提高系统的调节精度。Kp越大，系统的响应速度越快，系统的调节

3、精度越高，但易产生超调，甚至导致系统不稳定；Kp取值过小，则会降低调节精度，使响应速度缓慢，从而延长调节时间，使系统静态、动态特性变坏。 (2)积分作用系数Ki的作用是：消除系统的稳态误差。Ki越大，系统的稳态误差消除越快，但Ki过大，在响应过程的初期会产生积分饱和现象，从而引起响应过程的较大超调；若Ki过小，将使系统稳态误差难以消除，影响系统的调节精度。 (3)微分作用系数Kd的作用是：改善系统的动态特性。其作用主要是能反应偏差信号的变化趋势，并能在偏 差信号值变得太大之前，在系统中引入一个有效的早期修正信号，从而加快系统的动作速度，减少调节时 间。 2、模糊PID参数自整定控制器设计 2.

4、1控制器结构：由PID各个参数对系统的影响得到:(a)当误差|e|较大时，说明误差的绝对值较大，不论误差的变化趋势如何，都应该考虑控制器的Kp取较 大值，以提高响应的快速性；而为防止因为|ec|瞬时过大，kd应该取较小的值；为控制超调，ki也应该取 值很小。(b)当误差|e|在中等大小时，为保证系统的相应速度并控制超调，应减小Kp, Ki值应增大，Kd应适中。(c)当误差©较小时，为保证系统具有良好的稳态特性，应加大Kp、Ki的取值，同时为避免产生振荡，Kd的取值应该和16cl联系起来。模糊PID控制根据系统运行的不同状态，考虑Kp Ki、Kd三者的关联，根据工程经验设计模糊整定这三

5、个参数，选择输入语言变量为误差e和偏差变化率ec,语言变量值取NB,NM NSO, PS, PM, PB)七个模糊值；选择输出语言变量为©, Ki, Kd语言变量值也取NB, NM NSO, PS, PM, PB)七个模糊值，建立 Kp, Ki, Kd的模糊规则。2.3设计模糊控制器在MATLAB命令空间输入fuzzy,弹出对话框，然后进行设置，如图:S =S Editor zhinengtongzFiie Edit ViewzhmengkongzFIS Type：mamdaniecAnd methodmin<TCurrent VariableOr methodmaxNamek

6、mplicatiormin1TType庶ar虫七AggrcgcriionmaxTDefuzzificertio ncertroid<THedpinput-3 3CloseFIS Name：System '*zhinBngkcngz“； 2 inpifts, 3 outputs, and 49（Jew设置为两输入三输出结构。根据经验设置输入输出的隶属度函数都为:MembershipFunction Editor zhinengkongzRie Edit ViewFIS VariableskdReady设置完隶属度'函数后；按照上述三个模糊规则表格设置模糊规则如.图, HRu

7、le Editor： zhineriglcorigzFile Edit View Options1. W(e is rib) and - ec is nb 1 then (kp is pb)(ld is nlj)(kd is ps) (1)2. If (e is nbj and (ec is nm) then (kp is pbj(ki is nb)(kd is ns) (1)3. If (& is nbj and (ec is ns) then (kp is pmAlci is nm)(l£d is nbj (1)4. If (e is n 町 and (ec is zo)

8、 then (kp is pmAki isis nb (1)5. If (e is nb) and (ec is p$) then (kp is ps)(ki Is ns)(kd is nb) (1)6. If (e isnb)and (ec is pm)then (kp is zoXki is zo)(kd t£ nm)07. If (e is nb) and (ec is pb)thert (kp is zo)(ki Is zo)(kd is ps) (1)8. If e is nn) and (ec Is nb) then (kp Is Is nb)陋 is ias) (1)9

9、. If (e is nirr) and (ec Is nm) tnen (kp Isis nb 欣“is ns) (1)1D. If (e is nm) and (BC is ns) then (kp is pmXkl is nm)(kd is nb) 1 11. If (e is nm) and (BC is xoj then (kp is ps 丽 is ns)(kd is nm) 1)r Comedioin t w& And1!Delete ruleAdd rule Change ruleCloseFIS Name; zhinengAongzHelp”至磁盘，留给仿真调用，生成

10、的文件为:设置完毕后，保存文件zhi nen gk on gz.fis System Name='zh inengkon gz' Type='mamda ni' Versio n=2.0 Numl nputs=2 NumOutputs=3 NumRules=49 An dMethod='mi n' OrMethod='max' lmpMethod='mi n' AggMethod='max' DefuzzMethod=*ce ntroid'In put1 Name='e' R

11、an ge=-3 3 NumMFs=7 MF 仁nb':'zmf,-3-1 MF2=' nm'ftrimfJ-S -2 0 MF3=' ns'trimfJ-S -1 1 MF4=,zo,:'trimf,-2 0 2 MF5='ps,:'trimf,-1 1 3 MF6='pm':'trimf,0 2 3 MF7士pbTsmf,1 3Input2 Name='ec' Ran ge=-3 3 NumMFs=7 MF1=' nb'fzmfJ-S-l MF2=' nm

12、Ttrimf,卜3 -2 0 MF3=' ns'trimfJ-S -1 1 MF4='zo,:'trimf,-2 0 2 MF5='ps,:'trimf,-1 1 3 MF6=,pm':,trimf,0 2 3 MF7='pb,:'smf,1 3Outputl Name='kp' Ran ge=-3 3 NumMFs=7 MF1=' nb'fzmfJ-S-l MF2=' nmTtrimf,卜3 -2 0MF3rsTtrimf,卜3 -1 1 MF4=,zo,:'trimf,-

13、2 0 2 MF5='ps':'trimf,-1 1 3 MF6='pm':,trimf,0 2 3【LTcorJEWj-sullcoLLn o、corJEmEullCMLn tTeWENqullLLLwILl-WEnN ocoll 6 南 .pxll EeN Gndlno coL=ES.：.qdllZLW oCMo=Ef Edll9u_wcoL 工y.irsdllgLWCT0 e=Ew-:ozll寸 Ln HTcorJEWj-sullcoLn oQco=EE-;.-EUlleLW tTcoWENqullLLLnTSLWEnN coCOT 6 UBM&a

14、mp;lla EBN CM一 ndlno coL=ES.：-qdllZLW二(L)co寸T9 e 二eecogLDe 二(L)eco吗寸CM 二(L二 7 90e 二e8 L ZCMe 二(L)g 'Le 二 a) g 寸寸KL 二(L)OJ寸T9 L 二(L二 eglnL 二 a二 99 W L 二(L二 790L 二(L)COL Z 7 广 二a)g LZ二 ¥_n£ coL=ES?-qd.HZLW o。0=E_l-，Ed-H9LWco- L±ErsdllgLW n 0 e=E-s_:ozll寸 Ln2 7,344(1):13 1,6 1 4(1) :

15、 13 2,6 23 (1):13 3,632 (1):13 4,53 2(1):13 5,443(1):13 6,353 (1):13 7,354(1) : 14 1,624(1) : 14 2,623(1):14 3,533(1):14 4,443(1):14 5,353(1):14 6,26 3(1):14 7,264(1):15 1,524(1) : 15 2,534(1):15 3,444(1) : 15 4,354(1):15 5,354(1):15 6,264(1):15 7,274(1):16 1,547(1) : 16 2, 443 (1) : 16 3,355(1):16

16、4,255(1):16 5,265(1):16 6,275(1):16 7, 1 77(1) : 17 1,447(1) : 17 2,446(1):17 3,256(1):17 4,266(1):17 5,265(1):17 6, 1 75 (1) : 17 7, 1 77(1) : 13、系统仿真设控制对象为：G(s)=25/sa2+6s+25通过程序(含注释)仿真如下：clear all;cic;读取设定的FIS模糊控制器抽样时间为1MS构造系统 G(s)=25/sA2+6s+25a=readfis('zh inengkon gz'); %ts=O.OO1;%sys=tf

17、(25,1,6,25);dsys=c2d(sys,ts,'z'); %nu m,de n=tfdata(dsys,V); %液置系统初值离散化得到系数u1=0;u2=0;y1=0;y2=0;x=0,0,0;e1=0;ec 1=0;Kp0=0.3;Ki0=2;Kd0=1;fork=1:5000;%time(k)=k*ts;r(k)=1;%fpid=evalfis(e1 ,ec1 ,a); %Kp(k)=Kp0+fpid(1); %仿真5S输入阶跃信号 模糊推理 得到新的参数Ki(k)=KiO+fpid(2);Kd(k)=KdO+fpid(3);u(k)=Kp(k)*x(1)+Ki

18、(k)*x(2)+Kd(k)*x(3);%PID输出系统输出y(k)=-de n(2)*y1-de n(3)*y2+num(1 )*u(k)+num(2)*u1 + nu m(3)*u2;% xlabel(,time(s),);ylabel( 'Kd ); grid on;计算误差 更新数据对应KP 对应KI 对应KD重新计算ec 重新计算ee(k)=r(k)-y(k);%u2=u1;u1=u(k);%y2=yl;y 仁 y(k);x(1)=e(k);%x(2)=x(2)+e(k)*ts;%x(3)=e(k)-e1;%ec 仁 e(k)-e1;%e1=e(k);%end%丁印输出 figure(1);on;xlabel('time(s)');ylabel( fy');figure(2) plot(time,Kp,Y); xlabel('time(s)');ylabel( *Kp'); grid on;figure(3)plot(time,Ki,'b');grid on;xlabel('time(s),);ylabel( 'Ki');figure(4)plot(time,Kd；g');结果显示:输入输出对比图0.500511 522.5