《現(xiàn)代控制理論》課后習(xí)題答案

1、現(xiàn)代控制理論參考答案第一章答案1-1試求圖1-27系統(tǒng)的模擬結(jié)構(gòu)圖，并建立其狀態(tài)空間表達(dá)式。圖1-271-1試求圖1-27系統(tǒng)的模擬結(jié)構(gòu)圖，并建立其狀態(tài)空間表達(dá)式。圖1-27系統(tǒng)方塊結(jié)構(gòu)圖解：系統(tǒng)的模擬結(jié)構(gòu)圖如下：系統(tǒng)的狀態(tài)方程如下：x=x12Kx=bx2J32KK1Kx=pxnx+x+px3J3J4J5J61111x=x43x=Kx+KX51316KKKx=1xix+4U6K1K6Kppp令0(s)=y,則y=xi所以，系統(tǒng)的狀態(tài)空間表達(dá)式及輸出方程表達(dá)式為所以，系統(tǒng)的狀態(tài)空間表達(dá)式及輸出方程表達(dá)式為10 x1x2x3JK2-Px4X5X600K3-110-K100KnJ010001J000

2、0KpJ01KK1-i-Kpx+x4x5X600000Ki-Kpy=100001-2有電路如圖1-28所示。以電壓u(t)為輸入量，求以電感中的電流和電容上的電壓作為狀態(tài)變量的狀態(tài)方程，和以電阻r2上的電壓作為輸出量的輸出方程。既得R既得R11x=-jx一x+u1L1L3L111R1x=-x+x2L2L32211x=x+x3C1C2y=Rx22解：由圖，令=X1，1=,Uc=X3，輸出量y=R2X2Rx+LX+x=uTOC o 1-5 h z1113有電路原理可知：Lx2+Rx=x2223x=x+C123寫成矢量矩陣形式為：11010X10X55x1x2x3R1L1R2L21L010 x1x2

3、X31-4兩輸入uU2兩輸出人u1y2的系統(tǒng)其模擬結(jié)構(gòu)圖如圖i-30所示試求其狀態(tài)空間表達(dá)式和傳遞函數(shù)陣。u2圖1-30雙輸入-雙輸出系統(tǒng)模擬結(jié)構(gòu)圖bzhr1V+0i00X10i00X100_xb02+1X003X1A0bc0a601aa43u42解：系統(tǒng)的狀態(tài)空間表達(dá)式如下所示：aa2ii00a5x1x2x3a(sIA)=11s+a10a5000a6s1aa4300-1-00_0ab061s100aa0b432s1as+aW(s)=(siA)-1B=2iux100aW(W(s)=C(sI-A)-1B=10uys-100-1-00_as+a0a0216-10s-100aaa0b1543121-

4、5系統(tǒng)的動態(tài)特性由下列微分方程描述(2)節(jié)+5y+7y+3y=u+3u+2u列寫其相應(yīng)的狀態(tài)空間表達(dá)式，并畫出相應(yīng)的模擬結(jié)構(gòu)圖。X2=y，X3=X2=y，X3=y則有_010_X1x2001+0-3-7-5X31y，ux1xx1x2x3相應(yīng)的模擬結(jié)構(gòu)圖如下：1-6,試求出系統(tǒng)的約旦標(biāo)準(zhǔn)型的實現(xiàn)并畫出相應(yīng)的模擬結(jié)構(gòu)圖(2)1-6,試求出系統(tǒng)的約旦標(biāo)準(zhǔn)型的實現(xiàn)并畫出相應(yīng)的模擬結(jié)構(gòu)圖(2)已知系統(tǒng)傳遞函數(shù)W(s)=6(s+1)s(s+2)(s+3)2-1016(s+1)4333解：W(s)=+丄s(s+2)(s+3)2(s+3)2s+3s+2s1-300-41000-20313當(dāng)當(dāng)九1=1時,1-7

5、1-7給定下列狀態(tài)空間表達(dá)式01.x(si01.x(si-A)-1=1(s+3)(s+2)(s+1)(s+3-2(s+3)-s-5s+3s(s+3)s-100(s+1)(s+2)010_X1-2-30X2+1-11-3X321x2&3x12uX3（1）畫出其模擬結(jié)構(gòu)圖（2）求系統(tǒng)的傳遞函數(shù)解：(2)W(s)=(sI-A)=2s+30|sI-A=s(s+3)2+2(s+3)=(s+3)(s+2)(s+1)1r(s+3s+30_OW(s)-(si-A)-1B-2(s+3)s(s+3)01ux(s+3)(s+2)(s+1)-s-5s-1(s+1)(s+2)21(s+3)(s+2)(s+1)(s+_(

6、s+3)1.s1(s+3)(s+2)(s+1)_(s+3)1.s(s+3)(2s+1)(s+3)W(s)=C(si-A)-1B=louy1(s+3)(s+2)(s+1)(2s+1)(s+2)(s+1)1-8求下列矩陣的特征矢量-010-(3)A-302-12-7-6-10-解：A的特征方程|Xi-A-3九-2兀+62+11九+60127九+6解之得：九=一1,九=一2,九=一312310_p11P1102p21P21-7-6P31P3103-12解得：p2i-/1一Pn令Pn=1P11-1_P-121P1解得：p2i-/1一Pn令Pn=1P11-1_P-121P1Q11-1P1（或令P11=-

7、1得P1_010_PP1212302P-2P2222-12-7-6P32P32當(dāng)入1=-2時,1令P12=2解得：P22=_2P12,P32=2P12P12-2_P22-4P11P212（或令p12=t得P222-212_010_PP1313302P-3P2323-12-7-6P33P33p32當(dāng)九1=-3時,解得：p23=-3P13,P33二3P13令P13二1得1-9將下列狀態(tài)空間表達(dá)式化成約旦標(biāo)準(zhǔn)型(并聯(lián)分解)1&41102X1-12)3y1121201-2_X132X2+273X531311X1x2u解：A解：A的特征方程九-4|XlA|=1-1-12九2=（九一1）（九3）2=01九

8、一3九=3,九=11,2341-2_p11P11102p21=3P211-13P31P31當(dāng)入廣3時,令P1廣1解之得叮P3廣P11p11丁p21=1p1O111P141-2_p11p11丁當(dāng)九2=3時，102p21=3p21+11-13Lp31Lp3111解之得p12=p+1,P令p=112P241-2_p13p13當(dāng)當(dāng)3=1時，LL102p23=p231-131p33p332222=p328-58-5-3解之得P13=0，P23=2P33令P33=1得1100-12_=102T-1=11-210101-10-12318-1T-1B=LL11-2LL27=LL-5201-1531-34112

9、0_110314_CT=0111LL102203110111-10已知兩系統(tǒng)的傳遞函數(shù)分別為W(s)和W2(s)1s+40試求兩子系統(tǒng)串聯(lián)聯(lián)結(jié)和并聯(lián)連接時，系統(tǒng)的傳遞函數(shù)陣，并討論所得結(jié)果解：（1）串聯(lián)聯(lián)結(jié)1111s+31s+40s+10s+2s+1_s+1s+2_W(s)=W(s)W(s)=211s2+5s+7(s+1)(s+3)(s+2)(s+3)(s+4)11(s+1)2(s+1)(s+2)2）并聯(lián)聯(lián)結(jié)_11-11_s+1s+2+s+3s+4s+11八00_s+2_s+1_W(s)=W(s)土W(s)=11其中子系統(tǒng)1、2的傳遞函數(shù)陣分別為1-11（第其中子系統(tǒng)1、2的傳遞函數(shù)陣分別為W

10、(s)W(s)=1_1s1s+2W2（s）二0l求系統(tǒng)的閉環(huán)傳遞函數(shù)解：10W(s)W(s)121I+W/s)W_10W(s)W(s)121I+W/s)W_1ss+3s+2b+W(s)W(s)11=二112s+31ss+2s+1W(s)=tl+W(s)W(s)11W(s)=121s+3s+1s(s+3)s+2s+31ss+21s+1s_1s1s+211s+1s_s+2s(s+3)110s+1_s+3_s+1s+3(s+2)(s+1)01-111-11（第2版教材）已知如圖1-22所示的系統(tǒng)，其中子系統(tǒng)1、2的傳遞函數(shù)陣分別為0-0-2W(s)=1W2（s）二I+W(s)W(s)=11tl+W(

11、s)W(s)li=11s(s+1)s2+5s+2IlW(s)=1W2（s）二I+W(s)W(s)=11tl+W(s)W(s)li=11s(s+1)s2+5s+2Il+W(s)W(s)liW(s)=s(s+1)s2+5s+2-2s(s+1)(s+2)2s(s+2)s(s+2)s2+5s+22(s+2)s+1一一+s+1(s+l)2(3s+8)(s+2)2(s2+5s+2)s3+6s2+6s(s+2)(s2+5s+2)s+1s2+5s+2s+2s2+5s+21-12已知差分方程為y(k+2)+3y(k+1)+2y(k)=2u(k+1)+3u(k)試將其用離散狀態(tài)空間表達(dá)式表示，并使驅(qū)動函數(shù)u的系數(shù)

12、b（即控制列陣）為1）b=解法1：2z+3z2+3z+211+z+1z+2-1x(k+1)=0 x(k)+11u(k)求系統(tǒng)的閉環(huán)傳遞函數(shù)解：111111s+1s101s+1s2101_21s+2_s+2_11s+210sW(s)W(s)=11y(k)=&11x(k)解法2：11x(k+1)=x(k)12x(k+1)=-2x(k)-3x(k)+u212y(k)=3x(k)+2x(k)12x(k+1)=0-21-3x(k)+01u(k)y(k)=b2x(k)求T,使得T-iB=T-1AT=1101-2-3得T-1=01-4-5所以-1C=b2：0b-1所以，狀態(tài)空間表達(dá)式為11u(11u(k)z

13、(k+1)=51z(k)+51y(k)=b1z(k)第二章習(xí)題答案2-4用三種方法計算以下矩陣指數(shù)函數(shù)eAt。2)A=解：第一種方法：令解：第一種方法：令卜IA|=0九-1九-1-4-1九一1=0，即C-1)2-4=0。當(dāng)九=3當(dāng)九=3時,1特征矢量p1=_1p3p11=1141p213p21P11P21由A/1耳，得11當(dāng)當(dāng)21時特征矢量P：由AP21pP121241P22P2221p12P22i2P2，得可令P2=p12+P22一P124p+p1222一P2222T1eAt=211-1111_e3t+ete3tet0242244et1111e3t+ete3t+_24_2224e3t0第二種

14、方法，即拉氏反變換法rsiisiA=4s1sisiA1=1(s-3)(s+1)1s1ss1(s-3)(s+1),4,、(s-3)(s+1)1(s3)(s+1)s1(s3)(s+1)eAt=LeAt=L1r(siA)-111e3t+et22e3tet11e3tet4411e3t+et221r11)1r11+2js3s+1丿4js3s+1丿111r11)一+s3s+12js3s+1丿第三種方法，即凱萊哈密頓定理由第一種方法可知入-3，幾2=13-1e3t1et1113_-13_一e3t+et4e3t441et111-e3tet4L441414(13)10-(13)1eAt=e3t+e-t+e3t+

15、e-1144丿_01144丿4111e3t+e-t22e3t-e-t11e3te-14411e3t+e-t222-5下列矩陣是否滿足狀態(tài)轉(zhuǎn)移矩陣的條件，如果滿足，試求與之對應(yīng)的A陣。(t)=解：(3)因為4)2e-t-e-2te-t-e-2t(t)=解：(3)因為4)2e-t-e-2te-t-e-2t(0)=2e-2t-2e-t2e-2t-e-t(4)(t)=卩(e-1+e3t)2-e-1+e3t=I,所以該矩陣滿足狀態(tài)轉(zhuǎn)移矩陣的條件2e-t+2e-2t4e-2t+2e-1-0-2_t=0e-1+2e-2t4e-2t+e-t_1-3_01t=010所以該矩陣滿足狀態(tài)轉(zhuǎn)移矩陣的條件A=(&(t)

16、因為(0)=A=(&(t)13131et+e3te-t+e3t2244e-t+3e3t13et+e3t2201t=0t=01(-e-t+e3t(e-t+e3t2)2-6求下列狀態(tài)空間表達(dá)式的解：01101x+001y=(1,0)x初始狀態(tài)x(0)=H，輸入u(t)時單位階躍函數(shù)。解：s-1siA=0s(sI(sI-A)-1=-1s(t)=eAt=L-1(si-A)-111t1一一01因為u(t)=I(t)因為x(t)=Q(t)x(0)+jt(tT)BuG)dT0+jt+jt0tT1dT1t-丁+jt1tT_OdT01100111t+1-t21+21-1t12+t+12t+1y0 x=212+1

17、+12-9有系統(tǒng)如圖2.2所示試求離散化的狀態(tài)空間表達(dá)式。設(shè)采樣周期分別為Ts和1s,而ui和u2為分段常數(shù)。圖2.2系統(tǒng)結(jié)構(gòu)圖解：將此圖化成模擬結(jié)構(gòu)圖列出狀態(tài)方程x&=kux111x&=xu212y=x+2x21-10k0_ux=x+1100-1u11112y=2lXix2則離散時間狀態(tài)空間表達(dá)式為x(k+1)=G(T)x(k)+H(T)u(k)y(k)=cx(k)+Du(k)-10k0-2A=B=CT=100-11由G-10k0-2A=B=CT=100-11由G(T)=eAt和H(T)=JTeAtdtB得:0s+10e-Tr(si-a)-ieAt=L-1e-t0_k0一1-e-t0k0一k

18、(1-e-t)0_1-e-T1dt_0-1T-1+e-TT_0-1k(T-1+e-T)-T=L-11s000e-1x(k)+x(k+1)=H=JTeAtdt=JT0k(1-e-1)當(dāng)T=1時1-e-1u(k)ke-1-1y(k+1)=21x(k)當(dāng)T=0.1時x(k+當(dāng)T=0.1時x(k+1)=e-0.11-e-0.1x(k)+k(1-e-0.1)0k1-0.9)-0.1u(k)y(k+1)=21x(k)第三章習(xí)題3-1判斷下列系統(tǒng)的狀態(tài)能控性和能觀測性。系統(tǒng)中a,b,c,d的取值對能控性和能觀性是否有關(guān)，若有關(guān)，其取值條件如何？（1）系統(tǒng)如圖3.16所示：xuxu解：由圖可得：狀態(tài)空間表達(dá)式

19、為：圖3.16系統(tǒng)模擬結(jié)構(gòu)圖x=-ax+u11x=-解：由圖可得：狀態(tài)空間表達(dá)式為：圖3.16系統(tǒng)模擬結(jié)構(gòu)圖x=-ax+u11x=-bx22x=-cx+x+x=x+x-cx332x=x-dx43y=x3x1x2=x3x4-y=00-a000_x11丁0-b00 x02+11-c0 x03001-dx40uox1由于x、x、x與u無關(guān),234能觀的，為不能觀系統(tǒng)。因而狀態(tài)不能完全能控，為不能控系統(tǒng)。由于y只與x3有關(guān)，因而系統(tǒng)為不完全3）系統(tǒng)如下式：解：如狀態(tài)方程與輸出方程所示，A解：如狀態(tài)方程與輸出方程所示，A為約旦標(biāo)準(zhǔn)形。要使系統(tǒng)能控，控制矩陣b中相對于約旦塊的最后一行x-110_x-211

20、x=0-10 x+a022x00-2xb033u兀素不能為0,故有a豐0,b豐0。要使系統(tǒng)能觀，則C中對應(yīng)于約旦塊的第一列元素不全為0,故有CH0,d豐0。3-2時不變系統(tǒng)X=X=u-31-X+j1_1-311j1-X1-1y=試用兩種方法判別其能控性和能觀性。解：方法一：M=IbM=IbabL-2-2-2-2-31-1111_A=,B=,C=1-3111-1rankM=12,系統(tǒng)不能控。CCA1CCA11-2-41-1-24rankN=2,系統(tǒng)能觀。方法二：將系統(tǒng)化為約旦標(biāo)準(zhǔn)形。|xi|xiA=(X+31=0九=2,X=412則狀態(tài)矢量：AP=XPnP=11111AP22=AP22=XPnP

21、=2221-111T=11T=1-1T-1=T-1AT=-311T-1AT=-311-3111-1-20-411100T-1B=2110021111-11-1CT=2002T-陀中有全為零的行，系統(tǒng)不可控。CT中沒有全為0的列，系統(tǒng)可觀。3-3確定使下列系統(tǒng)為狀態(tài)完全能控和狀態(tài)完全能觀的待定常數(shù)a和Bii11112-1解：構(gòu)造能控陣：a+11a2要使系統(tǒng)完全能控，則a+1豐a，即a-a+1h01212構(gòu)造能觀陣：N=CN=CCA-1要使系統(tǒng)完全能觀，則1-aH-a，即aa+1h021123-4設(shè)系統(tǒng)的傳遞函數(shù)是y(s)_s+au(s)s3+10s2+27s+18當(dāng)a取何值時，系統(tǒng)將是不完全能控

22、或不完全能觀的？當(dāng)a取上述值時，求使系統(tǒng)的完全能控的狀態(tài)空間表達(dá)式。當(dāng)a取上述值時，求使系統(tǒng)的完全能觀的狀態(tài)空間表達(dá)式。解：(1)方法1:W(s)_凹_su(s)(s+1)(s+3)(s+6)系統(tǒng)能控且能觀的條件為W(s)沒有零極點(diǎn)對消。因此當(dāng)a=1，或a=3或a=6時，系統(tǒng)為不能控或不能觀。方法2：a-1a3a-6y(s)s+a106丄_丄u(s)(s丄1)(s丄3)(s丄6)s丄1s丄3s丄6九_1，九_3，九_6123-1-100-丁030X丄100-61系統(tǒng)能控且能觀的條件為矩陣C不存在全為0的列。因此當(dāng)a=1，或a=3或a=6時，系統(tǒng)為不能控或不能觀。(2)當(dāng)a=1，a=3或a=6時

23、，系統(tǒng)可化為能控標(biāo)準(zhǔn)I型uu-010-Ox=001x+0-18-27-101y=Da10 x(3)根據(jù)對偶原理，當(dāng)a=1,a=2或a=4時，系統(tǒng)的能觀標(biāo)準(zhǔn)II型為00-18-a10-27x+101-100u0lxy=t)3-6已知系統(tǒng)的微分方程為：+6y+11y+6y=6u試寫出其對偶系統(tǒng)的狀態(tài)空間表達(dá)式及其傳遞函數(shù)。解：a=6，a=11，a=6，a=3，b=601230系統(tǒng)的狀態(tài)空間表達(dá)式為010-x=001x+0-6-11-61y=D600 xu傳遞函數(shù)為s-10一-10s-10611s+610W(s)=C(sI-A)-1B=k0s3+6s2+11s+6其對偶系統(tǒng)的狀態(tài)空間表達(dá)式為：00-

24、6-610-11x+001-60u0lxy=D)傳遞函數(shù)為W(s)=s3-6s2一11s+63-9已知系統(tǒng)的傳遞函數(shù)為W(s)=斗s2+4s+3試求其能控標(biāo)準(zhǔn)型和能觀標(biāo)準(zhǔn)型。解:W(s)=s2+6s+8=1+2s+5s2+4s+3s2+4s+301-0_x+_-3-41系統(tǒng)的能控標(biāo)準(zhǔn)I型為uy=k2L+u能觀標(biāo)準(zhǔn)II能觀標(biāo)準(zhǔn)II型為u0100 x=-2-30X+1-11-32y=t001Xu0-35_X+1-42x=+u0-35_X+1-42x=+uy=t)11c3-10給定下列狀態(tài)空間方程，試判別其是否變換為能控和能觀標(biāo)準(zhǔn)型。M=bAbA2bL1-2-5-3711系統(tǒng)為不能控系統(tǒng)，不能變換為

25、能控標(biāo)準(zhǔn)型。c-001_N=CA=-1-1-3CA21-79rankM=23,rankN=3，系統(tǒng)為能觀系統(tǒng)，可以變換為能觀標(biāo)準(zhǔn)型。010-2-30，b=1-11-32解：A=c=bo111試將下列系統(tǒng)按能控性進(jìn)行分解12-廠010,b=00-4311)A=,C=1-11解：M=bAbA2b=-409rankM=23，系統(tǒng)不是完全能控的。0-100，R=Ab=0，R=123130=b=，其中R是任意的，只要滿足R滿秩。3c0-10301001得R-1=-100130c010構(gòu)造奇異變換陣R：cR1即R=cAA=R-1AR=cc3-121)_0-32-丁14-2b=R-1b=0001c0試將下列

26、系統(tǒng)按能觀性進(jìn)行結(jié)構(gòu)分解c=cR=12-1c12-1o010,b=00-431A=,C=1-11C_1-1則有N=CA=2-32CA24-74解：由已知得A=rankN=2v3,該系統(tǒng)不能觀12-12-1O010,b=00-431,C=1-11構(gòu)造非奇異變換矩陣R-1,0-1-30-101=cRox=h0氐y=3-13=cRox=h0氐y=3-13試將下列系統(tǒng)按能控性和能觀性進(jìn)行結(jié)構(gòu)分解100丁1)A=223,b=2-2012,C=11解：由已知得M=AAbAb221120126-2010-230%+2-7321uX=R-1AR張+R-1bu=000rankM=3，則系統(tǒng)能控c112_N=cA

27、=-125cA2-7411則系統(tǒng)能觀rankN=3,所以此系統(tǒng)為能控并且能觀系統(tǒng)77311取T=212c220126-2，則T-i=7c2-3141_2144-354002-10-5，B=T-1b=0c20140則A=c=cT3-14求下列傳遞函數(shù)陣的最小實現(xiàn)。c2=71323解：11a=1，B=，A=0011c1011-B=，C=，D=c_01c11cws+111-10-10000系統(tǒng)能控不能觀(s)=1)10-00C=CR=，D=c01000=1，所以最小實現(xiàn)為Am(s)=1)10-00C=CR=，D=c01000=1，所以最小實現(xiàn)為Am=11，100，D=1m00Cm/I八1/驗證：CU

28、i-A戶Bmmm=w(s)3-15設(shè)Z和Z是兩個能控且能觀的系統(tǒng)1201-0z：A=，b=11-3-411=-2，A2b=1，2C=12C=211z：21)試分析由Z和Z所組成的串聯(lián)系統(tǒng)的能控性和能觀性，并寫出其傳遞函數(shù);122)試分析由z和z所組成的并聯(lián)系統(tǒng)的能控性和能觀性，并寫出其傳遞函數(shù)。12解：1)Z1和Z2串聯(lián)111-，則R=01001取R-1=0-10-1R-1AR=，B=R-1B=000-10c01所以A=丫丫s+2當(dāng)為的輸出y是為的輸入u時，X=-2x+2x+x11223312_01_01o-0_-3-40 x+121-20y=00lx01-4M=bAbA2b=1-41301-

29、4則rankM=2v3，所以系統(tǒng)不完全能控。W(sW(s)=C(sI-A)-1B=(s+2)(s+3)(s+4)s2+7s+12當(dāng)2得輸出y是2的輸入u時221因為M=bAbA2b因為M=bAbA2b011-0_-3-41x+000-21u，y=21ox01-21-6-4rankN=23rankM=3則系統(tǒng)能控rankN=23c210_cA=-3-21cA2654因為N=則系統(tǒng)不能觀W(s)=C(si-A)-1B=s2+7s+12(2)工和工并聯(lián)12因為rankM=3，010-因為rankM=3，010-0_-3-40 x+100-21M=AAbAb2u，所以系統(tǒng)完全能控y=211x1-4-2

30、-413-4c211_N=cA=-3-2-2cA2654因為rankN=3，所以系統(tǒng)完全能觀w(w(s)=C(sI-A)-1B=f近2s+2+2八U1Xs+2)(s+3)現(xiàn)代控制理論第四章習(xí)題答案1判斷下列二次型函數(shù)的符號性質(zhì)：(1)Q(X)=X23X211x2+2XXXX2XX12312231313(2)v(x)=x2+4x2+x22xx6xx2xx13123解：(1)由已知得Q(x)=X3X1XX11x1X1X2X3x33x2X1111A=10，A1233=3710,A1602X33A2因此Q(X)不是正定的4-2已知二階系統(tǒng)的狀態(tài)方程x1x11TOC o 1-5

31、 h z_（aa&=1112x、aa.2122試確定系統(tǒng)在平衡狀態(tài)處大范圍漸進(jìn)穩(wěn)定的條件。解：方法（1）：要使系統(tǒng)在平衡狀態(tài)處大范圍漸進(jìn)穩(wěn)定，則要求滿足A的特征值均具有負(fù)實部。即：尢一a尢一aiia21a12九一a22=九2(a+a)九+aaaa112211221221有解，且解具有負(fù)實部。即：a+aaa112211221221方法（2）：系統(tǒng)的原點(diǎn)平衡狀態(tài)x=0為大范圍漸近穩(wěn)定，等價于AtP+PA=Q。e取Q=I，令P取Q=I，令P=P11P12P12P22，則帶入ATP+PA=Q，2a2a0_一P一-112111aa+aaP=0121122211202a2aP1112222211得到2a1

32、12a11a1202a21a+a11222a120a=4(a+a)(aaaa)豐0，則此方程組有唯一解。即211122112212212a22P=2(a+P=2(a+a1122)|A_(aIaI+a2+a22122a+aa)12222111(aa+aa)12222111IaI+a2+a21112其中detA=IAI=aaaa11221221要求P正定，則要求A+a2+a2A=P=212201112(a+a)|A1122(a+a)(a+a)2+(acc4(a+a)1122a)2210因此a+a011224-3試用lyapunov第二法確定下列系統(tǒng)原點(diǎn)的穩(wěn)定性1_x3解：(1)系統(tǒng)唯一的平衡狀態(tài)是x=0。選取Lyapunov函數(shù)為V(x)=x2+x20，則e12V(x)=2xxx+2xxx1122=2x(x+2x)+2x(2x3x)112212=2x2+6xx6x2112233=2(xx)2x20，則e12V(x)=2xxx+2xxx1122=2x(x+x)+2x(xx)112212=2x22x202221試確定平衡狀態(tài)的穩(wěn)定性。解：若采用克拉索夫斯基法，則依題意有：f(f(x)=2a(1+x)2xx21001a4ax3ax22Q(x)=Jt(x)+J(x)11a4ax3ax22202a8ax6ax2220