nmctrianglech.m

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

a1=0.5367,a2=0.7813,a3=0.2023,a4=0.6586,a5=0.2977,a6=0.4836,a7=0.3695,a8=0.1570,a9=0.4969,a10=0.4039,a11=0.7531,a12= 0.8922,
%a1=0.2734,a2=0.8602,a3=0.1945,a4=0.7977,a5=0.3789,a6=0.4578,a7=0.8406,a8=0.8688,a9=0.4914,a10=0.2250,a11=0.7813,a12=0.7125;
%a1=0.2688,a2=0.2930,a3=0.1133,a4=0.6664,a5=0.5234,a6=0.8469,a7=0.8063,a8=0.2102,a9=0.6008,a10=0.1734,a11=0.7117,a12=0.8773,
%a1=0.2914,a2=0.1453,a3=0.1398,a4=0.7969,a5=0.4609,a6=0.8781,a7=0.7664,a8=0.2383,a9=0.8445,a10=0.6039,a11=0.8094,a12=0.3648,
%a1=0.5,a2=0.5,a3=0.5,a4=0.5,a5=0.5,a6=0.5,a7=0.5,a8=0.5,a9=0.5,a10=0.5,a11=0.5,a12=0.5,
a=newfis('tm');
a=addvar(a,'input','e',[-4,4]);
a=addmf(a,'input',1,'NB','trimf',[-4,-4,-4+1/(1-a1)]);
a=addmf(a,'input',1,'NS','trimf',[-2-1/(1-a1),-2,-2+1/(1-a2)]);
a=addmf(a,'input',1,'ZO','trimf',[-1/(1-a2),0,1/(1-a3)]);
a=addmf(a,'input',1,'PS','trimf',[2-1/(1-a3),2,2+1/(1-a4)]);
a=addmf(a,'input',1,'PB','trimf',[4-1/(1-a4),4,4]);
a=addvar(a,'input','ec',[-4,4]);
a=addmf(a,'input',2,'NB','trimf',[-4,-4,-4+1/(1-a5)]);
a=addmf(a,'input',2,'NS','trimf',[-2-1/(1-a5),-2,-2+1/(1-a6)]);
a=addmf(a,'input',2,'ZO','trimf',[-1/(1-a6),0,1/(1-a7)]);
a=addmf(a,'input',2,'PS','trimf',[2-1/(1-a7),2,2+1/(1-a8)]);
a=addmf(a,'input',2,'PB','trimf',[4-1/(1-a8),4,4]);
a=addvar(a,'output','u',[-4,4]);
a=addmf(a,'output',1,'NB','trimf',[-4,-4,-4+1/(1-a9)]);
a=addmf(a,'output',1,'NS','trimf',[-2-1/(1-a9),-2,-2+1/(1-a10)]);
a=addmf(a,'output',1,'ZO','trimf',[-1/(1-a10),0,1/(1-a11)]);
a=addmf(a,'output',1,'PS','trimf',[2-1/(1-a11),2,2+1/(1-a12)]);
a=addmf(a,'output',1,'PB','trimf',[4-1/(1-a12),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('nmc');
y=ScopeData1.signals.values;
t=ScopeData1.time;
e=abs(50-y);
eval=trapz(t,e),