变论域模糊控制的S函数代码出错了
我做了一个变论域的模糊控制，可以仿真，但是模糊控制输出为0，不知道哪里出错了，请高人指点，谢谢！代码如下：（另外附件里是simulink仿真模型）
function=sfun_mohu(t,x,u,flag)
ke=3.00;
kc=3/60.00;
ku=30;
switch flag,
case 0
=mdlInitializeSizes(ke,kc,ku);
case 3
sys=mdlOutputs(t,x,u,ke,kc,ku);
case {1,2,4,9}
sys=[];
otherwise
error(['unhandled flag=',num2str(flag)]);
end
function=mdlInitializeSizes(ke,kc,ku)
sizes=simsizes; %read default control varibles
sizes.NumContStates=0; %no continuous states
sizes.NumDiscStates=0; %3 states and assume they are the P I D components
sizes.NumOutputs=1; %1 output
sizes.NumInputs=3; %3 inputs
sizes.DirFeedthrough=1; %input reflected directly in output
sizes.NumSampleTimes=1; %singal sampling period
sys=simsizes(sizes);
x0=[];
str=[];
ts=[-1 0];
function sys=mdlOutputs(t,x,u,ke,kc,ku)
m=1-0.99*exp(-0.5*u(1)^2);
n=1-0.99*exp(-0.5*u(1)^2);
p=1-0.99*exp(-0.5*u(3)^2);
a=newfis('vertical');
a=addvar(a,'input','e',[-4*m 4*m]);
a=addmf(a,'input',1,'NB','trimf',[-6,-4,-2]*m);
a=addmf(a,'input',1,'NS','trimf',[-4,-2,0]*m);
a=addmf(a,'input',1,'ZR','trimf',[-2,0,2]*m);
a=addmf(a,'input',1,'PS','trimf',*m);
a=addmf(a,'input',1,'PB','trimf',*m);
%
a=addvar(a,'input','de',[-4*n 4*n]);
a=addmf(a,'input',2,'NB','trimf',[-6,-4,-2]*n);
a=addmf(a,'input',2,'NS','trimf',[-4,-2,0]*n);
a=addmf(a,'input',2,'ZR','trimf',[-2,0,2]*n);
a=addmf(a,'input',2,'PS','trimf',*n);
a=addmf(a,'input',2,'PB','trimf',*n);
%
a=addvar(a,'output','u',[-4*p 4*p]);
a=addmf(a,'output',1,'NB','trimf',[-6,-4,-2]*p);
a=addmf(a,'output',1,'NS','trimf',[-4,-2,0]*p);
a=addmf(a,'output',1,'ZR','trimf',[-2,0,2]*p);
a=addmf(a,'output',1,'PS','trimf',*p);
a=addmf(a,'output',1,'PB','trimf',*p);
%
rr=[5 5 4 4 3
5 4 4 3 2
5 4 3 2 1
4 3 2 1 1
3 2 2 1 1];
r1=zeros(prod(size(rr)),3);k=1;
for i=1:size(rr,1)
for j=1:size(rr,2)
r1(k,:)=;
k=k+1;
end
end
=size(r1);
r2=ones(r,2);
ruleList=;
a=addrule(a,ruleList);
e1=ke*u(1);
de1=kc*u(2);
if e1>4
e1=4;%限制误差在基本论域
end
if e1<-4
e1=-4;
end
if de1>4
de1=4; %限制误差变化率
end
if de1<4
de1=-4;
end
in=;
sys=ku*evalfis(in,a); 不知道楼主的问题解决了没有倒数第二行是in=【】 还是rr[]
页:
[1]