matlab與模煳控制系統(tǒng)-資料下載頁
2025-05-10 12:22本頁面
【正文】 fuzdemos MATLAB應(yīng)用技術(shù) 46 6． 2． 3 MATLAB模糊邏輯工具箱命令函數(shù)應(yīng)用實例 設(shè)計典型二階環(huán)節(jié)： 20)( 2 ??? sssH的模糊控制器，使系統(tǒng)輸出盡快跟隨系統(tǒng)輸入。 MATLAB應(yīng)用技術(shù) 47 隸屬度函數(shù) 誤差隸屬度函數(shù) 誤差變化率隸屬度函數(shù) MATLAB應(yīng)用技術(shù) 48 輸出隸屬度函數(shù) MATLAB應(yīng)用技術(shù) 49 FC的模糊推量規(guī)則表 u e de NB NS ZR PS PB NB PB PB PS PS ZR NS PB PS PS ZR ZR ZR PS PS ZR ZR NS PS PS ZR ZR NS NS PB ZR ZR NS NS NB MATLAB應(yīng)用技術(shù) 50 MATLAB程序如下 num=20。 den=[ 1]。 [a1 b c d]=tf2ss(num,den)。 x=[0。0]。 T=。h=T。 N= R=*ones(1,N)。參考輸入 %定義輸入輸出變量與隸屬度函數(shù) a=newfis(39。Simple39。)。 a=addvar(a,39。input39。,39。e39。,[6 6])。 a=addmf(a,39。input39。,1,39。NB39。,39。trapmf39。,[6,6,5,3])。 a=addmf(a,39。input39。,1,39。NS39。,39。trapmf39。,[5,3,2,0])。 a=addmf(a,39。input39。,1,39。ZR39。,39。trimf39。,[2,0,2])。 a=addmf(a,39。input39。,1,39。PS39。,39。trapmf39。,[0,2,3,5])。 a=addmf(a,39。input39。,1,39。PB39。,39。trapmf39。,[3,5,6,6])。 a=addvar(a,39。input39。,39。de39。,[6,6])。 a=addmf(a,39。input39。,2,39。NB39。,39。trapmf39。,[6,6,5,3])。 a=addmf(a,39。input39。,2,39。NS39。,39。trapmf39。,[5,3,2,0])。 a=addmf(a,39。input39。,2,39。ZR39。,39。trimf39。,[2,0,2])。 a=addmf(a,39。input39。,2,39。PS39。,39。trapmf39。,[0,2,3,5])。 a=addmf(a,39。input39。,2,39。PB39。,39。trapmf39。,[3,5,6,6])。 a=addvar(a,39。output39。,39。u39。,[3,3])。 a=addmf(a,39。output39。,1,39。NB39。,39。trapmf39。,[3,3,2,1])。 a=addmf(a,39。output39。,1,39。NS39。,39。trimf39。,[2,1,0])。 a=addmf(a,39。output39。,1,39。ZR39。,39。trimf39。,[1,0,1])。 a=addmf(a,39。output39。,1,39。PS39。,39。trimf39。,[0,1,2])。 a=addmf(a,39。output39。,1,39。PB39。,39。trapmf39。,[1,2,3,3])。 %模糊規(guī)則矩陣 MATLAB應(yīng)用技術(shù) 51 %模糊規(guī)則矩陣 rr=[5 5 4 4 3 5 4 4 3 3 4 4 3 3 2 4 3 3 2 2 3 3 2 2 1]。 rr=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)。 MATLAB應(yīng)用技術(shù) 52 %模糊控制系統(tǒng)仿真 e=0。 de=0。 ke=30。 kd=20。 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 MATLAB應(yīng)用技術(shù) 53 %模糊推理，計算被控對象輸入 in=[e1 de1]。 u=ku*evalfis(in,a)。 uu(1,k)=u。 %控制作用于被控系統(tǒng)，計算系統(tǒng)輸出 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。 %計算系統(tǒng)輸出誤差及誤差變化率 e1=e。 e=yR(1,k)。 de=(ee1)/T。 end %模糊控制輸出曲線 kk=[1:N]*T。 figure(1)。 plot(kk,R,39。k39。,kk,yy,39。r39。)。 grid on MATLAB應(yīng)用技術(shù) 54 系統(tǒng)階躍響應(yīng)曲線
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1