%被控系统建模
num=20;
den=[8 6 1];
[a1,b,c,d]=tf2ss(num,den);%Transfer function to state-space conversion.
%系统参数
T=0.01;h=T;
N=500;
R=1.5*ones(1,N);%参考输入
uu=zeros(1,N);
YY=zeros(3,N);
kt=1;
for alpha=[0.25 0.5 0.75];
x=[0;0];
%创建FIS结构simple
a=newfis('simple');
%给FIS添加变量
a=addvar(a,'input','e',[-6,6]); %Parameter e
a=addvar(a,'input','ec',[-6,6]); %Parameter ec
a=addvar(a,'output','u',[-6,6]); %Parameter u
%定义输入和输出变量的隶属函数
a=addmf(a,'input',1,'NB','trapmf',[-6,-6,-5,-3]);
a=addmf(a,'input',1,'NS','trapmf',[-5,-3,-2,0]);
a=addmf(a,'input',1,'Z','trimf',[-2,0,2]);
a=addmf(a,'input',1,'PS','trapmf',[0,2,3,5]);
a=addmf(a,'input',1,'PB','trapmf',[3,5,6,6]);
a=addmf(a,'input',2,'NB','trapmf',[-6,-6,-5,-3]);
a=addmf(a,'input',2,'NS','trapmf',[-5,-3,-2,0]);
a=addmf(a,'input',2,'Z','trimf',[-2,0,2]);
a=addmf(a,'input',2,'PS','trapmf',[0,2,3,5]);
a=addmf(a,'input',2,'PB','trapmf',[3,5,6,6]);
a=addmf(a,'output',1,'NB','trapmf',[-6,-6,-5,-3]);
a=addmf(a,'output',1,'NS','trapmf',[-5,-3,-2,0]);
a=addmf(a,'output',1,'Z','trimf',[-2,0,2]);
a=addmf(a,'output',1,'PS','trapmf',[0,2,3,5]);
a=addmf(a,'output',1,'PB','trapmf',[3,5,6,6]);
%创建模糊规则库
for i=1:5
for j=1:5
rr(i,j)=round(alpha*i+(1-alpha)*j);%通过参数来修改规则
end
end
rr=6-rr;
r1=zeros(prod(size(rr)),3);k=1;
for i=1:size(rr,1)
for j=1:size(rr,2)
r1(k,:)=[i,j,rr(i,j)];
k=k+1;
end
end
[r,s]=size(r1);
r2=ones(r,2);
rulelist=[r1,r2];
a=addrule(a,rulelist);
%采用模糊控制器的系统仿真
e=0;de=0;
ke=30;kd=0.2;ku=1;
for k=1:N
%输入变量转换至邻域
e1=ke*e;
de1=kd*de;
if e1>=6
e1=6;
elseif e1<=-6
e1=-6;
end
if de1>=6
de1=6;
elseif de1<=-6
de1=-6;
end
%模糊推理,计算出被控对象的控制输入
in=[e1 de1];
u=ku*evalfis(in,a);
uu(1,k)=u;
%uu(1,k)=ku * evalfis([e1,de1],a);
%控制作用于被控系统,计算系统输出
k0=a1*x+b*u;
k1=a1*(x+h*k0/2)+b*u;
k2=a1*(x+h*k1/2)+b*u;
k3=a1*(x+h*k2)+b*u;
x=x+(k0+2*k1+2*k2+k3)*h/6;
y=c*x+d*u;
yy(kt,k)=y;
%计算系统输出误差及误差导数
e1=e;e=y-R(1,k);
de=(e-e1)/T;
end
kt=kt+1;
end
%控制输出曲线
kk=[1:N]*T;
figure(1);
plot(kk,R,'k',kk,yy(1,:),'k--',kk,yy(2,:),'r',kk,yy(3,:),'b')
xlabel('时间');ylabel('输出');
gtext('alpha=0.25');
gtext('alpha=0.5');
gtext('alpha=0.75');
评论0