1、第二章 控制系統(tǒng)的狀態(tài)空間分析在討論了狀態(tài)方程的描述、標(biāo)準(zhǔn)形和模型轉(zhuǎn)換后，本章將討論線性多變量系統(tǒng)的運動分析，包括線性狀態(tài)方程的求解和系統(tǒng)的能控與能觀性分析。2.1 定常系統(tǒng)狀態(tài)方程的解 線性系統(tǒng)的解 給定線性定常系統(tǒng)非齊次狀態(tài)方程為：(2.1.1）其中，且初始條件為。 將方程（2.1.1）寫為 在上式兩邊左乘e-At，可得 將上式由O積分到t，得故可求出其解為(2.1.2a)或(2.1.2b)式中為系統(tǒng)的狀態(tài)轉(zhuǎn)移矩陣。對于線性時變系統(tǒng)非齊次狀態(tài)方程， (2.1.3）類似可求出其解為(2.1.4)一般說來，線性時變系統(tǒng)的狀態(tài)轉(zhuǎn)移矩陣只能表示成一個無窮項之和，只有在特殊情況下，才能寫成矩陣指數(shù)函

2、數(shù)的形式。2.1.2 狀態(tài)轉(zhuǎn)移矩陣 1狀態(tài)轉(zhuǎn)移矩陣性質(zhì) 定義2.1 時變系統(tǒng)狀態(tài)轉(zhuǎn)移矩陣是滿足如下矩陣微分方程和初始條件(2.1.5)的解。下面不加證明地給出線性時變系統(tǒng)狀態(tài)轉(zhuǎn)移矩陣的幾個重要性質(zhì)：1、；2、；3、；4、當(dāng)A給定后， 唯一；5、計算時變系統(tǒng)狀態(tài)轉(zhuǎn)移矩陣的公式(2.1.6a) 上式一般不能寫成封閉形式，可按精度要求，用數(shù)值計算的方法取有限項近似。特別地，只有當(dāng)滿足即在矩陣乘法可交換的條件下，才可表示為如下矩陣指數(shù)函數(shù)形式(2.1.6b) 顯然，定常系統(tǒng)的狀態(tài)轉(zhuǎn)移矩陣不依賴于初始時刻，其性質(zhì)僅是上述時變系統(tǒng)的特例。 例21 試求如下線性定常系統(tǒng)的狀態(tài)轉(zhuǎn)移矩陣(t)和狀態(tài)轉(zhuǎn)移矩陣的逆

3、-1(t）。解 對于該系統(tǒng)，其狀態(tài)轉(zhuǎn)移矩陣由下式確定由于其逆矩陣為因此= 由于-1（t）=(-t)，故可求得狀態(tài)轉(zhuǎn)移矩陣的逆為 例2.2 求下列系統(tǒng)的時間響應(yīng)：式中，u(t)為t = 0時作用于系統(tǒng)的單位階躍函數(shù)，即u(t)=1(t)。解 對該系統(tǒng) 狀態(tài)轉(zhuǎn)移矩陣已在例2.1中求得，即因此，系統(tǒng)對單位階躍輸入的響應(yīng)為：或 如果初始狀態(tài)為零，即X(0)=0,可將X(t)簡化為2. 矩陣指數(shù)函數(shù)eAt的計算前已指出，狀態(tài)方程的解實質(zhì)上可歸結(jié)為計算狀態(tài)轉(zhuǎn)移矩陣，即矩陣指數(shù)函數(shù)eAt。如果給定矩陣A中所有元素的值，MATLAB將提供一種計算eAT的簡便方法，其中T為常數(shù)。 除了上述方法外，對eAt的計算

4、還有幾種分析方法可供使用。這里我們將介紹其中的四種計算方法。方法一：直接計算法(矩陣指數(shù)函數(shù))(2.1.7)可以證明，對所有常數(shù)矩陣A和有限的t值來說，這個無窮級數(shù)都是收斂的。方法二：對角線標(biāo)準(zhǔn)形與Jordan標(biāo)準(zhǔn)形法若可將矩陣A變換為對角線標(biāo)準(zhǔn)形，那么eAt可由下式給出式中，P是將A對角線化的非奇異線性變換矩陣。 類似地，若矩陣A可變換為Jordan標(biāo)準(zhǔn)形，則eAt可由下式確定出eAt = S e J t S 1(2.1.9) 例2.3 考慮如下矩陣A解 該矩陣的特征方程為因此，矩陣A有三個相重特征值=1。可以證明，矩陣A也將具有三重特征向量（即有兩個廣義特征向量）。易知，將矩陣A變換為Jo

5、rdan標(biāo)準(zhǔn)形的變換矩陣為 矩陣S的逆為于是注意到可得eAt = S e J t S 1即 方法三：拉氏變換法(2.1.10)為了求出eAt，關(guān)鍵是必須首先求出（sI-A）的逆。一般來說，當(dāng)系統(tǒng)矩陣A的階次較高時，可采用遞推算法。 例2.4 考慮如下矩陣A A試用前面介紹的兩種方法計算eAt 。解 方法一 由于A的特征值為0和-2（1=0，2= -2），故可求得所需的變換矩陣P為P =因此，由式(2.1.10）可得 方法二 由于可得因此方法四：化eAt為A的有限項法(Caley-Hamilton定理法)第四種是利用凱萊-哈密爾頓定理，化為A的有限項，然后通過求待定時間函數(shù)獲得的方法。凱萊-哈密

6、爾頓(Caley-Hamilton)定理在證明有關(guān)矩陣方程的定理或解決有關(guān)矩陣方程的問題時，凱萊-哈密爾頓定理是非常有用的。考慮n×n維矩陣A及其特征方程 凱萊-哈密爾頓定理指出，矩陣A滿足其自身的特征方程，即(2.1.11) 為了證明此定理，注意到（I-A）的伴隨矩陣adj(I-A)是的n -1次多項式，即式中，。由于可得 從上式可看出，A 和(i=1,2,，n)相乘的次序是可交換的。因此，如果(I-A）及其伴隨矩陣adj(I-A)中有一個為零，則其乘積為零。如果在上式中用A代替，顯然I-A為零。這樣即證明了凱萊-哈密爾頓定理。 最小多項式 按照凱萊-哈密爾頓定理，任一n×

7、;n維矩陣A滿足其自身的特征方程，然而特征方程不一定是A滿足的最小階次的純量方程。我們將矩陣A為其根的最小階次多項式稱為最小多項式，也就是說，定義n×n維矩陣A的最小多項式為最小階次的多項式()，即使得(A）= 0，或者 最小多項式在n×n維矩陣多項式的計算中起著重要作用。 假設(shè)的多項式d()是(I-A)的伴隨矩陣adj(I-A)的所有元素的最高公約式。可以證明，如果將d()的最高階次的系數(shù)選為1，則最小多項式()由下式給出：(2.1.12) 注意，n×n維矩陣A的最小多項式()可按下列步驟求出：1、根據(jù)伴隨矩陣adj(I-A)，寫出作為的因式分解多項式的adj(

8、I-A)的各元素；2、確定作為伴隨矩陣adj(I-A)各元素的最高公約式d()。選取d()的最高階次系數(shù)為1。如果不存在公約式，則d()=1；3、最小多項式()可由|I-A|除以d()得到。 設(shè)A的最小多項式階數(shù)為m。可以證明，采用賽爾維斯特內(nèi)插公式，通過求解行列式(2.1.13)即可求出。利用式(2.13）求解時，所得是以 (k=0,1,2,m-1)和 (i=1,2,3,，m)的形式表示的。此外，也可采用如下等價的方法。將式(2.13)按最后一行展開，容易得到(2.1.14）從而通過求解下列方程組：· ·(2.1.15)·可確定出(k=0,1,2，m-1)，進(jìn)而

9、代入式(2.14)即可求得。如果A為n×n維矩陣，且具有相異特征值，則所需確定的的個數(shù)為m=n，即有（2.1.16）如果A含有相重待征值，但其最小多項式有單根，則所需確定的的個數(shù)小于n.例2.5 考慮如下矩陣A試用化為A的有限項法計算。解 矩陣A的特征方程為可得相異特征值為1=0，2= -2。 由式(2.13），可得即 將上述行列式展開，可得或另一種可選用的方法是采用式(2.16）。首先，由確定待定時間函數(shù)和。由于1=0，2= -2，上述兩式變?yōu)?求解此方程組，可得因此，2.2能控性與能觀測性分析能控性和能觀測性的概念是由的。在用狀態(tài)空間法設(shè)計控制系統(tǒng)時，這兩個概念起到很重要的作用。

10、實際上，可控性和可觀性可以給出控制系統(tǒng)設(shè)計問題的完全解存在性的條件。如果所研究的系統(tǒng)不可控，那么控制系統(tǒng)設(shè)計問題的解不存在。可控和可觀測性決定了最優(yōu)控制問題解的存在性。例如，在極點配置問題中，狀態(tài)反饋的的存在性將由系統(tǒng)的能控性決定；在觀測器設(shè)計和最優(yōu)估計中，將涉及到系統(tǒng)的能觀測性條件。在本章中，我們的討論將限于線性系統(tǒng)。將首先給出能控性與能觀測性的定義，然后推導(dǎo)出判別系統(tǒng)能控和能觀測性的若干判據(jù)。 系統(tǒng)的能控性1 概述如果在一個有限的時間隔內(nèi)施加一個無約束的控制向量，使得系統(tǒng)由初始狀態(tài)x(to)轉(zhuǎn)移到任一狀態(tài)，則稱該系統(tǒng)在時刻to是能控的。如果系統(tǒng)的狀態(tài)x(to)在有限的時間間隔內(nèi)可由輸出的觀

11、測值確定，那么稱系統(tǒng)在時刻to是能觀測的。前已指出，在用狀態(tài)空間法設(shè)計控制系統(tǒng)時，這兩個概念起到非常重要的作用。實際上，雖然大多數(shù)物理系統(tǒng)是能控和能觀測的，然而其所對應(yīng)的數(shù)學(xué)模型可能不具有能控性和能觀測性。因此，必須了解系統(tǒng)在什么條件下是能控和能觀測的。上面給出了系統(tǒng)狀態(tài)能控與能觀測的定義，下面我們將首先推導(dǎo)狀態(tài)能控性的代數(shù)判據(jù)，然后給出狀態(tài)能控性的標(biāo)準(zhǔn)形判據(jù)。最后討論輸出能控性。2 定常系統(tǒng)狀態(tài)能控性的代數(shù)判據(jù)考慮線性連續(xù)時間系統(tǒng)：(2.2.1）其中，（單輸入），且初始條件為。如果施加一個無約束的控制信號，在有限的時間間隔tott1內(nèi)，使初始狀態(tài)轉(zhuǎn)移到任一終止?fàn)顟B(tài)，則稱由式（2.2.1）描述

12、的系統(tǒng)在t = to時為狀態(tài)(完全)能控的。如果每一個狀態(tài)都能控，則稱該系統(tǒng)為狀態(tài)(完全)能控的。下面我們將推導(dǎo)狀態(tài)能控的條件。不失一般性，設(shè)終止?fàn)顟B(tài)為狀態(tài)空間原點，并設(shè)初始時刻為零，即to=0。由上一章的內(nèi)容可知，式()的解為 利用狀態(tài)能控性的定義，可得或() 將寫為A的有限項的形式，即() 將式()代入式（2.2.2），可得()記則式()成為()如果系統(tǒng)是狀態(tài)能控的，那么給定任一初始狀態(tài)x(0)，都應(yīng)滿足式（2.2.5）。這就要求n×n維矩陣的秩為n。由此分析，可將狀態(tài)能控性的代數(shù)判據(jù)歸納為：當(dāng)且僅當(dāng)n×n維矩陣Q滿秩，即時，由式（2.2.1）確定的系統(tǒng)才是狀態(tài)能控的。

13、上述結(jié)論也可推廣到控制向量u為r維的情況。此時，如果系統(tǒng)的狀態(tài)方程為式中，那么可以證明，狀態(tài)能控性的條件為n×nr維矩陣的秩為n，或者說其中的n個列向量時線性無關(guān)的。通常，我們稱矩陣能控性矩陣。 例2.6 考慮由下式確定的系統(tǒng)： 由于即Q為奇異，所以該系統(tǒng)是狀態(tài)不能控的。 例2.7 考慮由下式確定的系統(tǒng)： 對于該情況，即Q為非奇異，因此系統(tǒng)是狀態(tài)能控的。3 狀態(tài)能控性條件的標(biāo)準(zhǔn)形判據(jù)關(guān)于定常系統(tǒng)能控性的判據(jù)很多。除了上述的代數(shù)判據(jù)外，本小節(jié)將給出一種相當(dāng)直觀的方法，這就是從標(biāo)準(zhǔn)形的角度給出的判據(jù)。 考慮如下的線性系統(tǒng)()式中，。 如果A的特征向量互不相同，則可找到一個非奇異線性變換矩

14、陣P，使得 注意，如果A的特征值相異，那么A的特征向量也互不相同；然而，反過來不成立。例如，具有相重特征值的n×n維實對稱矩陣也有可能有n個互不相同的特征向量。還應(yīng)注意，矩陣P的每一列是與i (i=1,2, ，n)有聯(lián)系的A的一個特征向量。 設(shè)x = P z() 將式()代入式()，可得 () 定義則可將式（2.2.8）重寫為 如果n×r維矩陣G 的任一行元素全為零，那么對應(yīng)的狀態(tài)變量就不能由任一來控制。由于狀態(tài)能控的條件是A的特征向量互異，因此當(dāng)且僅當(dāng)輸入矩陣沒有一行的所有元素均為零時，系統(tǒng)才是狀態(tài)能控的。在應(yīng)用狀態(tài)能控性的這一條件時，應(yīng)特別注意，必須將式（2.2.8）的

15、矩陣轉(zhuǎn)換成對角線形式。如果式（2.2.6）中的矩陣A不具有互異的特征向量，則不能將其化為對角線形式。在這種情況下，可將A化為Jordan標(biāo)準(zhǔn)形。例如，若A的特征值分別1,1,1,4,4,6,n,并且有n - 3個互異的特征向量，那么A的Jordan標(biāo)準(zhǔn)形為其中，在主對角線上的3×3和2×2子矩陣稱為Jordan塊。 假設(shè)能找到一個變換矩陣S，使得 如果利用x = S z()定義一個新的狀態(tài)向量z，將式（2.2.9）代入式（2.2.6）中，可得到()從而式（2.2.6）確定的系統(tǒng)的狀態(tài)能控性條件可表述為：當(dāng)且僅當(dāng)（1）式（2.2.10）中的矩陣J中沒有兩個Jordan塊與同一

16、特征值有關(guān)；（2）與每個Jordan塊最后一行相對應(yīng)的的任一行元素不全為零；（3）對應(yīng)于不同特征值的的每一行的元素不全為零時，則系統(tǒng)是狀態(tài)能控的。 例2.8 下列系統(tǒng)是狀態(tài)能控的： 下列系統(tǒng)是狀態(tài)不能控的： 4 用傳遞函數(shù)矩陣表達(dá)的狀態(tài)能控性條件狀態(tài)能控的條件也可用傳遞函數(shù)或傳遞矩陣描述。 狀態(tài)能控性的充要條件是在傳遞函數(shù)或傳遞函數(shù)矩陣中不出現(xiàn)相約現(xiàn)象。如果發(fā)生相約，那么在被約去的模態(tài)中，系統(tǒng)不能控。 例2.9 考慮下列傳遞函數(shù)：顯然，在此傳遞函數(shù)的分子和分母中存在可約的因子（s+2.5）(因此少了一階)。由于有相約因子，所以該系統(tǒng)狀態(tài)不能控。當(dāng)然，將該傳遞函數(shù)寫為狀態(tài)方程，可得到同樣的結(jié)論。

17、狀態(tài)方程為 由于即能控性矩陣的秩為1，所以可得到狀態(tài)不能控的同樣結(jié)論。5 輸出能控性在實際的控制系統(tǒng)設(shè)計中，需要控制的是輸出，而不是系統(tǒng)的狀態(tài)。對于控制系統(tǒng)的輸出，狀態(tài)能控性既不是必要的，也不是充分的。因此，有必要再定義輸出能控性。考慮下列狀態(tài)空間表達(dá)式所描述的線性定常系統(tǒng)式中，。如果能找到一個無約束的控制向量u(t)，在有限的時間間隔tott1內(nèi)，使任一給定的初始輸出y(to)轉(zhuǎn)移到任一最終輸出y(t1)，那么稱由式（2.2.11）和（2.2.12）所描述的系統(tǒng)為輸出能控的。可以證明，系統(tǒng)輸出能控的充要條件為：當(dāng)且僅當(dāng)m×(n+1)r維輸出能控性矩陣的秩為m時，由式（2.2.11）

18、和（2.2.12）所描述的系統(tǒng)為輸出能控的。注意，在式（2.2.12）中存在Du項，對確定輸出能控性是有幫助的。2.2.2 系統(tǒng)的能觀測性現(xiàn)在討論線性系統(tǒng)的能觀測性。考慮零輸入時的狀態(tài)空間表達(dá)式式中，。如果每一個狀態(tài)x(to)都可通過在有限時間間隔tott1內(nèi)，由y(t)觀測值確定，則稱系統(tǒng)為(完全)能觀測的。本節(jié)僅討論線性定常系統(tǒng)。不失一般性，設(shè)to=0。能觀測性的概念非常重要，這是由于在實際問題中，狀態(tài)反饋控制遇到的困難是一些狀態(tài)變量不易直接量測。因而在構(gòu)造控制器時，必須首先估計出不可量測的狀態(tài)變量。在“系統(tǒng)綜合”部分我們將指出，當(dāng)且僅當(dāng)系統(tǒng)是能觀測時，才能對系統(tǒng)狀態(tài)變量進(jìn)行觀測或估計。在

19、下面討論能觀測性條件時，我們將只考慮由式（2.2.13）和（2.2.14）給定的零輸入系統(tǒng)。這是因為，若采用如下狀態(tài)空間表達(dá)式則從而 由于矩陣A、B、C和D均為已知，u(t)也已知，所以上式右端的最后兩項為已知，因而它們可以從被量測值y(t)中消去。因此，為研究能觀測性的充要條件，只考慮式（2.2.13）和()所描述的零輸入系統(tǒng)就可以了。1. 定常系統(tǒng)狀態(tài)能觀測性的代數(shù)判據(jù) 考慮由式（3.13）和()所描述的線性定常系統(tǒng)。將其重寫為 易知，其輸出向量為 將寫為A的有限項的形式，即因而或() 顯然，如果系統(tǒng)是能觀測的，那么在0tt1時間間隔內(nèi)，給定輸出y(t)，就可由式（2.2.15）唯一地確定

20、出x(0)。可以證明，這就要求nm×n維能觀測性矩陣的秩為n。 由上述分析，我們可將能觀測的充要條件表述為：由式（2.2.13）和（2.2.14）所描述的線性定常系統(tǒng)，當(dāng)且僅當(dāng)n×nm維能觀測性矩陣的秩為n，即時，該系統(tǒng)才是能觀測的。例2.10 試判斷由式所描述的系統(tǒng)是否為能控和能觀測的。解 由于能控性矩陣的秩為2，即，故該系統(tǒng)是狀態(tài)能控的。 對于輸出能控性，可由系統(tǒng)輸出能控性矩陣的秩確定。由于的秩為1，即，故該系統(tǒng)是輸出能控的。為了檢驗?zāi)苡^測性條件，我們來驗算能觀測性矩陣的秩。由于 的秩為2，故此系統(tǒng)是能觀測的。2. 用傳遞函數(shù)矩陣表達(dá)的能觀測性條件類似地，能觀測性條件也

21、可用傳遞函數(shù)或傳遞函數(shù)矩陣表達(dá)。此時能觀測性的充要條件是：在傳遞函數(shù)或傳遞函數(shù)矩陣中不發(fā)生相約現(xiàn)象。如果存在相約，則約去的模態(tài)其輸出就不能觀測了。例2.11 證明下列系統(tǒng)是不能觀測的。式中解 由于能觀測性矩陣注意到即，故該系統(tǒng)是不能觀測的。 事實上，在該系統(tǒng)的傳遞函數(shù)中存在相約因子。由于X1(s)和U (s)之間的傳遞函數(shù)為 又Y (s)和X1(s)之間的傳遞函數(shù)為故Y(s)與U(s)之間的傳遞函數(shù)為顯然，分子、分母多項式中的因子（s+1）可以約去。這意味著，該系統(tǒng)是不能觀測的，或者說一些不為零的初始狀態(tài)x(0)不能由y(t)的量測值確定。3.說明當(dāng)且僅當(dāng)系統(tǒng)是狀態(tài)能控和能觀測時，其傳遞函數(shù)才

22、沒有相約因子。這意味著，可相約的傳遞函數(shù)不具有表征動態(tài)系統(tǒng)的所有信息。4. 狀態(tài)能觀測性條件的標(biāo)準(zhǔn)形判據(jù) 考慮由式（2.2.13）和（2.2.14）所描述的線性定常系統(tǒng)，將其重寫為 設(shè)非奇異線性變換矩陣P可將A化為對角線矩陣，式中，為對角線矩陣。定義式（2.2.16）和()可寫為如下對角線標(biāo)準(zhǔn)形因此或如果m×n維矩陣CP的任一列中都不含全為零的元素，那么系統(tǒng)是能觀測的。這是因為，如果CP的第i列含全為零的元素，則在輸出方程中將不出現(xiàn)狀態(tài)變量，因而不能由y上述判斷方法只適用于能將系統(tǒng)的狀態(tài)空間表達(dá)式（2.2.16）和()化為對角線標(biāo)準(zhǔn)形的情況。如果不能將式（2.2.16）和()變換為對角線標(biāo)準(zhǔn)形，則可利用一個合適的線性變換矩陣S，將其中的系統(tǒng)矩陣A變換為Jordan標(biāo)準(zhǔn)形。式中， J為Jordan標(biāo)準(zhǔn)形矩陣。 定義則式（2.2.16）和()可寫為如下Jordan標(biāo)準(zhǔn)形因此系統(tǒng)能觀測的充要條件為：（1） J中沒有兩個