北京交通大學(xué)現(xiàn)代控制理論matlab實驗報告(編輯修改稿)
2025-07-16 15:29 本頁面
【文章內(nèi)容簡介】 m(j)=1。 endend %計算各極點(diǎn)的重數(shù)Aj=diag(P)。for i=1:length(P)1 if Aj(i,i)==Aj(i+1,i+1) Aj(i,i+1)=1。 else Aj(i,i+1)=0。 endend %構(gòu)造系統(tǒng)矩陣AjB1=0。l=0。for j=1:length(m)l=l+m(j)。B1(l)=1。end Bj=B139。 %構(gòu)造輸入矩陣Bjn=1。l=m(1)。Cj(:,1:m(1))=rot90(R(1:m(1),:),3)。for k=2:length(m)n=l+1。l=l+m(k)。Cj(:,n:l)=rot90(R(n:l,:),3)。end %構(gòu)造輸出矩陣Cj if K==[ ] Dj=0。else Dj=K。 end %構(gòu)造直聯(lián)矩陣Dj Gj=ss(Aj,Bj,Cj,Dj)。（2）能控標(biāo)準(zhǔn)型函數(shù)ctrlts( )可將傳遞函數(shù)轉(zhuǎn)換為能控標(biāo)準(zhǔn)型。該函數(shù)的調(diào)用格式為： 其中num和den分別為傳遞函數(shù)的分子和分母多項式系數(shù)的行向量，為轉(zhuǎn)換得到的能控標(biāo)準(zhǔn)型。該函數(shù)的程序如下:function Gc=ctrlts(num,den) %將傳遞函數(shù)轉(zhuǎn)換為能控標(biāo)準(zhǔn)型m=length(num)1。n=length(den)1。if m==n[R,P,K]=residue(num,den)。num1=numK*den。A(n,:)=1*rot90(den(:,2:n+1),2)。A(1:n1,2:n)=eye(n1)。A(1:n1,1)=zeros(n1,1)。B=[zeros(n1,1)。1]。C=rot90(num1(:,2:n+1),2)。D=K。 else A(n,:)=1*rot90(den(:,2:n+1),2)。A(1:n1,2:n)=eye(n1)。A(1:n1,1)=zeros(n1,1)。B=[zeros(n1,1)。1]。C(:,1:m+1)=rot90(num,2)。C(:,m+2:n)=zeros(1,nm1)。 D=0。endGc=ss(A,B,C,D)。三 實驗步驟： 根據(jù)所給系統(tǒng)的已知條件，如傳遞函數(shù)、零極點(diǎn)模型或（A、B、C、D），實現(xiàn)狀態(tài)空間模型、傳遞函數(shù)模型、零極點(diǎn)增益模型之間的轉(zhuǎn)換，采用MATLAB的相關(guān)函數(shù)編寫m文件。 應(yīng)用系統(tǒng)建模工具，并聯(lián)、串聯(lián)、閉環(huán)、反饋等函數(shù)解
規(guī)章制度相關(guān)推薦
鄂ICP備17016276號-1