nmctrapch.m

求模糊控制器的梯型隶属函数优化后的最小目标值 请大家指教

a1=0.4352,a2=0.3867,a3=0.2672,a4=0.8844,m1=0.0070,m2=0.5352,m3=0.0531,m4=0.3141,a5=0.4578,a6=0.1344,a7=0.6844,a8=0.2648,m5=0.1641,m6=0.5406,m7=0.6391,m8=0.5648,a9=0.2289,a10=0.2109,a11=0.8547,a12=0.5102,m9=0.5242,m10=0.3672,m11=0.3414,m12=0.3578,
%a1=0.8195,a2=0.7359,a3=0.1844,a4=0.6914,m1=0.4453,m2=0.6320,m3=0.0289,m4=0.6969,a5=0.8664,a6=0.1984,a7=0.3867,a8=0.6828,m5=0.0766,m6=0.2125,m7=0.3891,m8=0.5328,a9=0.6875,a10=0.6937,a11=0.8203,a12=0.2055,m9=0.2531,m10=0.0305,m11=0.2008,m12=0.0898,
%a1=0.6141,a2=0.7539,a3=0.6875,a4=0.6172,m1=0.4836,m2=0.0617,m3=0.5938,m4=0.2188,a5=0.7477,a6=0.8992,a7=0.2594,a8=0.4422,m5=0.2414,m6=0.5102,m7=0.1203,m8=0.4125,a9=0.8430,a10=0.6164,a11=0.8352,a12=0.5016,m9=0.2766,m10=0.0563,m11=0.0063,m12=0.0867,
%a1=0.8977,a2=0.1008,a3=0.1320,a4=0.2039,m1=0.6297,m2=7.8125e-004,m3=0.7203,m4=0.2703,a5=0.3984,a6=0.6406,a7=0.1273,a8=0.1258,m5=0.7234,m6=0.5180,m7=0.0047,m8=0.0703,a9=0.4563,a10=0.5789,a11=0.1477,a12=0.1398,m9=0.0258,m10=0.1180,m11=0.4133,m12=0.1523,
%a1=0.6781,a2=0.3195,a3=0.8102,a4=0.6461,m1=0.0986,m2=0.0596,m3=0.6621,m4=0.4893,a5=0.1453,a6=0.1656,a7=0.6125,a8=0.8750,m5=0.7568,m6=0.2686,m7=0.4414,m8=0.1963,a9=0.8047,a10=0.8797,a11=0.6000,a12=0.5789,m9=0.2871,m10=0.1445,m11=0.4668,m12=0.3896,
%a1=0.3602,a2=0.7188,a3=0.2797,a4=0.3906,m1=0.4328,m2=0.4367,m3=0.0195,m4=0.0750,a5=0.3594,a6=0.3117,a7=0.5203,a8=0.1195,m5=0.4828,m6=0.2445,m7=0.7602,m8=0.3914,a9=0.5234,a10=0.2469,a11=0.8516,a12=0.8859,m9=0.5164,m10=0.4258,m11=0.5344,m12=0.1734,
%a1=0.3117,a2=0.7375,a3=0.3617,a4=0.3703,m1=0.2164,m2=0.2156,m3=0.0063,m4=0.4383,a5=0.1789,a6=0.6281,a7=0.5000,a8=0.4523,m5=0.4313,m6=0.5320,m7=0.7008,m8=0.2680,a9=0.6000,a10=0.1250,a11=0.8984,a12=0.5578,m9=0.7281,m10=0.6023,m11=0.0469,m12=0.4000,
%a1=0.7656,a2=0.1492,a3=0.3711,a4=0.5008,m1=0.2695,m2=0.7102,m3=0.0289,m4=0.3961,a5=0.4633,a6=0.6156,a7=0.1008,a8=0.5930,m5=0.1852,m6=0.0602,m7=0.6773,m8=0.0227,a9=0.8039,a10=0.1492,a11=0.8102,a12=0.3359,m9=0.6336,m10=0.0453,m11=0.1508,m12=0.4625,
%a1=0.7336,a2=0.2109,a3=0.3133,a4=0.4781,m1=0.8213,m2=0.0488,m3=0.9824,m4=0.0996,a5=0.1805,a6=0.7852,a7=0.5000,a8=0.4680,m5=0.9297,m6=0.2998,m7=0.1182,m8=0.4053,a9=0.4930,a10=0.6086,a11=0.3492,a12=0.5203,m9=0.1533,m10=0.3740,m11=0.6016,m12=0.5498,
%a1=0.5,a2=0.5,a3=0.5,a4=0.5,m1=0.5,m2=0.5,m3=0.5,m4=0.5,a5=0.5,a6=0.5,a7=0.5,a8=0.5,m5=0.5,m6=0.5,m7=0.5,m8=0.5,a9=0.5,a10=0.5,a11=0.5,a12=0.5,m9=0.5,m10=0.5,m11=0.5,m12=0.5,
a=newfis('tm');
a=addvar(a,'input','e',[-4,4]);
a=addmf(a,'input',1,'NB','trapmf',[-4,-4,-4+m1,-4+m1+(1-m1)/(1-a1)]);
a=addmf(a,'input',1,'NS','trapmf',[-2-m1-(1-m1)/(1-a1),-2-m1,-2+m2,-2+m2+(1-m2)/(1-a2)]);
a=addmf(a,'input',1,'ZO','trapmf',[-m2-(1-m2)/(1-a2),-m2,m3,m3+(1-m3)/(1-a3)]);
a=addmf(a,'input',1,'PS','trapmf',[2-m3-(1-m3)/(1-a3),2-m3,2+m4,2+m4+(1-m4)/(1-a4)]);
a=addmf(a,'input',1,'PB','trapmf',[4-m4-(1-m4)/(1-a4),4-m4,4,4]);
a=addvar(a,'input','ec',[-4,4]);
a=addmf(a,'input',2,'NB','trapmf',[-4,-4,-4+m5,-4+m5+(1-m5)/(1-a5)]);
a=addmf(a,'input',2,'NS','trapmf',[-2-m5-(1-m5)/(1-a5),-2-m5,-2+m6,-2+m6+(1-m6)/(1-a6)]);
a=addmf(a,'input',2,'ZO','trapmf',[-m6-(1-m6)/(1-a6),-m6,m7,m7+(1-m7)/(1-a7)]);
a=addmf(a,'input',2,'PS','trapmf',[2-m7-(1-m7)/(1-a7),2-m7,2+m8,2+m8+(1-m8)/(1-a8)]);
a=addmf(a,'input',2,'PB','trapmf',[4-m8-(1-m8)/(1-a8),4-m8,4,4]);
a=addvar(a,'output','u',[-4,4]);
a=addmf(a,'output',1,'NB','trapmf',[-4,-4,-4+m9,-4+m9+(1-m9)/(1-a9)]);
a=addmf(a,'output',1,'NS','trapmf',[-2-m9-(1-m9)/(1-a9),-2-m9,-2+m10,-2+m10+(1-m10)/(1-a10)]);
a=addmf(a,'output',1,'ZO','trapmf',[-m10-(1-m10)/(1-a10),-m10,m11,m11+(1-m11)/(1-a11)]);
a=addmf(a,'output',1,'PS','trapmf',[2-m11-(1-m11)/(1-a11),2-m11,2+m12,2+m12+(1-m12)/(1-a12)]);
a=addmf(a,'output',1,'PB','trapmf',[4-m12-(1-m12)/(1-a12),4-m12,4,4]);
f1=[1 1 1 1 1];
f2=[1 2 1 1 1];
f3=[1 3 2 1 1];
f4=[1 4 2 1 1];
f5=[1 5 3 1 1];
f6=[2 1 1 1 1];
f7=[2 2 2 1 1];
f8=[2 3 2 1 1];
f9=[2 4 3 1 1];
f10=[2 5 3 1 1];
f11=[3 1 2 1 1];
f12=[3 2 2 1 1];
f13=[3 3 3 1 1];
f14=[3 4 3 1 1];
f15=[3 5 4 1 1];
f16=[4 1 2 1 1];
f17=[4 2 3 1 1];
f18=[4 3 3 1 1];
f19=[4 4 4 1 1];
f20=[4 5 4 1 1];
f21=[5 1 3 1 1];
f22=[5 2 3 1 1];
f23=[5 3 4 1 1];
f24=[5 4 4 1 1];
f25=[5 5 5 1 1];
rulelist=[f1;f2;f3;f4;f5;f6;f7;f8;f9;f10;f11;f12;f13;f14;f15;f16;f17;f18;f19;f20;f21;f22;f23;f24;f25];
a=addrule(a,rulelist);
a=setfis(a,'DefuzzMethod','mom');
    sim('nmccc');
    y=ScopeData1.signals.values;
    t=ScopeData1.time;
    e=abs(50-y);
    eval2=trapz(t,e),