現(xiàn)代控制理論(第二章)

2.1線性定常齊次狀態(tài)方程的解(自由解)2.2矩陣指數(shù)函數(shù)——狀態(tài)轉(zhuǎn)移矩陣2.3線性定常系統(tǒng)非齊次方程的解2.4*線性時變系統(tǒng)的解2.5*離散時間系統(tǒng)狀態(tài)方程的解2.6*連續(xù)時間狀態(tài)空間表達(dá)式的離散化第二章控制系統(tǒng)狀態(tài)空間表達(dá)式的解現(xiàn)在是1頁\一共有55頁\編輯于星期三2.1線性定常齊次狀態(tài)方程的解(自由解)

所謂系統(tǒng)的自由解，是指系統(tǒng)輸入為零時，由初始狀態(tài)引起的自由運(yùn)動。此時，狀態(tài)方程為齊次微分方程：(1)若初始時刻

時的狀態(tài)給定為則式(1)有唯一確定解：(2)若初始時刻從開始，即則其解為：(3)

證明：

和標(biāo)量微分方程求解類似，先假設(shè)式(1)的解為的矢量冪級數(shù)形式現(xiàn)在是2頁\一共有55頁\編輯于星期三(4)代入式(1)得：(5)

既然式(4)是式(1)的解，則式(5)對任意時刻都成立，故的同次冪項(xiàng)的系數(shù)應(yīng)相等，有：現(xiàn)在是3頁\一共有55頁\編輯于星期三在式(4)中，令，可得：將以上結(jié)果代入式(4)，故得：(6)現(xiàn)在是4頁\一共有55頁\編輯于星期三

等式右邊括號內(nèi)的展開式是矩陣，它是一個矩陣指數(shù)函數(shù)，記為，即(7)于是式(6)可表示為：

再用代替即在代替的情況下，同樣可以證明式2)的正確性?，F(xiàn)在是5頁\一共有55頁\編輯于星期三2.2矩陣指數(shù)函數(shù)——狀態(tài)轉(zhuǎn)移矩陣2.2.1狀態(tài)轉(zhuǎn)移矩陣齊次微分方程(1)的自由解為：或該式反應(yīng)了狀態(tài)矢量由初始狀態(tài)到任意時刻的矢量變換關(guān)系，反應(yīng)了狀態(tài)矢量在空間隨時間轉(zhuǎn)移的規(guī)律，因此稱為狀態(tài)轉(zhuǎn)移矩陣?，F(xiàn)在是6頁\一共有55頁\編輯于星期三2.2矩陣指數(shù)函數(shù)——狀態(tài)轉(zhuǎn)移矩陣注：狀態(tài)矩陣一般不是常數(shù)，而是時間的函數(shù)

起始矢量可以任意取，系統(tǒng)求解區(qū)間可任意選定—狀態(tài)空間法的優(yōu)點(diǎn)

滿足初始狀態(tài)的解是：滿足初始狀態(tài)的解是：令：則有：線性定常系統(tǒng)的狀態(tài)轉(zhuǎn)移矩陣現(xiàn)在是7頁\一共有55頁\編輯于星期三2.性質(zhì)二或(2)3.性質(zhì)三或(3)1．性質(zhì)一這就是組合性質(zhì)，它意味著從轉(zhuǎn)移到0，再從0轉(zhuǎn)移到的組合。或(1)2.2.2狀態(tài)轉(zhuǎn)移矩陣(矩陣指數(shù)函數(shù))的基本性質(zhì)注：本性質(zhì)可用于判斷矩陣是否符合狀態(tài)轉(zhuǎn)移矩陣的條件現(xiàn)在是8頁\一共有55頁\編輯于星期三或(4)這個性質(zhì)說明，矩陣與A矩陣是可以交換的。注：本性質(zhì)還表明，由狀態(tài)轉(zhuǎn)移矩陣可反推A！5.性質(zhì)五

對于方陣A和B，當(dāng)且僅當(dāng)AB=BA時，有而當(dāng)AB≠BA是，則

這個性質(zhì)說明，除非距陣A與B是可交換的，它們各目的矩陣指數(shù)函數(shù)之積與其和的矩陣指數(shù)函數(shù)不等價。這與標(biāo)量指數(shù)函數(shù)的性質(zhì)是不同的。4.性質(zhì)四現(xiàn)在是9頁\一共有55頁\編輯于星期三1．若A為對角線矩陣，即(5)則(6)2.若A能夠通過非奇異變換予以對角線化，即2.2.3幾個特殊的矩陣指數(shù)函數(shù)現(xiàn)在是10頁\一共有55頁\編輯于星期三則(7)3.若A為約旦矩陣現(xiàn)在是11頁\一共有55頁\編輯于星期三則(8)4.若(9)現(xiàn)在是12頁\一共有55頁\編輯于星期三1.根據(jù)的定義直接計算2.變換A為約旦標(biāo)準(zhǔn)型（1）A特征根互異其中T是使A

變換為對角線矩陣的變換陣。由式（7），有：2.2.4

的計算編程，用計算機(jī)算，最終能得到收斂解。但很難得到解析解。例2-1現(xiàn)在是13頁\一共有55頁\編輯于星期三3.利用拉氏反變換法求(10)證明齊次微分方程兩邊取拉氏變換即故現(xiàn)在是14頁\一共有55頁\編輯于星期三4.應(yīng)用凱萊—哈密頓定理求對上式兩邊取拉氏反變換，從而得到齊次微分方程的解：（1）由凱萊—哈密頓定理，方陣A滿足其自身的特征方程，即所以有它是的線性組合。同理現(xiàn)在是15頁\一共有55頁\編輯于星期三以此類推，都可用線性表示。(2)在定義中，用上面的方法可以消去A的n及n以上的冪次項(xiàng)，即（11）（3）的計算公式現(xiàn)在是16頁\一共有55頁\編輯于星期三A的特征值互異時，則

證明根據(jù)A滿足其自身特征方程的定理，可知特征值和A

是可以互換的，因此，也必須滿足式（11），從而有：（12）現(xiàn)在是17頁\一共有55頁\編輯于星期三上式對求解，即得式（12）。A的特征值均相同，為時，則證明同上，有：(13)現(xiàn)在是18頁\一共有55頁\編輯于星期三上式對，求異數(shù)，有：再對求異數(shù)，有：重復(fù)以上步驟，最后有：由上面的n個方程，對求解，記得公式（13）?，F(xiàn)在是19頁\一共有55頁\編輯于星期三2）用標(biāo)準(zhǔn)型法求解特征值互異，轉(zhuǎn)化成對角標(biāo)準(zhǔn)型，且A為友矩陣特征值：例2-1，2-2，2-4：求以下矩陣A的狀態(tài)轉(zhuǎn)移矩陣[解]：1）直接算法（略）現(xiàn)在是20頁\一共有55頁\編輯于星期三3）用拉氏變換法求解

例2-6，利用凱萊-哈密頓定理—-----------------自學(xué)！例2-3與例2-7也請注意自學(xué)！現(xiàn)在是21頁\一共有55頁\編輯于星期三2.3線性定常系統(tǒng)非齊次方程的解

現(xiàn)在討論線性定常系統(tǒng)在控制作用作用下的強(qiáng)制運(yùn)動。此時狀態(tài)方程為非齊次矩陣微分方程：當(dāng)初始時刻

初始狀態(tài)時，其解為：式中，。（1）（2）當(dāng)初始時刻為，初始狀態(tài)為時，其解為：現(xiàn)在是22頁\一共有55頁\編輯于星期三式中，。（3）證明采用類似標(biāo)量微分方程求解的方法，將式（1）寫成：等式兩邊同左乘，得：即（4）現(xiàn)在是23頁\一共有55頁\編輯于星期三對式（4）在上間積分，有：整理后可得式（2）：同理，若對式（4）在上積分，即可證明式（3）。式（2）也可從拉氏變換法求得，對式（1）進(jìn)行拉氏變換，有：即現(xiàn)在是24頁\一共有55頁\編輯于星期三上式左乘，得：（5）注意式（5）等式右邊第二項(xiàng)，其中：兩個拉氏變換函數(shù)的積是一個卷積的拉氏變換，即以此代入式（5），并取拉氏反變換，即得：現(xiàn)在是25頁\一共有55頁\編輯于星期三

在特定控制作用下，如脈沖函數(shù)、階躍函數(shù)和斜坡函數(shù)的激勵下，則系統(tǒng)的解式（2）可以簡化為以下公式：1.脈沖響應(yīng)即當(dāng)時2.階躍響應(yīng)即當(dāng)時3.斜坡響應(yīng)即當(dāng)時（6）（7）（8）例2-8要求掌握！現(xiàn)在是26頁\一共有55頁\編輯于星期三例2-8：已知系統(tǒng)狀態(tài)方程中試求解該系統(tǒng)的單位階躍響應(yīng)。解法一：積分法現(xiàn)在是27頁\一共有55頁\編輯于星期三現(xiàn)在是28頁\一共有55頁\編輯于星期三例2-8：已知系統(tǒng)狀態(tài)方程中試求解該系統(tǒng)的單位階躍響應(yīng)。解法二：拉氏變換法現(xiàn)在是29頁\一共有55頁\編輯于星期三現(xiàn)在是30頁\一共有55頁\編輯于星期三2.4*線性時變系統(tǒng)的解2.4.1時變系統(tǒng)狀態(tài)方程解的特點(diǎn)為了討論時變系統(tǒng)狀態(tài)方程的求解方法，現(xiàn)在先討論一個標(biāo)量時變系統(tǒng)：采用分離變量法，將上式寫成：對上式兩邊積分得：（1）現(xiàn)在是31頁\一共有55頁\編輯于星期三因此（2）或者寫成：

仿照定常系統(tǒng)齊次狀態(tài)方程的求解公式，式(2)中的也可以表示為狀態(tài)轉(zhuǎn)移矩陣，不過這時狀態(tài)轉(zhuǎn)移矩陣不僅是時間t的函數(shù)，而且也是初始時刻t。的函數(shù)。故采用符號來表示這個二元函數(shù)：（3）于是式(2)可寫成：（4）現(xiàn)在是32頁\一共有55頁\編輯于星期三能否將式(3)這個關(guān)系式也推廣到矢量方程：

遺憾的是，只有當(dāng)滿足乘法可交換條件，上述關(guān)系才能成立?，F(xiàn)證明如下：使之有（5）

如果是齊次方程的解，那么必須滿足：（6）現(xiàn)在是33頁\一共有55頁\編輯于星期三把展開成冪級數(shù)：上式兩邊對時間取導(dǎo)數(shù)：(7)(8)(9)把式(7)兩邊左乘有：比較式(8)和式(9)，可以看出，要使現(xiàn)在是34頁\一共有55頁\編輯于星期三成立，其必要和充分條件是：（10）

即是乘法可交換的。但是，這個條件是很苛刻的一般是不成立的。從而時變系統(tǒng)的自由解，通常不能像定常系統(tǒng)那樣寫成一個封閉形式。2.4.2線性時變齊次矩陣微分方程的解

盡管線性時變系統(tǒng)的自由解不能像定常系統(tǒng)那樣寫成一個封閉的解析形式，但仍然能表示為狀態(tài)轉(zhuǎn)移的形式。對于齊次矩陣微方程：現(xiàn)在是35頁\一共有55頁\編輯于星期三（11）其解為：（12）

式中，類似于前述線性定常系統(tǒng)中的，它也是非奇異方陣，并滿足如下的矩陣微分方程和初始條件：（13）（14）證明將解式(12)代入式(11)，有現(xiàn)在是36頁\一共有55頁\編輯于星期三即又在解式(12)中令，有：即

這就證明了，滿足式(13)、式(14)的,按式(12)所求得的是齊次微分方程(11)的解。2.4.3狀態(tài)轉(zhuǎn)移矩陣基本性質(zhì)現(xiàn)在是37頁\一共有55頁\編輯于星期三與線性定常系統(tǒng)的轉(zhuǎn)移矩陣類似，同樣有：因?yàn)椋呵夜适?15)成立。2），見式（14）。（15）1）3）（16）現(xiàn)在是38頁\一共有55頁\編輯于星期三因?yàn)閺氖?14)和式(15)可得：或

那么無論右乘，或左乘，式(16)都成立，故是非奇異陣，其逆存在，且等于。4）見式（13）。在這里，一般是不能交換的。2.4.4線性時變系統(tǒng)非齊次狀態(tài)方程式的解現(xiàn)在是39頁\一共有55頁\編輯于星期三線性時變系統(tǒng)的非齊次狀態(tài)方程為：且的元素在時間區(qū)間內(nèi)分段連續(xù)，則其解為：（17）（18）

證明線性系統(tǒng)滿足疊加原理，故可將式(17)的解看成由初始狀態(tài)的轉(zhuǎn)移和控制作用激勵的狀態(tài)的轉(zhuǎn)移兩部分組成。即（19）代入式(17)，有：現(xiàn)在是40頁\一共有55頁\編輯于星期三即可知：在t?！玹區(qū)間積分，有：于是現(xiàn)在是41頁\一共有55頁\編輯于星期三

在式(19)中令，并注意到中，可知，這樣由上式即可得到式(18)。2.4.5狀態(tài)轉(zhuǎn)移矩陣的計算因?yàn)锳是常數(shù)矩陣，所以上式直接表示為：在定常系統(tǒng)中，齊次狀態(tài)方程的解是：式中，，只與有關(guān)?，F(xiàn)在是42頁\一共有55頁\編輯于星期三在時變系統(tǒng)中，齊次狀態(tài)方程的解，一般的表示為：前已證明，只有當(dāng)是可交換時，即（20）才有：在一般情況下

對于不滿足式(20)的時變系統(tǒng)，的計算，一般采用級數(shù)近似法，即現(xiàn)在是43頁\一共有55頁\編輯于星期三（21）

這個關(guān)系式的證明是十分簡單的，只需驗(yàn)證它滿足式(13)的矩陣方程和式(14)的起始條件即可?？芍?21)滿足式(13)和式(14)?，F(xiàn)在是44頁\一共有55頁\編輯于星期三2.5*離散時間系統(tǒng)狀態(tài)方程的解2.5.1遞推法線性定常離散時問控制系統(tǒng)的狀態(tài)方程為：這個一陣差分方程的解為：或（1）現(xiàn)在是45頁\一共有55頁\編輯于星期三即（2）2.5.2Z變換法對于線性定常離散系統(tǒng)的狀態(tài)方程，也可以來用Z變換法來求解。

設(shè)定常離散系統(tǒng)的狀態(tài)方程是：對上式兩端進(jìn)行Z變換，有：或現(xiàn)在是46頁\一共有55頁\編輯于星期三所以：對上式兩端取Z的反變換，得：(3)對式2)和式(3)比較，有：（4）（5）

如果要獲得采樣瞬時之問的狀態(tài)和輸出，只需在此采樣周期內(nèi)，即在內(nèi)，利用連續(xù)狀態(tài)方程解的表達(dá)式：現(xiàn)在是47頁\一共有55頁\編輯于星期三

為了突出地表示f的有效期在，可以令(這里0≤△≤1)于是上式變成：（6）

顯然，這個公式的形式和離散狀態(tài)方程是完全一致的，如果使△的值在0和1之間變動，那么便可獲得采樣瞬時之間全部的狀態(tài)和輸出信息。將式(2)和式(3)比較，有（7）（8）二者形式上雖有不同，但實(shí)際上是完全一樣的。現(xiàn)在是48頁\一共有55頁\編輯于星期三2.6*連續(xù)時間狀態(tài)空間表達(dá)式的離散化2.6.1離散化方法對于連續(xù)時間的狀態(tài)空間表達(dá)式：將其離散化之后．則得離散時間狀態(tài)空問表達(dá)式為：C和D則仍與式(1)中的一樣。（1）（2）式中（4）（3）現(xiàn)在是49頁\一共有55頁\編輯于星期三2.6.2近似離散化

在采樣周期T較小時，一般當(dāng)其為系統(tǒng)最小時間常數(shù)的l/10左右時，離散化的狀態(tài)方程可近似表示為：（5）也就是說：（6）（7）證明根據(jù)導(dǎo)數(shù)的定義