控制系統(tǒng)計算機(jī)仿真作業(yè)-閱讀頁
2025-07-22 14:28本頁面
【正文】 ,WIJ(k,2))=WIJ(k,3)。R=C*WA。Ab=Q\R。Y=Yt039。t=T0。for i=1:N for j=1:L1 K1=Ab*Y+b1*Y0。 K3=Ab*(Y+h*K2/2)+b1*Y0。 Y=Y+h*(K1+K2*2+K3*2+K4)/6。 t=[t,t(i)+h*L1]。,y39。r39。title(39。)。解：程序clear all。a=[1 0 1 0 1 1 0 1 1]39。c=[1 1 1 1 70 130 ]39。P=[a,b,c,d]。2 1 1。3 2 1。 4 8 1。6 5 1。7 6 1。9 7 1]。Y0=1。L1=10。Tf=3。A=P(:,1)。C=P(:,3)。m=length(WIJ(:,1))。W=zeros(n,n)。 else W(WIJ(k,1),WIJ(k,2))=WIJ(k,3)。 FI(i)=1。 FIJ(i)=h*h*C(i)/B(i)/2。FID(i)=0。 else end else FI(i)=exp(h*A(i)/B(i))。 FIJ(i)=h*C(i)/A(i)FIM(i)*B(i)/A(i)。FID(i)=0。 FID(i)=D(i)/B(i)。X=Y。Uk=zeros(n,1)。t=T0:h*L1:Tf。for k=1:N1 for l=1:L1 Ub=Uk。 Udot=(UkUb)/h。 X=FI39。.*Uk+FIJ39。 Y=FIC39。.*Uf。end[t39。]plot(t,y)24求下圖非線性系統(tǒng)的輸出響應(yīng)y(t)，并與無非線性環(huán)節(jié)情況進(jìn)行比較。0 1 20 0。10 1 1 0]。1 4 1。3 2 1。 Z=[0 0 0 0]。 h=。 n=4。 Tf=20。 Y0=10。 plot(t,y,39。) hold on Z=[4 0 0 0]。 sp4_4。g39。B=diag(P(:,2))。D=diag(P(:,4))。 W0=zeros(n,1)。 for k=1:m if(WIJ(k,2)==0)。 else W(WIJ(k,1),WIJ(k,2))=WIJ(k,3)。 end。 FI(i,i)=1。 FIJ(i,i)=h*h*C(i,i)/B(i,i)/2。FID(i,i)=0。 FID(i,i)=D(i,i)/B(i,i)。 FIM(i,i)=(1FI(i,i))*C(i,i)/A(i,i)。 FIC(i,i)=1。 if(D(i,i)~=0)。 FID(i,i)=D(i,i)/B(i,i)。X=Y。Uk=zeros(n,1)。 t=T0:h*L1:Tf。 for k=1:N1 for i=1:L1 Ub=Uk。 for i=1:n if(Z(i)~=0) if(Z(i)==1) Uk(i,i)=satu(Uk(i,i),S(i))。 end if(Z(i)==3) [Uk(i,i),Ubb(i,i)]=backlash(Ubb(i,i),Uk(i,i),Ub(i,i),s(i))。 Uf=2*UkUb。*X+FIM39。*Udot。 Y=FIC39。*Uf。n if(Z(i)~=0) if(Z(i)==4) Y(i,i)=satu(Y(i,i),S(i))。 end if(Z(i)==6) [Y(i,i),Ubb(i,i)]=backlash(Ubb(i,i),Y(i,i),Yb(i,i),s(i))。 end backlash函數(shù): function[Uc,Ubb]=backlash(Urb,Ur,Ucb,S1) if(UrUrb) if((UrS1)=Ucb) Uc=UrS1。 end else if(UrUrb) if((Ur+S1)=Ucb) Uc=Ur+S1。 end else Uc=Ucb。Satu函數(shù): function Uc=satu(Ur,S1) if(abs(Ur)=S1) if(Ur0) Uc=S1。 end else Uc=Ur。 else Uc=Ur+S1。 end25 已知系統(tǒng)，當(dāng)分別取0、1時，判別系統(tǒng)的能控性與能觀測性。A=[2 1 1。1 0 1]B=[1 1
法律信息相關(guān)推薦
鄂ICP備17016276號-1