MATLAB模糊控制仿真程序

1、個人收集整理僅供參考學(xué)習(xí)clear all;close all;ts=20;sys1=tf(1,60,1,'inputdelay',80);%生成或轉(zhuǎn)換傳遞函數(shù)模型dsys1=c2d(sys1,ts,'zoh');%將連續(xù)系統(tǒng)離散化num1,den1=tfdata(dsys1,'v');%傳統(tǒng)PID算法u_1=0;u_2=0;u_3=0;u_4=0;u_5=0;y_1=0;y_2=0;y_3=0;error_1=0;ei=0;for k=1:1:300time(k)=k*ts;rin(k)=40;yout(k)=-den1(2)*y_1+num1

2、(2)*u_5;error(k)=rin(k)-yout(k);ei=ei+error(k)*ts;kp=0.8;ki=0.005;kd=3.0;u(k)=kp*error(k)+kd*(error(k)-error_1)/ts+ki*ei;if u(k)>=110u(k)=110;endif u(k)<=-110u(k)=-110;end%返回PID參數(shù)u 5=u 4;u 4=u 3;u 3=u 2;u 2=u 1;u 1=u(k);y_3=y_2;y_2=y_1;y_1=yout(k);error_1=error(k);end%獲得傳遞函數(shù)模型數(shù)據(jù)%模*® PID控制

3、系統(tǒng)a=newfis('fuzzf);%創(chuàng)建新的模糊推理系統(tǒng)f1=1;a=addvar(a,'input','e',巧*f1,5*f1);%添加e的模糊語言變量a=addmf(a,'input',1,'NB','trimf, -7.5 -5 -2.5);%添加e的模糊語言變量的隸屬度函數(shù)(z型)a=addmf(a,'input',1,'NS','trimf, -5 -2.5 0);a=addmf(a,'input',1,'ZE','tr

4、imf,-2.5 0 2.5);a=addmf(a,'input',1,'PS','trimf,0 2.5 5);a=addmf(a,'input',1,'PB','smf,2.5 5 7.5);f2=1;a=addvar(a,'input','ec',-10*f2,10*f2);%添加ec的模糊語言變量a=addmf(a,'input'2'NB','trimf,-15 -10 -5);a=addmf(a,'input'2

5、9;NS','trimf,-10 -5 0);a=addmf(a,'input'2'Z','trimf, -5 -5.551e-017 5);a=addmf(a,'input'2'PS','trimf, 0 5 10);a=addmf(a,'input'2'PB','smf, 5 10 15);f3=1.5;a=addvar(a,'output','u',0*f3,220*f3);%添加u的模糊語言變量a=addmf(a,

6、9;output',1,'NB','trimf, -55 -4.441e-016 55);a=addmf(a,'output',1,'NS','trimf, 0 55 110);a=addmf(a,'output',1,'Z','trimf, 55 110 165);a=addmf(a,'output',1,'PS','trimf, 110 165 220);a=addmf(a,'output',1,'PB',&

7、#39;smf,165 220 275);rulelist=1 1 5 1 1;%編輯模糊規(guī)則1 2 5 1 1;1 3 5 1 1;1 4 5 1 1;1 5 5 1 1;2 1 5 1 1;2 2 4 1 1;2 3 4 1 1;2 4 3 1 1;2 5 2 1 1;3 1 5 1 1;3 2 4 1 1;3 3 3 1 1;3 4 2 1 1;3 5 2 1 1;4 1 4 1 1;4 2 3 1 1;4 3 2 1 1;4 4 2 1 1;4 5 1 1 1;5 1 1 1 1;5 2 1 1 1;5 3 1 1 1;5 4 1 1 1;5 5 1 1 1;a=addrule(a,r

8、ulelist);%添加模糊規(guī)則函數(shù)showrule(a)%顯示模糊規(guī)則函數(shù)a1=setfis(a,'DefuzzMethod','centroid');% 設(shè)置模糊系統(tǒng)特性writefis(a1,'fuzzf);% 保存模糊系統(tǒng)a2=readfis('fuzzf);%從磁盤讀出保存的模糊系統(tǒng)disp('fuzzy Controller table:e=-5,+5,ec=-10,+10');%顯示矩陣和數(shù)組內(nèi)容Ulist=zeros(5,5);% 全零矩陣for i=1:5for j=1:5e(i)=-4+i;ec(j)=-4+j

9、;%完成模糊推理計算Ulist(i,j)=evalfis(e(i),ec(j),a2);endendUlist=ceil(Ulist)sys1=tf(1,60,1,'inputdelay',80);dsys=c2d(sys1,ts,'zoh');num,den=tfdata(dsys,'v');U_1=0;U_2=0;U_3=0;U_4=0;U_5=0;Y_1=0;Y_2=0;Y_3=0;e_1=0.0;ec_1=0.0;Error_1=0;Ei=0for k=1:1:300time(k尸k*ts;rin(k)=40;U(k)=evalfis(e

10、_1,ec_1,a2);%完成模糊推理計算Yout(k尸-den(2)*Y_1+num(2)*U_5;Error(k尸rin(k)-Yout(k);ki=0.003; Ei=Ei+Error(k)*ts;U(k尸U(k)+ki*Ei;%返回參數(shù)U 5=U 4;U 4=U 3;U 3=U 2;U 2=U 1;U 1=U(k);Y_3=Y_2;Y_2=Y_1;Y_1=Yout(k);e_1=Error(k); ec_1=Error(k)-Error_1;Error_1=Error(k);end%朝正無窮方向取整figure(1); plotfis(a2);figure(2);plotmf(a,&#

11、39;input',1);figure(3);plotmf(a,'input',2);figure(4);plotmf(a,'output',1);figure;plot(time,rin,'b',time,yout,'r');xlabel('time(s)');ylabel('rin,yout');figure(2);plot(time,u,'r');xlabel('time(s)');ylabel('u');figure(3);plot(time,rin,'b',time,Yout,'r');xlabel('time(s)');ylabel('rin,yout');figure(4);plot(time,U,'r');xlab