基于遺傳算法的相關(guān)干涉儀測(cè)向方法*

楊 卓

(中國(guó)電子科技集團(tuán)公司第36研究所,浙江嘉興 314033)

基于遺傳算法的相關(guān)干涉儀測(cè)向方法*

楊 卓

(中國(guó)電子科技集團(tuán)公司第36研究所,浙江嘉興 314033)

相關(guān)干涉儀體制可以獲得未知輻射源電磁波的相位信息,借助相關(guān)算法解算出波達(dá)角,在工程上得到了廣泛應(yīng)用。但天線(xiàn)相位往往會(huì)受平臺(tái)影響發(fā)生畸變,并導(dǎo)致測(cè)向偏差。文中基于干涉儀體制,采用遺傳算法對(duì)測(cè)向系統(tǒng)的基線(xiàn)和權(quán)值進(jìn)行了優(yōu)化。結(jié)合典型實(shí)例,比較分析了優(yōu)化前后測(cè)向誤差的變化。證明該方法能夠顯著改善測(cè)向精度,提出的方法對(duì)機(jī)載測(cè)向具有實(shí)用意義。

相關(guān)干涉儀;測(cè)向;遺傳算法

0 引言

干涉儀測(cè)向作為無(wú)源偵察定位的一項(xiàng)關(guān)鍵技術(shù),因其測(cè)向靈敏度高、精度高和速度快的特點(diǎn)在軍事和民用領(lǐng)域都有著重要而廣泛的應(yīng)用。特別是近年來(lái),隨著我國(guó)航空航天平臺(tái)的井噴,基于空、天平臺(tái)的干涉儀測(cè)向得到了越來(lái)越多的應(yīng)用。

目前已經(jīng)有不少文獻(xiàn)針對(duì)干涉儀測(cè)向進(jìn)行了研究。文獻(xiàn)[1]對(duì)圓陣干涉儀解模糊性能進(jìn)行了系統(tǒng)分析,根據(jù)粗測(cè)入射角選擇基線(xiàn)組以實(shí)現(xiàn)高精度無(wú)模糊角度測(cè)量。文獻(xiàn)[2]分析了旋轉(zhuǎn)長(zhǎng)基線(xiàn)干涉儀測(cè)向的相位差觀(guān)測(cè)模型,進(jìn)而推導(dǎo)了無(wú)模糊測(cè)向條件。文獻(xiàn)[3]采用了虛擬陣列變換的方法,提高了入射信號(hào)的測(cè)向精度。文獻(xiàn)[4]提出一種基于時(shí)變長(zhǎng)基線(xiàn)虛擬2維干涉儀的偽線(xiàn)性最小二乘測(cè)向方法,測(cè)向精度接近克拉美-羅下限。文獻(xiàn)[5]深入探討了導(dǎo)致機(jī)載干涉儀測(cè)向誤差的外因和內(nèi)因。

上述文獻(xiàn)主要是對(duì)相關(guān)干涉儀的算法進(jìn)行了研究和改進(jìn),提高了測(cè)向精度,但計(jì)算量和復(fù)雜度有所提升,且未考慮安裝平臺(tái)對(duì)測(cè)向系統(tǒng)的影響。受平臺(tái)的安裝限制(特別是布局緊湊的中小型平臺(tái)),飛機(jī)機(jī)體、突出機(jī)體的任務(wù)設(shè)備很有可能對(duì)測(cè)向系統(tǒng)的接收相位產(chǎn)生影響,從而導(dǎo)致某些頻段發(fā)生測(cè)向模糊。由于很難定性定量評(píng)估飛機(jī)機(jī)體對(duì)相位的影響,一般是采取飛行校表的形式對(duì)測(cè)向模糊頻段進(jìn)行校正,費(fèi)時(shí)費(fèi)力。文中采用遺傳算法對(duì)干涉儀測(cè)向體制中的關(guān)鍵參量進(jìn)行優(yōu)化,得到最優(yōu)基線(xiàn)組合和權(quán)值,從而減少測(cè)向模糊現(xiàn)象。

1 基本原理

1.1 干涉儀測(cè)向原理

干涉儀測(cè)向的基本原理是利用相位匹配的多個(gè)信道,測(cè)量多個(gè)天線(xiàn)陣元接收到的來(lái)自同一個(gè)目標(biāo)輻射源的到達(dá)相位差[6-7]。基于實(shí)際工程需要,可以根據(jù)相位差與角度的函數(shù)關(guān)系計(jì)算出目標(biāo)輻射源的方位角度值;也可以將理論相位差與實(shí)際到達(dá)相位差進(jìn)行相關(guān)計(jì)算,以求得方位角度值。兩種方式并無(wú)本質(zhì)不同,文中以相關(guān)算法進(jìn)行分析。

以三天線(xiàn)為例,T1、T2、T3之間的距離(干涉儀基線(xiàn))為D12、D13、D23,若有波長(zhǎng)為λ的輻射源平面波由θ方向傳播而來(lái)(與天線(xiàn)軸法線(xiàn)夾角),它的相位差必為:

(1)

(2)

(3)

假設(shè)測(cè)向系統(tǒng)針對(duì)波長(zhǎng)為λ的輻射波,存儲(chǔ)了0°～360°范圍內(nèi)(假設(shè)等分為N份)的理論相位差組合[φ12,iφ13,iφ23,i],將實(shí)際相位差與理論相位差反復(fù)做相關(guān)計(jì)算,如式(4)、式(5)所示:

式中A12、A13、A23是相關(guān)運(yùn)算的權(quán)值。當(dāng)C(i)最小時(shí),β即是計(jì)算得到的入射角。

1.2 遺傳算法的應(yīng)用

基于此,擬采用遺傳算法對(duì)基線(xiàn)組合和權(quán)值進(jìn)行尋優(yōu)。遺傳算法具有較好的全局搜索性能,其優(yōu)化函數(shù)不受連續(xù)可微的約束,這兩個(gè)特點(diǎn)使其與復(fù)雜算法結(jié)合成為了可能。遺傳算法的基本概念不再贅述,具體和干涉儀相結(jié)合,主要分為5個(gè)過(guò)程[8],每個(gè)過(guò)程中需要注意的事項(xiàng)如下:

1)編碼和約束:首先對(duì)待優(yōu)化的基線(xiàn)和權(quán)值進(jìn)行編碼。對(duì)于m元陣來(lái)說(shuō),樣本個(gè)體由2m個(gè)變量組成,分別是{a1,a2,…,am,b1,b2,…,bm}。其中,b1,b2,…,bm代表基線(xiàn)配置,假設(shè)天線(xiàn)陣為7元陣,則兩兩天線(xiàn)之間有21種基線(xiàn)配置。因此b1,b2,…,b7的值為正整數(shù),且只可能在1到21之間取值;a1,a2,…,a7代表相應(yīng)的基線(xiàn)權(quán)重,只可能在0到1之間取值,且有如下約束條件:

a1+a2+a3+a4+a5+a6+a7=1

(6)

權(quán)重采用L位二進(jìn)制編碼表示,則每個(gè)樣本的染

色體可表示為:

(7)

2)初始群體生成:可以將部分有序的初始個(gè)體放入初始群體(例如基線(xiàn)對(duì)稱(chēng)分布,權(quán)值平均分布),提高進(jìn)化的收斂速度。

3)適應(yīng)度函數(shù)生成:適應(yīng)度函數(shù)標(biāo)明個(gè)體進(jìn)化的優(yōu)劣度。將解碼后的{a1,a2,…,a7,b1,b2,…,b7}代入1.1節(jié)式(4)、式(5)計(jì)算,經(jīng)全頻段、全方位統(tǒng)計(jì),并構(gòu)造如下適應(yīng)度函數(shù)(式中,β代表解算角度,DOA代表真實(shí)入射角):

(8)

(9)

4)進(jìn)化:主要包含“選擇”、“變異”、“交叉”等操作,不斷迭代產(chǎn)生下一代種群。

5)最終輸出:若判定進(jìn)化完成,則對(duì)最優(yōu)個(gè)體進(jìn)行解碼,輸出最終的基線(xiàn)和權(quán)值。

2 典型實(shí)例

如圖1所示,是某設(shè)備采用的不規(guī)則七元測(cè)向陣,對(duì)225～450 MHz進(jìn)行測(cè)向。其中1#、3#、4#、5#、6#天線(xiàn)分布在長(zhǎng)軸為1.2 m、短軸為0.9 m的橢圓和七單元等分線(xiàn)的交叉點(diǎn)上,2#、7#號(hào)天線(xiàn)相對(duì)橢圓有所偏移。

圖1 七元天線(xiàn)陣布局

根據(jù)實(shí)際飛行,在225～450 MHz頻段范圍內(nèi),0°～360°方位范圍內(nèi)共采集了8 015個(gè)點(diǎn)的相位組數(shù)據(jù)。

根據(jù)傳統(tǒng)的干涉儀分配方式,7條基線(xiàn)主要有3種分配方式:

a)基線(xiàn)配置為[1-2,2-3,3-4,4-5,5-6,6-7,1-7],權(quán)值均配為1/7;

b)基線(xiàn)配置為[1-3,2-4,3-5,4-6,5-7,1-6,2-7],權(quán)值均配為1/7;

c)基線(xiàn)配置為[1-4,2-5,3-6,4-7,1-5,2-6,3-7],權(quán)值均配為1/7。

認(rèn)為測(cè)向精度大于3°為測(cè)向出錯(cuò)點(diǎn),經(jīng)過(guò)相關(guān)計(jì)算,a)、b)、c)3種方式下,出錯(cuò)點(diǎn)分別是965個(gè)、1 057個(gè)、1 064個(gè),3種模式的測(cè)向出錯(cuò)點(diǎn)差異不大。測(cè)向精度色溫圖如圖2所示(白色代表測(cè)向出錯(cuò)點(diǎn)、黑色代表測(cè)向正常點(diǎn))。

圖2 優(yōu)化前測(cè)向精度色溫圖

將a)、b)、c) 3種模式作為初代種群中的一部分,進(jìn)行遺傳算法優(yōu)化,經(jīng)過(guò)10代后,進(jìn)化過(guò)程如圖3(a)所示。最優(yōu)值為:

d)基線(xiàn)配置:[1-3,3-4,4-5,5-6,1-6,2-7,2-6];權(quán)值:A13=0.118 5,A34=0.118 5,A45=0.118 5,A56=0.118 5,A16=0.118 5,A27=0.266 6,A26=0.141 1;出錯(cuò)點(diǎn)降為616。

圖3 優(yōu)化過(guò)程評(píng)估

將a)、b)、c)、d) 4種模式作為初代種群中的一部分,再次進(jìn)行遺傳算法優(yōu)化,經(jīng)過(guò)50代后,進(jìn)化過(guò)程如圖3(b)所示。最優(yōu)值為:

e)基線(xiàn)配置:[1-3,3-4,4-5,5-6,1-6,2-7,2-6];權(quán)值:A13=0.155 9,A34=0.108 3,A45=0.117 1,A56=0.135 8,A16=0.131 0,A27=0.243 6,A26=0.108 3;出錯(cuò)點(diǎn)降為593。

比較d)和e)的優(yōu)化結(jié)果,發(fā)現(xiàn)兩種模式下的基線(xiàn)配置一致,僅權(quán)值發(fā)生了微調(diào),其適應(yīng)度函數(shù)值也無(wú)大的變動(dòng),說(shuō)明其優(yōu)化迭代過(guò)程已趨穩(wěn)定。

圖4 優(yōu)化后測(cè)向精度色溫圖

優(yōu)化后的測(cè)向精度色溫圖如圖4所示(白色代表測(cè)向出錯(cuò)點(diǎn)、黑色代表測(cè)向正常點(diǎn)),可以看到,經(jīng)優(yōu)化后測(cè)向出錯(cuò)點(diǎn)明顯得到了改善。

3 結(jié)論

優(yōu)化后的基線(xiàn)配置如下圖5所示,1-3,3-4,4-5,5-6,1-6,2-7構(gòu)成了一個(gè)六邊對(duì)稱(chēng)結(jié)構(gòu)。其中2-7的權(quán)值明顯高于其他基線(xiàn),說(shuō)明其測(cè)向糾錯(cuò)作用較大。由于系統(tǒng)選擇了7條基線(xiàn)進(jìn)行計(jì)算,因此2-6作為補(bǔ)充基線(xiàn)進(jìn)入運(yùn)算。

圖5 優(yōu)化后測(cè)向基線(xiàn)配置

通過(guò)遺傳算法對(duì)干涉儀測(cè)向參數(shù)進(jìn)行優(yōu)化,得到了一組最佳基線(xiàn)和權(quán)值,測(cè)向出錯(cuò)點(diǎn)減少了40%以上。這種方法對(duì)于實(shí)際工程具有較大意義:

1)在過(guò)去,測(cè)向的基線(xiàn)和權(quán)值均是根據(jù)理論計(jì)算或者經(jīng)驗(yàn)進(jìn)行分配。結(jié)合實(shí)測(cè)相位值,文中給出了一種測(cè)向參數(shù)尋優(yōu)的策略,可以明顯降低測(cè)向出錯(cuò)概率。

2)在工程應(yīng)用上,通過(guò)暗室或者外場(chǎng)采集的數(shù)據(jù),即可對(duì)測(cè)向進(jìn)行尋優(yōu),得到定性定量結(jié)論,從而減少飛行校表的工作量。

3)下一步,在充分采集天線(xiàn)在各位置相位畸變的基礎(chǔ)上,甚至有可能對(duì)天線(xiàn)陣型進(jìn)行參數(shù)化,利用遺傳算法尋找最優(yōu)陣型,這也是我們今后可以努力的方向。

[1] 張亮, 徐振海, 熊子源, 等. 基于圓陣干涉儀的被動(dòng)導(dǎo)引頭寬帶測(cè)向方法 [J]. 系統(tǒng)工程與電子技術(shù), 2012, 34(3): 462-466.

[2] 張敏, 郭福成, 李騰, 等. 旋轉(zhuǎn)長(zhǎng)基線(xiàn)干涉儀測(cè)向方法及性能分析 [J]. 電子學(xué)報(bào), 2013, 41(12): 2422-2429.

[3] 吳奉微, 程婷, 賈可新, 等. 基于虛擬陣列變換的干涉儀測(cè)向算法 [J]. 現(xiàn)代雷達(dá), 2012, 34(3): 42-46.

[4] 張敏, 郭福成, 周一宇, 等. 時(shí)變長(zhǎng)基線(xiàn)2維干涉儀測(cè)向方法 [J]. 電子與信息學(xué)報(bào), 2013, 35(12): 2882-2888.

[5] 毛虎, 楊建波, 邱宏坤. 機(jī)載二維干涉儀測(cè)向誤差分析 [J]. 火力與指揮控制, 2011, 36(12): 67-71.

[6] WU Y W, RHODES S, SATORIUS E H. Direction of arrival estimation via extended phase interferometry [J]. IEEE Transaction on Aerospace and Electronic Systems, 1995, 31(1): 375-381.

[7] MURPHY S, EYRING P M. 2～18 GHz circular array interferometer DF antenna system [C]∥IEEE Antennas and Propagation Society, International Symposium, 1998, 4: 2332-2335.

[8] 楊卓, 薛正輝, 李偉明. 遺傳算法結(jié)合矩量法設(shè)計(jì)彈上小型寬帶天線(xiàn) [J]. 彈箭與制導(dǎo)學(xué)報(bào), 2009, 29(6): 231-234.

Direction Finding Method of Correlative Interferometer Based on Genetic Algorithm

YANG Zhuo

(No.36 Research Institute of China Electronics Technology Group Corporation, Zhejiang Jiaxing 314033, China)

Correlative interferometer could obtain unknown radiation source electromagnetic wave phase information, to calculate the arrival angle with the relevant algorithm, which had been widely used in engineering. However, the phase of the antenna was often distorted by platform impact, which leaded to the deviation of direction finding. Based on the interferometer system, the genetic algorithm was used to optimize the baseline and weight of the direction finding system. Combined with typical examples, the change of direction finding error before and after optimization was analyzed and compared. It was proved that this method could improve the direction finding precision, and the proposed method had practical significance fore airborne direction finding.

correlative interferometer; direction finding; genetic algorithm

2015-11-07

楊卓(1982-),男,貴州銅仁人,高級(jí)工程師,博士,研究方向:通信系統(tǒng)集成。

TN971.+4

A