實(shí)驗(yàn)連續(xù)系統(tǒng)的數(shù)字仿真

1、實(shí)驗(yàn)三連續(xù)系統(tǒng)的數(shù)字仿真【練習(xí)3.1】已知線性定常系統(tǒng)的狀態(tài)空間表達(dá)式為1122()()012()()()560xtxtutxtxt????????????????????????????????12()()12()xtytxt???????且初始狀態(tài)為零，試?yán)盟碾A龍格庫(kù)塔法求系統(tǒng)的單位階躍響應(yīng)。程序：r=1A=[0156]B=[20]C=[12]D=[0]Tf=input(仿真時(shí)間Tf=)h=input(計(jì)算步長(zhǎng)h=)x=zeros

2、(length(A)1)y=0t=0fi=1:TfhK1=AxBrK2=A(xhK12)BrK3=A(xhK22)BrK4=A(xhK3)Brx=xh(K12K22K3K4)6y=[yCx]t=[tt(i)h]endplot(ty)0123456789101.61.41.210.80.60.40.200.2【練習(xí)3.2】設(shè)單位反饋系統(tǒng)的開環(huán)傳遞函數(shù)為，試分別利用二階和四階4()(2)Gsss??龍格庫(kù)塔法求系統(tǒng)的單位階躍響應(yīng)。程序如下：

3、C=[C0]x=[xr]eAt=eye(size(A))AhA^2h^22A^3h^3(32)A^4h^4(432)fi=1:Tfhx=eAtxy=[yCx]t=[tt(i)h]endplot(ty)仿真時(shí)間Tf=5計(jì)算步長(zhǎng)h=0.100.511.522.533.544.55050100150200250300350【練習(xí)3.4】利用input函數(shù)修改課本P103例31所給的程序ex3_1.m，將其中的給定參數(shù)r、numo、deno、n

4、umh、denh利用鍵盤輸入，使其變?yōu)檫B續(xù)控制系統(tǒng)面向傳遞函數(shù)的通用數(shù)字仿真程序。仿真程序：r=input(r=)numo=input(numo=)deno=input(deno=)numh=input(?numh=)denh=input(?denh=)[numden]=feedback(numodenonumhdenh)[AbCd]=tf2ss(numden)Tf=input(Tf=)h=input(h=)x=zeros(length