控制系統(tǒng)計算機仿真作業(yè)-在線瀏覽
2024-08-17 14:28本頁面
【正文】 部分分式： [R,P,H]=residue(num,den)R = + P = + H = []（2） A=[ 5 。 1。2。0]B = 4 2 2 0 C=[0 2 0 2]C = 0 2 0 2 D=0D = 0零極點增益形式： [Z,P,K]=ss2zp(A,B,C,D)Z = + P = + + K = 傳遞函數(shù)形式： num=[0 4 14 22 15]num = 0 4 14 22 15 den=[1 4 ]den = 部分分式： [R，P，H]=residue(num,den)R = + P = + H = []12 用毆拉法matlab編程實現(xiàn)下列系統(tǒng)的輸出響應在上，時的數(shù)值解。t=0::1t = Columns 1 through 9 0 Columns 10 through 11 h=。 t=0::1。 y(1)=1。end plot(t,y,39。) hold on m=exp(1*t)m = Columns 1 through 9 Columns 10 through 11 plot(t,m,39。)13 用四階龍格－庫塔梯形法matlab編程實現(xiàn)12題的數(shù)值解，要求以圖形的方式通過與真值及毆拉法的比較，分析其精度。 y(1)=1。o39。r*39。g39。程序:disp(39。)h=1。for i=0:1:63 y=y+2^i。 運行結果：y= +019y= +01915編寫matlab的M函數(shù)，以實現(xiàn)。function sum=zuoye5(m)sum=0。 for j=1:i t=t*2。end二、控制系統(tǒng)分析21設典型閉環(huán)結構控制系統(tǒng)如下圖所示，當階躍輸入幅值時，(t)的響應。y(t)r(t)_2解：程序1：clear all。a=[ 1]。V=2。X0=[0 0 0 0]。T0=0。h=。a=a/a(1)。A=[rot90(rot90(eye(n1,n)))。B=[zeros(1,n1),1]39。C=[fliplr(b),zeros(1,nm1)]。X=X039。t=T0。for i=1:NK1=Ab*X+B*R。K3=Ab*(X+h*K2/2)+B*R。X=X+h*(K1+K2*2+K3*2+K4)/6。t=[t,t(i)+h]。,y39。b39。反饋系數(shù)為2時的階躍響應曲線39。 hold on。close all。b=[30 25]。b1=b。a1=a1/a1(1)。b=b/a(1)。m=length(b)。R=20。Tf=15。X10=[0 0 0 0 0]。X0=[0 0 0 0]。A1=a1(2:n1+1)。C1=[fliplr(b1),zeros(1,n1m1)]。A1=[rot90(rot90(eye(n11,n1)))。A=[rot90(rot90(eye(n1,n)))。B1=[zeros(1,n11),1]39。Ab=A1B1*C1*V。y1=0。N=round(TfT0)/h。K1=Ab*X1+B1*R。K3=Ab*(X1+h*K2/2)+B1*R。X1=X1+h*(K1+K2*2+K3*2+K4)/6。t=[t,t(i)+h]。y=0。for i=1:NK1=A*X+B*(Ry1(i))。K3=A*(X+h*K2/2)+B*(Ry1(i))。X=X+h*(K1+K2*2+K3*2+K4)/6。t=[t,t(i)+h]。,y39。r39。title(39。)。grid on22 下圖中，若各環(huán)節(jié)傳遞函數(shù)已知為：，， ；試列寫鏈接矩陣W、W0和非零元素陣WIJ，將程序sp4_2完善后，應用此程序求輸出的響應曲線。close all。b=[ 1 ]39。d=[0 0 0 0 0 0 0]39。WIJ=[1 0 1。2 9 1。4 3 1。5 4 1。6 7 。 8 6 1。n=9。Yt0=[0 0 0 0 0 0 0 0 0]。L1=10。Tf=3。A=diag(P(:,1))。C=diag(P(:,3))。m=length(WIJ(:,1))。W=zeros(n,n)。 else W(WIJ(k,1)
法律信息相關推薦
鄂ICP備17016276號-1