语言变量x,y,z的隶属度函数.

《语言变量x,y,z的隶属度函数.》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、语言变量X，Y，Z的隶属度函数.论域语言值-4-3-2-101234NB10.50000000NS00.510.500000ZO0000.510.5000PS000000.510.50PB00000000.51模糊控制规则表ZXNBNSZOPSPBYNBPBPBPSPSZONSPBPSPSZOZOZOPSPSZOZONSPSPSZOZONSNSPBZOZONSNSNB设计带有纯延迟的一阶惯性环节（假设T=6,=0.02）：G(s)=的模糊控制器，观察仿真结果。编程如下：%被控系统建模num=1;den=[6,1];[a1,b,c,d]=tf2ss(num,den);%传递函数转换到状态空间x=

2、[0];%系统参数T=0.01;h=T;td=0.02;N=1000;nd=td/T;%系统纯延迟R=ones(1,N);%参考输入%定义输入和输出变量及隶属度函数a=newfis('Simple');a=addvar(a,'input','e',[-44]);a=addmf(a,'input',1,'NB','trimf',[-4,-4,-2]);a=addmf(a,'input',1,'NS','trimf',[-4,-2,0]);a=addmf(a,'input',1,'ZO','trimf',[-2,0,2]);a=addmf(a,'input',1,'PS','trimf',[0,2

3、,4]);a=addmf(a,'input',1,'PB','trimf',[2,4,4]);a=addvar(a,'input','de',[-44]);a=addmf(a,'input',2,'NB','trimf',[-4,-4,-2]);a=addmf(a,'input',2,'NS','trimf',[-4,-2,0]);a=addmf(a,'input',2,'ZO','trimf',[-2,0,2]);a=addmf(a,'input',2,'PS','trimf',[0,2,4]);a=addmf(a,'input',2,'PB','trimf',[2,4,4]);a=addv

4、ar(a,'output','u',[-44]);a=addmf(a,'output',1,'NB','trimf',[-4,-4,-2]);a=addmf(a,'output',1,'NS','trimf',[-4,-2,0]);a=addmf(a,'output',1,'ZO','trimf',[-2,0,2]);a=addmf(a,'output',1,'PS','trimf',[0,2,4]);a=addmf(a,'output',1,'PB','trimf',[2,4,4]);%模糊规则矩阵rr=[5544354433443324332233221];r1=zeros(prod(si

5、ze(rr)),3);k=1;fori=1:size(rr,1)forj=1:size(rr,2)r1(k,:)=[i,j,rr(i,j)];k=k+1;endend[r,s]=size(r1);r2=ones(r,2);rulelist=[r1,r2];a=addrule(a,rulelist);%采用模糊控制器的二阶系统仿真e=0;de=0;ke=30;kd=5;ku=1;fork=1:N%输入变量变换至论域e1=ke*e;de1=kd*de;ife1>=4e1=4;elseife1<=-4e1=-4;endifde1>=4de1=4;elseifde1<=-4de1=-4;end%模糊推

6、理，计算出被控对象的控制输入in=[e1de1];u=ku*evalfis(in,a);uu(1,k)=u;ifk<=ndu=0;elseu=uu(1,k-nd);end%控制作用于被控系统，计算系统输出；利用龙格—库塔法进行系统仿真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(1,k)=y;%计算系统输出误差及误差导数e1=e;e=y-R(1,k);de=(e-e1)/T;end%典型二阶环节的模糊控制输出曲线kk=

7、[1:N]*T;figure(1);plot(kk,R,'k',kk,yy,'r');gridonxlabel('时间（秒）');ylabel('输出');