馮德龍,楊鎖昌，劉映希

(1.軍械工程學院 導彈工程系，河北 石家莊 050003；2.西安電子科技大學 雷達信號處理實驗室，陜西 西安 710071)

0 引言

導彈制導系統(tǒng)是一個具有隨機輸入的非線性時變系統(tǒng)，在此系統(tǒng)中含有非線性因素和一些隨機效應，這些非線性因素和隨機效應影響導彈制導精度。對具有隨機輸入的非線性時變系統(tǒng)進行統(tǒng)計分析，傳統(tǒng)的方法是采用Monte Carlo方法。但由于此方法存在計算量大、耗時過長的缺點，使得在一些對快速性能要求較高的環(huán)境無法使用該方法。協(xié)方差分析描述函數(shù)技術(covariance analysis describing equation technique，CADET)是在20世紀70年代由美國人A.Celb和R.S.Warren共同提出的對戰(zhàn)術飛行器制導系統(tǒng)性能進行統(tǒng)計分析的一種新方法。它解決了Monte Carlo方法快速性能較差的缺點[1]。協(xié)方差分析是一種用于分析具有隨機輸入的時變線性系統(tǒng)的計算機工具。應用這種方法，通過對非線性矩陣微分方程式進行直接積分，就可以產(chǎn)生作為一個時間函數(shù)的系統(tǒng)狀態(tài)向量協(xié)方差矩陣。應用這種技術，可以預估導彈制導系統(tǒng)的統(tǒng)計特性[2-4]。

1 CADET技術

1.1 協(xié)方差分析描述函數(shù)技術的基本原理

假設具有隨機輸入的連續(xù)非線性系統(tǒng)狀態(tài)方程為[5]

(1)

式中：f(x,t)為n維系統(tǒng)狀態(tài)變量x(t)的非線性向量函數(shù)。非線性時變系統(tǒng)方程如圖1所示。

圖1 非線性時變系統(tǒng)方程示意圖Fig.1 Schematic diagram of nonlinear time-varying system

非線性系統(tǒng)的傳播方程：

GT(t).

(2)

定義輔助矩陣N，使得

NP(t)=Ef(x,t)rT,

(3)

(4)

由于P(t)是正定的，故有

N=Ef(x,t)rTP-1.

(5)

根據(jù)上述統(tǒng)計線性化原理，令非線性函數(shù)向量f(x,t)可近似描述為x(t)的一個線性函數(shù)：

(6)

(7)

非線性時變系統(tǒng)均值方程和協(xié)方差傳播方程為

(8)

初始條件同上，利用方程(8)和初始條件式可直接解算系統(tǒng)變量的誤差統(tǒng)計特性。

1.2 高斯假設

目前，大多數(shù)協(xié)方差描述函數(shù)技術的研究都是采用假定系統(tǒng)狀態(tài)變量為高斯分布的辦法[6-7]。

(1) 高斯假設的重要結論

1)設系統(tǒng)的狀態(tài)變量都是正態(tài)的，則矢量非線性函數(shù)f(x,t)中的每一個標量非線性函數(shù)都可以單獨進行準線性化。

(9)

(10)

(2) 高斯假設的合理性

高斯分布的隨機過程通過非線性系統(tǒng)，其輸出量就變成非高斯分布的了[8]。即使這樣，高斯假設在工程上仍有它的實用價值，原因在于導彈制導系統(tǒng)的線性部分通常都具有較好的低通濾波作用。由中心極限定理可知，非高斯分布的隨機過程通過低通濾波后，將趨向于高斯隨機過程。由于系統(tǒng)線性部分的濾波作用，使得高斯假設基本成立。

(3) 根據(jù)上述理論假設[9]

1)x(t)服從聯(lián)合正態(tài)分布。

2) 非線性的隨機干擾以線性方式加入到系統(tǒng)當中，且服從正態(tài)分布。

2 一種改進比例制導律精度分析

2.1 制導回路模型

圖2 制導回路模型Fig.2 Guidance loop model

導彈自尋的末制導回路模型包括自動駕駛儀模型、導引頭模型、導彈目標相對運動模型、制導律和目標機動模型[10-11]，其相互關系如圖2所示。

2.2 制導律

本文對一種改進的比例導引律進行分析，其形式如下：

(11)

2.3 協(xié)方差分析

后面幾個字起了作用，花五奇這些年作為德公公貼身護衛(wèi)，早已習慣眉高眼低看人，老太醫(yī)往往品級不低，且能與后宮說上話，輕易不能得罪，于是他暫時放下兵器，不耐煩地問道：“老太醫(yī)？什么老太醫(yī)？”

圖3 制導回路Fig.3 Guidance loop

(12)

式中：

(13)

因為上述一階微分方程式是狀態(tài)的函數(shù)，它們可以寫成狀態(tài)空間表達如下：

(14)

定義：

則系統(tǒng)狀態(tài)方程可以表述為

(15)

式中：x(t)為系統(tǒng)的狀態(tài)；F(t)為系統(tǒng)動態(tài)特性矩陣；u(t)為具有譜密度矩陣Q(t)的白噪聲矢量，即

Q(t)=Eu(t)uT(t)，

(16)

(17)

其中協(xié)方差矩陣X(t)與x(t)的關系為

X(t)=Ex(t)xT(t).

(18)

如果隨機干擾具有零均值，則協(xié)方差矩陣的對角元素表示狀態(tài)變量的方差。協(xié)方差矩陣的非對角元素表示狀態(tài)變量之間的相關程度。

3 仿真分析

3.1 程序仿真

整個仿真程序都是采用Matlab編寫，仿真的流程圖如圖4所示。

圖4 CADET仿真流程圖Fig.4 Flow chart of CADET simulation

根據(jù)圖4所示的線性化末制導模型，導彈飛控系統(tǒng)為一階慣性環(huán)節(jié)且時間常數(shù)T為1 s，飛行時間tf為0～10 s，有效導航比為3，比例系數(shù)F為0.5，彈目相對速度為4 000 m/s，目標的加速度為5 m/s2。在制導系統(tǒng)中唯一的誤差源是目標機動特性為3g的階躍機動且開始時刻在整個飛行時間內(nèi)均勻分布。仿真結果如圖5～8所示。

圖5 比例導引律的協(xié)方差分析Fig.5 Covariance analysis of proportional guidance law

圖6 改進比例導引律的協(xié)方差分析Fig.6 Covariance analysis of improved proportional guidance law

圖7 比例導引律的導彈加速度變化Fig.7 Missile acceleration change of proportional guidance law

圖8 改進比例導引律的導彈加速度變化Fig.8 Missile acceleration change of improved proportional guidance law

3.2 結果分析

如圖8所示是整個飛行時間為10 s時，導彈與目標的相對距離的總的標準差。由圖5和圖6可以看出這種改進的比例導引律比之前的導引精度提高了大約3倍，所以這種改進是有效的、可行的。與伴隨技術不同，協(xié)方差分析方法在一次仿真中并不能提供所有不同飛行時間的脫靶量，但協(xié)方差分析方法可以得到在一次飛行時間內(nèi)所有時刻的相對距離y的信息，在飛行結束時，相對距離即為脫靶量，即σmiss=σy(tf)，而且也可以得到圖7和圖8所示的導彈加速度變化。

4 結束語

協(xié)方差分析描述函數(shù)技術是非線性時變系統(tǒng)統(tǒng)計分析的一種近似解析方法，它可以通過一次解析計算得到系統(tǒng)所有狀態(tài)變量的統(tǒng)計性能。該方法較傳統(tǒng)的蒙特卡羅試驗法要節(jié)省很多計算時間。本文對于一種改進的比例導引律，應用協(xié)方差分析法驗證了其可行性。協(xié)方差分析方法雖然也適用于分析隨機干擾的影響，但是只能獲得飛行末端時刻的脫靶量，對各干擾的影響需進行獨立仿真。

參考文獻：

[1] 畢開波，王曉東，劉智平.飛行器制導與控制及其MATLAB仿真技術[M].北京：國防工業(yè)出版社，2009：163-169.

BI Kai-bo, WANG Xiao-dong, LIU Zhi-ping. Vehicle Guidance and Control and Matlab Simulation [M]. Beijing: National Defense Industry Press, 2009：163-169.

[2] 郭慶，楊明，王子才.攔截器末制導段誤差模型的建立與仿真[J].計算機仿真，2007,24(3)：58-61.

GUO Qing, YANG Ming, WANG Zi-cai. Establishment and Simulation of Error Model for Interceptor’s Terminal Guidace Phase [J]. Computer Simulation, 2007, 24(3): 58-61.

[3] 蔣瑞民，周軍，郭建國. 導彈制導系統(tǒng)精度分析方法研究[J].計算機仿真，2011，28(5)：76-79.

JIANG Rui-min, ZHOU Jun, GUO Jian-guo. Research on Methods of Missile Guidace Precision Analysis [J]. Computer Simulation, 2011, 28(5): 76-79.

[4] 柴華，杏建軍，唐國軍.基于協(xié)方差分析描述函數(shù)技術的彈道導彈落點預報[J].導彈與航天運載技術，2011,4(4)：37-40.

CHAI Hua, XING Jian-jun, TANG Guo-jun. Ballistic Missile Impact Poit Prediction Based on CADET [J]. Missile and Space Vehicles, 2011, 4(4): 37-40.

[5] 周衛(wèi)文，梁曉庚，賈曉洪.自尋的導彈末制導系統(tǒng)仿真分析方法研究[J].計算機仿真，2012，29(8)：45-49.

ZHOU Wei-wen, LIANG Xiao-geng, JIA Xiao-hong. Study on Missile Homing Loop Performance Analysis Methods [J]. Computer Simulation, 2012, 29(8):45-49.

[6] YANUSHEVSKY R．Analysis and Design of Missile Guidance Systems in Frequency Domain [C]∥ 44th AIAA Space Sciences Meeting， Reno，Nevada，2006：821-825．

[7] SHNEYDOR N A．Missile Guidance and Pursuit[M]．Chichester:Horwood Publishing，1998．

[8] WANG Z M.A New Method for Estimating System Error of Guidance Instrument [J]. Science of China：Series E, 1998, 28(2): 160-167.

[9] 李海平，鐘瑞麟，魏岳，等.導彈末制導精度分析的協(xié)方差方法研究[J].戰(zhàn)術導彈技術，2004，6(1)：49-54.

LI Hai-ping, ZHONG Rui-lin, WEI Yue, et al. Research on Covariance Method of Accuracy Analysis in Missile Terminal Guidance [J]. Tactical Missile Technology, 2004, 6(1): 49-54.

[10] XU Jin-xiang. Key Techniques of Terminal Correction Mortar Projectiles [J]. Journal of Beijing Institute of Technology, 2007, 16(2)：168-172.

[11] FRORIEP R, JOOST D. Stocchastic Simulation Using Covarince Techniques: Modular Program Package for Nonlinear Missile Guidance [J]. Journal of Guidance and Control, 1984, 7(4): 509-512.

[12] Iain Matthews, Takahiro Ishikawa, Simon Baker. The Template Update Problem [C]∥IEEE Transactions on Pattern Analysis and Machine Intelligence, 2004, 6(6): 810-815.