模糊控制程序实例Word格式文档下载.docx

1、写成模糊推理句：if E=NB then U=PBif E=NS then U=PSif E=Z0 then U=Z0if E=PS then U=NSif E=PB then U=NB由上述的控制规则可得到模糊控制规则表，如表 5-4所示。表5-4模糊控制规则表EPSUZO模糊控制规则实际上是一组多重条件语句，它可以表示从误差论域 旦到控制量论域 U的模糊关系R。按着上述控制规则，可以得到该温度偏差与喷水阀门开度之间的模糊关系 R:R E U（NBe PBu）U（NSe PSu）U（Z0e Z0u）U（PSe NSu）U（PBe NBu ）计算模糊关系矩阵 R的子程序如F_Relation_

2、1.m 所示。%莫糊关系计算子程序F_Relation_1.cfunction R,mfe,mfu, ne,n u,Me=F_Relatio n_1%#!俞入模糊变量赋值表（表 5-3）#ne=7;%等级量e的个数nu=7;%等级量u的个数Me=0 0 0 0 0 0.5 1;0 0 0 0 1 0.5 0;0 0 0.5 1 0.5 0 0;0 0.5 1 0 0 0 0;1 0.5 0 0 0 0 0;Mu=Me;%#定义模糊变量及其语言值 仁PB,2=PS,3=O,4=NS,5=NB，并输入模糊控制规则表（表 5-4 ） #mfc=5;%模糊变量E的语言值个数，控制规则表列数E=5 4

3、3 2 1;UC=1 2 3 4 5;%#计算 R=EX U#R=zeros （ne,nu）;for i=1:mfciu=UC（i）;ie=E（i）;for k=1: nefor 1=1: nuif Me（ie,k）R（k,l）R（k,l）=Reu（k,l）;运行结果如下：1.0R 0由此，可以得到各种输入时，模糊控制器的输出:E 1，由表5-2可以查出单输入时： U EoR例如：单输入时，设e 1.2，根据前面的等级量划分原则，其等级量E PS 0 0 0 0 1 0.5 0。根据式（5-5）可以得到:U 0 0 0 0 1 0.5 0 00.5 0.5 1.0 0.5 0.5 0 0其模糊

4、决策子程序如 F_Deduce_1.m所示。%单输入时模糊决策子程序 F_Deduce_1.mfun ction FU=F_Deduce_1（fe,R, ne,nu）for j=1:if fe（j）EM;e1=EM;e ndif e1-EM;e1=-EM;j=e1+EM+1;Fi=1;for i=2:mu图5.11单输入模糊控制器的控制效果从上述的仿真结果可以看到，定值扰动时系统可以达到稳定，但是有很大的静差，不能满足 工程上的要求，究其原因是，模糊控制器实质上是一个具有继电器型非线性特性的控制器（如图 5.11中所示的u），没有积分作用，对于有自平衡对象一定会产生静差，而且系统极容易产生震荡

5、。 从图5.11就可看出，虽然设计的模糊控制器对定值扰动是稳定的，但对于内扰并不能使其稳定。非线性控制器的控制效果取决于各变量的论域及扰动量的大小 【，因此，模糊控制器的大范围工程应用还有许多问题需要研究。2、双输入模糊控制器的设计【例5.13】对于图5.10所示的系统，设计双输入模糊控制器，观察定值扰动和内部扰动的控制 效果。下面考虑两个输入的情况：设温度偏差e、偏差变化率ec及控制量u的实际论域：e ec u 1.5,1.5，选择它们的等级量论域分别为E 6, 5, 4, 3, 2,1, 0, 0, 1, 2, 3, 4, 5, 6EC6,5, 4,3,2, 1,0, 1, 2, 3, 4

6、, 5, 6U 7,2, 1,0, 1, 2, 3, 4, 5, 6, 7量化因子Ke ec2 64，Ku27 141.5（1.5） 3假设选取E的模糊变量词集为 NB,NM , NS, NO, PO, PS, PM , PB选取模糊变量E的赋值表如表5-5所示。表5-5 e的等级量与模糊量的关系-6-5-4-3-2-1-0+0+ 1+2+3+4+5+6PB0.10.40.8PM0.20.70.3PO0.6NONMNB选取EC、U的模糊变量词集为NB,NM , NS,ZO, PS, PM ,PB选取模糊变量EC的赋值表如表5-2所示，模糊变量U的赋值表如表5-6所示。 表5-6等级量U与模糊量

7、U的关系-7+7nm现将操作者在操作过程中遇到的各种出现的情况和相应的控制策略汇总为表 5-7。表5-7双输入时的模糊控制规则表nsnopopm下面说明建立模糊控制规则表的基本思想。首先考虑误差为负的情况，当误差（希望值减去温度 值）为负大时（说明温度高于希望值），若误差变化率也为负，这时误差有增大的趋势，为尽快消除已 有的负大误差并抑制误差变大，所以控制量的变化取正大（控制量增大，意味着喷水阀门开度增大， 喷水量增加，使得温度下降，下同） 。当误差为负而误差变化率为正时，系统本身已有减少误差的趋势，所以，为尽快消除误差而又不 超调，应取较小的控制量。当误差为负中时，控制量的变化应使误差尽快消

8、除，基于这种原则，控制量的变化选取同误差为 负大时相同。当误差为负小时，系统接近稳态，若误差变化微小时，选取控制量变化为正中，以抑制误差往负 方向变化；若误差变化为正时，系统本身有消除负小的误差的趋势，选取控制量变化为正小。上述选取控制量变化的原则：当误差大或较大时，选择控制量以尽快消除误差为主；当误差较小 时，选择控制量要注意防止超调，以系统稳定性为主要出发点。按着上述控制规则，可以得到温度偏差及偏差变化率与喷水阀门开度之间的模糊关系 R :R （E U）o（EC U ） E EC U（NBe NBec PBu）U（NMe NBec PBu ） U（ NSe NBec PMu ）UL L其中

9、角标E，EC U分别表示误差、误差变化率和控制量。设e、ec、u的等级量分别为 nfe、nfec、nfu，贝U R是一个nfe nfec行nfu列的矩阵。计算模糊关系矩阵R的子程序如F_Relation_2.m 所示。%莫糊关系计算子程序 F_Relation_2.cfun ctio n R, n,nE,n EC, nU, nfe,nfec,n fu,Me,Mec=fu nc（）%#3定义 E的模糊变量值仁PB,2=PM,3=PS,4=-O,5=+0,6=NS,7=NM,8=NB并输入表（5-5 ） # nE=8;%模糊变量E的语言值个数？控制表列数E=8 7 6 5 4 3 2 1;nfe=

10、14;%E的等级量个数Me=O 0 0 0 0 0 0 0 0 0 0.1 0.4 0.8 1.0;0 0 0 0 0 0 0 0 0 0.2 0.7 1.0 0.7 0.2;0 0 0 0 0 0 0 0.3 0.8 1.0 0.5 0.1 0 0;0 0 0 0 0 0 0 1.0 0.6 0.1 0 0 0 0;0 0 0 0 0.1 0.6 1.0 0 0 0 0 0 0 0;0 0 0.1 0.5 1.0 0.8 0.3 0 0 0 0 0 0 0;0.2 0.7 1.0 0.7 0.2 0 0 0 0 0 0 0 0 0;1.0 0.8 0.4 0.1 0 0 0 0 0 0 0

11、0 0 0;%#定义 EC的模糊变量值 1=PB,2=PM,3=PS,4=0,5=NS,6=NM,7=NB并输入表（5-2 ） #nEC=7;%模糊变量EC勺语言值个数？控制表行数EC=7 6 5 4 3 2 1;n fec=13;%EC的等级量个数Mec=0 0 0 0 0 0 0 0 0 0.1 0.4 0.8 1.0;0 0 0 0 0 0 0 0 0.2 0.7 1.0 0.7 0.2;0 0 0 0 0 0 0 0.9 1.0 0.7 0.2 0 0;0 0 0 0 0 0.5 1.0 0.5 0 0 0 0 0;0 0 0.2 0.7 1.0 0.9 0 0 0 0 0 0 0;0

12、.2 0.7 1.0 0.7 0.2 0 0 0 0 0 0 0 0;1.0 0.8 0.4 0.1 0 0 0 0 0 0 0 0 0;%#定义 U的模糊变量值 仁PB,2=PM,3=PS,4=0,5=NS,6=NM,7=NB并输入表（5-6 ） #nU=7;%模糊变量U的语言值个数U=7 6 5 4 3 2 1;nfu=15;%U的等级量个数Mu=0 0 0 0 0 0 0 0 0 0 0 0.1 0.4 0.8 1.0;0 0 0 0 0 0 0 0 0 0.2 0.7 1.0 0.7 0.2 0;0 0 0 0 0 0 0 0.4 1.0 0.8 0.4 0.1 0 0 0;0 0 0

13、 0 0 0 0.5 1.0 0.5 0 0 0 0 0 0;0 0 0 0.1 0.4 0.8 1.0 0.4 0 0 0 0 0 0 0;0 0.2 0.7 1.0 0.7 0.2 0 0 0 0 0 0 0 0 0;1.0 0.8 0.4 0.1 0 0 0 0 0 0 0 0 0 0 0;%#耕俞入模糊控制规则表（表 5-7）#nfc=8;%模糊变量E的语言值个数，模糊控制表列数 mfc=7;%模糊变量EC勺语言值个数，模糊控制表行数UC=1 1 2 2 2 3 4 4;1 1 2 2 2 3 4 4;1 1 2 3 3 4 5 6;1 2 3 4 4 5 6 7;2 2 3 4 5

14、5 6 7;2 3 4 5 6 6 7 7;3 4 5 6 6 7 7 7;%#计算 R= （ EX ECX U） #R=zeros（ nfe*nfec,n fu）; nfc%E X ECie=E（j）;iec=EC（i）; nfe nfecMec（iec,l） Reec（k,l）=Me（ie,k）; elseReec（k,l）=Mec（iec,l）; end%E xECX Uiu=UC（i,j）;n=0;n=n+1;for t=1: nfuif Reec（k,l）R（k,l） R（k,l）=Reecu（k,l）;然后双输入时，先根据 e、ec计算出相应的等级量 E、EC,再根据表（5-2）查

15、出相应的E, EC ,按式（5-6）计算出模糊控制器的输出。其模糊决策子程序如 F_Deduce_2.m所示。%双输入时模糊决策子程序 F_Deduce_2.mfunction FU=F_Deduce_2（fe,fec,R, n,n fe, nfec,nfu）%#5计算 ex EC#if fe（i）fec（j）feec（ n）=fe（i）;feec（ n）=fec（j）;%#5计算（E x EC）。R#if feec（i）R（i,l）fu（i）=feec（i）;fu（i）=R（i,l）;FU（l）=max（fu）;按着前面的设计过程，可以得到双输入模糊控制主程序，如 FC_MI_main.m所

16、示。%单输入单变量模糊控制仿真程序 FC_SI_main.m ，被控对象 W（s）=-1.1/（35s+1）A4）e0=0;%#设置e、ec、u的初始论域及其等级量论域，计算量化因子及比例因子 #em=2;EM=6;ecm=2;ECM=6;Kec=ECM/ecm;UM=7;um=2;U=-7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7;%#调用计算模糊关系矩阵 R的子程序F_Relation_2.m#R, n,nE,n EC, nU, nfe, nfec, nfu,Me,Mec=F_Relation_2; %# for k=1:%# 计算 FC 输岀 #%# 变量 e 模糊化处理 #% 计算输入e的确切量% 把e转换成等级量0;if （e1=0）&（ej=8;j=e1+EM+2;%得到e的模糊语言值mfE=Me（1,j）; nEif Me（i,j）mfE;Fi=i;mfE=Me（i,j）;fe=Me（Fi,:）;%# 变量 ec 模糊化处理 #ec=（e-eO）/DT; e0=e; 计算 ec 的确切量ec仁rou nd（Kec*ec）;% 把ec转换成等级量if ec1ECM;ec1=ECM;if ec1-ECM;e1