ex36.asv

运动补偿编码和检测

pc1=newfis('pc1');
pc1=addvar(pc1,'input','temperature',[0 30]);
pc1=addvar(pc1,'input','pressure',[0 3]);
pc1=addvar(pc1,'output','increment',[-10 10]);
pc1=addmf(pc1,'input',1,'cold','trapmf',[0 0 10 20]);
pc1=addmf(pc1,'input',1,'hot','trapmf',[10 20 30 30]);
pc1=addmf(pc1,'input',2,'normal','trimf',[0 1 2]);
pc1=addmf(pc1,'input',2,'high','trimf',[ 1 2 3 3]);
pc1=addmf(pc1,'output',1,'negative','trimf',[-10 -5 0]);
pc1=addmf(pc1,'output',1,'0','trimf',[-5 0 5]);
pc1=addmf(pc1,'output',1,'positive','trimf',[0 5 10]);
rulelist=[1 2 3 1 1;2 2 1 1 1;0 1 2 1 2];
pc1=addrule(pc1,rulelist);
evalfis([5 2],pc1)
t=5;P=[0.1:0.1:3]';
x=[t*ones(size(P)),P];
y=evalfis(x,pc1);