基于測控電磁波極化特征信息的飛行器姿態估計

孫紅梅*① 陳廣東② 張 弓①

?

基于測控電磁波極化特征信息的飛行器姿態估計

孫紅梅陳廣東張 弓

(南京航空航天大學電子信息工程學院 南京 210016)(南京航空航天大學無人機研究院 南京 210016)

該文應用正交解耦的姿態參數信號模型描述極化電磁波的3維波結構，將最大似然估計和MUSIC算法推廣應用于部分極化波參數估計，利用極化電磁波傳遞信號基站標定的姿態基準信息。飛行器通過單電磁矢量傳感器接收并處理信號，確定電磁波結構向量，計算得到地理坐標系下的飛行器姿態。此方法有別于需多點測量的三角計算法，只需一個地面基站信號和運動平臺上單一接收點，就可實現姿態感知，可替代姿態導航儀和航向導航儀。該文介紹了運動平臺姿態估計算法，為工程應用提供技術支持。

電磁矢量傳感器；極化狀態；波達方向；姿態導航；部分極化信號

1 引言

電磁波是雷達和通信等信息系統中承載信息的媒介，信號時延、多普勒頻率、空間到達角和極化狀態是其重要特征參量，能攜帶豐富的信息。在時域，電磁波的往返時間用來探測目標的距離；在頻域，多普勒頻移用來探測目標速度；在空域，波達角用來探測目標方位。因此對電磁信號的獲取和處理顯得尤為重要。目前提出的利用旋轉不變ESPRIT類算法研究波達角和極化角的聯合參量估計問題，所用陣列為均勻、矩形、稀疏電磁矢量傳感器陣列，它要求陣元間距離遠大于信號半波長。而MUSIC算法通過構造空域-極化域聯合參數估計譜，搜索MUSIC譜峰值位置來估計信號到達角和極化角，其對聯合譜的譜估計精度和分辨力這兩項性能指標要求很高。此外，國內外都研究在地球衛星導航系統接收端利用多點接收、三角計算的方法來感知姿態。相比之下，利用單一電磁矢量傳感器估計信號波達方向和極化參數的優勢在于：(1)適用于近場或遠場，寬帶或窄帶信號；(2)在空間域-極化域可分離多個極化信號；(3)陣列體積小，分辨率相對較高；(4)天線的工作方向覆蓋全空域，各方向信號接收性能一致，適宜機動平臺使用。

本文應用正交解耦的姿態參數信號模型描述極化電磁波的3維波結構，優化傳統的最大似然估計部分極化波參數估計方法，利用極化電磁波傳遞信號基站標定的姿態基準信息，避免了前面方法提出的一些標準要求，只需一個基站信號和飛行器上單一接收點，就可實現姿態感知，并用仿真實驗驗證了方法的有效性。

2 電磁矢量傳感器接收信號

(1)

主次極化分量導向矢量：

其中

空域導向因子：

極化域導向因子：

根據3維空間旋轉理論知，繞坐標軸的3次連續旋轉可以用1次旋轉替代，3個連續旋轉矩陣相乘可合并為1個旋轉矩陣：

(3)

可得：

3 信號參數最大似然估計

若要在噪聲環境下獲得高精度參數估計，可應用最大似然估計。假設式(1)中的在飛行器測控系統中為固定量，定義?？紤]到飛機與地面站距離不斷變化，飛機所處空間環境不斷變化，噪聲方差和接收信號功率也不斷變化。該信號模型快拍取樣的似然函數為：

似然函數對數為：

代入式(5)得：

圖2 為姿態導向矢量信號最大似然估計仿真結果，圖中信噪比0 dB，采樣快拍1000次，信號的極化橢圓率，極化度60%。圖2(a)和圖2(c)中接收到的信號1姿態參數分別為31°, 64°和17°，圖2(b)和2(d)中接收到的信號2姿態參數分別為2°,和13°，圖2(a)和圖2(b)為和時隨變化的3維圖形，由圖可見，在與信號源參數匹配位置出現清晰峰形。圖2(c)和圖2(d)分別為圖2(a)和圖2(b)的等高線圖，可方便觀察峰的位置。

4 MUSIC算法估計多信號參數

MUSIC算法正是利用信號子空間與噪聲子空間的正交性來估計信號DOA。利用噪聲子空間的所有向量構造MUSIC譜，空間方位譜中的峰值位置則對應信號的來波方向。MUSIC算法對天線陣列的形狀沒有要求，很大程度上提高了測向分辨率。估計多信號源參數時，基于特征值分解的MUSIC算法計算效率優于最大似然估計。將部分極化波唯一分解為兩個極化狀態正交的完全極化波的疊加，,分別為信號的主次極化分量導向矢量。具體步驟如下：

式中前后項分別為部分極化波主次極化分量零譜的倒數。通過式(7)求得極大值的位置來確定信號參數。

圖2 姿態導向矢量信號的最大似然估計結果(信號極化度60%，信噪比0 dB)

圖4為MUSIC算法估計多信號參數仿真結果，圖4中的兩個信號源極化度為60%，其余參數和采樣條件與圖2和圖3中分別出現的兩個信號相同，圖4采用姿態導向矢量信號計算MUSIC譜。由圖4(a)和圖4(b)可見，在信號參數匹配位置，出現了清晰的峰，表明采用姿態導向矢量信號計算MUSIC譜，能正確估計多信號參數。

圖5出現的兩個信號源與圖4相同，圖5采用空域-極化域導向矢量信號計算MUSIC譜。圖5(a)極化橢圓傾角與信號1匹配，由圖5(a)可見，在信號1參數匹配位置，出現了清晰的峰，同時存在信號2的干擾峰。圖5(b)極化橢圓傾角與信號2匹配，由圖5(b)可見，信號2的主峰旁出現寄生峰。

計算MUSIC譜估計多信號參數仿真結果表明，采用姿態導向矢量信號模型更有利于參數估計。

5 飛行器姿態估計

配置于地面的基站作為極化波的輻射源向飛機發射橢圓極化波，飛行器利用電磁矢量傳感器接收基站信號，估計出，得機身坐標系到波結構坐標系轉換矩陣；設基站測得的對飛機的仰角為，方位角為偏北(逆時針向)，通信傳遞給飛機，在地理坐標系下發射信號的極化參數為固定值，可直接存儲在飛機中，滿足遠場條件時，和也不變，也可直接存儲在飛機中。由得波結構坐標系到地理坐標系轉換矩陣為：

機身坐標系到地理坐標系轉換矩陣為：

(8)

圖4 姿態導向矢量信號的MUSIC算法(信號極化度60%，信噪比0 dB)

圖5 空域-極化域導向矢量信號的MUSIC算法(信號極化度60%，信噪比0 dB)

因機身坐標系中，飛機姿態矩陣為單位陣，所以式(8)就是地理坐標系下飛機姿態矩陣。該姿態矩陣已完整地表達了飛行器空間姿態，若將該姿態矩陣化為：

單電磁矢量傳感器可全向工作，即在飛行器的任何機動姿態下，電磁波姿態測量都將有效。與機載天線相比，地面基站的大型天線系統測得的參數精度高得多，仿真滿足遠場條件時，忽略誤差，僅考慮機載天線端測量參數的誤差，地面基站測得為。圖6為從單電磁矢量傳感器接收信號中估計飛行器姿態參數的仿真結果。地面站發射信號為320 MHz單色波。實驗仿真了在飛行器航向角旋轉90°，橫滾角上下擺動，俯仰角保持為0o的連續航姿變化過程。其間根據電磁波信號估計飛行器姿態角，作了36次姿態測量，進行每次估計的接收信號包含連續1024快拍，持續時間約，信號為部分極化波，極化橢圓率，連續快拍采樣頻率320 MHz。圖6(a)，圖6(b)，圖6(c)橫坐標為測量次數，縱坐標為橫滾角、俯仰角和航向角，坐標單位為°，圖中“”形線為理想角軌跡，“o”為本文最大似然估計方法估計得出的橫滾角、航向角和俯仰角，圖6(a)，圖6(b)，圖6(c)接收信號的極化度為60%,，信噪比10 dB。由圖6(a)，圖6(b)，圖6(c)可見本文最大似然估計方法估計得出的結果與理想值接近，在此條件下可以滿足普通飛行器控制的要求。圖6(d)，圖6(e)，圖6(f)為時，不同極化度情況下，前述仿真實驗的36次姿態角估計的標準差隨信噪比變化的曲線。圖6(d)，圖6(e)，圖6(f)橫坐標為接收信號的信噪比，縱坐標為36次實驗統計出的各姿態角參數估計標準差，單位為°。圖6(d)，圖6(e)，圖6(f)中信號的極化度分別為50%, 25%和75%，圖6(d)，圖6(e)，圖6(f)中“”，“×”和“?！本€分別為俯仰角、橫滾角和航向角參數估計標準差隨信噪比變化的曲線。圖6(d)，圖6(e)，圖6(f)可見各姿態角參數隨信噪比、極化度提高，估計誤差減小。估計精度在信噪比、極化度過小時，出現惡化。

6 結論

極化電磁波具備3維波結構，傳統的電磁波空域-極化域導向矢量信號模型割裂了自由度的空間姿態信息，不能對姿態參數進行整體優化，存在不足之處。本文提出的姿態導向矢量信號模型，姿態參數正交解耦，在估計信號參數時表現出更好的估計性能，為研究空域-極化域參數聯合估計問題提供了新思路。

部分極化波是電磁波存在的一般形式，本文將最大似然估計和MUSIC算法推廣應用于部分極化波參數估計。仿真數據結果可見，從單電磁矢量傳感器接收信號中估計飛行器姿態參數，精度達到飛行控制的要求，能作為自動駕駛儀航姿導航信息源。

附錄A 旋轉矩陣的指數表示

記式(A-1)為：

(A-2)

得：

由式(A-3)得：

(A-4)

繞坐標軸的3次連續旋轉矩陣相乘也可化為式(A-2)形式，從中計算出。

經過化簡正是式(A-1)。

附錄B 最大似然估計的分步算法[16,20]

同樣：

因此：

得：

因為：

所以：

代入式(A-1)得：

簡化為：

(A-7)

同理可得：

由此得到：

這樣：

(A-9)

根據式(A-9)：

因此：

根據

得：

得：

將各符號展開代入式(A-10)：

[1] Schmidt R O. Multiple emitter location and signal parameter Estimation[J]., 1986, 34(3): 276-280.

[2] Wong K T and Zoltowski M D. Closed-form direction finding and polarization estimation with arbitrarily spaced electromagnetic vector sensors at unknown locations[J]., 2000, 48(5): 671-681.

[3] Zoltowski M D and Wong K T. ESPRIT-based 2-D direction finding with a sparse uniform array of electromagnetic vector sensors[J]., 2000, 48(8): 2195-2204.

[4] Wong K T. Blind beamforming/geolocation for wideband- FFHs with unknown hop-sequences[J]., 2001, 37(1): 65-76.

[5] Roy R and Kailath T. ESPRIT-estimation of signal parameters via rotational invariance techniques[J].,, 1989, 37(7): 984-995.

[6] 徐振海. 極化敏感陣列信號處理的研究[D]. [博士論文], 國防科技大學, 2004.

Xu Zhen-hai. Study on the signal processing of polarized array[D]. [Ph.D. dissertation], National University of Defense Technology, 2004.

[7] Cho Am, Kang Young-shin,.. Altitude Integration of radar altimeter and GPS/INS for automatic takeoff and landing of a UAV[C]. 2011 11th International Conference on Control, Automation and Systems, KINTEX, Gyeonggi-do, 2011: 1429-1432.

[8] Wang Xiao-liang, Shao Xiao-wei, Gong De-ren,.. GPS/ VISNAV integrated relative navigation and attitude determination system for ultra-close spacecraft formation flying[J].2011, 22(2): 283-291.

[9] Wong K T and Zoltowski M D. Uni-vector-sensor ESPRIT for multisource azimuth, elevation, and polarization estimation[J]., 1997, 45(10): 1467-1474.

[10] Ko C C, Nehorai J, and Nehorai A. Separation and tracking of multiple broadband sources with one electromagnetic vector sensor[J]., 2002, 38(3): 1109-1116.

[11] Nehorai A and Paldi E. Vector-sensor array processing for electromagnetic source location[J]., 1994, 42(2): 376-398.

[12] Li Jian. On polarization estimation using a polarization sensitive array[C]. IEEE Proceedings Sixth Workshop on Satistical Signal and Array Processing, in Australia, 1992: 465-468.

[13] Wong K T. Direction finding/polarization estimation-dipole and/or loop triad(s)[J]., 2001, 37(2): 679-684.

[14] Shuster M D. A survey of attitude representations[J]., 1993, 41(4): 439-517.

[15] 張文典. 極化陣列信號參量估計方法研究[D]. [碩士論文], 吉林大學, 2007.

Zhang Wen-dian. Study on paremeters estimation of polarized array signal[D]. [Master dissertation], Jilin University, 2007.

[16] Marius P and Gershman A B. Maximum-likelihood direction-of-arrival estimation in the presence of unknown nonuniform noise[J]., 2001, 49(7): 1310-1324.

[17] 張小飛, 汪飛, 徐大專. 陣列信號處理的理論和應用[M]. 北京: 國防工業出版社, 2010: 41-88.

Zhang Xiao-fei, Wang Fei, and Xu Da-zhuan. The Thesis and Application of Array Signal Processing[M]. Beijing: National Defence Industry Press, 2010: 41-88.

[18] 鄧正隆. 慣性技術[M]. 哈爾濱: 哈爾濱工業大學出版社, 2006: 147-165.

Deng Zheng-long. Inertial Technology[M]. Harbin: Harbin Institute of Technology Press, 2006: 147-165.

[19] Hochwald B and Nehrai A. Concentrated Cramer-Rao bound expressions[J]., 1994, 40(2): 363-371.

[20] Stoica P and Nehorai A. MUSIC, maximum likelihood, and Cramer-Rao bound[J]., 1989, 37(5): 720-741.

Estimation of Aircraft Attitude/Heading Based on thePolarization of Electromagnetic Waves

Sun Hong-meiChen Guang-dongZhang Gong

(College of Electronic and Information Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China)(Unmanned Aerial Vehicle Research Department, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China)

This paper illustrates the full use of an orthogonal decoupling signal model to describe the 3-D construction of polarized electromagnetic waves. The maximum likelihood method and MUSIC algorithm were introduced into the parameter estimation of partially polarized waves. Attitude reference information for an aircraft that was calibrated by an information base station was delivered by polarized electromagnetic waves. Through receiving and processing the information of the aircraft, a single electromagnetic vector sensor can acquire the aircraft attitude in a geographic coordinate system. Differing from the triangle calculation method which relies on a multi-point measurement, the proposed method only needs a base station signal and a single receiver on the motion platform to realize the posture perception of aircraft. It can therefore serve as a substitute foraircraft attitude navigation and aircraft heading navigation. The introduction of the motion platform attitude estimation algorithm provides a technology support for engineering applications.

Elecrtomagnetic vector sensor; State of polarization; Direction Of Arrival (DOA); Attitude determination; Partially polarized signal

TN911

A

2095-283X(2013)04-0466-10

10.3724/SP.J.1300.2013.13087

2013-09-22收到，2013-12-05改回；2013-12-11網絡優先出版

國家自然科學基金(61271327, 61071163)資助課題

孫紅梅 876520283@qq.com

孫紅梅(1990-)，女，江蘇鹽城人，為南京航空航天大學碩士研究生，主要研究方向為通信與信息系統。

E-mail: 876520283@qq.com

陳廣東(1968-)，男，副研究員，研究方向為雷達與測控。

張 弓(1964-)，男，教授，博士生導師，主要研究方向為雷達信號處理、目標探測與識別、圖像分析與處理等。