控制系統(tǒng)計(jì)算機(jī)仿真作業(yè)(編輯修改稿)
2024-08-03 14:28 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 all。close all。a=[ 1]。b=[30 25]。a1=conv(a,[ 1])。b1=b。b1=b1/a1(1)。a1=a1/a1(1)。m1=length(b1)。b=b/a(1)。a=a/a(1)。m=length(b)。V=1。R=20。T0=0。Tf=15。h=。X10=[0 0 0 0 0]。n1=5。X0=[0 0 0 0]。n=4。A1=a1(2:n1+1)。A=a(2:n+1)。C1=[fliplr(b1),zeros(1,n1m1)]。C=[fliplr(b),zeros(1,nm)]。A1=[rot90(rot90(eye(n11,n1)))。fliplr(A1)]。A=[rot90(rot90(eye(n1,n)))。fliplr(A)]。B1=[zeros(1,n11),1]39。B=[zeros(1,n1),1]39。Ab=A1B1*C1*V。X1=X1039。y1=0。t=T0。N=round(TfT0)/h。for i=1:N。K1=Ab*X1+B1*R。K2=Ab*(X1+h*K1/2)+B1*R。K3=Ab*(X1+h*K2/2)+B1*R。K4=Ab*(X1+h*K3)+B1*R。X1=X1+h*(K1+K2*2+K3*2+K4)/6。y1=[y1,C1*X1]。t=[t,t(i)+h]。endX=X039。y=0。t=T0。for i=1:NK1=A*X+B*(Ry1(i))。K2=A*(X+h*K1/2)+B*(Ry1(i))。K3=A*(X+h*K2/2)+B*(Ry1(i))。K4=A*(X+h*K3)+B*(Ry1(i))。X=X+h*(K1+K2*2+K3*2+K4)/6。y=[y,C*X]。t=[t,t(i)+h]。end[t39。,y39。]plot(t,y,39。r39。)。title(39。反饋通道為慣性環(huán)節(jié)時(shí)的階躍響應(yīng)曲線39。)。hold on。grid on22 下圖中，若各環(huán)節(jié)傳遞函數(shù)已知為：，，，，， ；試列寫(xiě)鏈接矩陣W、W0和非零元素陣WIJ，將程序sp4_2完善后，應(yīng)用此程序求輸出的響應(yīng)曲線。G2(s)G5(s)G10(s)G4(s)y0_y7_G1(s)G3(s)G6(s)G7(s)G8(s)G9(s)解：程序clear all。close all。a=[1 0 1 0 1 1 0 1 1]39。b=[ 1 ]39。c=[1 1 1 1 70 130 ]39。d=[0 0 0 0 0 0 0]39。P=[a,b,c,d]。WIJ=[1 0 1。2 1 1。2 9 1。3 2 1。4 3 1。 4 8 1。5 4 1。6 5 1。6 7 。7 6 1。 8 6 1。9 7 1]。n=9。Y0=1。Yt0=[0 0 0 0 0 0 0 0 0]。h=。L1=10。T0=0。Tf=3。nout=7。A=diag(P(:,1))。B=diag(P(:,2))。C=diag(P(:,3))。D=diag(P(:,4))。m=length(WIJ(:,1))。W0=zeros(n,1)。W=zeros(n,n)。for k=1:m if(WIJ(k,2)==0) W0(WIJ(k,1))=WIJ(k,3)。 else W(WIJ(k,1),WIJ(k,2))=WIJ(k,3)。 endendQ=BD*W。R=C*WA。V1=C*W0。Ab=Q\R。b1=Q\V1。Y=Yt039。y=Y(nout)。t=T0。N=round((TfT0)/(h*L1))。for i=1:N for j=1:L1 K1=Ab*Y+b1*Y0。 K2=Ab*(Y+h*K1/2)+b1*Y0。 K3=Ab*(Y+h*K2/2)+b1*Y0。 K4=Ab*(Y+h*K3)+b1*Y0。 Y=Y+h*(K1+K2*2+K3*2+K4)/6。 end y=[y,Y(nout)]。 t=[t,t(i)+h*L1]。end[t39。,y39。]plot(t,y,39。r39。)。title(39。單位階躍響應(yīng)曲線39。)。grid on，并與四階龍格－庫(kù)塔法比較精度，同時(shí)分析兩種方法對(duì)仿真步長(zhǎng)選取的區(qū)別。解：程序clear all。close all。a=[1 0 1 0 1 1 0 1 1]39。b=[ 1 ]39。c=[1 1 1 1 70 130 ]39。d=[0 0 0 0 0 0 0]39。P=[a,b,c,d]。WIJ=[1 0 1。2 1 1。2 9 1。3 2 1。4 3 1。 4 8 1。5 4 1。6 5 1。6 7 。7 6 1。 8 6 1。9 7 1]。n=9。Y0=1。h=。L1=10。T0=0。Tf=3。nout=7。A=P(:,1)。B=P(:,2)。C=P(:,3)。D=P(:,4)。m=length(WIJ(:,1))。W0=zeros(n,1)。W=zeros(n,n)。for k=1:m if(WIJ(k,2)==0) W0(WIJ(k,1))=WIJ(k,
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1