導彈末端彈道擬合方法研究＊

袁震宇 鐘志通 張 林 張 宇

（海軍大連艦艇學院 大連 116018）

導彈末端彈道擬合方法研究＊

袁震宇 鐘志通 張 林 張 宇

（海軍大連艦艇學院 大連 116018）

導彈末端彈道估計問題可轉換為空間直線擬合問題，通過比較解決空間直線擬合問題的一般最小二乘算法、整體最小二乘算法和主成分分析法，得出使用主成分分析法進行空間直線擬合，理論嚴密，使用冪法對主成分分析法進行求解，計算簡單。算例分析結果表明了主成分分析法求解空間直線擬合問題的可行性和有效性。

末端彈道擬合；主成分分析法；冪法

Class NumberE927

1 引言

導彈末端彈道擬合是導彈作戰(zhàn)效能精確評估的一個重要環(huán)節(jié)。由于導彈末端彈道距離短，導彈速度高，可能的機動小，因此可以看成是一條直線。從而導彈末端彈道擬合可視為將導彈末端彈道測量點擬合成一條直線的過程，即空間直線擬合。

2 空間直線擬合問題

空間直線擬合問題［1～4］可以描述為：在三維空間給定n個點P1（x1，y1，z1），P2（x2，y2，z2），…，Pn（xn，yn，zn），設每個點在x，y，z三個方向上的測量誤差服從相同且相互獨立的正態(tài)分布，求空間直線L的方程：

其中P0（x0，y0，z0）為直線上已知一點，→v＝（X，Y，Z）為直線的方向向量，使得直線滿足構造的某種意義下的準則。

空間直線擬合相關研究較多，比較典型的方法有最小二乘算法、主成分分析算法等［2～3，5～7］。

3 最小二乘算法

對于直線方程（1），不妨設Y＝1（若Y＝0，則必然有X≠0或者Z≠0），則直線方程經過整理可以表示為

其中未知參數為x0，z0，X，Z，表示為方程組的形式如下：

寫成矩陣形式為

當有n個測量點時，將第i個點的坐標代入直線方程得：

聯(lián)立n個方程有

余下工作即是通過已知的觀測量xi、yi和zi，估計未知參數x0，z0，X，Z，進而求得直線方程。

1）一般最小二乘算法

類似平面直線擬合問題，將式（4）引入觀測誤差，表示為

根據最小二乘原理，在如下準則下：

解得：

化簡得：

通過簡單一個算例驗證：設四個測量點為（－10，－6，2）、（1，1，4.2）、（9，4.6，5.9）和（21，9.5，8.1），對以上四個空間測量點擬合空間直線，使用前述一般最小二乘算法，根據式（9）解得：

即直線方程為

將空間的測量點及擬合直線如圖1所示。

圖1 測量點及擬合直線

圖中“正方形”點為空間中待擬合的點，“圈形”點為空間點在各個坐標平面的投影，“正方形”點旁邊的直線為空間擬合直線，而各“圈形”點旁邊的直線則為空間擬合直線在各坐標平面上的投影，從圖1可以看出，在XZ平面上，直線有較好的擬合效果，而在XY平面和YZ平面上，擬合直線與對應的點存在明顯偏差。產生以上偏差的主要原因在于一般最小二乘算法使用的準則式（8）沒有考慮觀測點的縱坐標的誤差。

2）整體最小二乘算法

整體最小二乘算法［1～3，8］相比一般最小二乘算法，綜合考慮了坐標點的誤差，構造新的準則函數，如構造準則函數：

其中＾yi為觀測點縱坐標的估計值，yi－＾yi為觀測點縱坐標存在的誤差。

整體最小二乘算法相比經典最小二乘算法，理論上更加嚴密，參數估計結果為最優(yōu)解，但是求解過程更加復雜，不能直接獲得顯式的結果，必須通過迭代或者SVD算法等進行求解，求解過程較為繁瑣，使用性不佳［2，9～11］。

4 主成分分析算法

將三維空間中的點擬合成一條直線，其本質是降維，即將三維分布的點降為一維，即直線。依據主成分的思想，對三維空間的點集數據進行主成分分析，則其第一主成分可作為擬合直線的方向向量［1］。文獻［7］證明了在最小二乘準則下，空間直線必通過點集的數據中心，因此可以由點向式確定空間直線方程。據此給出直線擬合算法步驟如下：

1）將三維空間的n個點構造矩陣Xn×3，求出均值ˉX＝（ˉx，ˉy，ˉz），則擬合直線必然經過點集的數據中心（ˉx，ˉy，ˉz）。

2）求協(xié)方差矩陣為

3）使用冪法求協(xié)方差矩陣S最大特征值λ1對應的特征向量；

（1）任取非零向量s0，歸一化得＾s0；

（2）求si＝S·si－1，歸一化得＾si，i＝1，2…；

（3）判斷｜＾si－＾si－1｜是否小于給定的誤差限ε，“是”則＾si為所求最大特征值對應的特征向量→v＝（X，Y，Z），“否”則重復（2）、（3）兩個步驟。

由主成分的意義，各測量點到直線距離平方和為n（λ2＋λ3），其中λ2和λ3為S的其余兩個特征值，與式（10）物理意義相同［1］。

4）由點向式確定直線方程為

5 算例分析

表1為空間直線的一組實測數據，已知其中含有的模型誤差為δ＝±0.005。

表1 空間直線實測樣本數據

按照主成分分析算法步驟，計算擬合直線必然經過的數據中心點為

（ˉx，ˉy，ˉz）＝（7.5014，11.9991，16.5011）

通過冪法求得擬合直線的方向向量為

因此擬合直線的方程可以表示為

確定直線后可以容易地得到空間各點到擬合直線的距離，即擬合誤差Δi，具體值列于表1中。擬合總誤差∑Δi＝0.0307，誤差最大值0.0086，誤差最小值0.0009。與文獻［2，6］結果非常接近，如表2所示，說明本文方法的有效性。

表2 各種方法計算結果比較

6 結語

空間直線擬合不同于平面直線擬合，可以將直線擬合問題表示為簡單的線性形式，如果類似平面直線擬合進行簡化，則會因為模型簡化產生的誤差而使得擬合的直線存在較大的誤差；對于空間直線擬合使用整體最小二乘算法可以解決使用一般最小二乘方法產生模型誤差的問題，但是算法的求解過程較為繁瑣；本文采用主成分分析法解決空間直線擬合問題，理論嚴密，等價于整體最小二乘算法，求解過程中采用冪法求解協(xié)方差矩陣的最大特征值對應的特征向量，相比整體最小二乘算法計算簡單，算例分析表明，使用主成分分析法運算簡單，運算結果與其他文獻基本相同，證明了本文方法的可行性和有效性，因此在實際工程應用中可使用主成分分析法進行導彈末端彈道的估計。

［1］曹慧榮，王福昌.全最小二乘和全加權最小一乘準則下的空間直線擬合［J］.計算機應用與軟件，2010，27（9）：111－113.

［2］姚宜斌，黃書華，孔建，等.空間直線擬合的整體最小二乘算法［J］.武漢大學學報·信息科學版，2014，39（5）：571－574.

［3］王繼剛，周立，蔣廷臣.一種簡單的加權整體最小二乘直線擬合方法［J］.測繪通報，2014（4）：33－35.

［4］王偉峰，溫耐.空間直線擬合研究［J］.許昌學院學報，2010，29（5）：37－39.

［5］潘國榮，唐杭.特征分解與選權迭代在空間直線擬合中的應用［J］.東南大學學報（自然科學版），2013，43（SupⅡ）：250－255.

［6］郭際明，向巍，尹洪斌.空間直線擬合的無迭代算法［J］.測繪通報，2011（2）：24－26.

［7］Wuh Y，Chen Q，Yachidam.Face detection from color images using a fuzzy pattern matching method［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence，1999，21（6）：557－563.

［8］王炳杰，趙軍鵬，王春杰.基于三維最小二乘方法的空間直線度誤差評定［J］.北京航空航天大學學報，2014，40（10）：1477－1480.

［9］王媛，霍李，宮新宇，等.基于多點空間直線擬合的調炮精度檢測方法［J］.火力與指揮控制，2015，40（3）：130－133.

［10］王玥，趙振學，王江榮.基于空間直線擬合的古塔變形趨勢研究［J］.蘭州文理學院學報（自然科學版），2014，28（3）：32－36.

［11］陳天立，趙建勛，高揚.基于三維空間直線擬合的衛(wèi)星天線動態(tài)精密校準方法研究［J］.電子測量技術，2011，34（3）：99－102.

Missile Terminal Ballistic Curve FittingMethods

YUAN Zhenyu ZHONG Zhitong ZHANG Lin ZHANG Yu
（Dalian Naval Academy，Dalian 116018）

The missile terminal ballistic curve estimation can be transferred to be a linear fitting problem.Based on the comparison of general least square algorithm，total least square algorithm and principle component analysis，it is concluded that linear fitting using principle analysis is rigid reasoning，and solving the principle analysis by power method computation cost less.At the end，an example solving result indicates the feasibility and validity of the principle analysis method which is used to solve the 3－dimention linear fitting.

terminal ballistic curve fitting，principle component analysis，power method

E927DOI：10.3969／j.issn.1672－9730.2015.11.015

2015年5月7日，

2015年6月23日

學院科研發(fā)展基金資助。

袁震宇，男，碩士，講師，研究方向：艦載導彈戰(zhàn)斗使用、數據分析等。鐘志通，男，博士，副教授，研究方向：艦載導彈戰(zhàn)斗使用。張林，男，博士，講師，研究方向：艦載導彈戰(zhàn)斗使用。張宇，男，博士，講師，研究方向：數字圖像處理。