3第3章-控制系統(tǒng)數(shù)學(xué)模型及其轉(zhuǎn)換

3第3章-控制系統(tǒng)(kònɡzhìxìtǒnɡ)數(shù)學(xué)模型及其轉(zhuǎn)換第一頁，共79頁。仿真是對數(shù)學(xué)模型進行的試驗，建模是仿真的根底，系統(tǒng)模型化技術(shù)是系統(tǒng)仿真的核心(héxīn)。本章首先給出控制系統(tǒng)數(shù)學(xué)模型的分類，介紹控制系統(tǒng)常用的數(shù)學(xué)模型的描述形式，最后給出各種數(shù)學(xué)模型，如微分方程、傳遞函數(shù)、狀態(tài)空間等之間轉(zhuǎn)換和MATLAB實現(xiàn)。第二頁，共79頁。系統(tǒng)類型根據(jù)系統(tǒng)變量是時間連續(xù)函數(shù)還是時間離散函數(shù)，系統(tǒng)分為連續(xù)系統(tǒng)和離散系統(tǒng)。(1)連續(xù)系統(tǒng)——系統(tǒng)輸入、輸出信號(xìnhào)都是連續(xù)時間信號(xìnhào)?！惨话鉒、R、C電路〕(2)離散系統(tǒng)——系統(tǒng)輸入、輸出信號(xìnhào)都是離散時間信號(xìnhào)。〔數(shù)字計算機〕(3)混合系統(tǒng)——系統(tǒng)輸入、輸出信號(xìnhào)包含連續(xù)信號(xìnhào)和離散信號(xìnhào)。〔計算機控制系統(tǒng)〕連續(xù)時間系統(tǒng)的數(shù)學(xué)模型用微分方程描述。離散時間系統(tǒng)的數(shù)學(xué)模型用差分方程描述。第三頁，共79頁。二.線性和非線性系統(tǒng)(xìtǒng)根據(jù)輸入(shūrù)輸出關(guān)系是否同時滿足齊次性和疊加性，系統(tǒng)分為線性和非線性。假設(shè)系統(tǒng)在沒有外界信號作用之前處于靜止?fàn)顟B(tài)，在輸入(shūrù)信號u1和u2或a1u1和a2u2作用下，有式中，為任意實數(shù)，為輸入(shūrù)輸出之間函數(shù)關(guān)系。那么，該系統(tǒng)稱為線性系統(tǒng)，否那么是非線性系統(tǒng)。系統(tǒng)類型第四頁，共79頁。系統(tǒng)類型線性時變(shíbiàn)系統(tǒng)：非線性定常系統(tǒng)(xìtǒng)：線性定常系統(tǒng)(xìtǒng)：根據(jù)模型參數(shù)是否隨時間變化，線性系統(tǒng)又可細分為線性定常系統(tǒng)和線性時變系統(tǒng)。參數(shù)不隨時間變化的系統(tǒng)，稱為時不變系統(tǒng)或定常系統(tǒng)，否那么稱為時變系統(tǒng)。第五頁，共79頁。

系統(tǒng)類型根據(jù)系統(tǒng)輸入、輸出和內(nèi)部狀態(tài)呈現(xiàn)的規(guī)律，系統(tǒng)分為確定性系統(tǒng)與隨機性系統(tǒng)。輸入輸出之間函數(shù)關(guān)系能夠(nénggòu)用確定性模型描述的系統(tǒng)，稱為確定性系統(tǒng)，否那么稱為隨機系統(tǒng)(或不確定性系統(tǒng))。例如式中，

分別(fēnbié)為狀態(tài)變量和輸出變量，為噪聲(zàoshēng)。第六頁，共79頁。控制系統(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型一.連續(xù)(liánxù)系統(tǒng)式中，分別(fēnbié)為系統(tǒng)輸入量、輸出量，n為系統(tǒng)的階次，及各階導(dǎo)數(shù)的的初始值為為參數(shù)，均為實常數(shù)。輸出變量1．微分方程一個連續(xù)系統(tǒng)可以表示成高階微分方程，即第七頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型稍加整理(zhěnglǐ)，并記 2．傳遞函數(shù)假設(shè)系統(tǒng)的初始條件為零，那么(nàme)對微分方程兩邊取拉普拉斯變換后可得稱為系統(tǒng)的傳遞函數(shù)。第八頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型=AX+BU狀態(tài)方程Y=CX+DU輸出(shūchū)方程系統(tǒng)的動態(tài)特性是用由狀態(tài)變量構(gòu)成的一階微分方程組來描述。狀態(tài)空間表達式包括狀態(tài)方程和輸出(shūchū)方程。線性定常系統(tǒng)的狀態(tài)空間描述為3．狀態(tài)空間描述r維輸入向量

n維狀態(tài)向量m維輸出向量

第九頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型維輸出(shūchū)矩陣維直接(zhíjiē)傳遞矩陣維系統(tǒng)矩陣維輸入矩陣

第十頁，共79頁。控制系統(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

對于線性時變系統(tǒng)，系數(shù)(xìshù)矩陣A，B，C，D，均與時間t有關(guān)，狀態(tài)空間描述為第十一頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型1．差分(chàfēn)方程設(shè)系統(tǒng)差分(chàfēn)方程為引進(yǐnjìn)后移算子可得二.離散系統(tǒng)第十二頁，共79頁。控制系統(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型2．離散(lísàn)傳遞函數(shù)(Z傳函)假設(shè)系統(tǒng)的初始條件為零，即那么(nàme)得系統(tǒng)傳遞函數(shù)為

在初始條件為零時，與等價。第十三頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

3.離散狀態(tài)空間(kōngjiān)描述多變量離散狀態(tài)空間(kōngjiān)表達式第十四頁，共79頁。控制系統(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型三.MATLAB模型表示MATLAB的控制系統(tǒng)(kònɡzhìxìtǒnɡ)工具箱提供傳遞函數(shù)模型、零極點增益模型、狀態(tài)空間模型的生成函數(shù)。函數(shù)功能Sys=tf(num，den)生成傳遞函數(shù)模型Sys=zpk(z，p，k)生成零極點增益模型Sys=ss(a，b，c，d)生成狀態(tài)空間模型線性模型(móxíng)生成函數(shù)第十五頁，共79頁。控制系統(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

1.傳遞函數(shù)模型(móxíng)(transferfunctionmodel:TF)傳遞函數(shù)模型(móxíng)由分子和分母(fēnmǔ)多項式系數(shù)可以唯一確定傳遞函數(shù)。分子向量num=分母向量den=第十六頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

用命令tf()可建立(jiànlì)傳遞函數(shù)模型，或?qū)⒘銟O點增益模型和狀態(tài)空間模型化為傳遞函數(shù)模型。sys=tf(num，den)；%用于生成連續(xù)傳遞函數(shù)；sys=tf(num，den，Ts)；%用于生成離散傳遞函數(shù)；sys=tf(num，den，'Property1'，Value1，...，'PropertyN'，ValueN)；%用于生成具有LTI模型屬性的連續(xù)傳遞函數(shù)；sys=tf(num，den，Ts，‘Property1’，Value1，...，‘PropertyN’，ValueN)；%用于生成具有LTI模型屬性的離散傳遞函數(shù)；tfsys=tf(sys)；%用于將任意狀態(tài)空間模型SS或零極點增益模型ZPK的LTI對象sys轉(zhuǎn)換成傳遞函數(shù)形式第十七頁，共79頁。控制系統(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型給定(ɡěidìnɡ)SISO系統(tǒng)傳遞函數(shù)為使用(shǐyòng)MATLAB表示該傳遞函數(shù)num=[21];den=[341];sys1=tf(num,den)運行結(jié)果：Transferfunction:2s+1---------------3s^2+4s+1第十八頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

給定(ɡěidìnɡ)SISO系統(tǒng)傳遞函數(shù)為使用MATLAB表示(biǎoshì)該傳遞函數(shù)num=[1.322.5];den=[10.51.21];sys2=tf(num,den,'inputdelay',2)運行結(jié)果：Transferfunction:

exp(-2*s)*-------------------------s^3+0.5s^2+1.2s+1第十九頁，共79頁。控制系統(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

給定多輸入-多輸出(shūchū)MIMO系統(tǒng)MATLAB命令(mìnglìng)：num={11;2[11]};den={[11],[12];1,[12]};sys3=tf(num,den)第二十頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

結(jié)果(jiēguǒ)為：Transferfunctionfrominput1tooutput...1#1:-------s+1

#2:2

Transferfunctionfrominput2tooutput...1#1:-------s+2

s+1#2:-------s+2第二十一頁，共79頁。控制系統(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

MATLAB命令如下：num={[11],[10];1,2};%分子(fēnzǐ)den={[121],[102];[21],[11]};%分母sys4=tf(num,den,0.2)第二十二頁，共79頁。控制系統(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型Transferfunctionfrominput1tooutput...z+1#1:-------------z^2+2z+11#2:-------2z+1

Transferfunctionfrominput2tooutput...z#1:-------z^2+22#2:-----z+1第二十三頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

2.零極點增益模型(zero-pole-gainmodel:ZPK)零極點模型是傳遞函數(shù)的一種特殊(tèshū)形式離散系統(tǒng)連續(xù)(liánxù)系統(tǒng)第二十四頁，共79頁。控制系統(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型用命令zpk()可以建立零極點增益模型，或?qū)鬟f函數(shù)模型和狀態(tài)(zhuàngtài)空間模型變化為零極點增益模型。語法格式sys=zpk(z，p，k)sys=zpk(z，p，k，Ts)sys=zpk(z，p，k，'Property1'，Value1，...，'PropertyN'，ValueN)sys=zpk(z，p，k，Ts，'Property1'，Value1，...，'PropertyN'，ValueN)zsys=zpk(sys)zpk()函數(shù)調(diào)用方法與tf()一致。第二十五頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

給定一零極點(jídiǎn)增益模型使用MATLAB表示(biǎoshì)該傳遞函數(shù)。MATLAB命令如下：z={0;-0.5};p={0.3;[0.1-j,0.1+j]};k=[1;2];sys=zpk(z,p,k,[])第二十六頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

Zero/pole/gainfrominputtooutput...z#1:----------(z-0.3)2(z+0.5)#2:-----------------------(z^2-0.2z+1.01)Samplingtime:unspecified第二十七頁，共79頁。控制系統(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型將傳遞函數(shù)化為零極點增益(zēngyì)模型Matlab命令(mìnglìng)：num=[1.322.5];den=[10.51.21];sys2=tf(num,den,'inputdelay',2)sys6=zpk(sys2) 第二十八頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

Transferfunction:exp(-2*s)*------------------------------s^3+0.5s^2+1.2s+1Zero/pole/gain:1.3(s^2+1.538s+1.923)exp(-2*s)*---------------------------------------------(s+0.7307)(s^2-0.2307s+1.369)第二十九頁，共79頁。控制系統(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

3．狀態(tài)空間(kōngjiān)模型(state-spacemodel:SS)線性定常狀態(tài)空間(kōngjiān)模型描述為式中，為狀態(tài)(zhuàngtài)向量，為輸入向量，是輸出向量。第三十頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

用命令ss()可以建立狀態(tài)空間模型，或?qū)鬟f函數(shù)模型和零極點增益模型轉(zhuǎn)化為狀態(tài)空間模型。語法(yǔfǎ)調(diào)用格式:sys=ss(a，b，c，d)sys=ss(a，b，c，d，Ts)sys=ss(a，b，c，d，'Property1'，Value1，...，'PropertyN'，ValueN)sys=ss(a，b，c，d，Ts，'Property1'，Value1，...，'PropertyN'，ValueN)sys_ss=ss(sys)sys_ss=ss(sys，'minimal');%最小階實現(xiàn)ss函數(shù)的調(diào)用方法與tf()、zpk()一致。第三十一頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

狀態(tài)空間(kōngjiān)模型用MATLAB表示(biǎoshì)為A=[010;001;-5-20-1];B=[0;0;1];C=[100];D=0;sys=ss(A,B,C,D)第三十二頁，共79頁?？刂葡到y(tǒng)(kònɡzhìxìtǒnɡ)常用數(shù)學(xué)模型

結(jié)果(jiēguǒ)為：a=x1x2x3x1010x2001x3-5-20-1b=u1x10x20x31c=x1x2x3y1100d=u1y10Continuous-timemodel。第三十三頁，共79頁。利用(lìyòng)MATLAB實現(xiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

MATLAB實現(xiàn)模型(móxíng)轉(zhuǎn)換有兩種不同的方式。方式1：簡單的模型(móxíng)轉(zhuǎn)換首先生成任一指定的模型(móxíng)對象(tf，ss，zpk)，然后將該模型(móxíng)對象類作為輸入，調(diào)用欲轉(zhuǎn)換的模型(móxíng)函數(shù)即可。例如：欲將傳遞函數(shù)轉(zhuǎn)換為狀態(tài)空間模型(móxíng)sys=tf(num，den)；[a，b，c，d]=ss(sys)第三十四頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

方式2：直接調(diào)用模型轉(zhuǎn)換函數(shù)連續(xù)模型之間三種形式(xíngshì)的數(shù)學(xué)模型相互轉(zhuǎn)換函數(shù)類型包括tf2ss、ss2tf、zp2tf、tf2zp、ss2zp、zp2ss，共六個函數(shù)。第三十五頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換函數(shù)tf2ss用于將傳遞函數(shù)化成狀態(tài)(zhuàngtài)空間模型，調(diào)用格式如[A，B，C，D]=tf2ss(num，den);其中，輸入num，den分別為傳遞函數(shù)分子和分母多項式系數(shù)；輸出A，B，C，D為狀態(tài)(zhuàngtài)空間模型系數(shù)矩陣。第三十六頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

將傳遞函數(shù)化成狀態(tài)(zhuàngtài)空間表達式方式1：num={[023];[121]};%分子多項式系數(shù)den={[10.41];[10.41]};%分母多項式系數(shù)tfsys=tf(num,den);%生成傳遞函數(shù)sssys=ss(tfsys);%轉(zhuǎn)化為狀態(tài)(zhuàngtài)空間模型第三十七頁，共79頁。利用(lìyòng)MATLAB實現(xiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

結(jié)果(jiēguǒ)為：a=x1x2x220b=u1x12x20c=x1x2y20.80d=u1y10y21第三十八頁，共79頁。利用(lìyòng)MATLAB實現(xiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

方式(fāngshì)2：num=[023;121];den=[10.41];[A,B,C,D]=tf2ss(num,den)結(jié)果為：1.00000B=101.60000D=01第三十九頁，共79頁。利用(lìyòng)MATLAB實現(xiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換二.化傳遞函數(shù)為零極點(jídiǎn)增益模型函數(shù)tf2zp用于將傳遞函數(shù)化成零極點增益模型形式[z，p，k]=tf2zp(num，den)其中，num，den分別為傳遞函數(shù)的分子和分母多項式系數(shù)；z，p，k為零極點增益模型的零點(línɡdiǎn)、極點和增益向量。第四十頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換將傳遞函數(shù)化為零極點(jídiǎn)增益模型。方式(fāngshì)1：zpk(tf([21],[341]))Zero/pole/gain:0.66667(s+0.5)----------------(s+1)(s+0.3333)方式2：[z,p,k]=tf2zp([21],[341])z=p=k=因此零極點增益模型為第四十一頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換三.化零極點增益模型(móxíng)為狀態(tài)空間模型(móxíng)函數(shù)zp2ss用于將零極點(jídiǎn)增益模型化成狀態(tài)空間模型，zp2ss調(diào)用格式[A，B，C，D]=zp2ss(z，p，k)其中，輸入z，p，k為零極點增益模型的零點、極點和增益向量。輸出A，B，C，D為狀態(tài)空間模型系數(shù)矩陣。第四十二頁，共79頁。利用(lìyòng)MATLAB實現(xiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換零極點(jídiǎn)增益模型化成狀態(tài)空間表達式方式(fāngshì)1：z=[-2];p=[-1-3-3];k=2;zpksys=zpk(z,p,k);sssys=ss(zpksys)方式2:z=[-2];p=[-1-3-3];k=2;[A,B,C,D]=zp2ss(z,p,k)第四十三頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換四.化連續(xù)(liánxù)狀態(tài)方程為離散狀態(tài)方程設(shè)連續(xù)系統(tǒng)(xìtǒng)狀態(tài)空間方程為假設(shè)輸入端加上虛擬采樣開關(guān)和虛擬信號重構(gòu)器，輸出端加一個虛擬采樣開關(guān)，虛擬采樣周期為T，兩者同步。第四十四頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換1.采用零階保持器的離散狀態(tài)方程設(shè)kT及(k+1)T為兩個依次相連(xiānɡlián)的采樣瞬時，保持器為零階保持器，u(t)在kT及(k+1)T之間保持不變式(2)－式(1)乘以eAT，可得(1)(2)第四十五頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換離散(lísàn)的狀態(tài)空間表達式為式中，不變；為采樣(cǎiyànɡ)周期。第四十六頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換2.采用一階保持器的離散狀態(tài)方程保持器使為斜坡(xiépō)函數(shù)(梯形近似)，那么對離散(lísàn)的狀態(tài)空間表達式為式中，不變。第四十七頁，共79頁。利用(lìyòng)MATLAB實現(xiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

3.連續(xù)系統(tǒng)離散化的MATLAB方法c2d()可以將連續(xù)系統(tǒng)化成等價的離散化模型，格式如下：sysd=c2d(sys，Ts)sysd=c2d(sys，Ts，method)其中，sys為連續(xù)時間系統(tǒng)模型，Ts為采樣周期，單位為秒(s)。method定義離散化方法，method的取(1)‘zoh’——采用零階保持器。(2)‘foh’——采用一階保持器。(3)‘tustin’——采用雙線性(tustin)逼近(bījìn)方法。(4)‘matched’——采用SISO系統(tǒng)的零極點匹配法。默認時，method=‘zoh’。第四十八頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

采用一階保持器，采樣周期為0.5s，離散(lísàn)化以下狀態(tài)方程MATLAB程序(chéngxù)：sys=ss([01;0-2],[01]',[10],0);dss=c2d(sys,0.5,'foh')第四十九頁，共79頁。利用(lìyòng)MATLAB實現(xiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換結(jié)果(jiēguǒ)為：a=x1x2b=u1Discrete-timemodel

第五十頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

消去中間(zhōngjiān)項后，可得系統(tǒng)(xìtǒng)的傳遞函數(shù)矩陣

對方程式兩端分別取拉普拉斯變換(設(shè)初始條件為零)，得假設(shè)狀態(tài)空間表達式第五十一頁，共79頁。利用(lìyòng)MATLAB實現(xiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

將狀態(tài)空間方程化為傳遞函數(shù)的方法(fāngfǎ)有：1如果系統(tǒng)只有一個輸入，那么可采用[num，den]=ss2tf(A，B，C，D)或[num，den]=ss2tf(A，B，C，D，1)2對多輸入的系統(tǒng)，那么可采用[num，den]=ss2tf[A，B，C，D，iu]這里，iu是輸入的標(biāo)號。例如，如果系統(tǒng)有三個輸入(u1，u2，u3)，那么iu必須為1、2或3中的一個，其中1表示u1，2表示u2，3表示u3。第五十二頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換將狀態(tài)空間(kōngjiān)模型轉(zhuǎn)換為傳遞函數(shù)和零極點增益(zēngyì)模型第五十三頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換MATLAB命令(mìnglìng)：A=[010;001;-5-20-1];B=[0;0;1];C=[100];D=0;[num,den]=ss2tf(A,B,C,D); tf(num,den) [z,p,k]=ss2zp(A,B,C,D); zpk(z,p,k) 第五十四頁，共79頁。利用(lìyòng)MATLAB實現(xiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

Transferfunction:-1.11e-015s^2+1.421e-014s+1---------------------------------s^3+s^2+20s+5Zero/pole/gain:1----------------------------------(s+0.2524)(s^2+0.7476s+19.81)第五十五頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

多入單出系統(tǒng)狀態(tài)(zhuàngtài)空間表達式包括(bāokuò)兩個傳遞函數(shù)：Y1(s)/U1(s)、Y1(s)/U2(s)(當(dāng)考慮輸入u1時，可設(shè)u2為零。反之亦然)，試將其化成傳遞函數(shù)形式第五十六頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

MATLAB命令(mìnglìng)：A=[010;0 11;-5.008 25.1026-5.032471];B=[01;25.042;121.0053];C=[100];D=[00];[num,den]=ss2tf(A,B,C,D,1);Sys1=tf(num,den)傳遞函數(shù)Y1(s)/U1(s)為Transferfunction:-3.553e-015s^2+25.04s+247-------------------------------------------第五十七頁，共79頁。利用MATLAB實現(xiàn)(shíxiàn)數(shù)學(xué)模型之間的轉(zhuǎn)換

[num,den]=ss2tf(A,B,C,D,2);Sys2=tf(num,den)傳遞函數(shù)Y1(s)/U2(s)為Transferfunction:

------------------------------------------第五十八頁，共79頁。模型實現(xiàn)

所謂實現(xiàn)，就是根據(jù)描述系統(tǒng)輸入/輸出動態(tài)關(guān)系的運動方程式或傳遞函數(shù)建立(jiànlì)系統(tǒng)的狀態(tài)空間表達式。所求得的狀態(tài)空間表達式保持原傳遞函數(shù)的輸入/輸出關(guān)系，同時反映內(nèi)部動態(tài)變化。實現(xiàn)不是唯一的，會有無窮多個狀態(tài)空間表達式能夠獲得相同的輸入/輸出關(guān)系。第五十九頁，共79頁。模型實現(xiàn)并非任意的微分方程和傳遞函數(shù)都能求得其實現(xiàn)，實現(xiàn)存在(cúnzài)的條件是考慮單變量(biànliàng)線性定常系統(tǒng)傳遞