模糊控制實(shí)驗(yàn)_第1頁(yè)
模糊控制實(shí)驗(yàn)_第2頁(yè)
模糊控制實(shí)驗(yàn)_第3頁(yè)
模糊控制實(shí)驗(yàn)_第4頁(yè)
模糊控制實(shí)驗(yàn)_第5頁(yè)
已閱讀5頁(yè)，還剩10頁(yè)未讀， 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說(shuō)明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

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

2、模糊控制器的規(guī)則設(shè)計(jì)（5） 通過(guò)觀察器觀察規(guī)則情況在菜單view中的rules和surface選項(xiàng)分別對(duì)應(yīng)得是規(guī)則觀測(cè)器和曲面觀測(cè)器。 圖 3.7 規(guī)則觀測(cè)器圖 3.6 曲面觀測(cè)器 (6) 保存編輯好的FIS文件實(shí)驗(yàn)二：利用MATLAB軟件的M文件編輯器和實(shí)驗(yàn)一所生成的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ù)在水位

3、誤差E的論域?yàn)?1 1，誤差變化Ec的論域?yàn)?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 rule: 1x5 str

4、ucta=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實(shí)驗(yàn)三利用MATLAB軟件的M文件編輯器（也可選擇C語(yǔ)言）完成模糊控制查詢表的計(jì)算。注意：全部工作都由程序?qū)崿F(xiàn)，不利用FIS文件，不使用Fuzzy工具箱，包括模糊關(guān)系矩陣的求取及模糊推理、模糊判決。表1 偏差E的賦值表 eE

5、-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.710.70.3000

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

7、00000NS000.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.6,0.1,0,0,0,

8、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; 0,0,0,0,0,0,

9、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,0.7,1%7*13r

10、ulelist= 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ī)則表 %*計(jì)算C1A'* A = E(ii,:); %取A&#

11、39; 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) = Ra(j,i);elseA

12、_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計(jì)算完畢* %*計(jì)算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);elseRb(i,j) = B

13、(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)max = B_qux(j,

14、i);endendCiB(i) = max;end %*計(jì)算CiB'完畢* %*計(jì)算C1'=CiA'交CiB'*for i=1:13if CiA(i) >CiB(i)Ci(i) = CiB(i);elseCi(i) = CiA(i);endend %*計(jì)算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)行解模糊運(yùn)算,并映射為實(shí)際的控制量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

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說(shuō)明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽，若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間，僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請(qǐng)與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

評(píng)論

0/150

提交評(píng)論