模糊控制實驗演示教學(xué)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上 中南大學(xué)模糊控制課程實驗報告學(xué)生姓名： 彭雄威 指導(dǎo)教師： 袁 艷 學(xué) 院： 信息科學(xué)與工程學(xué)院 學(xué) 號： 專心-專注-專業(yè) 實驗一：本系統(tǒng)設(shè)計基于MATLAB圖形模糊推理系統(tǒng)，設(shè)計步驟如下：打開MATLAB，輸入指令fuzzy，打開模糊邏輯工具箱的圖形用戶界面窗口，新建一個Mamdani模糊推理系統(tǒng)。（1） 增加一個輸入變量，將輸入變量命名為E、Ec，將輸出變量命名U。這樣就建立了一個兩輸入單輸出的模糊推理系統(tǒng)。如圖1.1所示。 圖1.1增加一個輸入變量(2) 設(shè)計模糊化模塊：設(shè)計隸屬度函數(shù)論域范圍 圖3.2設(shè)計水位誤差E模塊 3.3設(shè)計水位誤差Ec模塊 圖3.4

2、設(shè)計水位輸出U模塊（4）模糊控制器的規(guī)則設(shè)計（5） 通過觀察器觀察規(guī)則情況在菜單view中的rules和surface選項分別對應(yīng)得是規(guī)則觀測器和曲面觀測器。 圖 3.7 規(guī)則觀測器圖 3.6 曲面觀測器 (6) 保存編輯好的FIS文件實驗二：利用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.FIS文件的模糊控制

3、查詢表，其中的數(shù)據(jù)在水位誤差E的論域為-1 1，誤差變化Ec的論域為-0.1 0.1內(nèi)可以任意取值。a=readfis(' fuzzf ')a = name: ' fuzzf ' type: 'mamdani' andMethod: 'min' orMethod: 'max' defuzzMethod: 'centroid' impMethod: 'min' aggMethod: 'max' input: 1x2 struct output: 1x1 struct r

4、ule: 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.4667 0.1363 -0.0014 -0.2271 0.4000 0.44120.1333實驗三利用MATLAB軟件的M文件編輯器（也可選擇C語言）完成模糊控制查詢表的計算。注意：全部工作都由程序?qū)崿F(xiàn)，不利用FIS文件，不使用Fuzzy工具箱，包括模糊關(guān)系矩陣的求取及模糊推理、模糊判決。表

5、1 偏差E的賦值表 eE-6-5-4-3-2-10123456NB10.80.70.40.40.10000000NM0.20.710.70.300000000NS00.10.30.710.70.2000000NO00000.10.61000000PO00000010.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

6、.710.70.3000000O00000.30.710.70.30000PS0000000.30.710.70.300PM000000000.30.710.70.3PB00000000000.30.71表3 推理語言規(guī)則表EcENBNMNSOPSPMPBNBPBPBPBPBPMOONMPBPBPBPBPMOONSPMPMPMPMONSNSNOPMPMPSONSNMNMPOPMPMPSONSNMNMPSPSPSONMNMNMNMPMOONMNSNSNSNSPBOONMNSNSNSNS表4 控制量U的賦值表U-6-5-4-3-2-10123456NB10.70.30000000000NM0.30

7、.710.70.300000000NS000.30.710.70.3000000O00000.30.710.70.30000PS0000000.30.710.70.300PM000000000.30.710.70.3PB00000000000.30.71程序如下：clc;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.

8、6,0.1,0,0,0,0; 0,0,0,0,0,0,0.2,0.7,1,0.7,0.3,0.1,0; 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.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;

9、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*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

10、,0.7,1%7*13rulelist= 1,1,1,1,2,4,4; 1,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; 3,3,4,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 %偏差率EC的模糊值控制 for ii=1:8 %偏差E的模糊規(guī)則控制 for jj=1:7 %偏差變化率EC的模糊規(guī)則控制A_rulelist = rulelist(ii,jj); %查模糊規(guī)則表 %*計算C1A'* A = E(

11、ii,:); %取A' C_A = U(A_rulelist,:); %取Ci for i=1:13 %求R1Afor j=1:13if(A(i) > C_A(1,j)Ra(i,j) = C_A(1,j);elseRa(i,j) = A(i);endendend AA = 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) = R

12、a(j,i);elseA_qux(j,i) = AA(j);endendend for i=1:13 %再取大，求和CiAmax = A_qux(1,i);for j=1:13if(max <A_qux(j,i)max = A_qux(j,i);endendCiA(i) = max;end %*CiA計算完畢* %*計算CiB'* B = EC(jj,:); %取B' ? C_B = U(A_rulelist,:); %取Ci ? for i=1:13 %求R1Bfor j=1:13if(B(1,i) > C_B(1,j)Rb(i,j) = C_B(1,j);els

13、eRb(i,j) = B(1,i);endendend BB = zeros(1,13); %取BB=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);endendend for i=1:13 %再取大，求和CiAmax = B_qux(1,i);for j=1:13if(max <B_qux(j,i)ma

14、x = B_qux(j,i);endendCiB(i) = max;end %*計算CiB'完畢* %*計算C1'=CiA'交CiB'*for i=1:13if CiA(i) >CiB(i)Ci(i) = CiB(i);elseCi(i) = CiA(i);endend %*計算C1'完畢* C(ii-1)*7+jj,:) = Ci; % 將Ci存到數(shù)組C中，C是56*13的數(shù)組endend %根據(jù)式Ui=C'=C1'并C2'并并3'for i=1:13max = C(1,i);for j=1:56if(max < C(j,i)max = C(j,i);endendUi(i) = max;end %用重心法進(jìn)行解模糊運算,并映射為實際的控制量sum_fenz = 0;sum_fenm = 0;for i=1:13sum_fenz = sum_fenz + (i-7)*Ui(i);sum_fenm = sum_fenm + Ui(i);endcore = sum_fenz/sum_fenm;U_control(iii