模糊控制实验

1、中南大学模糊控制课程实验报告学生姓名：彭雄威m指导教师：学院：信息科学与工程学院学号：114611167实验一：本系统设计基于MATLAB图形模糊推理系统，设计步骤如下：打开MATLAB, 输入指令fuzzy，打开模糊逻辑工具箱的图形用户界面窗口，新建一个Mamdani模糊推 理系统。增加一个输入变量，将输入变量命名为E、Ec,将输出变量命名U。这样就建立了 一个两输入单输出的模糊推理系统。如图1.1所示。图1.1增加一个输入变量设计模糊化模块：设计隶属度函数论域范围C lebership Function Editors 宠尤fFile Edit ViewFIS VariablesEZIM

2、e u o ftCCurrent VariableNameTypeRangeDisplay RangeinputCurrent Membership Function (click on MF tcNameTypeParamsHelpMEznr-3-1Ready lenberhip Ftmction Editor:File Edit Vi evinput variable 8output vriLjrWe u1NameNameCiisptav RangeCurrert Membership Function (click on IVF Id s图3-2设计水位误差E模块le Editor: f

3、uzzfFile Edit Vi Optionsif (e is N&)司nd 网度 N日)then (u N l旧1)It (e Is NB) and (ec fs NM) then fu 壕 MB) (1)h (e is NB) and (ec is NS) then (u is N的(1)If (e is NB；i and (ec is Z) then (d 135, If (b is NB) and (ec is PS) then (uts NS) (1)If (e Is NB) and (ec is PM) then (u is NS) (1)It (e is NB) and (ec

4、 is PB) then fu ts Z) (1)8. If (e is bM) and (ec is 阳)then (u is 浦)(1)9 lt(& is NM) and ec Is NMj then (u is NM) (1)W (e js NbO and (ec is NS) then (u is NW) (1)If (e 溶 N* and (ec is I) then (u is IMS) (1)(5)通过观察器观察规则情况在菜单view中的rules和surface选项分别对应得是规则观测器图3-6曲面观测器(6)保存编辑好的FIS文件实验二利用MATLAB软件的M文件编辑器和实验

5、一所生成的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.FIS文件的模糊控制查询表，其中的数据在水位误差E的论域为1 1,误差变化EC的论 域为.0.1 0.1内可以任意取值。a=readfis( fuzzf) a =name: fuzzftype: mamdaniandMethod: minorMethod: maxdefuzzMethod: centroidimpMe

6、thod: minaggMethod: maxinput: 1x2 structoutput: lxl structrule: 1x5 structa=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语言）完成模糊控制查询表的计算。注意：全部工作都由程序实现，不利用FI

7、S文件，不使用Fuzzy I具箱，包括模糊关系矩阵的 求取及模糊推理、模糊判决。表1偏差E的赋值表X-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的赋值表X-6-5-4-3-2-10123456NB10.70.30000000

8、000NM0.30.710.70.300000000NS000.30.710.70.3000000000000.30.710.70.30000PS0000000.30.710.70.300PM000000000.30.710.70.3PB00000000000.30.71表3推理语言规则表uNBNMNS0PSPMPBNBPBPBPBPBPM00NMPBPBPBPBPM00NSPMPMPMPM0NSNSNOPMPMPS0NSNMNMPOPMPMPS0NSNMNMPSPSPS0NMNMNMNMPM00NMNSNSNSNSPB00NMNSNSNSNS表4控制量U的赋值表9O_O_OoocnTo oo

9、 ou orHn oJ寸ozoSooooor*.6寸oooornIrnmoooor*.z6oCMooornIcnOoooz bz6oOoooC_A(l,j)Ra(iJ)= C_A(1J); else Ra(iJ) = A(i);endendendAA = zeros(l,13); %取 AA=A如 AA=1,0,0,0,0,0,0,0,0,0,0,0,0; %AA(l,iii) =E(iijii);AA(lJii)=l;for i=l:13% 先取小forj=l:13if(AA(j) Ra(j,i)A_qux(j,i) = Ra(jJ);elseA_qux(jJ) = AA(j);endend

10、endfor i=l:13%再取大求和CiAmax = A_qux(l,i);forj=l:13if(max Rb(jJB_qux(j,i) = Rb(jJ);elseB_qux(j,i) = BBj);endendendfori=l:13 %再取大，求和 CiA max = B_qux(lJ);forj=l:13if(max CiB(i)Ci(i) = CiB(i);elseCi(i)= CiA(i);endend%*计算 CT 完毕*C(ii-l)*7+jj,:) = Ci; % 将 Ci 存到数组 C 中，C 是 56*13 的数组 end end%根据式 ui=c=cr 并CZ并并3*for i=l:13max = C(1J);for j=l:56if(max 5.2353-5. 0000-4.9474-4.6923-4.5217-4.2581-4.2581-3.4706-3.4194-2.4667-2.1176-0.9600-0.5094实验四控制系统仿真实验。利用MATLAB的SIM