matlab計(jì)算模糊控制表

1、Matlab 計(jì)算模糊控制表實(shí)現(xiàn)程序如下：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,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,

2、0.1,0.4,0.7,0.8,1 %誤差的隸屬度函數(shù)表， 8*13 的矩陣eb=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 %誤差變化率的隸屬度函數(shù)表， 7*13 的矩陣；u=1,0

3、.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 %輸出控制量 u,7*13 的矩陣；rulelist=1,1,1,1,2,4,4;1,1,1,1,2,4,4;2,2,2,2,4,5,5;2,2,

4、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 %控制規(guī)則表，為 8*7=56 條規(guī)則；for ey=1:13 % 輸入變量e的13個(gè)取值：-6 6;for eby=1:13 %輸入變量eb的13個(gè)取值：-6 6;for ex=1:8% 輸入變量e的8個(gè)語(yǔ)言變量 NB-PB;for ebx=1:7 % 輸入變量eb的7個(gè)語(yǔ)言變量 NB-PB;arule=rulelist(ex,ebx); %取一條規(guī)則;C=u(arule,:);%取這條規(guī)則的C的隸屬度；代表取一行A=e(ex,:);%取這條規(guī)則的A的隸屬度；代

5、表取一行for i=1:13;for j=1:13;if (A(i)C(1,j);Ra(i,j)= C(1,j);elseRa(i,j)=A(i);endendend%算A C的蘊(yùn)含關(guān)系；AP=zeros(1,13);AP(1,ey)=e(ex,ey);for i=1:13;for j=1:13;if (AP(j)Ra(j,i);GDA(j,i)= Ra(j,i);elseGDA(j,i)=AP(j);endenden d%求合成運(yùn)算開(kāi)始的取小；for i=1:13;max=GDA(1,i);for j=1:13;if (maxC(1,j);Rb(i,j)= C(1,j);elseRb(i,j

6、)=B(i);end endend%算BC的蘊(yùn)含關(guān)系;BP=zeros(1,13);BP(1,eby)=eb(ebx,eby);for i=1:13;for j=1:13;if (BP(j)Rb(j,i);GDB(j,i)= Rb(j,i);elseGDB(j,i)=BP(j);endendend%求合成運(yùn)算開(kāi)始的取小;for i=1:13;max=GDB(1,i);for j=1:13;if (maxCPB(i)CP(i)=CPB(i);elseCP(i)=CPA(i);endend %第一個(gè)CP1計(jì)算結(jié)束；將來(lái)C是56*13的矩陣;CB(ex-1)*7+ebx,:)=CP;endendfo

7、r i=1:13;max=CB(1,i);for j=1:56;if(maxCB(j,i)max=CB(j,i);endendU(i)=max;end %56個(gè)CP整理結(jié)束；sumFZ=0;sumFM=0;for i=1:13;sumFZ=sumFZ+(i-7)*U(i); sumFM=sumFM+U(i);endcore=sumFZ/sumFM;UB(ey,eby)=core; % 一個(gè)控制表中的元素算出來(lái)了endendUB=round(UB*100)/100;% 保留兩位 display(UB);程序運(yùn)行結(jié)果如下：Coffirniand Wind w R XUE =Colunns 1 th

8、iough 3-5. 3500-B. 2400-5. 3500-5. 2；ao-5.3500-6,2400-1.69Q0-4. 2000-5.0000-4.9500-5.0000-1.9500-5.0000-4.9500-3. S600-3.7)00-4. 6900-4.5200-4,6900-4. 5200-4. 6900-4.5200-3.0500-2. 9300-4.2600-4.2600-4.2S00-4, 2500-4- 2600-4.2600-2.9300-2. 2900-4, mo-4. itiQQ-X7800-3, 750 Q-3,4700-3.4200-2+ 43QQ-l.

9、7?(J0-4. aooo4.0000-3.3600-3.0800-2. J 700-2.J2C0-L5O&0-1.0500-3.5900-3. 5500-5.9305-2. SOW-0. 9600-G.510000.5100-2.9200-3. 9200-2.3300-1.9100-0. 2600L.0500L GOOD2. 1200-L3I00-1.79Q0-0.57Q(J-0.3100Q.440QI. 79QQ2, 2QQ3200-l oaoo-1.0000D.2500o. 9ao1.42002.29002- 93004,2eaa-0.5800-0.64000.69001. 42001.

10、94002.93003. 05004. 5200-0. 2300-0.2100L 1200L7&002.36003.7100X SflOO4. 950000L29QQZ 00002.71004.26001. SBOO5, 2400Columns 9through 13-2.7100-2. 0000-1.290000-2,3500-1.7900-1 12006 240Q0,2300-I. 9400-L450O-0.6000o. 6ao0. 5S0C-1. 4200-0. 9400-0.2500l.OQOO1. 0000-0.4400-0.04000.16001.60001.63000.測(cè)r9lQQ2.33052h 92002.92000. 960D2.60002.93003. 55003. 59002.4：003.0BOO3.3A004. 00004.00003.47003.76003.78004. 00004. 00004.26004.26004.26004*2