模糊控制實(shí)驗(yàn)

1、中南大學(xué)模糊控制課程實(shí)驗(yàn)報告學(xué)生姓名：彭雄威m指導(dǎo)教師：學(xué)院：信息科學(xué)與工程學(xué)院學(xué)號：114611167實(shí)驗(yàn)一：本系統(tǒng)設(shè)計基于MATLAB圖形模糊推理系統(tǒng)，設(shè)計步驟如下：打開MATLAB, 輸入指令fuzzy，打開模糊邏輯工具箱的圖形用戶界面窗口，新建一個Mamdani模糊推 理系統(tǒng)。增加一個輸入變量，將輸入變量命名為E、Ec,將輸出變量命名U。這樣就建立了 一個兩輸入單輸出的模糊推理系統(tǒng)。如圖1.1所示。圖1.1增加一個輸入變量設(shè)計模糊化模塊：設(shè)計隸屬度函數(shù)論域范圍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設(shè)計水位誤差E模塊le Editor: f

3、uzzfFile Edit Vi Optionsif (e is N&)司nd 網(wǎng)度 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)通過觀察器觀察規(guī)則情況在菜單view中的rules和surface選項(xiàng)分別對應(yīng)得是規(guī)則觀測器圖3-6曲面觀測器(6)保存編輯好的FIS文件實(shí)驗(yàn)二利用MATLAB軟件的M文件編輯器和實(shí)驗(yàn)

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文件的模糊控制查詢表，其中的數(shù)據(jù)在水位誤差E的論域?yàn)? 1,誤差變化EC的論 域?yàn)?0.1 0.1內(nèi)可以任意取值。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實(shí)驗(yàn)三利用MATLAB軟件的M文件編輯器（也可選擇C語言）完成模糊控制查詢表的計算。注意：全部工作都由程序?qū)崿F(xiàn)，不利用FI

7、S文件，不使用Fuzzy I具箱，包括模糊關(guān)系矩陣的 求取及模糊推理、模糊判決。表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推理語言規(guī)則表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 存到數(shù)組 C 中，C 是 56*13 的數(shù)組 end end%根據(jù)式 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實(shí)驗(yàn)四控制系統(tǒng)仿真實(shí)驗(yàn)。利用MATLAB的SIM