【正文】 end y=[y,Y(nout)]。 Uf=2*UkUb。y=0。 FIM(i)=(1FI(i))*C(i)/A(i)。 endend for i=1:n if(A(i)==0)。B=P(:,2)。n=9。4 3 1。b=[ 1 ]39。)。 K4=Ab*(Y+h*K3)+b1*Y0。b1=Q\V1。W0=zeros(n,1)。T0=0。7 6 1。2 1 1。G2(s)G5(s)G10(s)G4(s)y0_y7_G1(s)G3(s)G6(s)G7(s)G8(s)G9(s)解：程序clear all。]plot(t,y,39。K2=A*(X+h*K1/2)+B*(Ry1(i))。K4=Ab*(X1+h*K3)+B1*R。X1=X1039。C=[fliplr(b),zeros(1,nm)]。h=。m1=length(b1)。 grid on程序2：clear all。end [t39。N=round(TfT0)/h。fliplr(A)]。n=4。進(jìn)一步考慮：當(dāng)反饋通道為時，以實現(xiàn)此時所對應(yīng)的y(t)。 end disp(y)。) lea=ym plot(t,lea,39。bo39。y(1)=1。 。自動化系《控制系統(tǒng)計算機仿真》課程組2017年3月《控制系統(tǒng)計算機仿真》上機報告一、Matlab基礎(chǔ)操作11用MATLAB語言求下列系統(tǒng)的狀態(tài)方程、傳遞函數(shù)、零極點增益和部分分式形式的模型參數(shù)，并分別寫出其相應(yīng)的數(shù)學(xué)模型表達(dá)式：（1）程序如下: num=[7,24,24]den=[10,35,50,24] [A,B,C,D]=tf2ss(num,den)系統(tǒng)的狀態(tài)方程：A = 0 0 0 0B = 1 0 0C = D = 0零極點增益形式：[Z,P,K]=tf2zp(num,den)Z = + P = + K = 部分分式： [R,P,H]=residue(num,den)R = + P = + H = []（2） A=[ 5 。t=0::1t = Columns 1 through 9 0 Columns 10 through 11 h=。) hold on m=exp(1*t)m = Columns 1 through 9 Columns 10 through 11 plot(t,m,39。r*39。for i=0:1:63 y=y+2^i。end二、控制系統(tǒng)分析21設(shè)典型閉環(huán)結(jié)構(gòu)控制系統(tǒng)如下圖所示，當(dāng)階躍輸入幅值時，(t)的響應(yīng)。X0=[0 0 0 0]。A=[rot90(rot90(eye(n1,n)))。t=T0。t=[t,t(i)+h]。 hold on。a1=a1/a1(1)。Tf=15。C1=[fliplr(b1),zeros(1,n1m1)]。Ab=A1B1*C1*V。K3=Ab*(X1+h*K2/2)+B1*R。for i=1:NK1=A*X+B*(Ry1(i))。,y39。grid on22 下圖中，若各環(huán)節(jié)傳遞函數(shù)已知為：，， ；試列寫鏈接矩陣W、W0和非零元素陣WIJ，將程序sp4_2完善后，應(yīng)用此程序求輸出的響應(yīng)曲線。WIJ=[1 0 1。6 7 。L1=10。m=length(WIJ(:,1))。Ab=Q\R。 K3=Ab*(Y+h*K2/2)+b1*Y0。r39。a=[1 0 1 0 1 1 0 1 1]39。3 2 1。9 7 1]。A=P(:,1)。 else W(WIJ(k,1),WIJ(k,2))=WIJ(k,3)。 else end else FI(i)=exp(h*A(i)/B(i))。X=Y。 Udot=(UkUb)/h。.*Uf。10 1 1 0]。 h=。 plot(t,y,39。B=diag(P(:,2))。 else W(WIJ(k,1),WIJ(k,2))=WIJ(k,3)。FID(i,i)=0。 if(D(i,i)~=0)。 t=T0:h*L1:Tf。 Uf=2*UkUb。*Uf。 end else if(UrUrb) if((Ur+S1)=Ucb) Uc=Ur+S1。 else Uc=U