控制系統(tǒng)計算機仿真作業(yè)(存儲版)
【正文】 k1 = k2 = k3 = k4 = y = Columns 1 through 9 Column 10 k1 = k2 = k3 = k4 = y = Columns 1 through 9 Columns 10 through 11 plot(t,y,39。y=39。format longfor i=0:m t=1。b=[30 25]。b=b/a(1)。Ab=AB*C*V。K4=Ab*(X+h*K3)+B*R。) title(39。a1=conv(a,[ 1])。V=1。n=4。fliplr(A)]。for i=1:N。endX=X039。y=[y,C*X]。反饋通道為慣性環(huán)節(jié)時的階躍響應(yīng)曲線39。c=[1 1 1 1 70 130 ]39。 4 8 1。Y0=1。B=diag(P(:,2))。 endendQ=BD*W。N=round((TfT0)/(h*L1))。end[t39。grid on，并與四階龍格－庫塔法比較精度，同時分析兩種方法對仿真步長選取的區(qū)別。WIJ=[1 0 1。6 7 。T0=0。W0=zeros(n,1)。 FIC(i)=1。 if(D(i)~=0) FIM(i)=(1FI(i))*D(i)/A(i) FIJ(i)=h*D(i)/A(i)FIM(i)*B(i)/A(i) FIC(i)=C(i)/D(i)A(i)/B(i)。N=length(t)。.*Udot。y(t)e(t)r(t)=1055P=[ 1 1。4 3 1]。 nout=4。 plot(t,y,39。W=zeros()。 FIM(i,i)=h*C(i,i)/B(i,i)。 FIJ(i,i)=h*C(i,i)/A(i,i)FIM(i,i)*B(i,i)/A(i,i)。y=0。 end if(Z(i)==2) Uk(i,i)=dead(Uk(i,i),S(i))。 Yb=Y。 end end end end y=[y,Y(nout)]。 else Uc=S1。1 0 1。a=1時num=[1]den=[10,27,18][A,B,C,D]=tf2ss(num,den)P=ctrb(A,B)N=rank(P)Q=obsv(A,C)N=rank(Q)A = 0B = 1 0C = 0 D = 0P = 0 N = 2Q = 0 0N = 2a=0時A = 0B = 1 0C = 0 0D = 0P = 0 N = 2Q = 0 0 0 0N = 0a=1時A = 0B = 1 0C = 0 D = 0P = 0 N = 2Q = 0 0N = 226 已知系統(tǒng)的狀態(tài)方程為，采用狀態(tài)反饋，將系統(tǒng)的幾點配置到3，求狀態(tài)反饋陣K。 end end Ubb=Ur。 end if(Z(i)==5) Y(i,i)=dead(Y(i,i),S(i))。*Uk+FIJ39。 Uk=W*Y+W0*Y0。 else end end end Y=zeros(n,1)。 else end else FI(i,i)=exp(h*A(i,i)/B(i,i))。 for i=1:n
法律信息相關(guān)推薦
鄂ICP備17016276號-1