狀態(tài)空間分析

關(guān)于狀態(tài)空間分析第一頁(yè)，共六十一頁(yè)，2022年，8月28日

由系統(tǒng)機(jī)理建立狀態(tài)空間表達(dá)式例9-2

力、彈簧、阻尼器組成的機(jī)械系統(tǒng)由牛頓定律其中，輸入為F(t),輸出為位移y(t).若已知初始位移和初始速度，則可唯一確定系統(tǒng)在輸入作用下的解。第二頁(yè)，共六十一頁(yè)，2022年，8月28日可選擇位移和速度作為狀態(tài)變量狀態(tài)方程為輸入系統(tǒng)的狀態(tài)空間表達(dá)式為第三頁(yè)，共六十一頁(yè)，2022年，8月28日例9-3

機(jī)械系統(tǒng)如圖所示，若不考慮重力的作用，試寫出以拉力F為輸入，質(zhì)量m1和m2的位移y1和y2為輸出的狀態(tài)空間表達(dá)式。解根據(jù)牛頓定律，分別列寫關(guān)于m1和m2的方程選擇4個(gè)互相獨(dú)立的狀態(tài)變量。第四頁(yè)，共六十一頁(yè)，2022年，8月28日第五頁(yè)，共六十一頁(yè)，2022年，8月28日9.3.3由微分方程建立狀態(tài)空間描述步驟：直接根據(jù)系統(tǒng)的物理機(jī)理建立相應(yīng)的微分(連續(xù)系統(tǒng)）或差分（離散系統(tǒng)）方程組。針對(duì)微分方程，定義一組狀態(tài)變量，建立狀態(tài)方程，并根據(jù)系統(tǒng)輸出和狀態(tài)之間的關(guān)系，建立系統(tǒng)的輸出方程。狀態(tài)變量的選取

1.狀態(tài)變量的選取是非唯一的。

2.選取方法

（1）可選取初始條件對(duì)應(yīng)的變量或與其相關(guān)的變量作為系統(tǒng)的狀態(tài)變量。

（2）可選取獨(dú)立儲(chǔ)能（或儲(chǔ)信息）元件的特征變量或與其相關(guān)的變量作為控制系統(tǒng)的狀態(tài)變量。（如電感電流i、電容電壓uc

、質(zhì)量m

的速度v

等）第六頁(yè)，共六十一頁(yè)，2022年，8月28日(1)

線性微分方程中不含有輸入函數(shù)導(dǎo)數(shù)項(xiàng)的系統(tǒng)的狀態(tài)空間表達(dá)式

設(shè)單輸入/單輸出的控制系統(tǒng)的動(dòng)態(tài)過(guò)程由下列n階微分方程來(lái)描述輸出信號(hào)及其各階導(dǎo)數(shù)u系統(tǒng)的輸入信號(hào)若已知初始條件及時(shí)刻的輸入信號(hào)u(t)，則系統(tǒng)在任何時(shí)刻的行為便可完全確定。第七頁(yè)，共六十一頁(yè)，2022年，8月28日狀態(tài)變量可取為可寫為第八頁(yè)，共六十一頁(yè)，2022年，8月28日狀態(tài)空間表達(dá)式C=[100…0]其中，各矩陣為,A=,B=第九頁(yè)，共六十一頁(yè)，2022年，8月28日由微分方程求狀態(tài)空間表達(dá)式的例題見書例9－3b根據(jù)上式繪制的狀態(tài)變量之間關(guān)系的方塊圖,也稱狀態(tài)變量圖,如圖所示，每個(gè)積分器的輸出都是對(duì)應(yīng)的一個(gè)狀態(tài)變量，狀態(tài)方程由積分器的輸入輸出關(guān)系確定，輸出方程在輸出端給出

：第十頁(yè)，共六十一頁(yè)，2022年，8月28日例9-4

設(shè)一控制系統(tǒng)的動(dòng)態(tài)過(guò)程用微分方程表示為u，y分別為系統(tǒng)的輸入和輸出信號(hào)，試求系統(tǒng)的狀態(tài)空間描述。解：選取狀態(tài)變量為寫成標(biāo)準(zhǔn)形式第十一頁(yè)，共六十一頁(yè)，2022年，8月28日(2)輸入信號(hào)包含導(dǎo)數(shù)項(xiàng)的n階微分方程系統(tǒng)的狀態(tài)空間描述控制系統(tǒng)由下列n階微分方程來(lái)描述若采用上述方法選取狀態(tài)變量，則可化成一階微分方程組方程組中，出現(xiàn)了u(t)的導(dǎo)數(shù)項(xiàng)，不能采用以上方法。第十二頁(yè)，共六十一頁(yè)，2022年，8月28日選擇狀態(tài)變量的原則：在包含狀態(tài)變量的n個(gè)一階微分方程構(gòu)成的系統(tǒng)狀態(tài)方程中，任何一個(gè)微分方程都不能含有作用函數(shù)/輸入量的導(dǎo)數(shù)項(xiàng)。

若輸入信號(hào)在時(shí)刻出現(xiàn)有限跳躍，如階躍函數(shù)，則在t0時(shí)刻出現(xiàn)脈沖函數(shù)，在t0時(shí)刻將是更高階的脈沖函數(shù)，狀態(tài)軌跡則在t0時(shí)刻產(chǎn)生無(wú)窮大跳躍。故此時(shí)不能按照前面的方法將系統(tǒng)的輸出量及其各階導(dǎo)數(shù)直接選作系統(tǒng)的狀態(tài)變量，因?yàn)檫@組狀態(tài)變量已經(jīng)不具備在已知系統(tǒng)輸入和初始狀態(tài)條件下完全確定系統(tǒng)未來(lái)運(yùn)動(dòng)狀態(tài)的特性。第十三頁(yè)，共六十一頁(yè)，2022年，8月28日選取狀態(tài)變量由此可得輸出方程n-1個(gè)狀態(tài)方程可寫成第十四頁(yè)，共六十一頁(yè)，2022年，8月28日對(duì)xn求導(dǎo)將y(n)代入第十五頁(yè)，共六十一頁(yè)，2022年，8月28日將上式中所有的輸出項(xiàng)以及輸出的導(dǎo)數(shù)項(xiàng)都用狀態(tài)和輸入的各階導(dǎo)數(shù)項(xiàng)表示有可寫成可以保證系統(tǒng)有唯一解。第十六頁(yè)，共六十一頁(yè)，2022年，8月28日

A=

B=C=[100…0]，將微分方程組改寫成矩陣向量形式第十七頁(yè)，共六十一頁(yè)，2022年，8月28日例9-5

設(shè)一控制系統(tǒng)的動(dòng)態(tài)方程用微分方程表示為

試求該控制系統(tǒng)的狀態(tài)空間描述。解：與下式比較可得選取狀態(tài)變量第十八頁(yè)，共六十一頁(yè)，2022年，8月28日狀態(tài)方程為系統(tǒng)狀態(tài)空間表達(dá)式為第十九頁(yè)，共六十一頁(yè)，2022年，8月28日一般形式：當(dāng)式中b0=0時(shí)，還可以按如下規(guī)則選擇另一組狀態(tài)變量。設(shè)第二十頁(yè)，共六十一頁(yè)，2022年，8月28日則得到輸出方程為第二十一頁(yè)，共六十一頁(yè)，2022年，8月28日繼續(xù)變換，則得到第二十二頁(yè)，共六十一頁(yè)，2022年，8月28日對(duì)x1求導(dǎo)，并將y

(n)用代入后整理得第二十三頁(yè)，共六十一頁(yè)，2022年，8月28日d=0第二十四頁(yè)，共六十一頁(yè)，2022年，8月28日例9-5（續(xù)）

設(shè)一控制系統(tǒng)的動(dòng)態(tài)方程用微分方程表示為

試求該控制系統(tǒng)的狀態(tài)空間描述。解：與下式比較可得選取狀態(tài)變量第二十五頁(yè)，共六十一頁(yè)，2022年，8月28日所以狀態(tài)空間表達(dá)式為對(duì)于一個(gè)給定的系統(tǒng)而言，狀態(tài)變量的選取并不是唯一的。第二十六頁(yè)，共六十一頁(yè)，2022年，8月28日例9-6

設(shè)有雙輸入－雙輸出的二階系統(tǒng)，其運(yùn)動(dòng)方程為

試寫出系統(tǒng)的狀態(tài)空間表達(dá)式解：按最高階導(dǎo)數(shù)項(xiàng)求解y1、y2的方法，對(duì)以上兩方程積分第二十七頁(yè)，共六十一頁(yè)，2022年，8月28日選取狀態(tài)變量由y1的表達(dá)式可得選取狀態(tài)變量則有由y2的表達(dá)式，有綜上第二十八頁(yè)，共六十一頁(yè)，2022年，8月28日將上述方程寫成矩陣方程形式，有系統(tǒng)輸出矩陣方程為第二十九頁(yè)，共六十一頁(yè)，2022年，8月28日

由狀態(tài)變量圖求狀態(tài)空間描述狀態(tài)變量圖描述系統(tǒng)狀態(tài)變量之間關(guān)系的圖，由積分環(huán)節(jié)、比例環(huán)節(jié)和相加符號(hào)組成。狀態(tài)變量圖的特點(diǎn)：每一個(gè)積分環(huán)節(jié)的輸出都代表系統(tǒng)的一個(gè)狀態(tài)變量第三十頁(yè)，共六十一頁(yè)，2022年，8月28日P276,例9-7：方框圖（結(jié)構(gòu)圖）狀態(tài)變量圖步驟：結(jié)構(gòu)圖由典型環(huán)節(jié)構(gòu)成；將一階慣性環(huán)節(jié)、二階震蕩環(huán)節(jié)分解成局部狀態(tài)圖；設(shè)定每一個(gè)積分環(huán)節(jié)的輸出為一個(gè)狀態(tài)變量，畫出狀態(tài)變量圖。第三十一頁(yè)，共六十一頁(yè)，2022年，8月28日例9-7

系統(tǒng)的閉環(huán)傳遞函數(shù)為試?yán)L制系統(tǒng)的狀態(tài)變量圖，并由圖列寫系統(tǒng)的狀態(tài)空間描述。解將系統(tǒng)的閉環(huán)傳遞函數(shù)改寫成令第三十二頁(yè)，共六十一頁(yè)，2022年，8月28日第三十三頁(yè)，共六十一頁(yè)，2022年，8月28日作業(yè)：P3489-3，9-6，9-9第三十四頁(yè)，共六十一頁(yè)，2022年，8月28日將傳遞函數(shù)轉(zhuǎn)換成狀態(tài)空間描述將狀態(tài)空間描述轉(zhuǎn)換成傳遞函數(shù)系統(tǒng)傳遞函數(shù)與狀態(tài)空間描述（1）傳遞函數(shù)轉(zhuǎn)換成狀態(tài)空間描述第三十五頁(yè)，共六十一頁(yè)，2022年，8月28日選取狀態(tài)變量第三十六頁(yè)，共六十一頁(yè)，2022年，8月28日若有，則輸出方程的前兩項(xiàng)為零。此法也可用來(lái)求取微分方程中含有輸入信號(hào)的導(dǎo)數(shù)時(shí)系統(tǒng)的狀態(tài)空間描述，注意與前面介紹的方法不同的是B和C不同。第三十七頁(yè)，共六十一頁(yè)，2022年，8月28日例9-8

設(shè)控制系統(tǒng)的傳遞函數(shù)為

試求該統(tǒng)的狀態(tài)空間描述。解狀態(tài)方程為輸出方程為第三十八頁(yè)，共六十一頁(yè)，2022年，8月28日

傳遞函數(shù)的串聯(lián)實(shí)現(xiàn)傳遞函數(shù)為兩多項(xiàng)式相除形式，分子多項(xiàng)式（Numerator）為分母多項(xiàng)式（Denominator）如果為G(s)的m個(gè)零點(diǎn)，為G（s）的n個(gè)極點(diǎn)，那么G(s)可以表示為：第三十九頁(yè)，共六十一頁(yè)，2022年，8月28日所以系統(tǒng)的實(shí)現(xiàn)可以由共n個(gè)環(huán)節(jié)串聯(lián)而成，如圖（a）所示。其結(jié)構(gòu)圖可以是如圖（b）中虛框表示。(a)對(duì)第一個(gè)環(huán)節(jié)，由于：第四十頁(yè)，共六十一頁(yè)，2022年，8月28日(b)m=n-1的情況我們令各個(gè)積分器的輸出為系統(tǒng)狀態(tài)變量，則得系統(tǒng)狀態(tài)方程為：第四十一頁(yè)，共六十一頁(yè)，2022年，8月28日第四十二頁(yè)，共六十一頁(yè)，2022年，8月28日傳遞函數(shù)的并聯(lián)實(shí)現(xiàn)系統(tǒng)傳遞函數(shù)為系統(tǒng)的特征方程。當(dāng)Den(s)=0有n個(gè)不等的特征根（）G(s)可以分解為n個(gè)分式之和，即：，稱作系統(tǒng)對(duì)應(yīng)極點(diǎn)pi的留數(shù)。

其中其中，第四十三頁(yè)，共六十一頁(yè)，2022年，8月28日上式可以用如圖所示的并聯(lián)方式實(shí)現(xiàn)。(a)(b)并聯(lián)實(shí)現(xiàn)(無(wú)重根)第四十四頁(yè)，共六十一頁(yè)，2022年，8月28日從圖（b）我們可得系統(tǒng)的狀態(tài)方程：輸出方程為：寫成矢量形式為：請(qǐng)注意，這里的系統(tǒng)矩陣A為一標(biāo)準(zhǔn)的對(duì)角型。第四十五頁(yè)，共六十一頁(yè)，2022年，8月28日

如果狀態(tài)變量選擇為第四十六頁(yè)，共六十一頁(yè)，2022年，8月28日那么系統(tǒng)輸出則為同樣，經(jīng)過(guò)反拉氏變換并展成矩陣形式有對(duì)角陣標(biāo)準(zhǔn)型

第四十七頁(yè)，共六十一頁(yè)，2022年，8月28日當(dāng)上述G(s)的分母Den(s)=0有重根時(shí)，不失一般性，假設(shè)：即為q重根，其它為單根。這時(shí)G(S)可以分解為：其中：

i=1,2,…,qj=q+1,q+2,…,n第四十八頁(yè)，共六十一頁(yè)，2022年，8月28日并聯(lián)實(shí)現(xiàn)（有重根）第四十九頁(yè)，共六十一頁(yè)，2022年，8月28日取圖中每個(gè)積分器輸出為狀態(tài)變量，則有：第五十頁(yè)，共六十一頁(yè)，2022年，8月28日注意這里的A為一約當(dāng)標(biāo)準(zhǔn)型。第五十一頁(yè)，共六十一頁(yè)，2022年，8月28日例

求下列傳遞函數(shù)的并聯(lián)實(shí)現(xiàn)解：分母各項(xiàng)多項(xiàng)式分解可得系統(tǒng)并聯(lián)實(shí)現(xiàn)的動(dòng)態(tài)方程為：第五十二頁(yè)，共六十一頁(yè)，2022年，8月28日系統(tǒng)動(dòng)態(tài)方程和系統(tǒng)傳遞函數(shù)（陣）都是控制系統(tǒng)兩種經(jīng)常使用的數(shù)學(xué)模型。動(dòng)態(tài)方程不但體現(xiàn)了系統(tǒng)輸入輸出的關(guān)系，而且還清楚地表達(dá)了系統(tǒng)內(nèi)部狀態(tài)變量的關(guān)系。傳遞函數(shù)只體現(xiàn)了系統(tǒng)輸入與輸出的關(guān)系。從傳遞函數(shù)到動(dòng)態(tài)方程是個(gè)系統(tǒng)實(shí)現(xiàn)的問(wèn)題，這是一個(gè)比較復(fù)雜的并且是非唯一的過(guò)程。但從動(dòng)態(tài)方程到傳遞函數(shù)（陣）卻是一個(gè)唯一的、比較簡(jiǎn)單的過(guò)程。(2)將狀態(tài)空間描述轉(zhuǎn)換成傳遞函數(shù)第五十三頁(yè)，共六十一頁(yè)，2022年，8月28日單輸入/單輸出系統(tǒng)的狀態(tài)空間描述轉(zhuǎn)換成傳遞函數(shù)設(shè)一單輸入/單輸出系統(tǒng)的狀態(tài)空間描述為D為標(biāo)量設(shè)初始條件為零，對(duì)上式進(jìn)行拉氏變換可得系統(tǒng)的傳遞函數(shù)為第五十四頁(yè)，共六十一頁(yè)，2022年，8月28日例9-9

設(shè)系統(tǒng)的狀態(tài)空間描述為

求系統(tǒng)的傳遞函數(shù)。解狀態(tài)空間描述的[A,B,C]分別為第五十五頁(yè)，共六十一頁(yè)，2022年，8月28日系統(tǒng)的傳遞函數(shù)為第五十六頁(yè)，共六十一頁(yè)，2022年，8月28日

多輸入/多輸出系統(tǒng)的狀態(tài)空間描述轉(zhuǎn)換成傳遞函數(shù)矩陣

系統(tǒng)的輸出信號(hào)

系統(tǒng)的輸入信號(hào)

系統(tǒng)的狀態(tài)變量第五十七頁(yè)，共六十一頁(yè)，2022年，8月28日與單輸入/單輸出系統(tǒng)的狀態(tài)方程形式相同，僅是矩陣B,C,D的維數(shù)不同。將狀態(tài)方程進(jìn)行拉氏變換，得－－－傳