自動控制原理第2章

第二章線性系統(tǒng)的數(shù)學描述

2.1線性系統(tǒng)的時域數(shù)學模型

2.2傳遞函數(shù)

2.3結構圖

2.4信號流圖及梅遜公式

2.5線形系統(tǒng)的狀態(tài)變量描述（了解）

2.6線性定常系統(tǒng)數(shù)學模型的MATLAB實現(xiàn)

小結為什么要建立系統(tǒng)模型？物理模型：一個理想的物理系統(tǒng)數(shù)學模型：物理模型的數(shù)學描述調(diào)查研究，抓本質(zhì)、主流;忽略非本質(zhì)和次要因數(shù)求解觀察線性微分方程性能指標傳遞函數(shù)時間響應頻率響應拉氏變換拉氏反變換估算估算計算傅氏變換S=jω頻率特性建?！⑵鸨容^簡單又能反映實際物理過程的模型。

建模的線性化問題

兩種基本方法：機理分析法和實驗辨識法。2.1線性系統(tǒng)的時間函數(shù)描述

控制系統(tǒng)中的輸出量和輸入量通常都是時間t的函數(shù)。給定量和擾動量稱為輸入量；被控量稱作輸出量。很多常見的元件或系統(tǒng)的輸出量和輸入量之間的關系都可以用一個微分方程表示，方程中含有輸出量、輸入量及它們各自對時間的導數(shù)或積分。這種微分方程又稱為動態(tài)方程或運動方程。微分方程的階數(shù)一般是指方程中最高導數(shù)項的階數(shù)，又稱為系統(tǒng)的階數(shù)。對于單輸入、單輸出線性定常系統(tǒng)，采用下列微分方程來描述：式中，r(t)是系統(tǒng)的輸入信號，c(t)是系統(tǒng)的輸出信號；c(n)(t)為c(t)對時間t的n階導數(shù)；ai(i=1，2，…，n)和bj(j=0，1，…，m)是由系統(tǒng)的結構參數(shù)決定的系數(shù)。

(2.1)c(t)=H(t)r(t)r(t):系統(tǒng)輸入c(t)：系統(tǒng)響應H(t)：算子

commonly，列寫控制系統(tǒng)運動方程的步驟是:1）分析系統(tǒng)工作原理及各變量之間的關系，找出系統(tǒng)的輸入量和輸出量；所謂標準形式是指在系統(tǒng)運動方程中將輸入變量及其導數(shù)置于等號的右邊，將輸出變量及其導數(shù)置于等號左邊，等號兩邊的導數(shù)項均按降冪排列，并且將系數(shù)規(guī)劃為反映系統(tǒng)動態(tài)特性的參數(shù)，如時間常數(shù)、阻尼系數(shù)等。2）根據(jù)系統(tǒng)運動特性的基本定律，一般從系統(tǒng)的輸入端開始依次寫出各元件的運動方程，在列寫元件運動方程時，需要考慮相接元件間的相互作用；3）由組成系統(tǒng)各元件的運動方程中，消去中間變量，求取只含有系統(tǒng)輸入和輸出變量及其各階導數(shù)的方程，并將其化為標準形式。2.1.1電氣系統(tǒng)

電氣系統(tǒng)中最常見的裝置是由電阻、電容、運算放大器等元件組成的電路，又稱電氣網(wǎng)絡。我們將電阻、電感和電容等本身不含有電源的器件稱為無源器件，而將運算放大器這樣本身包含電源的器件稱為有源器件。僅由無源器件構成的電氣網(wǎng)絡稱為無源網(wǎng)絡；如果電氣網(wǎng)絡中含有有源器件或電源，就稱之為有源網(wǎng)絡。

圖

2.1RLC無源網(wǎng)絡

例2.1圖2.1是由電阻R、電感L和電容C組成的無源網(wǎng)絡，試列寫以ui(t)為輸入量，以uo(t)為輸出量的網(wǎng)絡微分方程。消去中間變量i(t)，可得描述該無源網(wǎng)絡輸入輸出關系的微分方程

(2.2)解：設回路電流為i(t)，由基爾霍夫電壓定律可寫出回路方程為上式也可以寫為

(2.3)其中，T1=L/R，T2=RC。方程(2.2)和(2.3)就是所求的微分方程。這是一個典型的二階線性常系數(shù)微分方程，對應的系統(tǒng)稱為二階線性定常系統(tǒng)。

例2.2

圖2.2是一個由理想運算放大器組成的電容負反饋電路，電壓ui(t)和uo(t)分別表示輸入量和輸出量，試確定這個電路的微分方程式。圖

2.2電容負反饋電路

解：理想運算放大器正、反相輸入端電位相同，且輸入電流為零。根據(jù)基爾霍夫電流定律有？

解理想運算放大器正、反相輸入端電位相同，且輸入電流為零。根據(jù)基爾霍夫電流定律有

整理后得

(2.4)或為

(2.5)其中，T=RC為時間常數(shù)。方程(2.4)和(2.5)就是該系統(tǒng)的微分方程，這是一個一階系統(tǒng)。

圖2.3機械阻尼器示意圖k是彈簧系數(shù)，m是運動部件質(zhì)量，μ是阻尼器的阻尼系數(shù)；外力f(t)是系統(tǒng)的輸入量，位移y(t)是系統(tǒng)的輸出量。例2.3(2.6)2.1.2機械系統(tǒng)

例2.3

試確定圖2.3系統(tǒng)的微分方程。

解根據(jù)牛頓運動定律，運動部件在外力作用下克服彈簧拉力ky(t)、阻尼器的阻力 ，將產(chǎn)生加速度力所以系統(tǒng)的運動方程為或?qū)懗?/p>

(2.7)這也是一個二階線性常微分方程。比較表達式(2.7)和(2.3)可以發(fā)現(xiàn)，兩個不同的物理系統(tǒng)具有相同形式的運動方程，即具有相同的數(shù)學模型。

例2.4圖2.4表示一個單擺系統(tǒng)，輸入量為零(不加外力)，輸出量為擺幅θ(t)。擺錘的質(zhì)量為M，擺桿長度為l，阻尼系數(shù)為μ，重力加速度為g。試建立系統(tǒng)的運動方程。圖2.4單擺運動示意圖解對于圖2.4所示的單擺系統(tǒng)，根據(jù)牛頓運動定律可以直接推出如下系統(tǒng)運動方程：顯然方程(2.8)是一個二階的非線性微分方程(因為含有sinθ)，但是在擺幅較小的情況下，單擺運動方程可以認為是線性的，對應的微分方程為

(2.8)(2.9)

在工程實際中，大多數(shù)系統(tǒng)是非線性的。比如，彈簧的剛度與其形變有關系，因此彈簧系數(shù)k實際上是其位移的函數(shù)，而并非常數(shù)；電阻、電容和電感等參數(shù)值與周圍的環(huán)境(溫度、濕度、壓力等)及流經(jīng)它們的電流有關，也并非常值；電動機本身的摩擦、死區(qū)等非線性因素會使其運動方程復雜化而成為非線性方程。非線性系統(tǒng)的分析一般比線性系統(tǒng)復雜。但是當控制系統(tǒng)在圍繞平衡點附近的小范圍內(nèi)動作時，通常采用泰勒級數(shù)展開的方法，可將非線性系統(tǒng)線性化為平衡點附近的線性系統(tǒng)，從而使問題簡化。卷積單位脈沖響應：對輸入r(t)分解，脈沖寬度，則發(fā)生在時刻脈沖強度為的脈沖即可表示為則當取為在區(qū)間內(nèi)的變量時，則輸入r(t)的近似表示為：則響應可近似表示為：當時，響應為：2.2傳

遞

函

數(shù)

2.2.1拉氏變換

1.拉氏變換的定義若將實變量t的函數(shù)f(t)乘上指數(shù)函數(shù)e-st(其中s=σ+jω是一個復數(shù))，并且在［0，+∞］上對t積分，就可以得到一個新的函數(shù)F(s)，稱F(s)為f(t)的拉氏變換，并用符號L［f(t)］表示。

(2.10)上式就是拉氏變換的定義式。從這個定義可以看出，拉氏變換將原來的實變量函數(shù)f(t)轉(zhuǎn)化為復變量函數(shù)F(s)。

通常將F(s)稱作f(t)的象函數(shù)，將f(t)稱作F(s)的原函數(shù)。

2.拉氏變換的基本定理(2.11)函數(shù)放大k倍的拉氏變換等于該函數(shù)拉氏變換的k倍，即

(2.12)1)線性定理兩個函數(shù)和的拉氏變換，等于每個函數(shù)拉氏變換的和，即2)微分定理如果初始條件

成立，則有

(2.13)

3)積分定理一個函數(shù)積分后再取拉氏變換等于這個函數(shù)的拉氏變換除以復參數(shù)s，即

重復運用式(2.14)可以推出

(2.14)(2.15)

4)初值定理函數(shù)f(t)在t=0時的函數(shù)值可以通過f(t)的拉氏變換F(s)乘以s取s→∞時的極限而得到，即

(2.16)

5)終值定理函數(shù)f(t)在t→+∞時的函數(shù)值(即穩(wěn)定值)可以通過f(t)的拉氏變換F(s)乘以s取s→0時的極限而得到，即

(2.17)2.2.2傳遞函數(shù)的定義和特點

1.傳遞函數(shù)的定義線性定常系統(tǒng)的傳遞函數(shù)，定義為零初始條件下，式中，r(t)和c(t)分別是系統(tǒng)的輸入信號和輸出信號；系統(tǒng)輸出量的拉氏變換與輸入量的拉氏變換之比。設線性定常系統(tǒng)由下面的n階線性常微分方程描述：

c(n)(t)為c(t)對時間t的n階導數(shù)；ai(i=0，1，…，n)和bj(j=0，1，…，m)是由系統(tǒng)的結構參數(shù)決定的常系數(shù)。如果r(t)和c(t)及其各階導數(shù)在t=0時的值均為零，即滿足如下的零初始條件：則根據(jù)拉氏變換的定義和性質(zhì)，對式(2.18)進行拉氏變換，并令C(s)=L［c(t)］，R(s)=L［r(t)］，可得

由傳遞函數(shù)的定義可得系統(tǒng)的傳遞函數(shù)為

式中，M(s)和N(s)分別稱為傳遞函數(shù)G(s)的分子多項式和分母多項式。

例2.5

試確定圖2.1所示的RLC無源網(wǎng)絡系統(tǒng)的傳遞函數(shù)。解

由例2.1可知，RLC無源網(wǎng)絡系統(tǒng)的微分方程為

在零初始條件下，對上述方程中各項求拉氏變換，并令Uo(s)=L［uo(t)］，

Ui(s)=L［ui(t)］，可得復頻域的代數(shù)方程

(LCs2+RCs+1)Uo(s)=Ui(s)所以系統(tǒng)的傳遞函數(shù)為

例2.6

試確定如圖2.2所示的運算放大器電路的傳遞函數(shù)。

解由例2.2可知，運算放大器電路系統(tǒng)的微分方程為

在零初始條件下，對上述方程中各項求拉氏變換，

得

所以，系統(tǒng)的傳遞函數(shù)為

(2.21)

例2.7

試確定如圖2.3所示的機械阻尼系統(tǒng)的傳遞函數(shù)。

解

由例2.3可知，該系統(tǒng)的運動方程為

在零初始條件下，對上式進行拉氏變換，

可得

(ms2+μs+k)Y(s)=F(s)所以系統(tǒng)的傳遞函數(shù)為

(2.22)

2.傳遞函數(shù)的特點

(1)傳遞函數(shù)的概念適用于線性定常系統(tǒng)，傳遞函數(shù)的結構和各項系數(shù)(包括常數(shù)項)完全取決于系統(tǒng)本身結構，因此，它是系統(tǒng)的動態(tài)數(shù)學模型，而與輸入信號的具體形式和大小無關，也不反映系統(tǒng)的任何內(nèi)部信息。因此可以用圖2.5的方塊圖來表示一個具有傳遞函數(shù)G(s)的線性系統(tǒng)。該圖說明，系統(tǒng)輸入量和輸出量的因果關系可以用傳遞函數(shù)聯(lián)系起來。

圖2.5傳遞函數(shù)的圖示但是同一個系統(tǒng)若選擇不同的量作為輸入量和輸出量，所得到的傳遞函數(shù)可能不同。所以談到傳遞函數(shù)，必須指明輸入量和輸出量。傳遞函數(shù)的概念主要適用于單輸入、單輸出的情況。若系統(tǒng)有多個輸入信號，在求傳遞函數(shù)時，除了指定的輸入量以外，其它輸入量(包括常值輸入量)一概視為零；對于多輸入、多輸出線性定常系統(tǒng)，求取不同輸入和輸出之間的傳遞函數(shù)將得到系統(tǒng)的傳遞函數(shù)矩陣。

(2)傳遞函數(shù)是在零初始條件下定義的。控制系統(tǒng)的零初始條件有兩層含義：一是指輸入量在t≥0時才起作用；二是指輸入量加于系統(tǒng)之前，系統(tǒng)處于穩(wěn)定工作狀態(tài)。

(3)傳遞函數(shù)是復變量s的有理真分式函數(shù)，具有復變函數(shù)的所有性質(zhì)；并且理論分析和實驗都指出，對于實際的物理系統(tǒng)和元件而言，輸入量和它所引起的響應(輸出量)之間的傳遞函數(shù)，分子多項式M(s)的階次m總是小于分母多項式N(s)的階次n，即m＜n。這個結論可以看作是客觀物理世界的基本屬性。它反映了一個基本事實：一個物理系統(tǒng)的輸出不可能立即復現(xiàn)輸入信號，只有經(jīng)過一段時間后，輸出量才能達到輸入量所要求的數(shù)值。對于具體的控制元件和系統(tǒng)，我們總是可以找到形成上述事實的原因。例如對于機械系統(tǒng)，由于物體都有質(zhì)量，物體受到外力和外力矩作用時都要產(chǎn)生形變，相互接觸并存在相對運動的物體之間總是存在摩擦，這些都是造成機械裝置傳遞函數(shù)分母階次高于分子階次的原因。電氣網(wǎng)絡中，由運算放大器組成的電壓放大器，如果考慮其中潛在的電容和電感，輸出電壓和輸入電壓間的傳遞函數(shù)，分子多項式的階次一定低于分母多項式的階次。

(4)傳遞函數(shù)與線性常微分方程一一對應。傳遞函數(shù)分子多項式系數(shù)和分母多項式系數(shù)，分別與相應微分方程的右端及左端微分算符多項式系數(shù)相對應。所以，將微分方程的算符d/dt用復數(shù)s置換便可以得到傳遞函數(shù)；反之，將傳遞函數(shù)中的復數(shù)s用算符d/dt置換便可以得到微分方程。例如，由傳遞函數(shù)

可得s的代數(shù)方程

(a0s2+a1s+a2)C(s)=(b1s+b2)R(s)用算符d/dt置換復數(shù)s，便得到相應的微分方程

(5)傳遞函數(shù)不能反映系統(tǒng)或元件的學科屬性和物理性質(zhì)。物理性質(zhì)和學科類別截然不同的系統(tǒng)可能具有完全相同的傳遞函數(shù)。例如，RLC電路和機械阻尼系統(tǒng)的傳遞函數(shù)在適當?shù)膮?shù)代換后可以具有相同的形式，但是兩者屬于完全不同的學科領域。另一方面，研究某一種傳遞函數(shù)所得到的結論，可以適用于具有這種傳遞函數(shù)的各種系統(tǒng)，不管它們的學科類別和工作機理如何不同。這就極大地提高了控制工作者的效率。

(6)傳遞函數(shù)除具有式(2.19)表示的分子、分母多項式形式外，還具有如下兩種常見形式：(2.23)(2.24)表達式(2.23)和(2.24)分別稱為傳遞函數(shù)的零極點形式和時間常數(shù)形式。式(2.23)的特點是每個一次因子項中s的系數(shù)為1。M(s)=0和N(s)=0的根zi(i=1，2，…，m)和pj(j=1，2，…，n)分別稱為傳遞函數(shù)的零點和極點，k稱為傳遞函數(shù)的增益或根軌跡增益。由于M(s)和N(s)的系數(shù)均為實數(shù)，因此零極點是實數(shù)或共軛復數(shù)。式(2.24)的特點是各個因式的常數(shù)項均為1，τi(i=1，2，…，m)和Tj(j=1，2，…，n)為系統(tǒng)中各環(huán)節(jié)的時間常數(shù)，K為系統(tǒng)的放大倍數(shù)。

(7)令系統(tǒng)的傳遞函數(shù)分母等于零，所得方程稱為特征方程，即N(s)=0。特征方程的根稱為特征根，也就是系統(tǒng)的極點。

非線性數(shù)學模型的線性化泰勒級數(shù)展開條件：函數(shù)在給定區(qū)域內(nèi)各階導數(shù)存在，在給定工作點的鄰域內(nèi)。注意事項：1工作點；2變換范圍；3存在不連續(xù)的情況；4線性化后得到的是增量方程。2.2.3典型環(huán)節(jié)傳遞函數(shù)

任何一個復雜系統(tǒng)都是由有限個典型環(huán)節(jié)組合而成的。典型環(huán)節(jié)通常分為以下六種。特點：輸入輸出量成比例，無失真和時間延遲。實例：電子放大器，齒輪，電阻(電位器)，感應式變送器等。

(2.25)(2.26)1.比例環(huán)節(jié)比例環(huán)節(jié)又稱放大環(huán)節(jié)，該環(huán)節(jié)的運動方程和相對應的傳遞函數(shù)分別為式中K為增益。

2.慣性環(huán)節(jié)慣性環(huán)節(jié)又稱非周期環(huán)節(jié)，該環(huán)節(jié)的運動方程和相對應的傳遞函數(shù)分別為

(2.27)(2.28)式中T為時間常數(shù)，K為比例系數(shù)。特點：含一個儲能元件，對突變的輸入，其輸出不能立即復現(xiàn)，輸出無振蕩。實例：直流伺服電動機的勵磁回路。3.純微分環(huán)節(jié)純微分環(huán)節(jié)常簡稱為微分環(huán)節(jié)，其運動方程和傳遞函數(shù)分別為

(2.29)(2.30)特點：輸出量正比輸入量變化的速度，能預示輸入信號的變化趨勢。實例：實際中沒有純粹的微分環(huán)節(jié)，它總是與其他環(huán)節(jié)并存的。實際中可實現(xiàn)的微分環(huán)節(jié)都具有一定的慣性，其傳遞函數(shù)如下：(2.31)4.積分環(huán)節(jié)積分環(huán)節(jié)的動態(tài)方程和傳遞函數(shù)分別為

(2.32)(2.33)特點：輸出量與輸入量的積分成正比例，當輸入消失，輸出具有記憶功能。實例：電動機角速度與角度間的傳遞函數(shù)，模擬計算機中的積分器等。

寫出傳遞函數(shù)列方程：I1I2ui(t)uo(t)RC作拉氏變換得：5.振蕩環(huán)節(jié)振蕩環(huán)節(jié)的運動方程和傳遞函數(shù)分別為

(2.34)(2.35)式中為振蕩環(huán)節(jié)的阻尼比，為時間常數(shù)，ωn為系統(tǒng)的自然振蕩角頻率(無阻尼自振角頻率)，并且有

特點：環(huán)節(jié)中有兩個獨立的儲能元件，并可進行能量交換，其輸出出現(xiàn)振蕩。實例：RLC電路的輸出與輸入電壓間的傳遞函數(shù)，以及機械阻尼系統(tǒng)的傳遞函數(shù)。

振蕩環(huán)節(jié)的單位階躍響應6.純時間延時環(huán)節(jié)延時環(huán)節(jié)的動態(tài)方程和傳遞函數(shù)分別為

(2.36)(2.37)特點：輸出量能準確復現(xiàn)輸入量，但要延遲一固定的時間間隔τ。實例：管道壓力、

流量等物理量的控制，其數(shù)學模型就包含有延遲環(huán)節(jié)。

式中τ為該環(huán)節(jié)的延遲時間。實例滯后環(huán)節(jié)的單位階躍響應具有傳輸滯后的系統(tǒng)2.3結

構

圖2.3.1結構圖的組成與繪制

1.結構圖的組成

(1)結構圖的每一元件用標有傳遞函數(shù)的方框表示，方框外面帶箭頭的線段表示這個環(huán)節(jié)的輸入信號(箭頭指向方框)和輸出信號(箭頭離開方框)，其方向表示信號傳遞方向。箭頭處標有代表信號物理量的符號字母，如圖2.6所示。

圖

2.6元件的結構圖

(2)然后把系統(tǒng)中所有元件都用上述方框形式表示，按系統(tǒng)輸入信號經(jīng)過各元件的先后次序，依次將代表各元件的方塊用連接線連接起來。顯然，前后兩方塊連接時，前面方塊輸出信號必為后面方塊的輸入信號。(3)對于閉環(huán)系統(tǒng)，需引入兩個新符號，分別稱為相加點和分支點。需要注意的是，分支點無論從一條信號流線或一個分支點引出多少條信號流線，它們都代表一個信號，即原始大小的信號。圖

2.7結構圖的相加點和分支點

圖2.7(a)是系統(tǒng)的比較元件，表示兩個以上信號的代數(shù)運算。箭頭指向的信號流線表示它的輸入信號，箭頭離開它的信號流線表示它的輸出信號，附近的＋、－號表示信號之間的運算關系是相加還是相減。在框圖中，可以從一條信號流線上引出另一條或幾條信號流線，而信號引出的位置稱為分支點或引出點如圖2.7(b)所示。連接方式串聯(lián)前一環(huán)節(jié)的輸出作為后一環(huán)節(jié)的輸入，等效環(huán)節(jié)為各個串聯(lián)環(huán)節(jié)的傳遞函數(shù)積。并聯(lián)各環(huán)節(jié)的輸入信號相同，輸出信號相加（或相減），等效環(huán)節(jié)為各個環(huán)節(jié)的相加。反饋連接將系統(tǒng)或環(huán)節(jié)的輸出信號反饋到輸入端，與輸入信號進行比較。

2.結構圖的繪制

(1)列寫系統(tǒng)的微分方程組，并求出其對應的拉氏變換方程組。

(2)從輸出量開始寫，以系統(tǒng)輸出量作為第一個方程左邊的量。

(3)每個方程左邊只有一個量。從第二個方程開始，每個方程左邊的量是前面方程右邊的中間變量。列寫方程時盡量用已出現(xiàn)過的量。

(4)輸入量至少要在一個方程的右邊出現(xiàn)；除輸入量外，在方程右邊出現(xiàn)過的中間變量一定要在某個方程的左邊出現(xiàn)。

(5)按照上述整理后拉氏變換方程組的順序，從輸出端開始繪制系統(tǒng)的結構圖。

例2.8

在圖2.8(a)中，電壓u1(t)、u2(t)分別為輸入量和輸出量，繪制系統(tǒng)的結構圖。

圖

2.8RC濾波電路結構圖

解:

電氣網(wǎng)絡采用電路理論中“運算阻抗”的概念和方法，不列寫微分方程就可以方便地求出相應的傳遞函數(shù)。電阻R的運算阻抗就是電阻R本身，電感L的運算阻抗是Ls，電容C的運算阻抗是1/(Cs)，其中s是拉氏變換的復參量。把電路中的電阻R、電感L和電容C全換成運算阻抗，把電流i(t)和電壓u(t)全換成相應的拉氏變換式I(s)和U(s)，把運算阻抗當作普通電阻。這樣從形式上看，在零初始條件下，電路中的運算阻抗和電流、電壓的拉氏變換式之間的關系滿足各種電路定律，如歐姆定律、基爾霍夫定律。從而采用普通的電路定律，經(jīng)過簡單的代數(shù)運算就可求解I(s)和U(s)及相應的傳遞函數(shù)。采用運算阻抗的方法又稱運算法，相應的電路圖稱為運算電路。圖2.8(a)對應的運算電路如圖2.8(b)所示。設中間變量I1(s)、I2(s)和U3(s)。從輸出量U2(s)開始按上述步驟列寫系統(tǒng)方程式：按照上述方程的順序，從輸出量開始繪制系統(tǒng)的結構圖，其繪制結果如圖2.8(c)所示(注意這是一個還沒有經(jīng)過簡化的系統(tǒng)結構圖)。值得注意的是，一個系統(tǒng)可以具有不同的結構圖，但由結構圖得到的輸出和輸入信號的關系都是相同的。

2.3.2閉環(huán)系統(tǒng)的結構圖

E(s)和B(s)分別為偏差信號和反饋信號的拉氏變換，H(s)為閉環(huán)系統(tǒng)中的反饋傳遞函數(shù)。圖2.9閉環(huán)系統(tǒng)結構圖2.3.2閉環(huán)系統(tǒng)的結構圖

一個閉環(huán)負反饋系統(tǒng)通常用圖2.9所示的結構圖來表示。

輸出量C(s)反饋到相加點，并且在相加點與參考輸入量R(s)進行比較。

圖中各信號之間的關系為

C(s)=G(s)E(s)E(s)=R(s)-B(s)B(s)=H(s)C(s)反饋信號B(s)與偏差信號E(s)之比，叫做開環(huán)傳遞函數(shù)，即

輸出量C(s)和偏差信號E(s)之比，叫做前向傳遞函數(shù)，即

如果反饋傳遞函數(shù)等于1，那么開環(huán)傳遞函數(shù)和前向傳遞函數(shù)相同，并稱這時的閉環(huán)反饋系統(tǒng)為單位反饋系統(tǒng)。從圖2.9可以推出系統(tǒng)輸出量C(s)和輸入量R(s)之間的關系，具體推導如下：C(s)=G(s)E(s)E(s)=R(s)-B(s)=R(s)-H(s)C(s)

消去E(s)可得

C(s)=G(s)［R(s)-H(s)C(s)］

所以有

(2.38)上式就是系統(tǒng)輸出量C(s)和輸入量R(s)之間的傳遞函數(shù)，稱為閉環(huán)傳遞函數(shù)。這個傳遞函數(shù)將閉環(huán)系統(tǒng)的動態(tài)特性與前向通道環(huán)節(jié)和反饋通道環(huán)節(jié)的動態(tài)特性聯(lián)系在一起。由方程(2.38)可得

可見，閉環(huán)系統(tǒng)的輸出量取決于閉環(huán)傳遞函數(shù)和輸入量的性質(zhì)。2.3.3擾動作用下的閉環(huán)系統(tǒng)

實際的系統(tǒng)經(jīng)常會受到外界擾動的干擾，通常擾動作用下的閉環(huán)系統(tǒng)的結構圖可由圖2.10表示。從圖2.10可知，這個系統(tǒng)存在兩個輸入量，即參考輸入量R(s)和擾動量N(s)。

圖

2.10擾動作用下的閉環(huán)系統(tǒng)結構圖

根據(jù)線性系統(tǒng)滿足疊加性原理的性質(zhì)，可以先對每一個輸入量單獨地進行處理，然后將每個輸入量單獨作用時相應的輸出量進行疊加，就可得到系統(tǒng)的總輸出量。對于圖2.10所示的系統(tǒng)，研究擾動量N(s)對系統(tǒng)的影響時，可以假設參考輸入信號R(s)=0，經(jīng)過簡單的推導可以得出系統(tǒng)對擾動的響應CN(s)為

所以，系統(tǒng)輸出對擾動的傳遞函數(shù)ΦN(s)=CN(s)/N(s)為

(2.39)同樣在分析系統(tǒng)對參考輸入的響應時，可以假設擾動量N(s)=0，這時系統(tǒng)對參考輸入量R(s)的響應CR(s)為

所以，系統(tǒng)輸出對參考輸入的傳遞函數(shù)Φ(s)=CR(s)/R(s)為

(2.40)根據(jù)線性系統(tǒng)的疊加原理可知，參考輸入量R(s)和擾動量N(s)同時作用于系統(tǒng)時，系統(tǒng)的響應(總輸出)C(s)為

2.3.4結構圖的簡化和變換規(guī)則

利用結構圖分析和設計系統(tǒng)時，常常要對結構圖進行簡化和變換。對結構圖進行簡化和變換的基本原則是等效原則，即對結構圖任何部分進行變換時，變換前后該部分的輸入量、輸出量及其相互之間的數(shù)學關系應保持不變。1.串聯(lián)環(huán)節(jié)的簡化幾個環(huán)節(jié)的結構圖首尾連接，前一個結構圖的輸出是后一個結構圖的輸入，稱這種結構為串聯(lián)環(huán)節(jié)。圖2.11(a)是三個環(huán)節(jié)串聯(lián)的結構。根據(jù)結構圖可知：以下是根據(jù)等效原則給出的幾條結構圖的變換規(guī)則。消去中間變量X1(s)和X2(s)得

所以此系統(tǒng)的等效傳遞函數(shù)為

圖

2.11三個環(huán)節(jié)串聯(lián)

上述結論可以推廣到任意個環(huán)節(jié)的串聯(lián)，即n個環(huán)節(jié)(每個環(huán)節(jié)的傳遞函數(shù)為Gi(s)，i=1，2，…，n)串聯(lián)的等效傳遞函數(shù)等于n個傳遞函數(shù)相乘。

G(s)=G1(s)G2(s)…Gn(s)(2.42)

2.并聯(lián)環(huán)節(jié)的簡化兩個或多個環(huán)節(jié)具有同一個輸入信號，而以各自環(huán)節(jié)輸出信號的代數(shù)和作為總的輸出信號，這種結構稱為并聯(lián)環(huán)節(jié)。圖2.12(a)表示三個環(huán)節(jié)并聯(lián)的結構，根據(jù)結構圖可知：圖

2.12三個環(huán)節(jié)并聯(lián)

所以，整個系統(tǒng)的等效傳遞函數(shù)為

(2.43)

3.反饋回路的簡化圖

2.13基本反饋回路的簡化

（b）±圖2.13(a)表示一個基本的反饋回路。根據(jù)2.3.2節(jié)的分析和得到的閉環(huán)傳遞函數(shù)形式可以推出(2.44)所以，圖2.13(a)所表示的反饋系統(tǒng)結構圖可簡化為圖2.13(b)。

4.相加點和分支點的移動在結構圖的變換中經(jīng)常要求改變相加點和分支點的位置。一般包括相加點前移、相加點后移、分支點前移和分支點后移四種基本情況，以及相鄰相加點和相鄰分支點之間的移動。1)相加點前移圖2.14(a)和圖2.14(b)分別表示相加點前移變換前后的系統(tǒng)結構圖。圖

2.14相加點前移

可以看出兩圖具有如下相同的輸入、

輸出關系：2)相加點后移

圖

2.15相加點后移

可以看出兩圖具有如下相同的輸入、

輸出關系：

3)分支點前移

圖

2.16分支點前移

可以看出兩圖具有如下相同的輸入、

輸出關系：

C(s)=R(s)G(s)

4)分支點后移圖

2.17分支點后移

可以看出兩圖具有如下相同的輸入、

輸出關系：

5)相鄰相加點之間的移動D=A±B±C=A±C±B

并且，這個結論對于多個相鄰的相加點也適用。

如圖2.18所示，相鄰相加點之間可以互換位置而不改變該結構輸入和輸出信號之間的關系。圖

2.18相加點之間的移動

6)相鄰分支點之間的移動從一個信號流線上無論分出多少條信號線，它們都代表同一個信號。所以在一條信號流線上的各分支點之間可以隨意改變位置，不必作任何其他改動(如圖2.19所示)。

圖

2.19相鄰分支點的移動

例2.9

試簡化圖2.20系統(tǒng)的結構圖，并求系統(tǒng)的傳遞函數(shù)C(s)/R(s)。

圖

2.20系統(tǒng)結構圖

解在圖2.20中，如果不移動相加點或分支點的位置就無法進行結構圖的等效運算。采用以下步驟簡化原圖：①利用分支點后移規(guī)則，將G3(s)和G4(s)之間的分支點移到G4(s)方框的輸出端(注意不宜前移)，變換結果如圖2.21(a)所示；②將G3(s)、G4(s)和H3(s)組成的內(nèi)反饋回路簡化(如圖2.21(b)所示)，其等效傳遞函數(shù)為③再將G2(s)、G34(s)、H2(s)和1/G4(s)組成的內(nèi)反饋回路簡化(見圖2.21(c))。

其等效傳遞函數(shù)為

④將G1(s)、G23(s)和H1(s)組成的反饋回路簡化便求得系統(tǒng)的傳遞函數(shù)為

圖

2.21系統(tǒng)結構圖簡化

應當指出，在結構圖簡化過程中，兩個相鄰的相加點和分支點不能輕易交換?？傊鶕?jù)實際系統(tǒng)中各環(huán)節(jié)(子系統(tǒng))的結構圖和信息流向，可建立系統(tǒng)的結構圖。在確定系統(tǒng)的輸入量和輸出量后，經(jīng)過對系統(tǒng)結構圖的簡化和運算，就能求出系統(tǒng)的傳遞函數(shù)。

這是經(jīng)典控制理論中利用傳遞函數(shù)來建立系統(tǒng)數(shù)學模型的基本方法。

圖

2.22信號流圖2.4信

號

流

圖

1.信號流圖中的術語下面結合圖2.23介紹信號流圖的有關術語。輸入節(jié)點(源)

僅具有輸出支路的節(jié)點。如圖2.23中的x1。

圖

2.23信號流圖

前向通道回路不接觸回路輸出節(jié)點(阱)

僅有輸入支路的節(jié)點。有時信號流圖中沒有一個節(jié)點是僅具有輸入支路的。我們只要定義信號流圖中任一變量為輸出變量，然后從該節(jié)點變量引出一條增益為1的支路，即可形成一輸出節(jié)點，如圖2.23中的x6。

混合節(jié)點既有輸入支路又有輸出支路的節(jié)點。如圖2.23中的x2，x3，x4，x5。

通道沿支路箭頭方向而穿過各相連支路的途徑。如果通道與任一節(jié)點相交不多于一次，就叫開通道。如果通道的終點就是起點，并且與任何其他節(jié)點相交不多于一次，就稱作閉通道。

前向通道如果從輸入節(jié)點(源)到輸出節(jié)點(阱)的通道上，通過任何節(jié)點不多于一次，則該通道叫前向通道。如

前向通道增益前向通道上各支路增益之乘積，用Pk表示?；芈菲瘘c和終點在同一節(jié)點，并與其它節(jié)點相遇僅一次的通路，也就是閉合通道。以下是圖2.23中的一些回路：

回路增益回路中所有支路的乘積，用La表示。

不接觸回路回路之間沒有公共節(jié)點時，這種回路叫做不接觸回路。在信號流圖中，可以有兩個或兩個以上的不接觸回路。

例如：

就是不接觸回路的例子。上述定義可以類推到系統(tǒng)的結構圖中，從而采用梅遜公式(后面將介紹)求取由結構圖表示的系統(tǒng)的閉環(huán)傳遞函數(shù)。

2.信號流圖的性質(zhì)

(1)信號流圖適用于線性系統(tǒng)。

(2)支路表示一個信號對另一個信號的函數(shù)關系，信號只能沿支路上的箭頭指向傳遞。

(3)在節(jié)點上可以把所有輸入支路的信號疊加，并把相加后的信號送到所有的輸出支路。

(4)具有輸入和輸出支路的混合節(jié)點，通過增加一個具有單位增益的支路可以把它作為輸出節(jié)點來處理。

(5)對于一個給定的系統(tǒng)，信號流圖不是唯一的。由于描述同一個系統(tǒng)的方程可以表示為不同的形式，因此可以畫出不同的信號流程圖。信號流程圖的化簡串聯(lián)支路的總傳輸?shù)扔诟髦穫鬏斨e。并聯(lián)支路的總傳輸?shù)扔诟髦穫鬏斨??；旌瞎?jié)點可以通過移動支路的方法消去?；丨h(huán)可以根據(jù)反饋連接的規(guī)則化為等效支路。用梅遜公式可以直接求信號流圖從輸入節(jié)點到輸出節(jié)點的增益，其表達式為

(2.45)

P——系統(tǒng)總增益(對于控制系統(tǒng)的結構圖而言，就是輸入到輸出的傳遞函數(shù))；

k——前向通道數(shù)目；

Pk——第k條前向通道的增益；

3.梅遜公式☆

Δ——信號流圖的特征式，它是信號流圖所表示的方程組的系數(shù)矩陣的行列式。在同一個信號流圖中不論求圖中任何一對節(jié)點之間的增益，其分母總是Δ，變化的只是其分子。它可以通過下面的表達式計算：其中，∑L(1)——所有不同回路增益乘積之和；∑L(2)——所有任意兩個互不接觸回路增益乘積之和；∑L(3)——所有任意三個互不接觸回路增益乘積之和；∑L(m)——所有任意m個不接觸回路增益乘積之和；

Δk

——信號流圖中除去與第k條前向通道Pk相接觸的支路和節(jié)點后余下的信號流圖的特征式，稱為Pk的余因式。

Pk

——第k條前向通道的增益；例2.10根據(jù)梅遜公式求下圖的信號流圖總輸出

解：此系統(tǒng)有六個回環(huán)，即ab、cd、ef、gh、ij和kfdb，因此：∑L1=ab+cd+

ef+gh+ij+kfdb兩個互不接觸的回環(huán)有七種組合，即abef、abgh、abij、cdgh、cdij、efij及kfdbij，有∑L2=abef+abgh+abij+cdgh+cdij+efij+kfdbij三個互不接觸的回環(huán)有ab、ef和ij，有∑L3=abefij特征式：△=1-∑L1+∑L2-∑L3

從源節(jié)點到匯節(jié)點有兩條前向通道。一條為acegi，它與所有的回環(huán)均有接觸，因此：

P1=acegi，△1=1另一條前向通道為kgi它不與回環(huán)cd接觸，所以

P2=kgi

，△2=1-cd將以上結果代入得到總傳輸T=（P1△1+P2△2

）/（1-∑L1+∑L2-∑L3）

=[acegi+kgi（1-cd）]/[1-（ab+cd+

ef+gh+ij+kfdb）+（abef+abgh+abij+cdgh+cdij+efij+kfdbij）-abefij]

例2.11系統(tǒng)的方塊圖如圖2.24所示，試用梅遜公式求系統(tǒng)的傳遞函數(shù)C(s)/R(s)。

圖2.24系統(tǒng)的方塊圖有三個獨立回路：

解

從圖中可以看出，該框圖只有一個前向通道，其增益為

沒有兩個及兩個以上的互相獨立回路。所以，特征式Δ為連接輸入節(jié)點和輸出節(jié)點的前向通道的余因式Δ1，可以通過除去與該通道接觸的回路的方法得到。因為通道P1與三個回路都接觸，所以有

Δ1=1

因此，輸入量R(s)和輸出量C(s)之間的總增益或閉環(huán)傳遞函數(shù)為

例2.11

利用梅遜公式確定圖2.8(c)所表示的系統(tǒng)的傳遞函數(shù)Φ(s)=U2(s)/U1(s)及ΦE(s)=E(s)/U1(s)。(自己練習)

解

該圖有三個反饋回路：回路1和回路3不接觸，所以有

以U2(s)作為輸出信號時，該系統(tǒng)只有一條前向通道。

且有

該前向通道與各回路都有接觸，所以

Δ1=1故

以E(s)作為輸出信號時，該系統(tǒng)也只有一條前向通道。

且

P1=1這條前向通道與回路1相接觸，故

所以

總之，當求解系統(tǒng)的傳遞函數(shù)時，簡單的系統(tǒng)可以直接用結構圖運算，既清楚又方便；復雜的系統(tǒng)可以將其看作信號流圖后，再利用梅遜公式計算。

需要強調(diào)的是，在利用梅遜公式時，要考慮周到，不能遺漏任何應當計算的回路和前向通道。

2.5線性系統(tǒng)的狀態(tài)空間描述（了解）

2.5.1狀態(tài)空間描述的基本概念

狀態(tài)指系統(tǒng)的運動狀態(tài)。設想有一個質(zhì)點作直線運動，這個系統(tǒng)的狀態(tài)就是質(zhì)點每一個時刻的位置和速度。

狀態(tài)變量指足以完全表征系統(tǒng)運動狀態(tài)的最小個數(shù)的一組變量。若知道這些變量在任何初始時刻t0的值和t≥t0時系統(tǒng)所加的輸入函數(shù)，便可完全確定在任何t＞t0時刻的狀態(tài)。

一個用n階微分方程描述的系統(tǒng)，有n個獨立變量，當這n個獨立變量的時間響應都求得時，系統(tǒng)的運動狀態(tài)也就完全被揭示了。因此可以說，系統(tǒng)的狀態(tài)變量就是n階系統(tǒng)的n個獨立變量。需要指出，對同一個系統(tǒng)，選取哪些變量作為狀態(tài)變量并不是唯一的，但這些變量必須是互相獨立的，且個數(shù)等于微分方程的階數(shù)。對于一般物理系統(tǒng)，微分方程的階數(shù)唯一地取決于系統(tǒng)中獨立儲能元件的個數(shù)。因此，系統(tǒng)狀態(tài)變量的個數(shù)又可以說等于系統(tǒng)中獨立儲能元件的個數(shù)。

狀態(tài)向量如果n個狀態(tài)變量用x1(t)，

x2(t)，…，xn(t)表示，并把這些狀態(tài)變量看作是向量x(t)的分量，則向量x(t)稱為狀態(tài)向量。記為

或

狀態(tài)空間以狀態(tài)變量x1(t)，x2(t)，…，xn(t)為坐標軸構成的n維空間。系統(tǒng)在任意時刻的狀態(tài)x(t)都可用狀態(tài)空間中的一個點來表示。已知初始時刻t0的狀態(tài)x(t0)，可得到狀態(tài)空間中的一個初始點。隨著時間的推移，x(t)將在狀態(tài)空間中描繪出一條軌跡，稱為狀態(tài)軌跡。

狀態(tài)方程描述系統(tǒng)的狀態(tài)變量與系統(tǒng)輸入量之間關系的一階微分方程組，稱為系統(tǒng)的狀態(tài)方程。

輸出方程描述系統(tǒng)輸出量與狀態(tài)變量間的函數(shù)關系式，稱為系統(tǒng)的輸出方程。

狀態(tài)空間表達式狀態(tài)方程與輸出方程總合起來，就構成對一個系統(tǒng)動態(tài)的完整描述，稱之為狀態(tài)空間表達式。

通常，對于單變量系統(tǒng)，狀態(tài)方程習慣寫成如下形式：(2.46)輸出方程為

(2.47)寫成矩陣向量形式為

(2.48)式中=x［x1，x2，…，xn］T表示n維狀態(tài)向量；A、B、C、d分別表示系統(tǒng)內(nèi)部狀態(tài)的系數(shù)矩陣(系統(tǒng)矩陣)、輸入對狀態(tài)作用的輸入矩陣、輸出與狀態(tài)關系的輸出矩陣、直接聯(lián)系輸入量與輸出量的直接傳遞函數(shù)(或稱前饋系數(shù))。

推廣到p輸入、q輸出的系統(tǒng)，其狀態(tài)空間表達式為

(i=1，2，…，n)(2.49)(j=1，2，…，q)(2.50)寫成矩陣向量形式為

(2.51)式中x和A同單變量系統(tǒng)。

表示p維輸入向量；表示輸入矩陣；表示q維輸出向量；表示輸出矩陣；表示直接傳遞函數(shù)矩陣。

上述系統(tǒng)可簡稱為系統(tǒng)(A，B，C，D)。

圖

2.25狀態(tài)空間表達式的框圖

例2.12設有一個RLC電路(如圖2.26所示)，試確定它的狀態(tài)空間表達式。圖2.26RLC電路

例2.12設有一個RLC電路(如圖2.26所示)，試確定它的狀態(tài)空間表達式。

解根據(jù)電學原理，可得

(2.52)選取u2(t)和i(t)作為系統(tǒng)的狀態(tài)變量，可把式(2.52)表示為如下兩個一階微分方程：(2.53)取狀態(tài)變量x1=u2(t)，x2=i(t)，輸出y=u2=x1，則系統(tǒng)的狀態(tài)空間表達式為

(2.54)(2.55)寫成矩陣向量形式為

(2.56)(2.57)圖2.3機械阻尼器示意圖k是彈簧系數(shù)，m是運動部件質(zhì)量，μ是阻尼器的阻尼系數(shù)；外力f(t)是系統(tǒng)的輸入量，位移y(t)是系統(tǒng)的輸出量。試建立其狀態(tài)空間模型。

例2.13

對圖2.3所示的機械阻尼系統(tǒng)，試建立其狀態(tài)空間模型。

解根據(jù)系統(tǒng)的運動方程式(2.6)，選取狀態(tài)變量x1(t)=y(t)、x2(t)=y(t)，令輸入為系統(tǒng)所受的外力u(t)=f(t)，輸出為質(zhì)量塊的位移y(t)，經(jīng)過簡單的推導可知系統(tǒng)的狀態(tài)空間表達式為

.(2.58)(2.59)可將式(2.58)和式(2.59)寫成標準形式

式中，2.5.2狀態(tài)空間表達式與傳遞函數(shù)的關系

同一系統(tǒng)的兩種不同模型(傳遞函數(shù)和狀態(tài)空間表達式)之間存在著內(nèi)在的聯(lián)系，并且可以互相轉(zhuǎn)化。以下是對單輸入、單輸出系統(tǒng)的討論。(2.60)該系統(tǒng)在狀態(tài)空間可表示為

設要研究的系統(tǒng)的傳遞函數(shù)為(2.61)式中為狀態(tài)向量，u，y分別為輸入量和輸出量。在零初始條件假設下，方程

所以有

其中I為單位矩陣。用(sI-A).-1乘上式兩邊，有

(2.62)(2.63)(2.64)因此有

根據(jù)傳遞函數(shù)的定義可知，系統(tǒng)的傳遞函數(shù)與狀態(tài)空間描述之間的關系為

對于單輸入單輸出系統(tǒng)，采用式(2.65)可以求出系統(tǒng)的傳遞函數(shù)，對于多輸入多輸出系統(tǒng)，用式(2.65)求出的是用于描述多變量系統(tǒng)輸入輸出關系的傳遞函數(shù)矩陣。

(2.65)

例2.14

試根據(jù)例2.13中系統(tǒng)的狀態(tài)空間描述求系統(tǒng)的傳遞函數(shù)。解將例2.13中系統(tǒng)狀態(tài)空間表達式的系數(shù)矩陣(A，B，C，D)代入方程(2.65)可得因為

所以系統(tǒng)的傳遞函數(shù)為

2.5.3狀態(tài)空間表達式的建立

線性定常系統(tǒng)的狀態(tài)空間表達式也可以由系統(tǒng)的微分方程或傳遞函數(shù)來建立。(2.66)

情形一線性微分方程中不含輸入的導數(shù)項，傳遞函數(shù)沒有零點?？紤]下列n階系統(tǒng)：其對應的傳遞函數(shù)為

(2.67)如果已知y(0)，y(0)，…，y(n.1)(0)和t≥0時的輸入量u(t)，就可以完全確定系統(tǒng)未來的行為。因此，可以將y(t)，y(t)，…，y(n.1)(t)作為系統(tǒng)的狀態(tài)變量，即

..則方程(2.66)可表示為下列微分方程組

即

(2.68)式中

輸出量可由下式確定

即

(2.69)式中，所以，當線性微分方程中不含輸入的導數(shù)項，傳遞函數(shù)沒有零點時，系統(tǒng)的一種狀態(tài)空間表達式可由方程(2.68)和(2.69)給出。

在這種情況下，系統(tǒng)的微分方程和相應的傳遞函數(shù)為

(2.70)(2.71)此時，如果再取y(t)，y(t)，…，y(n-1)(t)作為系統(tǒng)的狀態(tài)變量，則構成的n個一階微分方程不能唯一確定系統(tǒng)的狀態(tài)，其主要原因是

.情形二線性微分方程含有輸入的導數(shù)項，傳遞函數(shù)有零點。中含有輸入的導數(shù)項。因此，狀態(tài)變量的設置，必須能消除狀態(tài)方程中u的導數(shù)項。其中一種選取方法是構造如下n個狀態(tài)變量(還有其他的選取方法)：從而，可以得到由如下這組狀態(tài)變量描述的系統(tǒng)的狀態(tài)方程：

寫成向量.矩陣形式，則系統(tǒng)的狀態(tài)方程和輸出方程為

式中，例

2.15

已知某控制系統(tǒng)的運動方程為

其中，u(t)和y(t)分別為系統(tǒng)的輸入和輸出。選取狀態(tài)變量x1=y，x2=，寫出系統(tǒng)的狀態(tài)空間表達式。

解

由已知條件容易得到

將其寫成向量矩陣形式，并考慮到x1=y，可得系統(tǒng)狀態(tài)空間表達式為

例

2.16

已知控制系統(tǒng)的微分方程為

試寫出其狀態(tài)空間表達式。

解

這是一個含輸入導數(shù)項的運動方程，由已知條件可知

狀態(tài)變量按照情形二選取，則

于是得到系統(tǒng)的狀態(tài)空間表達式為

2.6線性定常系統(tǒng)數(shù)學模型的MATLAB實現(xiàn)

1.MATLAB建立系統(tǒng)數(shù)學模型的方法解輸入上述傳遞函數(shù)的MATLAB程序如下：

%ex.2.17num=[12,24,12,20];%分母den=[2,4,6,2,2];%分子G=tf(num,den)%傳遞函數(shù)例2.17

若給定系統(tǒng)的傳遞函數(shù)為試用MATLAB語句表示該傳遞函數(shù)。程序第一行是注釋語句，不執(zhí)行；第二、三行分別按降冪順序輸入給定傳遞函數(shù)的分子和分母多項式的系數(shù)；第四行建立系統(tǒng)的傳遞函數(shù)模型。運行結果顯示