基于隨機魯棒分析的輸出反饋特征結(jié)構(gòu)配置

金憲哲 吳森堂

(北京航空航天大學 自動化科學與電氣工程學院,北京 100191)

基于隨機魯棒分析的輸出反饋特征結(jié)構(gòu)配置

金憲哲 吳森堂

(北京航空航天大學 自動化科學與電氣工程學院,北京 100191)

針對輸出反饋特征結(jié)構(gòu)配置在參數(shù)不確定性系統(tǒng)設計中的魯棒性問題,提出一種基于隨機魯棒分析的輸出反饋特征結(jié)構(gòu)配置優(yōu)化方法.該方法通過隨機魯棒分析準確度量了閉環(huán)系統(tǒng)的魯棒性,確立控制系統(tǒng)設計要求與待設計參數(shù)間的直接聯(lián)系,并運用優(yōu)化技術實現(xiàn)閉環(huán)系統(tǒng)穩(wěn)定性與性能間的折衷,最大化控制系統(tǒng)的魯棒性.通過在某高超聲速飛行器橫側(cè)向解耦控制系統(tǒng)設計上的應用,驗證了該方法的有效性.

特征結(jié)構(gòu)配置;隨機魯棒分析;魯棒控制;反饋控制

特征結(jié)構(gòu)配置(EA,Eigenstructure Assignment)是一種通過極點配置,使閉環(huán)特征值及其相應的特征向量滿足指定的設計要求,以獲得期望的閉環(huán)系統(tǒng)響應的現(xiàn)代控制設計方法,其廣泛應用在飛行控制系統(tǒng)的靜態(tài)解耦設計上[1-3].由于實際系統(tǒng)中不是所有狀態(tài)都是可測的,當采用輸出反饋特征結(jié)構(gòu)配置時,因系統(tǒng)信息的缺失,在設計過程中無法獲知未配置特征值的變化情況,設計結(jié)果可能導致系統(tǒng)是不穩(wěn)定的.同時,外界干擾等因素引起的參數(shù)不確定性也會影響閉環(huán)系統(tǒng)的穩(wěn)定性.另外,運用 EA設計出的定常增益輸出反饋矩陣能否使系統(tǒng)的特征結(jié)構(gòu)在整個不確定性參數(shù)變化范圍內(nèi)滿足設計要求,也是需要重點關注的問題,這在進行解耦控制設計時尤為重要.目前,魯棒特征結(jié)構(gòu)配置已經(jīng)取得了許多研究成果[3-7],但在系統(tǒng)參數(shù)不確定性信息的利用、魯棒性優(yōu)化以及工程實用性等方面仍需要更進一步的探討.

魯棒控制的重要元素之一是分析控制系統(tǒng)魯棒性的技術,本文中引入隨機魯棒分析[8-9](SRA,Stochastic Robustness Analysis)對控制系統(tǒng)的魯棒性進行度量.SRA是一種用于評估線性時不變系統(tǒng)魯棒性的數(shù)值方法,通過對閉環(huán)系統(tǒng)特征值以及性能的蒙特卡羅估計,獲得不確定性參數(shù)統(tǒng)計意義下的閉環(huán)系統(tǒng)不穩(wěn)定概率和性能不滿足概率,作為閉環(huán)系統(tǒng)的魯棒穩(wěn)定性度量和魯棒性能度量.使用不穩(wěn)定概率和性能不滿足概率量化系統(tǒng)魯棒性的一個優(yōu)點是:它能夠明確地指出系統(tǒng)可允許多大的參數(shù)不確定性而不失魯棒性.SRA中二項置信區(qū)間的應用保證了魯棒性度量的精確性.

本文針對輸出反饋特征結(jié)構(gòu)配置在參數(shù)不確定性系統(tǒng)設計中的魯棒控制問題,提出一種基于SRA的特征結(jié)構(gòu)配置方法——SRA-EA.SRA-EA在設計過程中充分考慮了設計要求以及系統(tǒng)參數(shù)不確定性的信息,同時兼顧特征值及相應特征向量的變化情況,最大化地實現(xiàn)了控制系統(tǒng)的魯棒性.

1 特征結(jié)構(gòu)配置

考慮線性時不變系統(tǒng):

式中,x∈Rn,u∈ Rm,y∈ Rr和 q∈Rl分別為系統(tǒng)的狀態(tài)向量、輸入向量、輸出向量和系統(tǒng)參數(shù)向量;A(q),B(q)和 C(q)分別為系統(tǒng)矩陣、輸入矩陣和輸出矩陣,它們是關于系統(tǒng)參數(shù)向量 q的任意函數(shù).假設 q有已知的或估計的概率密度函數(shù)f(q),統(tǒng)計地描述了系統(tǒng)的參數(shù)不確定性.

對某一確定的 q0,采用輸出反饋

進行特征結(jié)構(gòu)配置,K∈Rm×r為增益輸出反饋矩陣,閉環(huán)特征值 λi(i=1,2,…,r)與相應的特征向量 vi(i=1,2,…,r)有如下關系:

式中,I為單位矩陣.由于缺少足夠的系統(tǒng)信息,輸出反饋特征結(jié)構(gòu)配置一般無法將 vi精確配置到期望的特征向量,實際得到的特征向量只是一種最可能的結(jié)果,它是由(λiI-A(q0))-1·B(q0)的列所張成子空間上的投影.有關的計算在文獻[10]中有詳細描述.

為便于計算,通過線性變換 x=Tx～,其中

Q為使 rank(T)=n的任意矩陣,將原系統(tǒng)

變換為

原系統(tǒng)的特征值與變換后系統(tǒng)的特征值一致,而相應的特征向量之間有如下關系:

由文獻[10],有

根據(jù)以上算法得到的 K,雖能夠精確配置閉環(huán)系統(tǒng)的 r個特征值,但缺少其余 n-r個特征值的信息而無法確定閉環(huán)系統(tǒng)是穩(wěn)定的.這種情況在系統(tǒng)存在參數(shù)不確定性,即 q按 f(q)指定的規(guī)律變化時,會變得更為復雜,因為 K只是在 q0的條件下得到的.此外,閉環(huán)系統(tǒng)特征向量也會因參數(shù)不確定性而偏離期望特征向量,閉環(huán)系統(tǒng)可能無法達到期望的性能.

2 SRA-EA

針對上述問題,本文在輸出反饋特征結(jié)構(gòu)配置過程中,運用 SRA設計對控制系統(tǒng)的魯棒穩(wěn)定性和魯棒性能進行分析,并在此基礎上對 K進行尋優(yōu),最大化閉環(huán)系統(tǒng)的魯棒性.

2.1 待設計參數(shù)向量的選取

若確定了特征向量各元素的構(gòu)成形式,期望閉環(huán)特征值與 K之間存在著某種映射關系.這里將需要配置的閉環(huán)特征值 λi(i=1,2,…,r)作為待設計參數(shù)向量 h,h∈Rr,即

期望閉環(huán)特征值的可允許變化范圍構(gòu)成了 h的設計參數(shù)空間.

K可根據(jù) h和標稱系統(tǒng)由上述算法得到.一般地,將 q變化范圍內(nèi)的 A(q)、B(q)和 C(q)的均值所構(gòu)成的系統(tǒng)設為標稱系統(tǒng).

2.2 魯棒穩(wěn)定性度量

對于某一確定的 K,若輸出反饋配置的閉環(huán)特征值 λi(i=1,2,…,r)在整個參數(shù)不確定性變化范圍內(nèi)滿足設計要求,同時,未配置的特征值λi(i=r+1,r+2,…,n-r)均具有負實部,則認為控制系統(tǒng)是魯棒穩(wěn)定的.

由于每次關于穩(wěn)定性的蒙特卡羅估計只會產(chǎn)生兩種結(jié)果:穩(wěn)定或不穩(wěn)定,因此,每次的蒙特卡羅估計可用二項隨機變量 X1表示.若閉環(huán)系統(tǒng)在一次蒙特卡羅估計中是魯棒穩(wěn)定的,則令 X1=0;否則,X1=1.閉環(huán)系統(tǒng)的不穩(wěn)定概率估值可近似為[9]

式中,J為蒙特卡羅估計次數(shù).J越大,p1越接近真實的不穩(wěn)定概率.

2.3 魯棒性能度量

若將閉環(huán)系統(tǒng)的性能指標轉(zhuǎn)化成二項決策(binomial decision)的形式,就可用類似地的方法得到閉環(huán)系統(tǒng)的性能不滿足概率估值,以描述控制系統(tǒng)的魯棒性能.一般地,系統(tǒng)的性能很難僅用一個指標描述其優(yōu)劣,常用的性能指標有調(diào)節(jié)時間、超調(diào)量以及特定應用要求(如解耦設計中特征結(jié)構(gòu)所要滿足的條件)等.控制系統(tǒng)的魯棒性能由各項魯棒性能度量綜合描述.

本文采用確定的性能可接受邊界作為此項性能的判據(jù),若閉環(huán)系統(tǒng)的某項性能在可接受邊界內(nèi),則相應的二項隨機變量 Xi=0(i≥2);否則,Xi=1.因此,閉環(huán)系統(tǒng)該項指標的性能不滿足概率估值為

2.4 二項置信區(qū)間的應用

由于 pi都是估計值,它接近真實不穩(wěn)定概率或性能不滿足概率 p-i的精度需要通過二項置信區(qū)間來保證.

由文獻[9]可知,關于不穩(wěn)定概率(或性能不滿足概率)P的置信水平為 1-α的置信區(qū)間為

文獻[9]還同時給出了在給定準確性系數(shù) ε和置信系數(shù) α情況下所需最少的 J,從而解決了蒙特卡羅估計內(nèi)在的計算問題.

判斷待設計參數(shù)向量 h1所對應的控制系統(tǒng)是否比 h2所對應的控制系統(tǒng)更加魯棒,可通過計算它們的不穩(wěn)定概率或性能不滿足概率之差來確定.設 P1為 h1所對應不穩(wěn)定概率或性能不滿足概率,P2為 h2所對應不穩(wěn)定概率或性能不滿足概率,并假設 P1＞P2,則兩個概率之差 ΔP?P1-P2的置信區(qū)間可根據(jù) P1和 P2的獨立置信區(qū)間給出[8]:

ΔP的置信區(qū)間可由下式給出:

若 ΔP的置信區(qū)間包含了 0,說明兩個控制系統(tǒng)在各自的估計次數(shù)和置信系數(shù)情況下,魯棒性上沒有顯著差異,即認為兩個控制系統(tǒng)具有相同的魯棒穩(wěn)定性或魯棒性能;否則,概率越小的控制系統(tǒng)其魯棒性越強.

2.5 魯棒性優(yōu)化設計

魯棒特征結(jié)構(gòu)配置的目的是設計常值反饋矩陣,使閉環(huán)系統(tǒng)的特征值及相應特征向量中指定元素不會因系統(tǒng)參數(shù)不確定性的影響而超出設計要求,從而保證系統(tǒng)的魯棒穩(wěn)定性和魯棒性能.

因此,魯棒特征結(jié)構(gòu)配置的問題可以描述為:在 h的設計參數(shù)空間內(nèi)尋找 h*,使得不穩(wěn)定概率以及各項性能不滿足概率最小.但是,魯棒控制系統(tǒng)設計的一個不可避免的問題是,系統(tǒng)的魯棒穩(wěn)定性與某些性能指標之間存在著某種折衷,對于此類折衷問題,可應用優(yōu)化技術實現(xiàn).

設優(yōu)化代價函數(shù)為

式中,M為關注的控制系統(tǒng)設計指標的個數(shù);pi為各項設計指標的不滿足概率;wi為 pi的加權,通過調(diào)整 wi可實現(xiàn)穩(wěn)定性與各項性能指標間的折衷.

判斷是否成功實現(xiàn)魯棒特征結(jié)構(gòu)配置設計,取決于待解決的實際問題.一般只要滿足下列條件之一即可:

1)pi最小;

2)在滿足可靠性要求的前提下,pi足夠小.

若獲得的 h*無法滿足設計要求,則可以根據(jù)pi的值修改相應的設計指標,或改變 h的設計空間范圍,重新進行設計.

3 設計算例

某高超聲速飛行器橫側(cè)向非線性狀態(tài)方程為

式中,β,φ,p,r,δa,δr分別為側(cè)滑角、滾轉(zhuǎn)角 、滾轉(zhuǎn)角速度、偏航角速度、升降副翼偏角和方向舵偏角 ;T,Y,α,θ,q,m,V分別為推力、側(cè)力 、迎角 、俯仰角、俯仰角速度、質(zhì)量和飛行速度;Iz分別為滾轉(zhuǎn)力矩、偏航力矩以及飛行器各軸的轉(zhuǎn)動慣量;Gya為重力在氣流坐標系 y軸上的投影;以上各參數(shù)遵循英美坐標系定義.Tδa,kδa,Tδr,kδr分別為升降副翼和方向舵的時間常數(shù)和增益;δac,δrc分別為升降副翼偏轉(zhuǎn)指令和方向舵偏轉(zhuǎn)指令.

應用 SRA-EA方法進行橫側(cè)向飛行控制系統(tǒng)設計,要求在飛行高度H∈[18000,25000]m,V∈[1700,2400]m/s,質(zhì)量 m∈ [121000,180000]kg的飛行包線內(nèi)采用輸出反饋特征結(jié)構(gòu)配置實現(xiàn)飛行器荷蘭滾模態(tài)與滾轉(zhuǎn)模態(tài)之間的解耦.輸出反饋為 β,φ,p,r.

設 H,V和 m按均勻分布進行攝動,待設計參數(shù)向量 h=[λ1λ2λ3λ4]中各元素分別對應閉環(huán)荷蘭滾模態(tài)(β,r)和滾轉(zhuǎn)模態(tài)(φ,p).h的設計參數(shù)空間為

各模態(tài)所對應的期望特征向量為

式中,×表示特征向量中未指定的分量.

為實現(xiàn)模態(tài)間的解耦,在應用 SRA-EA設計時的性能指標要求:期望閉環(huán)特征向量中指定為1或 -1的分量在整個飛行包線內(nèi)變化量不超過5%,分量為 0的變化不超過 ±0.05.

置信水平取 1-α=0.95,準確性系數(shù) ε=0.0001,可確定 J的最小值為 29955,增加額外的估計次數(shù)可進一步提高估計精度,最終設定 J=30000.不穩(wěn)定概率以及性能不滿足概率的期望值均為 1.2×10-3.

不穩(wěn)定概率的加權 w1=10,性能不滿足概率的加權 wi=1(i≥2).

用遺傳算法對 h進行尋優(yōu),最終得

高超聲速飛行器橫側(cè)向解耦控制系統(tǒng)在初始條件為 β=1°時的非線性隨機時間響應曲線如圖1所示.

圖 1 初始條件為 β=1°時的非線性隨機響應

初始條件為 φ=20°時的非線性隨機時間響應曲線如圖 2所示.

圖 2 初始條件為 φ=20°時的非線性隨機響應

仿真結(jié)果表明:所設計的控制系統(tǒng)能夠在給定的飛行包線內(nèi)實現(xiàn)荷蘭滾模態(tài)與滾轉(zhuǎn)模態(tài)間的解耦.

4 結(jié) 論

SRA-EA有效解決輸出反饋特征結(jié)構(gòu)配置在參數(shù)不確定性系統(tǒng)設計中存在的魯棒性問題.該方法充分利用了系統(tǒng)參數(shù)不確定性信息,設計結(jié)果對控制系統(tǒng)魯棒性的影響以及各設計指標間的折衷直接體現(xiàn)在設計中,便于設計人員了解和改進設計過程,有效避免產(chǎn)生過度保守的設計,具有較高的工程實用價值.

References)

[1]劉小剛,吳梅,安錦文.基于特征結(jié)構(gòu)配置/定量反饋理論的飛行控制系統(tǒng)設計研究[J].彈箭與制導學報,2006,26(1):340-343 Liu Xiaogang,Wu Mei,An Jinwen.Application of EA based QFT in design of lateral flight control system[J].Journal of Projectiles,Rockets,Missiles and Guidance,2006,26(1):340-343(in Chinese)

[2]李帆,周鳳岐,周軍.導彈基于特征結(jié)構(gòu)配置的輸出反饋解耦控制[J].彈箭與制導學報,2001,21(3):5-8 Li Fan,Zhou Fengqi,Zhou Jun.Decoupling control of missile by output feedback eigenstructure assignment[J].Journal of Projectiles,Rockets,Missiles and Guidance,2001,21(3):5-8(in Chinese)

[3]Seo Young Bong,Choi JaeWeon.Eigenstructure assignment considering probability of instability with flight control application[J].International Journal of Control,Automation,and Systems,2007,5(6):607-613

[4]王德軍,李元春.線性系統(tǒng)的魯棒特征結(jié)構(gòu)配置[J].吉林大學學報:信息科學版,2004,22(5):476-480 Wang Dejun,Li Yuanchun.Rob usteig enstru cture assign-ment for a class of linear systems[J].Journal of Jilin University:Information Science Edition,2004,22(5):476-480(in Chinese)

[5]潘常春,陳欣.基于魯棒特征結(jié)構(gòu)配置的無人機直接側(cè)力控制[J].飛行力學,2004,22(3):84-87 Pan Changchun,Chen Xin.UAV direct side force controllaw design with robust eigenstructure assignment[J].Flight Dynamics,2004,22(3):84-87(in Chinese)

[6]黃玲,段廣仁,于海華.廣義二階動力學系統(tǒng)的魯棒特征結(jié)構(gòu)配置[J].控制理論與應用,2009,26(3):238-242 Huang Ling,Duan Guangren,Yu Haihua.Robust eigen structure assignment in descriptor second-order linear systems[J].Control Theory＆Application,2009,26(3):238-242(in Chinese)

[7]Seo Young Bong,Choi Jae Weon,Man Hyung Lee.Eigen structure assignment for LTI systems with stochastic parameter variations[C]//Proceedings of the 2000 American Control Conference.Chicago:ACC,2000:3812-3816

[8]Ryan Laura Ray.Stochastic robustness of linear multivariable control systems:towards comprehens-ive robustness analysis[D].Princeton:Department of Mechanical and Aerospace Engineering,Princeton University,1991

[9]吳森堂.飛航導彈制導控制系統(tǒng)隨機魯棒分析與設計[M].北京:國防工業(yè)出版社,2010:8-23 Wu Sentang.Stochastic robustness analysis and design for guidance and control system of winged missile[M].Beijing:National Defense Industry Press,2010:8-23(in Chinese)

[10]Stevens Brian L,Lewis Frank L.Aircraft control and simulation[M].New York:Wiley,1992:347-353

(編 輯:劉登敏)

Output feedback eigen structure as signment based on stochastic robustness analysis

Jin Xianzhe Wu Sentang

(School of Automation Science and Electrical Engineering,Beijing University of Aeronautics and Astronautics,Beijing 100191,China)

An approach for output feedback eigenstructure assignment(EA)of linear multi variable system with parameter uncertainty was presented based on stochastic robustness analysis(SRA),which was used to solve the problem of the eigen structure robustness.With SRA,the robustness measures of closed-loop control system could be achieved accurately.And the direct relationship between the control system design specifications and the design parameters was established.The optimization algorithm was employed to maximize the control system robustness through trading off bet ween the stability and the performances.At last,SRA-EA was applied to the decoupled design for lateral control system of the hypersonic vehicle.The numerical simulations demonstrated the effectiveness of the presented approach.

eigenstructure assignment;robustness;robust control;feedback control

TP 13

A

1001-5965(2011)04-0487-05

2010-10-19

十一五國防基礎科研資助項目(A 212006×××)

金憲哲 (1973-),男,黑龍江佳木斯人,博士生,jinxianzhe_first@sina.com.