基于LMI的連續(xù)時(shí)間線性時(shí)滯系統(tǒng)有限頻模型降階

杜鑫，范培兵，劉夫偉

（1.上海大學(xué)機(jī)電工程與自動(dòng)化學(xué)院，上海 200072；

2.華為技術(shù)有限公司，上海 200040）

基于LMI的連續(xù)時(shí)間線性時(shí)滯系統(tǒng)有限頻模型降階

杜鑫1，范培兵1，劉夫偉2

（1.上海大學(xué)機(jī)電工程與自動(dòng)化學(xué)院，上海 200072；

2.華為技術(shù)有限公司，上海 200040）

針對(duì)連續(xù)時(shí)間線性時(shí)滯系統(tǒng),研究了在系統(tǒng)工作頻率范圍為已知有限區(qū)間情形下的有限頻模型降階問題.通過(guò)引入有限頻域內(nèi)誤差傳遞函數(shù)的最大奇異值函數(shù)作為指標(biāo)函數(shù),對(duì)模型逼近性能進(jìn)行了刻畫,進(jìn)而結(jié)合一些基礎(chǔ)性的矩陣不等式技術(shù)和線性時(shí)滯系統(tǒng)性能進(jìn)行分析,得到了保持降階模型穩(wěn)定性的有限頻模型逼近性能優(yōu)化設(shè)計(jì)條件,這些條件以線性矩陣不等式（linear matrix inequalities,LMIs）的形式表示,易于檢驗(yàn)和數(shù)值求解.最后,算例驗(yàn)證了結(jié)果的有效性.

模型降階；有限頻域；線性時(shí)滯系統(tǒng)；線性矩陣不等式

建立一個(gè)描述系統(tǒng)動(dòng)力學(xué)的數(shù)學(xué)模型是對(duì)許多工程實(shí)際系統(tǒng)進(jìn)行仿真、分析及綜合利用的必要前提.對(duì)于一些儲(chǔ)能元件個(gè)數(shù)較多的動(dòng)態(tài)系統(tǒng)如電路系統(tǒng)、機(jī)械系統(tǒng)等來(lái)說(shuō)，通常通過(guò)機(jī)理分析等方法得到的系統(tǒng)精確數(shù)學(xué)模型的階數(shù)較高，從而給后續(xù)的系統(tǒng)仿真、分析帶來(lái)了數(shù)值計(jì)算及物理實(shí)現(xiàn)上的困難.針對(duì)這類問題，一個(gè)廣泛采用的解決方案是對(duì)原模型進(jìn)行模型降階［1-3］，即通過(guò)數(shù)學(xué)方法求取一個(gè)盡可能逼近原系統(tǒng)在工作頻率范圍內(nèi)的動(dòng)態(tài)響應(yīng)特性的降階模型.在過(guò)去的30多年時(shí)間里，模型降階問題受到了包括應(yīng)用數(shù)學(xué)、控制理論等領(lǐng)域的廣泛關(guān)注.常見的模型降階方法有基于矩匹配的有理插值方法［1,3］、基于奇異值分解的平衡截?cái)喾椒ǎ?-3］及各類數(shù)值優(yōu)化方法［4］等.

時(shí)滯是一類常見的工程現(xiàn)象，常出現(xiàn)在網(wǎng)絡(luò)控制系統(tǒng)、傳輸線系統(tǒng)、化工過(guò)程系統(tǒng)等［5-6］實(shí)際工程系統(tǒng)中.在過(guò)去20年時(shí)間中對(duì)時(shí)滯系統(tǒng)各類問題的研究得到了廣泛的關(guān)注［7-16］，其中時(shí)滯系統(tǒng)的模型降階問題有著重要的工程應(yīng)用背景，如RLC電路內(nèi)部鏈接傳輸線系統(tǒng)的建模與仿真等.因此關(guān)于時(shí)滯系統(tǒng)的模型降階研究已成為時(shí)滯系統(tǒng)研究領(lǐng)域的一個(gè)重要分支，相應(yīng)地出現(xiàn)了一些結(jié)果和報(bào)道，如文獻(xiàn)［7］給出了基于矩匹配法的時(shí)滯系統(tǒng)模型降階算法，這類算法使得降階模型在給定插值頻率點(diǎn)處及其鄰域內(nèi)都有很好的逼近效果，但是對(duì)逼近誤差進(jìn)行分析估計(jì)比較困難；文獻(xiàn)［8］從標(biāo)準(zhǔn)線性系統(tǒng)格拉姆矩陣的定義出發(fā)，給出了一類時(shí)滯相關(guān)的基于位置格拉姆矩陣的平衡截?cái)嗄Ｐ徒惦A方法；文獻(xiàn)［9-10］從格拉姆矩陣的數(shù)值計(jì)算方法出發(fā)，給出了一類快速的基于估計(jì)格拉姆矩陣的平衡截?cái)嗄Ｐ徒惦A方法，但這兩類時(shí)滯系統(tǒng)的平衡截?cái)喾椒ǘ紵o(wú)法提供標(biāo)準(zhǔn)線性系統(tǒng)平衡截?cái)嗬碚撝械哪Ｐ捅平`差來(lái)估計(jì)上界，也無(wú)法保持原系統(tǒng)的穩(wěn)定性.除這兩類主要應(yīng)用于大規(guī)模系統(tǒng)快速模型降階的方法之外，針對(duì)相對(duì)系統(tǒng)模型階數(shù)不太高的模型降階則主要考慮采用數(shù)值優(yōu)化方法［11-16］，該方法大多基于近年來(lái)在各類控制系統(tǒng)相關(guān)問題研究中廣泛采用的線性矩陣不等式優(yōu)化方法［17］.但是，現(xiàn)有的結(jié)果中均采用全頻范圍內(nèi)的逼近性能指標(biāo)，而該指標(biāo)只能對(duì)全頻域內(nèi)的逼近誤差進(jìn)行優(yōu)化［12-16］.另外，許多實(shí)際工程系統(tǒng)的工作頻率往往為一個(gè)已知的有限頻率區(qū)間［18-21］，現(xiàn)有的結(jié)果在逼近誤差上就會(huì)顯得比較保守，特別是在頻率范圍較小的情況下，基于全頻性能指標(biāo)而得到的工作頻率范圍內(nèi)的逼近性能往往是不夠理想的.

本研究考慮已知工作頻率范圍情形下的連續(xù)時(shí)間線性時(shí)滯系統(tǒng)的模型降階問題，首先引入一個(gè)有限頻誤差系統(tǒng)傳遞函數(shù)的最大奇異值函數(shù)來(lái)作為模型在工作頻率范圍內(nèi)逼近精度的衡量指標(biāo)，然后結(jié)合投影引理、Schur補(bǔ)引理等基本的不等式工具以及一些已有的線性時(shí)滯系統(tǒng)有限頻性能分析結(jié)果［22-24］，推導(dǎo)出相應(yīng)的可優(yōu)化有限頻逼近性能的線性矩陣不等式設(shè)計(jì)條件，同時(shí)還給出了保持降階模型穩(wěn)定性的設(shè)計(jì)條件，進(jìn)而將所考慮的模型降階問題轉(zhuǎn)化為一組基于線性矩陣不等式（linear matrix inequalities，LMIs）的優(yōu)化問題，可方便地通過(guò)Matlab中的LMIs求解工具箱進(jìn)行相應(yīng)的數(shù)值求解.最后，通過(guò)數(shù)值算例驗(yàn)證了本方法的有效性.

1 問題描述

考慮如下方程描述的連續(xù)時(shí)間線性時(shí)滯系統(tǒng)：

式中，x（t）∈Rn為系統(tǒng)狀態(tài)向量，u（t）∈Rl為系統(tǒng)輸入向量，y（t）∈Rm是系統(tǒng)輸出向量，A，Ad，B，C，Cd，D為系統(tǒng)參數(shù)矩陣，d＞0為未知的定常時(shí)延，假設(shè)其上界為已知常數(shù)該系統(tǒng)在頻域內(nèi)的表現(xiàn)形式為

式中，U（jω），X（jω）和Y（jω）分別為信號(hào)u（t），x（t）和y（t）經(jīng)傅里葉變換后的頻域表達(dá)形式.系統(tǒng)輸入輸出信號(hào)之間的關(guān)系可進(jìn)一步表達(dá)為Y（jω）=G（jω）U（jω），其中G（jω）為系統(tǒng)的傳遞函數(shù)，

本研究討論由方程（1）或（2）描述的線性時(shí)滯系統(tǒng)有限頻模型降階問題，即求取一個(gè)階數(shù)低于給定系統(tǒng)的降階系統(tǒng)：

式中，xr（t）∈Rr為降階系統(tǒng)的低維狀態(tài)向量，u（t）∈Rl為和原系統(tǒng)相同的輸入向量，yr（t）∈Rm為降階系統(tǒng)的輸出向量，Ar,Adr,Br,Cr,Cdr,Dr為待求的參數(shù)矩陣.同樣，降階系統(tǒng)也可以在頻域上進(jìn)行描述，其輸入輸出的頻域關(guān)系為Yr（jω）=Gr（jω）U（jω），其中Gr（jω）為降階系統(tǒng)的傳遞函數(shù)，

且該降階系統(tǒng)與原系統(tǒng)在給定頻率范圍內(nèi)的誤差逼近性能滿足

以上文章并非寫出履歷的事實(shí)，而是從“追求真理”的角度整理自己的思想問題。而從書名可知，為使佛教相關(guān)人士注意，此文章包含偏差及強(qiáng)調(diào)。敘述自誕生至開始修習(xí)洋學(xué)的前半部分相當(dāng)于長(zhǎng)岡時(shí)期的事，后半部分則相當(dāng)于東京大學(xué)時(shí)期的事。

來(lái)作為模型逼近性能的刻畫指標(biāo)，通過(guò)設(shè)計(jì)相應(yīng)的優(yōu)化算法使得該指標(biāo)函數(shù)極小，進(jìn)而給出相應(yīng)的降階系統(tǒng)模型.

對(duì)此，我是不敢茍同的。有次，我實(shí)在沒忍住，跟她分享了自己的想法：“你總覺得自己是個(gè)姑娘，所以不該為生活奔波，理應(yīng)坐享其成。可是生活從來(lái)不會(huì)因?yàn)槟闶枪媚锞蜁?huì)對(duì)你格外開恩。”

引理2（投影引理） 給定Γ,Λ,Φ，存在矩陣F滿足ΓFΛT+（ΓFΛT）T+Φ＜0的充分必要條件是

2 主要結(jié)果

首先給出滿足有限頻域內(nèi)性能指標(biāo)的降階系統(tǒng)模型參數(shù)矩陣設(shè)計(jì)條件，然后給出保持降階模型穩(wěn)定性的設(shè)計(jì)條件，根據(jù)得出的這兩個(gè)設(shè)計(jì)條件，降階模型的參數(shù)矩陣可以方便地通過(guò)求解一組線性矩陣不等式優(yōu)化問題得出.

“直開臨杪”之“杪”，《全宋詞》有一缺字符，而宋刻本有。依律“凡間謫墮”前缺一字，《全宋詞》有一空方(缺字符)，而宋刻本連書不空。

定理1 考慮方程（1）給出的連續(xù)時(shí)間線性時(shí)滯系統(tǒng)，其中系統(tǒng)的工作頻率范圍ω∈［ω1,ω2］，如果存在對(duì)稱正定矩陣Q11,Q22,Z11,Z22，對(duì)稱矩陣P11,P22,P11,P22以及矩陣使得下列線性矩陣不等式成立：

一是檢測(cè)+宣傳：營(yíng)造良好氛圍。一要堅(jiān)持強(qiáng)化謀劃布局，提升食品安全宣傳資源整合能力、創(chuàng)新線上線下宣傳模式，形成以上帶下、以下促上，上下共同發(fā)力，整體推進(jìn)。二要打破部門間信息壁壘，發(fā)揮各部門便民信息引流作用，將檢測(cè)宣傳內(nèi)容與便民信息打造升級(jí)為涵蓋全方位監(jiān)管工作的綜合性信息平臺(tái)。

那么可給出待求解的降階系統(tǒng)參數(shù)矩陣：

本研究重點(diǎn)討論在已知系統(tǒng)工作頻率范圍為有限區(qū)間ω∈［ω1,ω2］時(shí)求取最佳的降階模型參數(shù)矩陣，使得降階系統(tǒng)和原給定系統(tǒng)在已知的工作頻率范圍內(nèi)具有最好的逼近性能.為此，這里選取誤差傳遞函數(shù)在已知工作頻率區(qū)間上的最大奇異值

2項(xiàng)研究[18，21]報(bào)道了中藥冷熱交替浸泡治療腦卒中后肩手綜合征的BI評(píng)分結(jié)果。各研究間統(tǒng)計(jì)學(xué)異質(zhì)性較大(P=0.11，I2=61%)，采用隨機(jī)效應(yīng)模型進(jìn)行分析。Meta分析結(jié)果顯示：試驗(yàn)組BI評(píng)分高于對(duì)照組，差異有統(tǒng)計(jì)學(xué)意義[SMD=2.68，95%CI(1.94,3.43)，P<0.000 01]，見圖4。

證明 由給定時(shí)滯系統(tǒng)模型及待求的降階時(shí)滯系統(tǒng)模型，可得出二者之間的誤差系統(tǒng)為

式中，ye（t）=y（t）-yr（t）為模型逼近誤差，xe（t）=［x（t）Txr（t）T］T為誤差系統(tǒng)的增廣狀態(tài)向量.誤差系統(tǒng)的參數(shù)矩陣分別為

肌膚補(bǔ)水除了使用補(bǔ)水保濕的孕婦護(hù)膚品外，最健康的方式就是通過(guò)飲食補(bǔ)水來(lái)潤(rùn)澤肌膚，再輔助孕期專用的護(hù)膚品，加上規(guī)律的飲水，就可以幫助孕媽們甩掉缺水肌膚，重拾水潤(rùn)嫩白肌，也讓冬日不再那么干燥。

結(jié)合誤差系統(tǒng)的參數(shù)矩陣形式，可將定理1中的不等式（7）整理為如下緊致形式：

其中

由式（12），根據(jù)引理2（即投影引理）可知

其中M1i,M2i分別表示相應(yīng)的決策變量集合?1,?2中不同的決策矩陣變量.由定理1和2可知

式中，“?”表示相應(yīng)的矩陣對(duì)稱項(xiàng)，Tr=［Ir×r0r×（n-r）］T,ωc=0.5（ω1+ω2），以及

注1 結(jié)合定理1和2中給出的設(shè)計(jì)條件，可以令有限頻域內(nèi)逼近誤差上界γ為優(yōu)化目標(biāo)，其余涉及的矩陣變量為決策變量，進(jìn)而將本有限頻模型降階問題轉(zhuǎn)化為下列線性矩陣不等式優(yōu)化問題：

由上式以及奇異值的定義可得

即原模型和降階模型之間滿足式（9）所描述的有限頻模型逼近誤差上界，證畢.

3、費(fèi)用與人員費(fèi)用。長(zhǎng)期以來(lái)，科研經(jīng)費(fèi)管理執(zhí)行以收付實(shí)現(xiàn)制為基礎(chǔ)的行政事業(yè)單位會(huì)計(jì)，沒有費(fèi)用的概念。黨的十八屆三中全會(huì)指出，“財(cái)政是國(guó)家治理的基礎(chǔ)和重要支柱”，要“建立權(quán)責(zé)發(fā)生制的政府綜合財(cái)務(wù)報(bào)告制度”。

在許多實(shí)際應(yīng)用中，通常給定的系統(tǒng)都是滿足穩(wěn)定性要求的，對(duì)模型降階問題而言也往往是要求降階系統(tǒng)保持原系統(tǒng)的穩(wěn)定性特征，因此在定理1之外還需要附加能使降階模型穩(wěn)定的條件.

定理2 考慮式（1）給出的連續(xù)時(shí)間線性時(shí)滯系統(tǒng)以及滿足定理1線性矩陣不等式條件的降階系統(tǒng)（8），如果同時(shí)存在對(duì)稱矩陣使得定理1中的矩陣變量同時(shí)滿足

則按照式（8）給出的降階系統(tǒng)為穩(wěn)定的線性時(shí)滯系統(tǒng).

證明 根據(jù)降階模型參數(shù)矩陣的求解式（8），可將式（18）重新整理為

對(duì)式（19）應(yīng)用引理2（投影引理）可得

由式（20）及文獻(xiàn)［23-24］中給出的關(guān)于線性時(shí)滯系統(tǒng)穩(wěn)定性的判據(jù)，可知滿足不等式（18）的降階系統(tǒng)穩(wěn)定.證畢.

基于此，在將PBL融入文檢課教學(xué)時(shí)，教師應(yīng)當(dāng)注意以下幾點(diǎn)：①在實(shí)施PBL教學(xué)法之前，教師應(yīng)做好充分的引導(dǎo)工作，對(duì)PBL的優(yōu)勢(shì)、特點(diǎn)、流程等進(jìn)行充分地講解，讓學(xué)生積極接受并快速適應(yīng)新的教學(xué)模式。②在教學(xué)過(guò)程中，教師要注意學(xué)生的個(gè)體差異，對(duì)基礎(chǔ)薄弱或自學(xué)能力差的學(xué)生加強(qiáng)引導(dǎo)。③項(xiàng)目數(shù)量不宜過(guò)多，以免學(xué)習(xí)和反思不夠深入。④教師要積極跟進(jìn)本學(xué)科最新的研究動(dòng)態(tài)，不斷拓展自身的知識(shí)結(jié)構(gòu)，提升協(xié)調(diào)能力和管理能力，以便更好地勝任PBL教學(xué)。

注意，一方面Q＞0,Z＞0，另一方面（jω-jω1）（jω-jω2）＞0對(duì)任意ω∈［ω1,ω2］均成立，此外根據(jù)文獻(xiàn)［18］中的分析結(jié)果可知因此由式（15）可得

該優(yōu)化問題很容易通過(guò)Matlab的LMIs工具箱得到求解.

注2 定理1給出了當(dāng)系統(tǒng)工作頻率范圍在一個(gè)已知的有限頻率區(qū)間下的降階模型參數(shù)設(shè)計(jì)條件，在實(shí)際工作中還有可能出現(xiàn)系統(tǒng)工作頻率范圍為一個(gè)高頻區(qū)間的情形（即|ω|> ωh）.對(duì)這類問題可以采用和定理1類似的推導(dǎo)方法，結(jié)合文獻(xiàn)［22-23］中給出的高頻區(qū)間線性時(shí)滯系統(tǒng)分析結(jié)果來(lái)給出相應(yīng)的設(shè)計(jì)條件.

2.2 2型糖尿病 我國(guó)2型糖尿病的發(fā)病率逐年增高，這與經(jīng)濟(jì)增長(zhǎng)、人民飲食習(xí)慣改變密切相關(guān)。基因組學(xué)研究表明，長(zhǎng)期高脂飲食會(huì)增加厚壁菌門與擬桿菌的比值，并增加變形菌的數(shù)量，導(dǎo)致胰島素抵抗和肥胖，從而增加糖尿病的發(fā)病風(fēng)險(xiǎn)[58]。腸道菌群的改變?nèi)缂俳z酵母屬、鏈球菌屬、埃希菌屬及腸球菌屬的減少會(huì)降低5-羥色胺的合成，埃希桿菌屬、芽孢桿菌屬以及酵母菌屬的減少會(huì)影響多巴胺和(或)去甲腎上腺素的產(chǎn)生，乳桿菌和雙歧桿菌減少會(huì)影響丁酸的產(chǎn)生[59]。5-羥色胺、多巴胺等可以改善肥胖，而丁酸可以改善胰島素抵抗[60]并激活腸道腸異生[61]，改善血糖[62]，從而治療2型糖尿病。

注3 對(duì)一些具有某些結(jié)構(gòu)特性的線性時(shí)滯系統(tǒng)（如輸出時(shí)滯矩陣Cd=0），往往希望降階模型也具備類似的結(jié)構(gòu)特征，因此在進(jìn)行優(yōu)化計(jì)算時(shí)可以對(duì)相應(yīng)的矩陣變量作相應(yīng)的結(jié)構(gòu)限制來(lái)解決（如要保持降階模型的輸出時(shí)滯矩陣為零的狀態(tài)，可以令決策變量

面向未來(lái)，海爾廚電繼續(xù)完善“智慧大成套”戰(zhàn)略的步伐不會(huì)停歇，致力于通過(guò)科技、模式、服務(wù)等多維度的創(chuàng)新，為消費(fèi)者提供更加舒適、便捷的智慧廚房應(yīng)用體驗(yàn)，也為廚電行業(yè)的轉(zhuǎn)型升級(jí)提供參考方向。據(jù)悉，2019年1季度，海爾廚電將再推出5套成套廚電，全方位滿足用戶的專屬個(gè)性化定制需求。

6.dabusu qami arun-a bui ɡebel ，dalai tenɡkis ba na arun-a要問鹽是哪里來(lái)的，是從大海和湖泊里來(lái)的)

3 算例

考慮在文獻(xiàn)［14］中給出的一個(gè)連續(xù)時(shí)間線性時(shí)滯系統(tǒng)（1），其參數(shù)矩陣為

式中，延時(shí)時(shí)長(zhǎng)d=1.6.在本例中假設(shè)系統(tǒng)的工作頻率范圍ω∈［-0.5,0.5］.注意到系統(tǒng)參數(shù)矩陣Cd=0，為保持這樣一個(gè)限制結(jié)構(gòu)，這里采用由定理1和2共同給出的線性矩陣不等式設(shè)計(jì)條件來(lái)進(jìn)行降階模型參數(shù)矩陣的優(yōu)化計(jì)算，同時(shí)令其中的矩陣變量其Matlab的計(jì)算結(jié)果為

原系統(tǒng)、降階系統(tǒng)以及誤差系統(tǒng)在給定工作頻率范圍內(nèi)的傳遞函數(shù)最大奇異值曲線如圖1所示.相應(yīng)的工作頻率范圍內(nèi)的誤差滿足σmax（G（jω）-Gr（jω））＜0.245 6，?ω∈［-0.5,0.5］.注意到模型降階的本質(zhì)需求為通過(guò)降階模型來(lái)模擬原系統(tǒng)在給定輸入信號(hào)激勵(lì)下的輸出響應(yīng)，為進(jìn)一步說(shuō)明本研究在針對(duì)輸入信號(hào)為有限頻信號(hào)時(shí)的有效性和優(yōu)越性，下面給出在兩類輸入信號(hào)頻率的帶寬均不超過(guò)給定頻率范圍情況下的時(shí)域仿真結(jié)果.首先考慮當(dāng)輸入為常見的階躍信號(hào)（頻率為零）的情況，采用下式給出的階躍信號(hào)

圖1 原系統(tǒng)、降階系統(tǒng)以及誤差系統(tǒng)在工作頻率范圍內(nèi)的傳遞函數(shù)的最大奇異值Fig.1 Maximum singular value of transfer functions of original system,reduced-order system and error system

在Matlab中進(jìn)行仿真可以得到如圖2所示的仿真結(jié)果.從圖2可以看出，在相同的階躍信號(hào)激勵(lì)下，原系統(tǒng)的輸出信號(hào)對(duì)降階系統(tǒng)的輸出信號(hào)具有良好的逼近效果.

4)自愈合效果最佳組合是微膠囊摻量 0.8%、微膠囊粒徑200 μm、養(yǎng)護(hù)齡期7 d、養(yǎng)護(hù)溫度40 ℃。

圖2 原系統(tǒng)和降階系統(tǒng)在階躍輸入信號(hào)u1（t）下的輸出響應(yīng)Fig.2 Output response of the original system and the reduced-order system with step input signals u1（t）

接下來(lái)考慮一個(gè)能量有界的三角形輸入信號(hào)：

容易驗(yàn)證，該信號(hào)的主要頻率分量屬于給定的頻率區(qū)間.在該信號(hào)激勵(lì)下，原系統(tǒng)和降階系統(tǒng)的輸出響應(yīng)如圖3所示.

圖3中給出的原系統(tǒng)輸出信號(hào)和降階系統(tǒng)輸出信號(hào)具有非常相似的動(dòng)態(tài)特性，從時(shí)域上說(shuō)明降階系統(tǒng)和原系統(tǒng)在有限頻輸入信號(hào)下的良好逼近性能.圖2和3的結(jié)果表明本方法不失為一類有效的線性時(shí)滯系統(tǒng)有限頻模型降階問題的解決方案.

圖3 原系統(tǒng)和降階系統(tǒng)在三角形輸入信號(hào)u2（t）下的輸出信號(hào)響應(yīng)Fig.3 Output signals of original system and reduced system with sinusoidal input signal u2（t）

4 結(jié)束語(yǔ)

本研究針對(duì)時(shí)間連續(xù)線性時(shí)滯系統(tǒng)，分析了在系統(tǒng)工作頻率范圍為已知有限區(qū)間情形下的模型降階問題.通過(guò)引入合適的用于刻畫有限頻域內(nèi)模型逼近性能的指標(biāo)函數(shù)，以及結(jié)合矩陣不等式技術(shù)和線性時(shí)滯系統(tǒng)的性能分析方法，推導(dǎo)出了用于求解降階系統(tǒng)模型參數(shù)矩陣的逼近性能優(yōu)化設(shè)計(jì)條件，并給出了保持系統(tǒng)穩(wěn)定的設(shè)計(jì)條件.數(shù)值算例驗(yàn)證了所提出方法的有效性，特別是在系統(tǒng)工作頻率范圍較小的情況下，所得降階模型通常可以很好地逼近原模型在工作頻率范圍內(nèi)的時(shí)頻域動(dòng)態(tài)特性.

［1］ROZZA G.Reduced order methods for modeling and computational reduction［M］.Berlin：Springer,2014：25-60.

［2］PETER B,MICHAEL H,MATEN T,et al.Model reduction for circuit simulation［M］.Berlin：Springer,2011：3-25.

［3］ANTOULAS A C,SORENSEN D C.Approximation of large-scale dynamical systems：an overview［J］.Applied Mathematics and Computer Science,2011,11（5）：1093-1122.

［4］EBIHARA Y,TOMOMICHI H.On H∞model reduction using LMIs［J］.IEEE Transactions on Automatic Control,2004,49（7）：1187-1191.

［5］ZHONG Q C.Robust control of time-delay systems［M］.Berlin：Springer,2006：68-90.

［6］張仲華,孟慶勛,鎖要紅.一類具有病菌信息交流機(jī)制的時(shí)滯模型的穩(wěn)定性［J］.上海大學(xué)學(xué)報(bào)（自然科學(xué)版）,2013,19（3）：308-314.

［7］BEATTIE C,GUGERCIN S.Interpolatory projection methods for structure-preserving model reduction［J］.Systems&Control Letters,2009,58（3）：225-232.

［8］ELIAS J,TOBIAS D,WIM M.Model reduction of time-delay systems using position balancing and delay Lyapunov equations［J］.Mathematics of Control,Signals,and Systems,2013,25（2）：147-166.

［9］WANG X,WANG Q,ZHANG Z,et al.Balanced truncation for time-delay systems via approximate gramians［C］//Proceedings of the 16th Asia and South Pacific Design Automation Conference. 2011：55-60.

［10］WANG X,ZHANG Z,WANG Q,et al.Gramian-based model order reduction of parameterized time-delay systems［J］.International Journal of Circuit Theory and Applications,2014,42（7）：687-706.

［11］ZHOU P,XIANG B,FU J,et al.Model approximation of multiple delay transfer function models using multiple-point step response fitting［J］.International Journal of Control,Automation and Systems,2012,10（4）：180-185.

［12］XU S,LAM J,HUANG S,et al.H∞model reduction for linear time-delay systems：continuoustime case［J］.International Journal of Control,2001,74（11）：1062-1074.

［13］XU S,LAM J,ZOU Y.New results on delay-dependent robust H∞control for systems with time-varying delays［J］.Automatica,2006,42（2）：343-348.

［14］GAO H,LAM J,WANG C,et al.Hankel norm approximation of linear systems with time-varying delay continuous and discrete cases［J］.International Journal of Control,2004,77（17）：1503-1520.

［15］WU L,ZHENG W X.Weighted H∞model reduction for linear switched systems with time-varying delay［J］.Automatica,2009,45（1）：186-193.

［16］LU H,ZHOU W,XU Y,et al.Time-delay dependent H∞model simplification for singular systems with Markovian jumping parameters［J］.Optimal Control Applications and Methods,2011, 32（4）：379-395.

［17］鐘金標(biāo),杜鑫,朱訓(xùn)林.一類離散奇異Markov跳變系統(tǒng)的鎮(zhèn)定性［J］.上海大學(xué)學(xué)報(bào)（自然科學(xué)版）, 2014,20（4）：513-520.

［18］DU X,YANG G H.H∞model reduction of linear continuous-time systems over finite-frequency interval［J］.IET Control Theory&Applications,2010,4（3）：499-508.

［19］DING D W,DU X,LI X.Finite-frequency model reduction of two-dimensional digital filters［J］. IEEE Transactions on Automatic Control,2014,6（6）：1624-1629.

［20］SHI X,DING D W,LI X,et al.Model reduction of discrete-time switched linear systems over finite-frequency ranges［J］.Nonlinear Dynamics,2013,71（1/2）：361-370.

［21］LI X,YU C,GAO H.Frequency-limited H-infinity model reduction for positive systems［J］.IEEE Transactions on Automatic Control,2014,60（4）：1093-1098.

［22］ZHANG X N,YANG G H.Performance analysis for multi-delay systems in finite frequency domains［J］.International Journal of Robust and Nonlinear Control,2012,22（8）：933-944.

［23］ZHANG X N,YANG G H.Delay-dependent state feedback control with small gain conditions in finite frequency domains［J］.International Journal of Systems Science,2011,42（3）：369-375.

［24］沈喆,邢俊雷,楊光紅,等.時(shí)滯依賴觀測(cè)器的線性系統(tǒng)故障檢測(cè)［J］.控制理論與應(yīng)用, 2013,30（5）：592-596.

本文彩色版可登陸本刊網(wǎng)站查詢：http：//www.journal.shu.edu.cn

LMI based approach for finite-frequency model order reduction of continuous-time linear time-delayed systems

DU Xin1，F(xiàn)AN Peibing1，LIU Fuwei2
（1.School of Mechatronics Engineering and Automation，Shanghai University，Shanghai 200072，China；
2.Huawei Technology Co.，Ltd.，Shanghai 200040，China）

Model order reduction of continuous-time linear time-delayed systems over limited frequency intervals is discussed in this paper.The approximation performance is characterized by introducing an index associated with the finite-frequency maximum singular value of the error transfer function.With the aid of some fundamental matrix inequality techniques,sufficient criterion for stability of the reduced-order model and optimizing finite-frequency approximation error is derived.The model order reduction problems can be tackled by solving the corresponding linear matrix inequalities（LMIs）based optimization problems.A numerical example is given to show effectiveness of the proposed technique.

model order reduction；finite-frequency；linear time-delayed system；linear matrix inequality（LMI）

TP 273

A

1007-2861（2016）04-0408-13

10.3969/j.issn.1007-2861.2015.04.014

2014-12-02

國(guó)家自然科學(xué)基金資助項(xiàng)目（61304143，61174085）

杜鑫（1983—），男，博士，研究方向?yàn)槟Ｐ徒惦A、控制系統(tǒng)優(yōu)化設(shè)計(jì)等. E-mail：duxin@shu.edu.cn