控制系統(tǒng)的計(jì)算機(jī)輔助設(shè)計(jì)1_08

1、1 本章內(nèi)容本章內(nèi)容(1)利用MATLAB實(shí)現(xiàn)串聯(lián)頻率校正串聯(lián)頻率校正的三種方法；(2)利用MATLAB實(shí)現(xiàn)系統(tǒng)狀態(tài)反饋系統(tǒng)狀態(tài)反饋的兩種方法；(3)利用MATLAB實(shí)現(xiàn)系統(tǒng)狀態(tài)觀測(cè)器系統(tǒng)狀態(tài)觀測(cè)器的兩種方法；(4)利用MATLAB實(shí)現(xiàn)帶狀態(tài)觀測(cè)器的狀態(tài)反饋系統(tǒng)帶狀態(tài)觀測(cè)器的狀態(tài)反饋系統(tǒng)；(5)利用MATLAB實(shí)現(xiàn)系統(tǒng)的解耦；(6)利用MATLAB實(shí)現(xiàn)狀態(tài)反饋的線性二次型最優(yōu)控制器 的設(shè)計(jì)；(7) 利用MATLAB實(shí)現(xiàn)輸出反饋的線性二次型的最優(yōu)控制。 第第8章章 控制系統(tǒng)的計(jì)算機(jī)輔助設(shè)計(jì)控制系統(tǒng)的計(jì)算機(jī)輔助設(shè)計(jì)2并聯(lián)校正并聯(lián)校正 E(s) C(s) R(s) + - G(s) 固有特性 E(s

2、) C(s) R(s) Gc(s) + - G0(s) + - 固有系統(tǒng) 校正裝置 )()(1)()(00sGsGsGsGc局部閉環(huán)傳函局部閉環(huán)傳函 優(yōu)點(diǎn)：高靈敏度優(yōu)點(diǎn)：高靈敏度 高穩(wěn)定度高穩(wěn)定度3)(L20-60-40-2060400.0010.110.01-80-40-20-40820) 105. 005. 0)(1125. 0)(110() 1100(10)(22223sssssssGK420lgK0dB0-180L()()cLgcg5 控制系統(tǒng)的設(shè)計(jì)，就是在系統(tǒng)中引入控制系統(tǒng)的設(shè)計(jì)，就是在系統(tǒng)中引入適當(dāng)?shù)沫h(huán)節(jié)，用以對(duì)原有系統(tǒng)的某些性能適當(dāng)?shù)沫h(huán)節(jié)，用以對(duì)原有系統(tǒng)的某些性能進(jìn)行校正，使之達(dá)

3、到理想的效果，故進(jìn)行校正，使之達(dá)到理想的效果，故又稱又稱為系統(tǒng)的校正為系統(tǒng)的校正，下面介紹幾種常用的系統(tǒng)，下面介紹幾種常用的系統(tǒng)校正方法的計(jì)算機(jī)輔助設(shè)計(jì)實(shí)現(xiàn)。校正方法的計(jì)算機(jī)輔助設(shè)計(jì)實(shí)現(xiàn)。68.1 頻率法的串聯(lián)校正方法頻率法的串聯(lián)校正方法 應(yīng)用頻率法對(duì)系統(tǒng)進(jìn)行校正，其目的是改應(yīng)用頻率法對(duì)系統(tǒng)進(jìn)行校正，其目的是改變系統(tǒng)的頻率特性形狀，使校正后的系統(tǒng)頻變系統(tǒng)的頻率特性形狀，使校正后的系統(tǒng)頻率特性具有合適的低頻、中頻和高頻特性以率特性具有合適的低頻、中頻和高頻特性以及足夠的穩(wěn)定裕量，從而滿足所要求的性能及足夠的穩(wěn)定裕量，從而滿足所要求的性能指標(biāo)。指標(biāo)?？刂葡到y(tǒng)中控制系統(tǒng)中常用的串聯(lián)校正裝置常用的串聯(lián)

4、校正裝置是帶有是帶有單單零點(diǎn)與單極點(diǎn)零點(diǎn)與單極點(diǎn)的濾波器，若其的濾波器，若其零點(diǎn)比極點(diǎn)更零點(diǎn)比極點(diǎn)更靠近原點(diǎn)，靠近原點(diǎn)，則稱之為則稱之為超前校正超前校正，否則稱之為，否則稱之為滯后校正。滯后校正。78.1.1 基于頻率響應(yīng)法的串聯(lián)超前校正基于頻率響應(yīng)法的串聯(lián)超前校正1.超前校正裝置的特性超前校正裝置的特性設(shè)設(shè)超前超前校正裝置的傳遞函數(shù)為校正裝置的傳遞函數(shù)為其頻率特性為其頻率特性為) 1(11cTsTssG）（) 1(j1j1jcTTG）（jT1T1圖8-1 超前校正裝置的零極點(diǎn)圖o8(1) 極坐標(biāo)圖極坐標(biāo)圖 超前校正裝置的極坐標(biāo)圖如超前校正裝置的極坐標(biāo)圖如圖圖8-2所示。所示。 當(dāng)當(dāng)0變化時(shí)，

5、變化時(shí)，Gc(j)的相位角的相位角 0，Gc(j)的軌跡的軌跡為一半圓，由圖可得超前校正為一半圓，由圖可得超前校正的最大超前相位角的最大超前相位角m為為 (8-3) 令令 可得對(duì)應(yīng)于最大相位角可得對(duì)應(yīng)于最大相位角m時(shí)時(shí)的頻率的頻率m為為11sin1m0d)(dT1mj01)0(圖8-2 超前校正裝置的極坐標(biāo)圖m9(2) 對(duì)數(shù)坐標(biāo)圖對(duì)數(shù)坐標(biāo)圖 超前校正裝置的對(duì)數(shù)坐標(biāo)圖如超前校正裝置的對(duì)數(shù)坐標(biāo)圖如圖圖8-3所示。所示。 當(dāng)當(dāng) 由此可見(jiàn)，由此可見(jiàn)，超前校正裝置超前校正裝置是一是一個(gè)個(gè)高通濾波器高通濾波器(高頻通過(guò)，低高頻通過(guò)，低頻被衰減頻被衰減),它主要能使系統(tǒng)的它主要能使系統(tǒng)的瞬態(tài)響應(yīng)得到顯著改善

6、瞬態(tài)響應(yīng)得到顯著改善，而穩(wěn)，而穩(wěn)態(tài)精度的提高則較小。態(tài)精度的提高則較小。越大，越大，微分作用越強(qiáng)微分作用越強(qiáng)，從而超調(diào)量和，從而超調(diào)量和過(guò)渡過(guò)程時(shí)間等也越小。過(guò)渡過(guò)程時(shí)間等也越小。0|j|lg20, 0c）（Glg20|j|lg20,c）（Glg10|j|lg20,cm）（GT1T10lg10lg20L( )dB0mm圖8-3 超前校正裝置的對(duì)數(shù)坐標(biāo)圖m102.串聯(lián)超前校正方法串聯(lián)超前校正方法 超前校正超前校正裝置的主要作用是通過(guò)其裝置的主要作用是通過(guò)其相位相位超前超前效應(yīng)來(lái)改變頻率響應(yīng)曲線的形狀，效應(yīng)來(lái)改變頻率響應(yīng)曲線的形狀，產(chǎn)產(chǎn)生足夠大的相位超前角生足夠大的相位超前角,以補(bǔ)償原來(lái)系統(tǒng)中以補(bǔ)

7、償原來(lái)系統(tǒng)中元件造成的過(guò)大的相位滯后。因此校正時(shí)元件造成的過(guò)大的相位滯后。因此校正時(shí)應(yīng)使校正裝置的應(yīng)使校正裝置的最大超前相位角出現(xiàn)在最大超前相位角出現(xiàn)在校校正后正后系統(tǒng)的開(kāi)環(huán)剪切頻率系統(tǒng)的開(kāi)環(huán)剪切頻率(幅頻特性的交接幅頻特性的交接頻率頻率)c處處。11利用頻率法設(shè)計(jì)超前校正裝置的步驟利用頻率法設(shè)計(jì)超前校正裝置的步驟:(1)根據(jù)性能指標(biāo)對(duì)根據(jù)性能指標(biāo)對(duì)穩(wěn)態(tài)誤差穩(wěn)態(tài)誤差系數(shù)的要求系數(shù)的要求,確定確定開(kāi)環(huán)增益開(kāi)環(huán)增益k；(2)利用確定的開(kāi)環(huán)增益利用確定的開(kāi)環(huán)增益k，畫出，畫出未校正未校正系統(tǒng)的系統(tǒng)的Bode圖，圖，并求出其并求出其相位裕量相位裕量r0和幅值裕量和幅值裕量kg；(3)確定為使相位裕量達(dá)

8、到要求值，所需增加的超前相位確定為使相位裕量達(dá)到要求值，所需增加的超前相位角角，即，即r -r0+.式中式中 r為要求的相位裕量，為要求的相位裕量，是考是考慮到系統(tǒng)增加串聯(lián)超前校正裝置后系統(tǒng)的剪切頻率要向慮到系統(tǒng)增加串聯(lián)超前校正裝置后系統(tǒng)的剪切頻率要向右移而右移而附加的相位角附加的相位角，一般取，一般取=5 15 ；(4)令超前校正裝置的最大超前相位角令超前校正裝置的最大超前相位角=，則由下式，則由下式可可求出求出校正裝置的參數(shù)校正裝置的參數(shù)；ccsin1sin112(5)若將校正裝置的最大超前相位角處的頻率若將校正裝置的最大超前相位角處的頻率作為作為校正后系統(tǒng)的剪切頻率校正后系統(tǒng)的剪切頻率，

9、則有，則有即即或或由此可見(jiàn)，由此可見(jiàn)，未校正系統(tǒng)未校正系統(tǒng)的的對(duì)數(shù)幅頻對(duì)數(shù)幅頻特性幅值等于特性幅值等于-20lg（-10lg ）時(shí)的頻率即為時(shí)的頻率即為c；0| )(j)(j|lg20coccGGlg20| )(j|lg200| )(j|lg20lg20ococGG1| )(j|coG13(6) 根據(jù)根據(jù)=c，利用下式，利用下式求參數(shù)求參數(shù)T(7) 畫出校正后系統(tǒng)的畫出校正后系統(tǒng)的Bode圖，圖，檢驗(yàn)檢驗(yàn)性能指性能指標(biāo)是否已全部達(dá)到要求，若不滿足要求，標(biāo)是否已全部達(dá)到要求，若不滿足要求，可增大可增大值，從第三步起重新計(jì)算。值，從第三步起重新計(jì)算。c1T14例例8-1 設(shè)有一單位反饋系統(tǒng)，其開(kāi)環(huán)

10、傳遞函數(shù)為設(shè)有一單位反饋系統(tǒng)，其開(kāi)環(huán)傳遞函數(shù)為要求系統(tǒng)的穩(wěn)態(tài)速度誤差系數(shù)要求系統(tǒng)的穩(wěn)態(tài)速度誤差系數(shù)kv=20(1/s),相位裕量相位裕量 r50度度,幅值裕量幅值裕量kg10dB,試確定串聯(lián)校正裝置試確定串聯(lián)校正裝置解解 根據(jù)根據(jù)可求出可求出k =40，即，即 根據(jù)串聯(lián)超前校正的設(shè)計(jì)步驟，可編寫以下根據(jù)串聯(lián)超前校正的設(shè)計(jì)步驟，可編寫以下m文件。文件。ex8_1.m202)2(lim)(lim000kssksssGkssv)2(40)(0sssG)2()(0ssksG15執(zhí)行后可得如下結(jié)果及圖執(zhí)行后可得如下結(jié)果及圖8-4所示曲線所示曲線num/den = 0.22541 s + 1 - 0.05

11、3537 s + 1num/den = 9.0165 s + 40 - 0.053537 s3 + 1.1071 s2 + 2 s 校正前：幅值裕量校正前：幅值裕量=Inf dB,相位裕量相位裕量=17.9642 校正后：幅值裕量校正后：幅值裕量=Inf dB,相位裕量相位裕量=50.7196 16 圖圖8-4 超前校正裝置及校正前后系統(tǒng)的伯德圖超前校正裝置及校正前后系統(tǒng)的伯德圖178.1.2 基于頻率響應(yīng)法的串聯(lián)滯后校正基于頻率響應(yīng)法的串聯(lián)滯后校正 1.滯后校正裝置的特性滯后校正裝置的特性 設(shè)滯后校正裝置的傳遞設(shè)滯后校正裝置的傳遞函數(shù)為函數(shù)為其頻率特性為其頻率特性為) 1(11cTsTssG

12、）（) 1(j1j1(jcTTG）圖8-5 滯后校正裝置的零極點(diǎn)圖jT1T1o18（1）極坐標(biāo)圖）極坐標(biāo)圖 滯后校正裝置的極坐標(biāo)圖滯后校正裝置的極坐標(biāo)圖如圖如圖8-6所示。由圖可知，所示。由圖可知，當(dāng)當(dāng)=0變化時(shí)，變化時(shí)，Gc(j)的相位角的相位角0，Gc(j)的根的根軌跡為一半圓。軌跡為一半圓。 同理可求得同理可求得最大滯后相最大滯后相位角位角和對(duì)應(yīng)的頻率和對(duì)應(yīng)的頻率分別為分別為11sin1mT1m11)(0j0m圖8-6 滯后校正裝置的極坐標(biāo)圖19（2）對(duì)數(shù)坐標(biāo)圖）對(duì)數(shù)坐標(biāo)圖滯后校正裝置的對(duì)數(shù)坐標(biāo)圖如圖滯后校正裝置的對(duì)數(shù)坐標(biāo)圖如圖8-7所示所示由此可見(jiàn)，由此可見(jiàn)，滯后校正滯后校正裝置是一個(gè)

13、裝置是一個(gè)低低通濾波器通濾波器（低頻通過(guò)，高頻被衰（低頻通過(guò)，高頻被衰減），且減），且越大，越大，高頻衰減高頻衰減越厲越厲害，抗高頻干擾性能越好，但使響害，抗高頻干擾性能越好，但使響應(yīng)速度變慢，故滯后校正能使穩(wěn)態(tài)應(yīng)速度變慢，故滯后校正能使穩(wěn)態(tài)得到顯著提高，但瞬態(tài)響應(yīng)時(shí)間卻得到顯著提高，但瞬態(tài)響應(yīng)時(shí)間卻隨之而增加，隨之而增加，越大，積分作用越越大，積分作用越強(qiáng)，穩(wěn)態(tài)誤差越小強(qiáng)，穩(wěn)態(tài)誤差越小。0|lg20, 0c）（jGlg20|lg20,c）（jGlg10|lg20,cm）（jGlg20lg10T1T1mmdB00)(L圖8-7 滯后校正裝置的極坐標(biāo)圖m202.串聯(lián)滯后校正方法串聯(lián)滯后校正方法

14、滯后校正裝置的主要作用是在滯后校正裝置的主要作用是在高頻段造成幅值衰減高頻段造成幅值衰減，降低降低系統(tǒng)的系統(tǒng)的剪切頻率剪切頻率，以便能使系統(tǒng)獲得充分的相位裕，以便能使系統(tǒng)獲得充分的相位裕量，但應(yīng)同時(shí)保證系統(tǒng)在新的剪切頻率附近的相頻特性量，但應(yīng)同時(shí)保證系統(tǒng)在新的剪切頻率附近的相頻特性曲線變化不大。曲線變化不大。利用頻率法設(shè)計(jì)滯后校正裝置的步驟利用頻率法設(shè)計(jì)滯后校正裝置的步驟:（1）根據(jù)性能指標(biāo)對(duì)）根據(jù)性能指標(biāo)對(duì)穩(wěn)態(tài)誤差穩(wěn)態(tài)誤差系數(shù)的要求，確定系數(shù)的要求，確定開(kāi)環(huán)增開(kāi)環(huán)增益益k；（2）利用已確定的開(kāi)環(huán)增益）利用已確定的開(kāi)環(huán)增益k，畫出未校正系統(tǒng)的，畫出未校正系統(tǒng)的Bode圖，圖，并求出其相位裕量并

15、求出其相位裕量r0和幅值裕量和幅值裕量kg；21（3）如）如未校正系統(tǒng)未校正系統(tǒng)的相位和幅值裕量不滿足要求，的相位和幅值裕量不滿足要求，尋尋找一新的剪切頻率找一新的剪切頻率c，在，在c處開(kāi)環(huán)傳遞函數(shù)的相位角處開(kāi)環(huán)傳遞函數(shù)的相位角應(yīng)滿足下式應(yīng)滿足下式 Go(jc)= -180+ r + 式中式中 r為要求的相角裕量，為要求的相角裕量，是為補(bǔ)償滯后校正裝置的是為補(bǔ)償滯后校正裝置的相位滯后而附加的相位角，一般取相位滯后而附加的相位角，一般取=5 12 ； (4) 為使為使滯后校正裝置對(duì)系統(tǒng)的滯后校正裝置對(duì)系統(tǒng)的相位滯后影響較小相位滯后影響較?。ㄒ唬ㄒ话阆拗圃诎阆拗圃? 12 ），），m應(yīng)遠(yuǎn)離應(yīng)遠(yuǎn)離c

16、，一般取滯后校正裝，一般取滯后校正裝置的第一個(gè)交接頻率：置的第一個(gè)交接頻率：1=1/T=(1/51/10)c（即（即mc），此時(shí)有），此時(shí)有|Gc(jc)|=-20lg。1取得愈小，取得愈小，對(duì)系統(tǒng)的相位裕量影響愈小，但太小則校正裝置的時(shí)對(duì)系統(tǒng)的相位裕量影響愈小，但太小則校正裝置的時(shí)間常數(shù)間常數(shù)T將很大，這也是不允許的；將很大，這也是不允許的；22(5)確定使校正后系統(tǒng)的幅值曲線在新的剪切頻確定使校正后系統(tǒng)的幅值曲線在新的剪切頻率率c處下降到處下降到0dB所需的衰減量所需的衰減量20lg|Go(jc)|，并根據(jù)，并根據(jù)20lg |Go(jc)Gc(jc)|20lg |Go(jc)|-20lg=

17、0即即 =|Go(jc)|求求出校正裝置的參數(shù)出校正裝置的參數(shù)； (6) 畫出校正后系統(tǒng)的畫出校正后系統(tǒng)的Bode圖，圖，檢驗(yàn)檢驗(yàn)性能指標(biāo)性能指標(biāo)是否已全部達(dá)到要求，若不滿足要求，可增是否已全部達(dá)到要求，若不滿足要求，可增大大值，從第三步起重新計(jì)算。值，從第三步起重新計(jì)算。23例例8-2 設(shè)有一單位負(fù)反饋系統(tǒng)的開(kāi)環(huán)傳遞函數(shù)為設(shè)有一單位負(fù)反饋系統(tǒng)的開(kāi)環(huán)傳遞函數(shù)為 要求系統(tǒng)的穩(wěn)態(tài)速度誤差系數(shù)要求系統(tǒng)的穩(wěn)態(tài)速度誤差系數(shù)kv=5(1/s)，相位裕量，相位裕量r400，幅值裕量，幅值裕量kg10dB，試確定串聯(lián)校正裝置。，試確定串聯(lián)校正裝置。解解 根據(jù)根據(jù)可求出可求出k=5，即，即) 125. 0)(1

18、()(0sssksG5) 125. 0)(1(lim)(lim000ksssksssGkssv) 125. 0)(1(5)(0ssssG24根據(jù)串聯(lián)滯后校正的設(shè)計(jì)步驟，可編寫以下根據(jù)串聯(lián)滯后校正的設(shè)計(jì)步驟，可編寫以下m文件。文件。 ex8_2.m25執(zhí)行后可得如下結(jié)果及圖執(zhí)行后可得如下結(jié)果及圖8-8所示曲線。所示曲線。num/den = 8.3842s + 1 - 59.7135s + 1num/den = 41.9208s+5 - 14.9284s4+74.8918s3+60.9635s2+s校正前：幅值裕量校正前：幅值裕量= -3.8573e-015 dB，相位裕量，相位裕量=7.3342

19、e-006 校正后：幅值裕量校正后：幅值裕量=15.8574dB，相位裕量，相位裕量=40.6552 26圖圖8-8 滯后校正裝置及校正前后系統(tǒng)的伯德圖滯后校正裝置及校正前后系統(tǒng)的伯德圖278.1.3 基于頻率響應(yīng)法的串聯(lián)滯后基于頻率響應(yīng)法的串聯(lián)滯后-超前校正超前校正1.滯后滯后-超前校正裝置的特性超前校正裝置的特性設(shè)滯后設(shè)滯后-超前校正裝置的傳遞函數(shù)為超前校正裝置的傳遞函數(shù)為上式等號(hào)右邊的上式等號(hào)右邊的第一項(xiàng)第一項(xiàng)產(chǎn)生產(chǎn)生超前超前網(wǎng)絡(luò)的作網(wǎng)絡(luò)的作用，而用，而第二項(xiàng)第二項(xiàng)產(chǎn)生產(chǎn)生滯后滯后網(wǎng)絡(luò)的作用。網(wǎng)絡(luò)的作用。) 1,(1111G122211cTTsTsTsTsTs）（28(1) 極坐標(biāo)圖極坐

20、標(biāo)圖 滯后滯后-超前校正裝置的極超前校正裝置的極坐標(biāo)圖如圖坐標(biāo)圖如圖8-9所示。所示。 由圖可知，當(dāng)角頻率由圖可知，當(dāng)角頻率在在00之間變化時(shí)之間變化時(shí), 滯后滯后-超前校正裝置起著超前校正裝置起著相位相位滯后校正滯后校正的作用；當(dāng)?shù)淖饔?；?dāng)在在0之間變化時(shí)之間變化時(shí),它起它起著著超前校正超前校正的作用，對(duì)應(yīng)的作用，對(duì)應(yīng)相位角為零的頻率相位角為零的頻率0為為2101TT010mIeR0圖8-9 滯后-超前校正裝置的極坐標(biāo)圖29（2）對(duì)數(shù)坐標(biāo)圖）對(duì)數(shù)坐標(biāo)圖 滯后滯后-超前校正超前校正裝置的對(duì)數(shù)坐裝置的對(duì)數(shù)坐標(biāo)圖如圖標(biāo)圖如圖8-10所示。從圖可所示。從圖可清楚看出，當(dāng)清楚看出，當(dāng)00時(shí)滯后時(shí)滯后-

21、超前校正裝置超前校正裝置起著起著相位滯后相位滯后校正校正的作用；的作用；當(dāng)當(dāng)0A=-2 -1 1;1 0 1;-1 0 1; b=1;1;1; rc=rank(ctrb(A,b); p=-1,-2,-3; K=acker(A,b,p)結(jié)果顯示結(jié)果顯示K= -1 2 443 對(duì)于對(duì)于多變量系統(tǒng)的極點(diǎn)配置多變量系統(tǒng)的極點(diǎn)配置，MATLAB控制系統(tǒng)工具箱中也給出了函數(shù)控制系統(tǒng)工具箱中也給出了函數(shù)place( )，其調(diào)用格式為其調(diào)用格式為K=place(A,B,P)例例8-5 已知系統(tǒng)的狀態(tài)方程為已知系統(tǒng)的狀態(tài)方程為 求使?fàn)顟B(tài)反饋系統(tǒng)的閉環(huán)極點(diǎn)為求使?fàn)顟B(tài)反饋系統(tǒng)的閉環(huán)極點(diǎn)為-2,-3,(-1j3)/2

22、的狀態(tài)反饋陣的狀態(tài)反饋陣K。uxx331134029514102033821310140044解解 MATLAB程序?yàn)槌绦驗(yàn)?ex8_5.m執(zhí)行后得執(zhí)行后得K= 32.5923 65.6844 58.8332 46.6557 32.5923 65.6844 58.8332 46.6557 55.4594 111.8348 103.6800 81.0239 55.4594 111.8348 103.6800 81.0239452.部分極點(diǎn)配置部分極點(diǎn)配置 在一些特定的應(yīng)用中，有時(shí)沒(méi)有必要在一些特定的應(yīng)用中，有時(shí)沒(méi)有必要去對(duì)所有的極點(diǎn)進(jìn)行重新配置，而只需去對(duì)所有的極點(diǎn)進(jìn)行重新配置，而只需對(duì)其中若干

23、個(gè)極點(diǎn)進(jìn)行配置，使得其他對(duì)其中若干個(gè)極點(diǎn)進(jìn)行配置，使得其他極點(diǎn)保持原來(lái)的值，例如若系統(tǒng)開(kāi)環(huán)模極點(diǎn)保持原來(lái)的值，例如若系統(tǒng)開(kāi)環(huán)模型是不穩(wěn)定的，則型是不穩(wěn)定的，則可以將那些不穩(wěn)定的可以將那些不穩(wěn)定的極點(diǎn)配置成穩(wěn)定的值極點(diǎn)配置成穩(wěn)定的值，而不去改變那些而不去改變那些原本穩(wěn)定的極點(diǎn)原本穩(wěn)定的極點(diǎn)。作這樣配置的。作這樣配置的前提前提條條件是件是原系統(tǒng)沒(méi)有重極點(diǎn)原系統(tǒng)沒(méi)有重極點(diǎn)，這就能保證由，這就能保證由系統(tǒng)特征向量構(gòu)成的矩陣是非奇異的。系統(tǒng)特征向量構(gòu)成的矩陣是非奇異的。46 假設(shè)假設(shè)xi為對(duì)應(yīng)于為對(duì)應(yīng)于i的特征向量，即的特征向量，即A xi =i xi，這樣可，這樣可以對(duì)各個(gè)特征值構(gòu)造特征向量矩陣以對(duì)各

24、個(gè)特征值構(gòu)造特征向量矩陣X=x1,x2,xn，由，由前面的假設(shè)可知前面的假設(shè)可知X矩陣為非奇異的，故可以得出其逆陣矩陣為非奇異的，故可以得出其逆陣T=X-1，且令，且令T的第的第i個(gè)行向量為個(gè)行向量為Ti，且想把，且想把i配置到配置到i的的位置，則可以定義變量位置，則可以定義變量ri=（i-i）/bi，其中，其中bi為向量為向量Tb的第的第i個(gè)分量，這時(shí)配置全部的極點(diǎn)，則可以得出狀態(tài)個(gè)分量，這時(shí)配置全部的極點(diǎn)，則可以得出狀態(tài)反饋陣反饋陣 特別地，特別地，若不想對(duì)哪個(gè)極點(diǎn)進(jìn)行重新配置，則可以將若不想對(duì)哪個(gè)極點(diǎn)進(jìn)行重新配置，則可以將對(duì)應(yīng)的項(xiàng)從上面的求和式子中刪除就可以得出相應(yīng)的狀對(duì)應(yīng)的項(xiàng)從上面的求

25、和式子中刪除就可以得出相應(yīng)的狀態(tài)反饋陣態(tài)反饋陣，它能按指定的方式進(jìn)行極點(diǎn)配置。，它能按指定的方式進(jìn)行極點(diǎn)配置。niiiTrK147例例8-6對(duì)于例對(duì)于例8-4所示系統(tǒng)，實(shí)際上只有一所示系統(tǒng)，實(shí)際上只有一個(gè)不穩(wěn)定的極點(diǎn)個(gè)不穩(wěn)定的極點(diǎn)1,若僅將此極點(diǎn)配置到若僅將此極點(diǎn)配置到-5，試采用，試采用部分極點(diǎn)配置方法部分極點(diǎn)配置方法對(duì)其進(jìn)行。對(duì)其進(jìn)行。解解 MATLAB程序?yàn)槌绦驗(yàn)閑x8_6.m執(zhí)行后得執(zhí)行后得K= 1.5000 -1.5000 -6.0000488.2.2 狀態(tài)觀測(cè)器狀態(tài)觀測(cè)器1.全維狀態(tài)觀測(cè)器的設(shè)計(jì)全維狀態(tài)觀測(cè)器的設(shè)計(jì) 極點(diǎn)配置是基于狀態(tài)反饋，因此狀態(tài)極點(diǎn)配置是基于狀態(tài)反饋，因此狀態(tài)x

26、必須必須可量測(cè)，可量測(cè)，當(dāng)狀態(tài)不能量測(cè)時(shí)，則應(yīng)設(shè)計(jì)狀態(tài)當(dāng)狀態(tài)不能量測(cè)時(shí)，則應(yīng)設(shè)計(jì)狀態(tài)觀測(cè)器來(lái)估計(jì)狀態(tài)觀測(cè)器來(lái)估計(jì)狀態(tài)。對(duì)于系統(tǒng)對(duì)于系統(tǒng)若系統(tǒng)完全能觀測(cè)若系統(tǒng)完全能觀測(cè)，則可構(gòu)造如圖，則可構(gòu)造如圖8-12所示的狀所示的狀態(tài)觀測(cè)器。態(tài)觀測(cè)器。CxyBuAxx 49uyBACBACL-+-圖8-12 帶狀態(tài)觀測(cè)器的系統(tǒng)-+xx 50由上圖可得觀測(cè)器的狀態(tài)方程為由上圖可得觀測(cè)器的狀態(tài)方程為即即 其特征多項(xiàng)式為其特征多項(xiàng)式為 f（s）=|sI-(A-LC)| 由于工程上要求能比較快速的逼近由于工程上要求能比較快速的逼近x，只要只要調(diào)整反饋陣調(diào)整反饋陣L，觀測(cè)器的極點(diǎn)就可以任意配置，觀測(cè)器的極點(diǎn)就可以任意

27、配置達(dá)到要求的性能達(dá)到要求的性能，所以，所以，觀測(cè)器的設(shè)計(jì)與狀態(tài)反饋觀測(cè)器的設(shè)計(jì)與狀態(tài)反饋極點(diǎn)配置的設(shè)計(jì)類似極點(diǎn)配置的設(shè)計(jì)類似。LyxLCBuxAxLyBuxLCAx)(51 假 定 單 變 量 系 統(tǒng) 所 要 求 的假 定 單 變 量 系 統(tǒng) 所 要 求 的 n 個(gè) 觀 測(cè) 器 的 極 點(diǎn) 為個(gè) 觀 測(cè) 器 的 極 點(diǎn) 為1,2,n，則可求出，則可求出期望的狀態(tài)觀測(cè)器的特征方程期望的狀態(tài)觀測(cè)器的特征方程為為 f*(s)=(-1)(-2)(-n)=sn + a1sn-1 + + an這時(shí)可求得反饋陣這時(shí)可求得反饋陣L為為 式中式中 ， , f*（A）是將系統(tǒng)期望的觀測(cè)器）是將系統(tǒng)期望的觀測(cè)器特

28、征方程特征方程中中s換成系統(tǒng)矩陣換成系統(tǒng)矩陣A后的矩陣多項(xiàng)式。后的矩陣多項(xiàng)式。100)(10*VAfL10nCACACV52 利用對(duì)偶原理，可使設(shè)計(jì)問(wèn)題大為簡(jiǎn)化利用對(duì)偶原理，可使設(shè)計(jì)問(wèn)題大為簡(jiǎn)化，求解，求解過(guò)程如下：過(guò)程如下：首先構(gòu)造首先構(gòu)造系統(tǒng)式系統(tǒng)式(8-14)的的對(duì)偶系統(tǒng)對(duì)偶系統(tǒng) (8-15)然后，根據(jù)下式可然后，根據(jù)下式可求得狀態(tài)觀測(cè)器的反饋陣求得狀態(tài)觀測(cè)器的反饋陣L。 LT=acker(AT,CT,P)或或 LT=place(AT,CT,P)其中其中 P為給定的極點(diǎn)，為給定的極點(diǎn)，L為狀態(tài)觀測(cè)器的反饋陣。為狀態(tài)觀測(cè)器的反饋陣。zBwnCAzTTT53例例8-7 已知開(kāi)環(huán)系統(tǒng)已知開(kāi)環(huán)系

29、統(tǒng)其中其中設(shè)計(jì)全維狀態(tài)觀測(cè)器設(shè)計(jì)全維狀態(tài)觀測(cè)器，使觀測(cè)器的閉環(huán)極點(diǎn)為，使觀測(cè)器的閉環(huán)極點(diǎn)為 -2j23,-5。CxybuAxx 001,100,6116100010CbA54解解 為求出狀態(tài)觀測(cè)器的反饋陣為求出狀態(tài)觀測(cè)器的反饋陣L，先為原，先為原系統(tǒng)系統(tǒng)構(gòu)造一對(duì)偶系統(tǒng)構(gòu)造一對(duì)偶系統(tǒng)。然后采用極點(diǎn)配置方法然后采用極點(diǎn)配置方法對(duì)對(duì)偶系統(tǒng)進(jìn)行閉環(huán)對(duì)對(duì)偶系統(tǒng)進(jìn)行閉環(huán)極點(diǎn)位置的配置極點(diǎn)位置的配置，得到反饋陣，得到反饋陣K，從而可由，從而可由對(duì)對(duì)偶原理偶原理得到原系統(tǒng)的狀態(tài)觀測(cè)器的反饋陣得到原系統(tǒng)的狀態(tài)觀測(cè)器的反饋陣L。MATLAB程序?yàn)槌绦驗(yàn)閑x8_7.mzBwnCAzTTT55執(zhí)行后得執(zhí)行后得The

30、Rank of Obstrabilaty Matrixr0 = 3 3L = 3.0000 3.0000 7.0000 7.0000 -1.0000 -1.0000 由于由于rankr0=3，所以系統(tǒng)能觀測(cè)，因此可設(shè)計(jì)全，所以系統(tǒng)能觀測(cè)，因此可設(shè)計(jì)全維狀態(tài)觀測(cè)器。維狀態(tài)觀測(cè)器。562.2.降維觀測(cè)器的設(shè)計(jì)降維觀測(cè)器的設(shè)計(jì) 前面所討論的狀態(tài)觀測(cè)器的維數(shù)和被控系統(tǒng)的前面所討論的狀態(tài)觀測(cè)器的維數(shù)和被控系統(tǒng)的維數(shù)相同，故稱為全維觀測(cè)器，實(shí)際上系統(tǒng)的輸維數(shù)相同，故稱為全維觀測(cè)器，實(shí)際上系統(tǒng)的輸出出y總是能夠觀測(cè)的。因此，可以總是能夠觀測(cè)的。因此，可以利用系統(tǒng)的輸利用系統(tǒng)的輸出量出量y來(lái)直接產(chǎn)生部分狀態(tài)變

31、量來(lái)直接產(chǎn)生部分狀態(tài)變量，從而降低觀測(cè)，從而降低觀測(cè)器的維數(shù)。假設(shè)系統(tǒng)是完全能觀測(cè)器，若狀態(tài)器的維數(shù)。假設(shè)系統(tǒng)是完全能觀測(cè)器，若狀態(tài)x為為n維，輸出維，輸出y為為m維，由于維，由于y是可量測(cè)的，因此是可量測(cè)的，因此只需對(duì)只需對(duì)n-m個(gè)狀態(tài)進(jìn)行觀測(cè)個(gè)狀態(tài)進(jìn)行觀測(cè)，也就是說(shuō)用，也就是說(shuō)用(n-m)維維的狀態(tài)觀測(cè)器可以代替全維觀測(cè)器，這樣觀測(cè)器的狀態(tài)觀測(cè)器可以代替全維觀測(cè)器，這樣觀測(cè)器的結(jié)構(gòu)可以大大簡(jiǎn)化。的結(jié)構(gòu)可以大大簡(jiǎn)化。5758Matlab程序?yàn)椋篹x8_8598.2.3 帶狀態(tài)觀測(cè)器的狀態(tài)反饋系統(tǒng)帶狀態(tài)觀測(cè)器的狀態(tài)反饋系統(tǒng) 狀態(tài)觀測(cè)器狀態(tài)觀測(cè)器解決了受控系統(tǒng)的狀態(tài)重解決了受控系統(tǒng)的狀態(tài)重構(gòu)問(wèn)題

32、，構(gòu)問(wèn)題，為為那些狀態(tài)變量不能直接量測(cè)得那些狀態(tài)變量不能直接量測(cè)得到的系統(tǒng)到的系統(tǒng)實(shí)現(xiàn)狀態(tài)反饋創(chuàng)造了條件實(shí)現(xiàn)狀態(tài)反饋創(chuàng)造了條件。帶狀。帶狀態(tài)觀測(cè)器的狀態(tài)反饋系統(tǒng)由三部分組成，態(tài)觀測(cè)器的狀態(tài)反饋系統(tǒng)由三部分組成，即即原系統(tǒng)原系統(tǒng)、觀測(cè)器觀測(cè)器和和控制器控制器，圖，圖8-13是一是一個(gè)帶有全維觀測(cè)器的狀態(tài)反饋系統(tǒng)。個(gè)帶有全維觀測(cè)器的狀態(tài)反饋系統(tǒng)。60uryBACBACL-+-圖8-13 帶狀態(tài)觀測(cè)器的狀態(tài)反饋系統(tǒng)-+Kxx +-61設(shè)能控能觀測(cè)的受控系統(tǒng)為設(shè)能控能觀測(cè)的受控系統(tǒng)為 （8-21）狀態(tài)反饋控制律為狀態(tài)反饋控制律為 （8-22）狀態(tài)觀測(cè)器方程為狀態(tài)觀測(cè)器方程為 （8-23）由以上三式可得

33、由以上三式可得閉環(huán)系統(tǒng)的狀態(tài)空間表達(dá)式閉環(huán)系統(tǒng)的狀態(tài)空間表達(dá)式為為CxyBuAxx xKruLyBuxLCAx)(CxyBrxBKLCALCxxBrxBKAxx)(62 可以證明，由觀測(cè)器構(gòu)成的狀態(tài)反饋閉環(huán)系可以證明，由觀測(cè)器構(gòu)成的狀態(tài)反饋閉環(huán)系統(tǒng)，其統(tǒng)，其特征多項(xiàng)式特征多項(xiàng)式等于等于狀態(tài)反饋部分狀態(tài)反饋部分的特征多項(xiàng)的特征多項(xiàng)式式|sI-(A-BK)|和和觀測(cè)器部分觀測(cè)器部分的特征多項(xiàng)式的特征多項(xiàng)式|sI-(A-LC)|的的乘積乘積，而且兩者相互獨(dú)立。因此，只要系，而且兩者相互獨(dú)立。因此，只要系統(tǒng)統(tǒng)0(A,B,C)能控能觀測(cè)，則系統(tǒng)的狀態(tài)反饋陣能控能觀測(cè)，則系統(tǒng)的狀態(tài)反饋陣K和觀測(cè)器反饋陣和

34、觀測(cè)器反饋陣L可分別根據(jù)各自的要求，可分別根據(jù)各自的要求，獨(dú)立獨(dú)立進(jìn)行配置進(jìn)行配置，這種性質(zhì)被稱為，這種性質(zhì)被稱為分離特性分離特性。 同理，用降維觀測(cè)器構(gòu)成的反饋系統(tǒng)也具有分同理，用降維觀測(cè)器構(gòu)成的反饋系統(tǒng)也具有分離特性離特性63例例8-9 已知開(kāi)環(huán)系統(tǒng)已知開(kāi)環(huán)系統(tǒng)設(shè)計(jì)狀態(tài)反饋使閉環(huán)極點(diǎn)為設(shè)計(jì)狀態(tài)反饋使閉環(huán)極點(diǎn)為 -1.8j2.4，而且狀態(tài)不，而且狀態(tài)不可量測(cè)，因此設(shè)計(jì)狀態(tài)觀測(cè)器使其閉環(huán)極點(diǎn)為可量測(cè)，因此設(shè)計(jì)狀態(tài)觀測(cè)器使其閉環(huán)極點(diǎn)為-8,-8。解解 狀態(tài)反饋和狀態(tài)觀測(cè)器的設(shè)計(jì)分開(kāi)進(jìn)行狀態(tài)反饋和狀態(tài)觀測(cè)器的設(shè)計(jì)分開(kāi)進(jìn)行，狀態(tài)觀，狀態(tài)觀測(cè)器的設(shè)計(jì)借助于對(duì)偶原理。在設(shè)計(jì)之前，應(yīng)先判別測(cè)器的設(shè)計(jì)借助于

35、對(duì)偶原理。在設(shè)計(jì)之前，應(yīng)先判別系統(tǒng)的能控性和能觀測(cè)性，系統(tǒng)的能控性和能觀測(cè)性，MATLAB的程序的程序 ex8_9.mxyuxx011006 .201064執(zhí)行后得執(zhí)行后得The rank of Controllability Matrixrc = 2The rank of Observability Matrixro = 2K = 29.6000 3.6000L = 16.0000 84.6000658.2.4 離散系統(tǒng)的極點(diǎn)配置和狀態(tài)觀測(cè)器的設(shè)計(jì)離散系統(tǒng)的極點(diǎn)配置和狀態(tài)觀測(cè)器的設(shè)計(jì) 離散離散系統(tǒng)的系統(tǒng)的極點(diǎn)配置極點(diǎn)配置和和狀態(tài)觀測(cè)器狀態(tài)觀測(cè)器的設(shè)計(jì)的的設(shè)計(jì)的求解過(guò)程與連續(xù)系統(tǒng)基本相同，在求

36、解過(guò)程與連續(xù)系統(tǒng)基本相同，在MATLAB中，中，可直接采用工具箱中的可直接采用工具箱中的place( )和和acker( )函數(shù)函數(shù)進(jìn)行設(shè)計(jì)，這里不在贅述。進(jìn)行設(shè)計(jì)，這里不在贅述。66Matlab程序?yàn)槌绦驗(yàn)閑x8_10.m678.2.5 系統(tǒng)解耦系統(tǒng)解耦 在在多變量系統(tǒng)多變量系統(tǒng)中，如果中，如果傳遞函數(shù)陣不是傳遞函數(shù)陣不是對(duì)角矩陣對(duì)角矩陣，則不同的輸入與輸出之間存在，則不同的輸入與輸出之間存在著耦合，即第著耦合，即第i輸入不但會(huì)對(duì)第輸入不但會(huì)對(duì)第i輸出有影響，輸出有影響，而且還會(huì)影響到其他的輸出，就而且還會(huì)影響到其他的輸出，就給控制系給控制系統(tǒng)的設(shè)計(jì)造成了很大的麻煩統(tǒng)的設(shè)計(jì)造成了很大的麻煩，

37、故在多變量，故在多變量控制系統(tǒng)的設(shè)計(jì)中就出現(xiàn)了解耦控制方法。控制系統(tǒng)的設(shè)計(jì)中就出現(xiàn)了解耦控制方法。68假設(shè)控制系統(tǒng)的狀態(tài)空間表達(dá)式為假設(shè)控制系統(tǒng)的狀態(tài)空間表達(dá)式為 (8-25) 其中其中 A:nn;B:nr;C:mn;D:mr引入狀態(tài)反饋引入狀態(tài)反饋 (8-26)其中其中 R為為r1參考輸入向量，在解耦控制中實(shí)際還應(yīng)參考輸入向量，在解耦控制中實(shí)際還應(yīng)要求要求r=m，亦即系統(tǒng)的輸入個(gè)數(shù)等于輸出個(gè)數(shù)，這時(shí)，亦即系統(tǒng)的輸入個(gè)數(shù)等于輸出個(gè)數(shù)，這時(shí)閉環(huán)系統(tǒng)的傳遞函數(shù)矩陣可以寫成閉環(huán)系統(tǒng)的傳遞函數(shù)矩陣可以寫成KxHRuHDBBKAsIDKCsRsYsG)()()()(1DuCxyBuAxx 69 若閉環(huán)系

38、統(tǒng)的若閉環(huán)系統(tǒng)的mr矩陣矩陣G(s)為對(duì)角的非奇異矩陣，則為對(duì)角的非奇異矩陣，則稱該系統(tǒng)是動(dòng)態(tài)解耦的系統(tǒng)稱該系統(tǒng)是動(dòng)態(tài)解耦的系統(tǒng)，若，若G(0)為對(duì)角非奇異矩為對(duì)角非奇異矩陣，且系統(tǒng)為穩(wěn)定的，則稱該系統(tǒng)是靜態(tài)解耦的。陣，且系統(tǒng)為穩(wěn)定的，則稱該系統(tǒng)是靜態(tài)解耦的。 在給定的控制結(jié)構(gòu)下，若系統(tǒng)的在給定的控制結(jié)構(gòu)下，若系統(tǒng)的D矩陣為矩陣為0，則，則閉環(huán)傳閉環(huán)傳遞函數(shù)陣遞函數(shù)陣G(s)可以簡(jiǎn)化成可以簡(jiǎn)化成 (8-28)BHBKAsICsRsYsG1)()()()(70 由上式可見(jiàn)，若由上式可見(jiàn)，若H矩陣為奇異矩陣，則矩陣為奇異矩陣，則G(s)矩陣必矩陣必為奇異的，所以為奇異的，所以為使得系統(tǒng)可以解耦，首

39、先應(yīng)該要求為使得系統(tǒng)可以解耦，首先應(yīng)該要求H為非奇異矩陣為非奇異矩陣。對(duì)于給定系統(tǒng)對(duì)于給定系統(tǒng),狀態(tài)方程可以狀態(tài)方程可以寫成為可控標(biāo)準(zhǔn)型寫成為可控標(biāo)準(zhǔn)型，故，故其中其中HsCBaBBKACCBsBKAsIsGnnn)(|1)(211.|011asasBKAsInnn71 首先這里將給出首先這里將給出能解耦的條件能解耦的條件：可以證明，若按下面：可以證明，若按下面方法生成的矩陣方法生成的矩陣B*為非奇異的為非奇異的，若取，若取H=(B*)-1，則由，則由前面給出的控制格式得出的系統(tǒng)能解耦原系統(tǒng)。前面給出的控制格式得出的系統(tǒng)能解耦原系統(tǒng)。 （8-29） 式中式中 C,C, Cm為為C矩陣的行向量矩

40、陣的行向量，參數(shù)，參數(shù)d,d,dm是在保證是在保證B為非奇異的前提下任選區(qū)間為非奇異的前提下任選區(qū)間0,n-1上的整數(shù)。若確定了上的整數(shù)。若確定了di參數(shù)參數(shù)，則可以直接獲，則可以直接獲得得解耦矩陣解耦矩陣BACBACBACBmdmdd21211111mdmdACACHK72例例8-11 對(duì)如下系統(tǒng)進(jìn)行解耦對(duì)如下系統(tǒng)進(jìn)行解耦解解 MATLAB程序?yàn)槌绦驗(yàn)镋xample8_11.mxyuxx11100133320130002000173執(zhí)行后可得執(zhí)行后可得H = 1.0000 0 1.0000 0 -1.3333 0.3333 -1.3333 0.3333K = -1.0000 0 0 -1.0

41、000 0 0 1.6667 1.3333 3.0000 1.6667 1.3333 3.0000n1 = 0 1.0000 -0.0000 -0.0000 0 1.0000 -0.0000 -0.0000 0 0.0000 0.0000 0.0000 0 0.0000 0.0000 0.0000d1 = 1.0000 -0.0000 -0.0000 0 1.0000 -0.0000 -0.0000 0n2 = 0 0 0 00 0 0 0 0 0.0000 1.0000 0 0 0.0000 1.0000 0d2 = 1.0000 -0.0000 -0.0000 0 1.0000 -0.00

42、00 -0.0000 074 亦即系統(tǒng)亦即系統(tǒng)解耦后解耦后的傳遞函數(shù)陣為的傳遞函數(shù)陣為 解耦控制系統(tǒng)的目的是將原模型變換成解耦的模解耦控制系統(tǒng)的目的是將原模型變換成解耦的模型，而并型，而并不必去考慮變換之后的響應(yīng)品質(zhì)不必去考慮變換之后的響應(yīng)品質(zhì)，因?yàn)轫?，因?yàn)轫憫?yīng)品質(zhì)這類問(wèn)題可以在解耦之后應(yīng)品質(zhì)這類問(wèn)題可以在解耦之后按照單變量系統(tǒng)進(jìn)按照單變量系統(tǒng)進(jìn)行設(shè)計(jì)補(bǔ)償行設(shè)計(jì)補(bǔ)償，單回路的設(shè)計(jì)當(dāng)然可以采用單變量系，單回路的設(shè)計(jì)當(dāng)然可以采用單變量系統(tǒng)的各種方法，例如可以采用超前滯后補(bǔ)償，統(tǒng)的各種方法，例如可以采用超前滯后補(bǔ)償，PI設(shè)設(shè)計(jì)以及計(jì)以及PID設(shè)計(jì)等，并能保證這樣設(shè)計(jì)出來(lái)的控制設(shè)計(jì)等，并能保證這樣設(shè)計(jì)

43、出來(lái)的控制器不會(huì)去影響其他回路。器不會(huì)去影響其他回路。1001001)(223ssssssG758.2.6 狀態(tài)估計(jì)器或觀測(cè)器狀態(tài)估計(jì)器或觀測(cè)器假設(shè)控制系統(tǒng)的狀態(tài)空間表達(dá)式為假設(shè)控制系統(tǒng)的狀態(tài)空間表達(dá)式為函數(shù)函數(shù)estim( )將生成下述狀態(tài)和輸出估計(jì)器將生成下述狀態(tài)和輸出估計(jì)器DuCxyBuAxx xICxyxCyLxAx)( 76 在在MATLAB中，函數(shù)中，函數(shù)estim( )的調(diào)用格式如下的調(diào)用格式如下est=estim(A,B,C,D,L) 其中其中 A,B,C,D為系統(tǒng)系數(shù)矩陣，為系統(tǒng)系數(shù)矩陣，L為為狀態(tài)估計(jì)增狀態(tài)估計(jì)增益矩陣益矩陣。狀態(tài)估計(jì)增益矩陣。狀態(tài)估計(jì)增益矩陣L可由極點(diǎn)配置

44、函數(shù)可由極點(diǎn)配置函數(shù)place( )形成形成，或者由，或者由Kalman濾波函數(shù)濾波函數(shù)kalman生成。利用以上命令可生成給定增益矩陣生成。利用以上命令可生成給定增益矩陣L下的下的狀態(tài)空間模型狀態(tài)空間模型A,B,C,D的輸出估計(jì)器的輸出估計(jì)器est。7778例例8-13 利用例利用例8-7所得的狀態(tài)觀測(cè)器的反饋陣所得的狀態(tài)觀測(cè)器的反饋陣L，求，求其系統(tǒng)的狀態(tài)估計(jì)器。其系統(tǒng)的狀態(tài)估計(jì)器。解解 MATLAB程序?yàn)槌绦驗(yàn)锳=0 1 0;0 0 1;-6 -11 -6;b=0;0;1;C=1 0 0; L=3;7;-1est=estim(A,b,C,0,L)執(zhí)行后得執(zhí)行后得est = -3.0000

45、 1.0000 0 -3.0000 1.0000 0 -7.0000 0 1.0000 -7.0000 0 1.0000 -5.0000 -11.0000 -6.0000 -5.0000 -11.0000 -6.0000798.2.7 系統(tǒng)控制器系統(tǒng)控制器假設(shè)控制系統(tǒng)的狀態(tài)空間表達(dá)式為假設(shè)控制系統(tǒng)的狀態(tài)空間表達(dá)式為利用函數(shù)利用函數(shù)reg( )可生成下述控制器可生成下述控制器DuCxyBuAxx xKuLyxKLDBLCAx)( 80 在在MATLAB中，函數(shù)中，函數(shù)reg( )的調(diào)用格式為的調(diào)用格式為est=reg(A,B,C,D,K,L) 其中其中 A,B,C,D為系統(tǒng)系數(shù)矩陣，為系統(tǒng)系數(shù)矩

46、陣，K為狀為狀態(tài)反饋增益矩陣，態(tài)反饋增益矩陣，L為狀態(tài)估計(jì)增益矩陣。為狀態(tài)估計(jì)增益矩陣。利用以上命令可生成給定狀態(tài)反饋增益利用以上命令可生成給定狀態(tài)反饋增益矩陣矩陣K及狀態(tài)估計(jì)增益矩陣及狀態(tài)估計(jì)增益矩陣L下的狀態(tài)空下的狀態(tài)空間模型間模型A,B,C,D的控制器的控制器est。假定系統(tǒng)。假定系統(tǒng)的所有輸出可測(cè)。的所有輸出可測(cè)。81例例8-14 利用例利用例8-7所得的狀態(tài)觀測(cè)器的反饋陣所得的狀態(tài)觀測(cè)器的反饋陣L，求，求其系統(tǒng)的控制器。假設(shè)狀態(tài)反饋陣其系統(tǒng)的控制器。假設(shè)狀態(tài)反饋陣K=-1 2 4。解解 MATLAB程序?yàn)槌绦驗(yàn)锳=0 1 0;0 0 1;-6 -11 -6;b=0;0;1;C=1 0

47、 0;K=-1 2 4;L=3;7;-1est=reg(A,b,C,0,K,L)執(zhí)行后得執(zhí)行后得 est = -3 1 0 -3 1 0 -7 0 1 -7 0 1 -4 -13 -10 -4 -13 -10828.3 最優(yōu)控制系統(tǒng)設(shè)計(jì)最優(yōu)控制系統(tǒng)設(shè)計(jì) MATLAB控制系統(tǒng)工具箱中也提控制系統(tǒng)工具箱中也提供了很多函數(shù)用來(lái)進(jìn)行系統(tǒng)的最優(yōu)供了很多函數(shù)用來(lái)進(jìn)行系統(tǒng)的最優(yōu)控制設(shè)計(jì)，相關(guān)函數(shù)如表控制設(shè)計(jì)，相關(guān)函數(shù)如表8-3所示。所示。83848.3.1 狀態(tài)反饋的狀態(tài)反饋的線性二次型最優(yōu)控制線性二次型最優(yōu)控制設(shè)線性定常系統(tǒng)的狀態(tài)空間表達(dá)式為設(shè)線性定常系統(tǒng)的狀態(tài)空間表達(dá)式為 (8-31)式中式中 A:nn

48、;B:nr;C:mn并設(shè)目標(biāo)函數(shù)為并設(shè)目標(biāo)函數(shù)為二次型性能指標(biāo)二次型性能指標(biāo) (8-32) 式中式中 Q(t)為為nn半正定實(shí)對(duì)稱矩陣，半正定實(shí)對(duì)稱矩陣，R(t)為為rr正定實(shí)對(duì)稱矩陣。一般情況下，假定這兩個(gè)矩陣正定實(shí)對(duì)稱矩陣。一般情況下，假定這兩個(gè)矩陣為定常矩陣，它們分別決定了系統(tǒng)暫態(tài)誤差與控為定常矩陣，它們分別決定了系統(tǒng)暫態(tài)誤差與控制能量消耗之間的相對(duì)重要性。制能量消耗之間的相對(duì)重要性。S為對(duì)稱半正定終為對(duì)稱半正定終端的加權(quán)陣，它為常數(shù)。端的加權(quán)陣，它為常數(shù)。)()()()()(tCxtytButAxtx f0d)()()()()()(21)()(21TTffTttttutRtutxtQt

49、xtSxtxJ85 當(dāng)當(dāng)x(tf)值固定時(shí)，則為終端控制問(wèn)題，特別是值固定時(shí)，則為終端控制問(wèn)題，特別是當(dāng)當(dāng)x(tf)0 時(shí)，則為調(diào)節(jié)器問(wèn)題；當(dāng)時(shí)，則為調(diào)節(jié)器問(wèn)題；當(dāng)t0 , tf均固均固定時(shí)，則為暫態(tài)過(guò)程最優(yōu)控制。定時(shí)，則為暫態(tài)過(guò)程最優(yōu)控制。 最優(yōu)控制問(wèn)題是為給定的線性系統(tǒng)式（最優(yōu)控制問(wèn)題是為給定的線性系統(tǒng)式（8-31）尋找一個(gè)最優(yōu)控制律尋找一個(gè)最優(yōu)控制律u*(t)，使系統(tǒng)從初始狀態(tài)，使系統(tǒng)從初始狀態(tài)x(t0)轉(zhuǎn)移到終端狀態(tài)轉(zhuǎn)移到終端狀態(tài)x(tf)，且滿足性能指標(biāo)式，且滿足性能指標(biāo)式（8-32）最小。它可以用變分法、極大值原理）最小。它可以用變分法、極大值原理和動(dòng)態(tài)規(guī)劃等三種方法中的任一種求解。這里和動(dòng)態(tài)規(guī)劃等三種方法中的任一種求解。這里我們采用極大值原理求解我們采用極大值原理求解u*(t)。86 MATLAB的控制系統(tǒng)工具箱中也提供了完整的解決的控制系統(tǒng)工具箱中也提供了完整的解決線性二次型最優(yōu)控制的函數(shù)，其中命令線性二次型最優(yōu)控制的函數(shù)，其中命令lqr( )和和lqry( )可以直接求解可以直接求解二次型調(diào)節(jié)器二次型調(diào)節(jié)器問(wèn)題及相關(guān)的問(wèn)題及相關(guān)的Riccati方程，它