模糊控制试验

1、F南大学模糊控制课程实验报告彭雄威学生姓名: 指导教师:18 / 14院:信息科学与工程学院t=r.号 :1146111671.1所示。实验一：本系统设计基于MATLAB图形模糊推理系统，设计步骤如下：打开MATLAB，输入指令fuzzy，打开模糊逻辑工具箱的图形用户界面窗口，新建一 个Mamdani模糊推理系统。(1)增加一个输入变量，将输入变量命名为E、Ec,将输出变量命名U。这样就建立了一个两输入单输出的模糊推理系统。如图图1.1增加一个输入变量(2)设计模糊化模块：设计隶属度函数论域范围寸 Be wb er ship Funct ion Edit or; fuz ziewbershij

2、) Funct ion Edit or; fuzzfJ Bewbership Funct ion Edit or; fuz图3.2设计水位误差E模块3.3设计水位误差Ec模块图3.4设计水位输出U模块(4)模糊控制器的规则设计> Rule Editoe: fuzzfPile E时itOptions2 IfCe is NG) a闻 Cec is NM) then (u is WB (1)3 If (fi is MeXrwJ Cec is NS) then (u is NM) (1)4 If Ce is we) a网 Cec is Z theru is NMi) 0】 i.lfCe is N

3、B) srwi Cec is PS) then (w is NS) Cl) 6 If Ce is MB)也网 Cec is PMD then ( is W(1) 7.ifCe is 血urtfA Cec is pB) then (u is Z)(1) S.lf Ce is Mvl) arT 佃c is MB) then (u is NB) (1>9. If e is M町 arid (ec is NM)thein (u is NM) 0)10. If (e is NM) and (ec is NS) then (u is NM)11. If (e is NM) and (ec is 2

4、then (u is NSj e isandec EThenu isNM NS ZPSPMrlMBNMNS ZPS PMrwk not notr CornectianWfeisht:O or© arwle = D3411 ITi*69晶1112131415Ifi17 IB ia 2D 21 22 旳 24Jb 26 2? je 29FIS 忖sme: fuzzf Melp Cl05e I1Delete ruleAdd rue I Change fuie (5)通过观察器观察规则情况在菜单view中的rules和surface选项分别对应得是规则观测器和曲面观测器。Kil l p Vi

5、 ever = "FiiTTfFileu-0haut n QPlot pcintE：QU1比池 1 iFtT njhf jrjn'jvn| up Eh1 CIOS(6 I图3.7规则观测器? Suxrace Viewer: fuzzfPile E时itOptionsK (input)X grids:15R 亡f. Iriput:42-2-4=0-2e” Y (input):Y grids:ec15k Z (outputXEvaluateHsip Gloss M*R&a血图3.6曲面观测器（6）保存编辑好的FIS文件Save FIS保存在（X）: I work&

6、 c? 1011 first, fiz I d iwi hlh bl llnnlll rwi II 百 gifuzif fii WBMimiBIHBIII HI IlTGizri f . £i s保存II*. £i a实验二：利用 MATLAB 软件的 M 文件编辑器和实验一所生成的 fuzzf.FIS 文件，在 M 文件编辑器中输入： a=readfis('fuzzf');evalfis(-0.5,-0.07;-0.5,0;-0.5,0.07; 0, -0.07;0,0;0,0.07;0.5,-0.07;0.5,0;0.5,0.07,a)便可得fuzzf.

7、FIS文件的模糊控制查询表，其中的数据在水位误差E的论域为-1 1,误差变化EC的论域为-0.10.1 内可以任意取值。a=readfis(' fuzzf ')name: ' fuzzf 'type: 'mamdani' andMethod: 'min'orMethod: 'max'defuzzMethod: 'Centroid'impMethod: 'min'aggMethod: 'max'input: 1x2 struCtoutput: 1x1 struCtrul

8、e: 1x5 struCt a=readfis(' fuzzf ');evalfis(-0.5,-0.07;-0.5,0;-0.5,0.07; 0, -0.07;0,0;0,0.07;0.5,-0.07;0.5,0;0.5,0.07,a) ans =-0.2000-0.4444-0.46670.1363-0.0014-0.22710.40000.44120.1333实验三利用MATLAB软件的M文件编辑器（也可选择 C语言）完成模糊控制查询表的计算。注意：全部工作都由程序实现，不利用FIS文件，不使用Fuzzy工具箱，包括模糊关系矩阵的求取及模糊推理、模糊判决。表1偏差E的赋值

9、表-6-5-4-3-2-10123456NB10.80.70.40.40.10000000NM0.20.710.70.300000000NS00.10.30.710.70.2000000NO00000.10.61000000P000000010.60.10000PS0000000.20.710.70.30.10PM000000000.20.710.70.3PB000000000.10.40.70.81表2偏差变化率EC的赋值表eceC-6-5-4-3-2-10123456NB10.70.30000000000NM0.30.710.70.300000000NS000.30.710.70.3000

10、000000000.30.710.70.30000PS0000000.30.710.70.300PM000000000.30.710.70.3PB00000000000.30.71表3推理语言规则表7NBNMNS0PSPMPBNBPBPBPBPBPM00NMPBPBPBPBPM00NSPMPMPMPM0NSNSNOPMPMPS0NSNMNMPOPMPMPS0NSNMNMPSPSPS0NMNMNMNMPM00NMNSNSNSNSPB00NMNSNSNSNS表4控制量U的赋值表U-6-5-4-3-2-10123456NB10.70.30000000000NM0.30.710.70.30000000

11、0NS000.30.710.70.3000000O00000.30.710.70.30000PS0000000.30.710.70.300PM000000000.30.710.70.3PB00000000000.30.71程序如下:cic;clear;E = 1,0.8,0.7,0.4,0.1,0,0,0,0,0,0,0,0;0.2,0.7,1,0.7,0.3,0,0,0,0,0,0,0,0;0,0.1,0.3,0.7,1,0.7,0.2,0,0,0,0,0,0;0,0,0,0,0.1,0.6,1,0,0,0,0,0,0;0,0,0,0,0,0,1,0.6,0.1,0,0,0,0;0,0,0,

12、0,0,0,0.2,0.7,1,0.7,0.3,0.1,0;0,0,0,0,0,0,0,0,020.7,1,0.7,0.3;0,0,0,0,0,0,0,0,0.1,0.4,0.7,0.8,1%8*13EC =1,0.7,0.3,0,0,0,0,0,0,0,0,0,0;0.3,0.7,1,0.7,0.3,0,0,0,0,0,0,0,0;0,0,0.3,0.7,1,0.7,0.3,0,0,0,0,0,0;0,0,0,0,0.3,0.7,1,0.7,0.3,0,0,0,0;0,0,0,0,0,0,0.3,0.7,1,0.7,0.3,0,0;0,0,0,0,0,0,0,0,0.3,0.7,1,0.7,

13、0.3;0,0,0,0,0,0,0,0,0,0,0.3,0.7,1%7*13U=1,0.7,0.3,0,0,0,0,0,0,0,0,0,0;0.3,0.7,1,0.7,0.3,0,0,0,0,0,0,0,0;0,0,0.3,0.7,1,0.7,0.3,0,0,0,0,0,0;0,0,0,0,0.3,0.7,1,0.7,0.3,0,0,0,0;0,0,0,0,0,0,0.3,0.7,1,0.7,0.3,0,0;0,0,0,0,0,0,0,0,0.3,0.7,1,0.7,0.3;0,0,0,0,0,0,0,0,0,0,0.3,0.7,1%7*13 rulelist= 1,1,1,1,2,4,4;1

14、,1,1,1,2,4,4;2,2,2,2,4,5,5;2,2,3,4,5,6,6;2,2,3,4,5,6,6;334,6,6,6,6;4,4,6,7,7,7,7;4,4,6,7,7,7,7%8*7for iii=1:13%偏差E的模糊值控制for jjj=1:13 for ii=1:8 for jj=1:7A_rulelist = rulelist(ii,jj); % 查模糊规则表 %* 计算 C1A'* A = E(ii,:);%取 A'C_A = U(A_rulelist,:);%取 Cifor i=1:13%求 R1Afor j=1:13 if(A(i) > C_A

15、(1,j) Ra(i,j) = C_A(1,j);%偏差率EC的模糊值控制%偏差E的模糊规则控制%偏差变化率EC的模糊规则控制elseRa(i,j) = A(i);endendendAA = zeros(1,13); %取 AA=A' 如 AA=1,0,0,0,0,0,0,0,0,0,0,0,0; %AA(1,iii) = E(ii,iii);AA(1,iii) = 1;for i=1:13 for j=1:13%先取小if(AA(j) > Ra(j,i) A_qux(j,i) = Ra(j,i);elseA_qux(j,i) = AA(j);endendendendfor i=

16、1:13max = A_qux(1,i);for j=1:13if(max <A_qux(j,i) max = A_qux(j,i); end%再取大，求和 CiAend CiA(i) = max;%*CiA%* 计算B = EC(jj,:);C_B = U(A_rulelist,:);计算完毕 *CiB'*%取 B'?%取 Ci?for i=1:13%求 R1Bfor j=1:13if(B(1,i) > C_B(1,j)elseend endendRb(i,j) = C_B(1,j);Rb(i,j) = B(1,i);BB = zeros(1,13); %取 BB

17、=B' 如 BB=1,0,0,0,0,0,0,0,0,0,0,0,0; %BB(1,jjj) = EC(jj,jjj);BB(1,jjj) = 1;for i=1:13 for j=1:13%先取小if(BB(j) >Rb(j,i) B_qux(j,i) = Rb(j,i);elseB_qux(j,i) = BB(j);endendendendfor i=1:13max = B_qux(1,i); for j=1:13 if(max <B_qux(j,i)max = B_qux(j,i); end end CiB(i) = max;%再取大，求和 CiA%* 计算 CiB&

18、#39;完毕 *%* 计算 C1'=CiA'交 CiB'* for i=1:13if CiA(i) >CiB(i)Ci(i) = CiB(i);elseCi(i) = CiA(i);endendend%* 算 C1'完毕 *C(ii-1)*7+jj,:) = Ci; % 将 Ci 存到数组 C 中，C是 56*13 的数组end%根据式Ui=C'=C1'并C2'并 for i=1:13 max = C(1,i);for j=1:56并 3'if(max < C(j,i)max = C(j,i);endendUi(i)

19、= max;end%用重心法进行解模糊运算，并映射为实际的控制量 sum_fe nz = 0;sum_fe nm = 0;for i=1:13sum_fe nz = sum_fe nz + (i-7)*Ui(i);sum_fe nm = sum_fe nm + Ui(i);endcore = sum_fe nz/sum_fe nm;U_c on trol(iii,jjj) = core;end endroun d(U_co ntrol.*10000)/10000;%对控制量保留小数点后两位U_con trol =U_c on trol 计算结果：U control =jyju-b. 2 泊,-

20、b.JfcUU-b.-b.JbJU-b.JJSJ-4.-U. Mb创-N 11 JT. JUUUjQ-SlOMO-s.orco-4.舸 W-D.C050-t附划-ISS?-3.3571-1.T9F50.2414-4-保 3-4. b21?-4世卫3-4.bai/-4.t.!?23-ZSJ2-l.S3rb-1.1211-t.fcy?bO.tJliiu.tsda-d. ssai-d.疋自IT一址飢-4.2581-4.2551-i. 25-2.9302-2. i9t9-l.d211-0.335-o 一制附1.tD03I. CQOCl-4.umo-4.ULOQ-3.?686-3.47 阳一乳5羽-2.4318-1. TAS'b-0. 44(0-0.13ei4C.lfiOJj.eooji.esie-Q+ OCdQ7+范汕-3h0769-2h i65?-E* IL