石　榮，李　瀟，劉　暢

(電子信息控制重點(diǎn)實(shí)驗(yàn)室,　成都 610036)

?

·信號處理·

基于矢量合成的相干信號干涉儀測向模型

石榮，李瀟，劉暢

(電子信息控制重點(diǎn)實(shí)驗(yàn)室,成都 610036)

基于通道間相位差測量的傳統(tǒng)干涉儀測向模型在對相干多信號進(jìn)行測向時會產(chǎn)生較大的誤差，甚至測向失效，在分析其失效原因的基礎(chǔ)上，利用相干多信號矢量合成方法，根據(jù)干涉儀基線幾何結(jié)構(gòu)，構(gòu)建了新的相干信號干涉儀測向模型，并詳細(xì)討論了模型的求解過程。該模型不僅能夠?qū)ο喔啥嘈盘枌?shí)施測向，而且通過對求解結(jié)果的分析可推斷產(chǎn)生相干效應(yīng)的信號個數(shù)，也能實(shí)現(xiàn)對傳統(tǒng)單信號的測向，即具有較好的向下兼容性。最后通過仿真驗(yàn)證了模型的有效性，從而為干涉儀測向理論的發(fā)展和抗干擾應(yīng)用提供了重要參考。

干涉儀；相干多信號；矢量合成；信號個數(shù)判斷；測向模型；向下兼容性

0　引　言

采用干涉儀對射頻信號的來波方向進(jìn)行測量在射電天文、電子對抗和頻譜監(jiān)測等領(lǐng)域獲得了十分廣泛的應(yīng)用，這其中所使用的測向模型大多以干涉儀各通道間的相位差測量為基礎(chǔ)，從而使得在同一個工作頻率上瞬時只能對一個信號實(shí)施測向。如果有兩個或兩個以上的同頻復(fù)制的信號同時到達(dá)干涉儀，它們相互之間就會產(chǎn)生干涉，通常稱之為相干信號。傳統(tǒng)干涉儀測向模型在對相干信號進(jìn)行測向時，通常會產(chǎn)生較大的誤差，甚至完全失效。利用此現(xiàn)象文獻(xiàn)[1]提出對干涉儀的測向過程實(shí)施相干干擾，從而達(dá)到擾亂電子偵察系統(tǒng)的目的。到目前為止雖然有大量的文獻(xiàn)針對陣列天線對相干信號的測向問題進(jìn)行過研究與討論，所提出的方法主要包括：基于空間平滑的各種MU-SIC多重信號分類算法[2-5]、基于各種去相關(guān)的ESPRIT旋轉(zhuǎn)不變子空間算法[6-9]、基于DFT離散傅里葉變換及變參考陣元特征法[10]、基于虛擬陣列的算法[11]、基于壓縮感知的方法[12]等。但是，上述文獻(xiàn)中所提出的方法如何在干涉儀測向中進(jìn)行應(yīng)用并沒有相關(guān)的分析與探討。

針對這一問題，本文分析了傳統(tǒng)干涉儀測向模型在對相干信號測向時產(chǎn)生失效的原因，在此基礎(chǔ)上利用同時到達(dá)干涉儀的多個相干信號具有相同頻率這一特點(diǎn)，采用多信號矢量合成方法，建立了新的理論模型，通過對其求解結(jié)果的分析，可以對同時到達(dá)干涉儀的相干信號的個數(shù)進(jìn)行推斷。如果判定的結(jié)果是當(dāng)前到達(dá)干涉儀的信號個數(shù)只有一個時，該模型所得測向結(jié)果與傳統(tǒng)針對單信號的干涉儀測向模型所得結(jié)果完全一致，所以該模型具有較好的向下兼容性。相關(guān)的理論分析和仿真驗(yàn)證詳細(xì)闡述如下。

1　傳統(tǒng)干涉儀測向模型

假設(shè)在所關(guān)注的時段內(nèi)干涉儀只接收到一個信號S1(t)，其波長為λ，從與天線視軸夾角為θ1的方向到達(dá)測向天線A和B，兩天線之間的距離為d，如圖1所示。

圖1　單基線干涉儀測向

于是干涉儀兩個接收通道中的信號分別為

SA,1(t)=a1·sin(2πfct+φ1)+nA(t)

(1)

SB,1(t)=a1·sin(2πfct+φ1+φAB,1)+nB(t)

(2)

式中：fc為信號的載波頻率；φ1為初始相位；a1為信號幅度；φAB，1為信號S1(t)到達(dá)兩天線的路徑不同而引起的相位差；nA(t)與nB(t)分別是兩個接收通道引入的噪聲。在一定信噪比條件下，由式(1)、式(2)可求解出φAB，1，于是信號的來波方向θ1為

(3)

由式(3)可知，在單基線干涉儀測向應(yīng)用中為了避免相位差測量過程中的相位模糊，要求d≤λ/2。如果需要進(jìn)一步提高測向精度，可使用多基線干涉儀，即用長基線獲得較高精度的相位差測量值，用短基線來解相位差模糊。但無論干涉儀的長短基線如何設(shè)計(jì)，傳統(tǒng)求解方法仍是通過通道間信號的相位差測量來實(shí)現(xiàn)的。

2　在相干信號條件下傳統(tǒng)模型的失效性分析

對于單頻正弦波信號而言，當(dāng)兩個信號具有相同頻率時，這就成為一對相干信號。如果相干信號S1(t)與S2(t)同時到達(dá)干涉儀的測向天線A和B，且其來波方向分別為θ1和θ2，如圖2所示。

圖2　干涉儀對兩個同時到達(dá)信號的測向

此情況下干涉儀兩接收通道中的信號分別表示為

SA，2(t)=a1·sin(2πfct+φ1)+

a2·sin(2πfct+φ2)+nA(t)

(4)

SB,2=a1·sin(2πfct+φ1+φAB,1)+

a2·sin(2πfct+φ2+φAB,2)+nB(t)

(5)

式中：φ2為信號S2(t)的初始相位；a2為其幅度；φAB,2為信號S2(t)到達(dá)兩天線的路徑不同而引起的相位差。由式(3)可知φAB,1和φAB,2由下式所決定

φAB,1=2πdsin(θ1)/λ

(6)

φAB,2=2πdsin(θ2)/λ

(7)

實(shí)際上這兩個相干信號在同一個接收通道內(nèi)可以合成一個穩(wěn)定的同頻信號，即

SA,2(t)=a1·γA·sin(2πfct+φA)+nA(t)

(8)

SB,2(t)=a1·γB·sin(2πfct+φB)+nB(t)

(9)

式中：γA、γB分別為A、B通道中信號相干合成后的幅度變化系數(shù)；φA、φB分別為信號相干合成后的相位，上述四個參數(shù)可由式(4)、式(5)計(jì)算如下

(10)

(11)

(12)

(13)

由上可見，在兩個相干信號合成之后干涉儀兩個通道所接收到的合成信號的幅度在大多數(shù)情況下是不相等的，即γA≠γB；另一方面，干涉儀兩個通道所接收到的合成信號之間的相位差φAB,T為

φAB,T=φA-φB

(14)

在兩個相干信號到達(dá)角不相同時，相位差φAB,T既不等于φAB,1，也不等于φAB,2。顯然利用φAB,T來推算出的來波方向角既不等于θ1，也不等于θ2。雖然上述分析是以單基線干涉儀為基礎(chǔ)進(jìn)行的，但對于多基線干涉儀來說，上述分析流程同樣適用，所以傳統(tǒng)干涉儀測向方法對兩個相干信號的測向錯誤在理論上就是不可避免的。

3　基于矢量合成的相干多信號干涉儀測向模型

實(shí)際上由式(4)、式(5)所表達(dá)的相干信號的時域信號模型可知，相干信號具有相同的載頻fc，如果將時域信號模型轉(zhuǎn)換為復(fù)平面上的矢量信號模型，雖然每一個信號矢量都以相同的角速度2πfc繞原點(diǎn)旋轉(zhuǎn)，但各個信號矢量的相對位置關(guān)系是保持不變的，在復(fù)平面上以信號的幅度和相位為基礎(chǔ)，將產(chǎn)生相干的N個信號Si,c表示成復(fù)矢量形式如下

Si,c=ai·exp(jφi) ，i=1,2,…,N

(15)

式中：ai表示第i個信號的幅度；φi表示第i個信號到達(dá)干涉儀中參考單元天線時的相位，如圖3所示。

圖3　第i個信號到達(dá)干涉儀的方位關(guān)系圖示

圖3中的1#天線為參考天線單元，θi是第i個信號Si(t)的來波方向與干涉儀法向之間的夾角，d1k表示1#天線單元與k#天線單元之間的距離。于是所有相干信號在1#天線處所合成的綜合信號矢量S1,z可以表示為

(16)

同理，所有相干信號在k#天線(k=2, 3, …,M，M為干涉儀單元天線的總數(shù))處所合成的綜合信號矢量Sk,z可以表示為

(17)

式(17)中的附加相位項(xiàng)是由于信號到達(dá)不同天線單元時的延遲所產(chǎn)生的。為了形成統(tǒng)一的數(shù)學(xué)表達(dá)式，可從形式上定義1#天線與其自身的距離d11=0，于是基于矢量合成的相干多信號干涉儀測向理論模型可統(tǒng)一表示為

k=1,2,…,M

(18)

式中：nk為第k#天線測量時所引入的噪聲矢量。

4　測向模型的求解及相干信號個數(shù)的判斷

在式(18)所表達(dá)的模型中待求解的未知數(shù)的個數(shù)為3N；而方程的個數(shù)為2M，在不考慮通道噪聲影響的情況下，滿足如下條件時該模型才有解

2M≥3N

(19)

(20)

式中：‖·‖p表示p范數(shù)，在p=2時，式(20)則表示均方誤差最小，diff(k)表示第k#天線單元對最終解所引入的噪聲影響項(xiàng)，即

(21)

顯然滿足式(20)的最優(yōu)解可通過對整個解空間進(jìn)行優(yōu)化搜索的方法來得到。但上述求解過程中，還隱含了一個假設(shè)，即已知同時到達(dá)干涉儀的相干信號的個數(shù)N。實(shí)際上在求解式(20)之前，N是一個未知數(shù)，因?yàn)槲覀儾⒉恢浪玫降臏y量數(shù)據(jù)中究竟有多少個相干信號。所以在求解該式之前，需要利用式(19)對N的上限值Nmax做一個估計(jì)

Nmax≤?2M/3」

(22)

式中：?X」表示不超過X的最大整數(shù)。首先，假設(shè)有Nmax個相干信號同時到達(dá)干涉儀；然后，對式(20)進(jìn)行求解，即可得到Nmax個信號參數(shù)解集合{ai,φi,θi|i=1,2,…,Nmax}。在此基礎(chǔ)上，對上述Nmax個解按照對應(yīng)的信號幅度ai的大小由大至小排序如下

aj1≥aj2≥…≥ajNmax

(23)

如果在上述信號個數(shù)推斷過程中，直接得到的結(jié)果是Nth=1，這說明此時干涉儀所接收到的是單信號。在此條件下既可按照傳統(tǒng)的通道間相位差測量方式來獲得該信號的來波方向，也可以通過求解式(20)來獲得該信號的來波方向，在此情況下兩種方法所得到的解應(yīng)具有一致性。這也同時說明，新模型相對于傳統(tǒng)模型來講，更具有普適性，同時對傳統(tǒng)模型也具有向下兼容性。

5　仿真驗(yàn)證

在如下仿真過程中采用(幅度，初相，來波方向)的三元組形式來進(jìn)行信號參數(shù)和求解結(jié)果的描述，其中幅度為相對值，初相以到達(dá)干涉儀第一個單元天線位置處為參考，來波方向以與干涉儀天線陣的法向夾角來表達(dá)。

1)干涉儀對兩個相干信號的測向

仿真條件：一維干涉儀由5個單元天線組成，它們之間的間距分別為0.05m，0.20m，0.80m，3.00m，兩個信號的載波頻率均為3GHz，其信號參數(shù)分別以三元組形式表示為(1.0，40°，50.1°)，(1.6，110°，19.9°)，如圖4所示。

圖4　干涉儀及兩相干信號來波方向示意圖

按照前文所提出的信號矢量表示方法，以1#天線作為參考天線，可以測量得到在五個單元天線處的綜合信號矢量如圖5所示。

圖5　干涉儀五個天線所測得的雙信號合成矢量圖示

圖5中各信號矢量均添加了標(biāo)準(zhǔn)差為0.02的高斯白噪聲。根據(jù)式(22)，在此條件下最多可實(shí)現(xiàn)三個相干信號的測向。于是首先假設(shè)三個不同方向的同頻信號同時到達(dá)，采用式(20)進(jìn)行求解。由于整個解空間比較大，所以在求解搜索策略上采用了由粗至精的兩步搜索法，粗搜索的步進(jìn)較大，幅度步進(jìn)為0.25，相位步進(jìn)為30°，來波方向步進(jìn)為5°。通過粗搜索可得到幾個候選結(jié)果，其中一個表示為(1.5，120°，20°)；(1，30°，50°)；(0.25，60°，15°)。在幾個候選結(jié)果附近區(qū)域都做精搜索，精搜索的幅度步進(jìn)為0.1，相位步進(jìn)為5.0°，來波方向步進(jìn)為0.2°。通過精搜索所得到的結(jié)果為(1.6，115.0°，20.0°)；(1.1，35.0°，50.2°)；(0.1，55.0°，15°)。

上述結(jié)果是在三個同時到達(dá)的相干信號假設(shè)條件下得到的，其信號幅度值分別為1.6，1.1，0.1。前兩個數(shù)值明顯大于第三個，且第三個值趨近于零，由此判斷：所接收到的相干信號的個數(shù)為2。在此條件下，采用兩個信號的合成模型再對式(20)進(jìn)行求解，最終可得這兩個相干信號的來波方向分別為19.88°，50.02°；誤差分別為-0.02°，-0.08°。

2)干涉儀對單信號的測向

除了信號個數(shù)減少至一個之外，即信號參數(shù)為(1.0，40.0°，50.1°)，其他仿真條件同前。同樣可測量得到五個綜合信號矢量如圖6所示。

圖6　干涉儀五個天線所測得的單信號矢量圖示

由圖6可見，五個信號矢量的幅度幾乎全部一樣，這也滿足單信號到達(dá)干涉儀時的特征。同樣通過解空間搜索，得到的結(jié)果以三元組形式表示為(1.0，40.0°，50.0°)，(0，0，0)，(0，0，0)。顯然此時只有一個信號到達(dá)干涉儀天線陣。在一個信號條件下，再對式(20)進(jìn)行求解，從而可得該信號的來波方向?yàn)?0.09°，誤差為-0.01°。

6　結(jié)束語

本文針對干涉儀應(yīng)用中的相干多信號測向問題，提出了一個新的基于矢量合成的理論模型，該模型消除了傳統(tǒng)的基于通道間相位差計(jì)算而帶來的相干信號情況下的測向失效條件，通過解空間搜索優(yōu)化方法可直接求解出各個相干信號的來波方向。在求解過程中，不僅可實(shí)現(xiàn)對相干信號個數(shù)的推斷，而且還能得到產(chǎn)生相干的各個信號的幅度大小和相對相位值。另一方面，該模型對單信號測向同樣適用，所以相對于傳統(tǒng)方法來說具有向下兼容性。相關(guān)的仿真結(jié)果也印證了理論分析的有效性。雖然本文對相干信號源建模時以單頻信號為對象展開討論，實(shí)際上對于調(diào)制信號，如相位編碼、調(diào)頻信號等，在對信號幅值與初相參數(shù)測量之后，也可采用此方法進(jìn)行測向。綜上所述，本文從理論和仿真上對基于矢量合成的相干信號干涉儀測向方法進(jìn)行了分析，后續(xù)還需繼續(xù)開展應(yīng)用方面的研究工作，以解決工程中相關(guān)的應(yīng)用性問題。

[1]楊軍佳, 畢大平, 陳慧. 相位法測向系統(tǒng)相干干擾效果分析[J]. 電子信息對抗技術(shù), 2012, 27(6): 47-49.

YANG Junjia, BI Daping, CHEN Hui. Coherent jamming effect analysis on interferometer[J]. Electronic Information Warfare Technology, 2012, 27(6): 47- 49.

[2]司錫才，甄佳奇，那振宇. 二維相干信號測向估計(jì)方法[J]. 系統(tǒng)工程與電子技術(shù), 2009, 31(8): 1982-1984.

SI Xicai, ZHEN Jiaqi, NA Zhenyu. Approach to 2-D angles direction estimation for coherent signals[J]. Systems Engineering and Electronics, 2009, 31(8): 1982-1984.

[3]肖雷，管振輝，楊春華. 基于MUSIC算法對相干信號DOA估計(jì)的研究[J]. 艦船電子對抗, 2010, 33(6): 114-117.

XIAO Lei, GUAN Zhenhui, YANG Chunhua. Research into the DOA estimation of coherent signals based on MUISC algorithm[J]. Shipboard Electronic Countermeasure, 2010, 33(6): 114-117.

[4]吳小強(qiáng)，刁鳴，繆善林. 基于內(nèi)插陣列變換的相干源測向方法研究[J]. 應(yīng)用科技, 2007, 34(2): 1-4.

WU Xiaoqiang, DIAO Ming, MIAO Shanlin. Study on DOA estimation of coherent sources based on interpolated array transformation[J]. Applied Science and Technology, 2007, 34(2): 1-4.

[5]張聰，胡謀法，盧煥章. 基于虛擬陣列空間平滑的相干信號DOA估計(jì)[J]. 電子學(xué)報, 2010, 38(4): 929-933.

ZHANG Cong, HU Moufa, LU Huanzhang. Virtual array-based spatial smoothing method for direction finding of coherent signals[J]. Acta Electronica Sicinca, 2010, 38(4): 929-933.

[6]何勁，嵇瑋瑋，劉衛(wèi)平，等. 基于線性tripole陣列的相干信號角度和極化估計(jì)與跟蹤[J]. 上海航天, 2015, 32(1): 17- 21.

HE Jin, JI Weiwei, LIU Weiping, et al. Angle-polarization estimation & tracking for coherent signals using a linear tripole array[J]. Aerospace Shanghai. 2015, 32(1): 17-21.

[7]甄佳奇，司錫才，王桐，等. 基于虛擬陣列的相干信號快速測向算法[J]. 吉林大學(xué)學(xué)報(工學(xué)版), 2010, 40(3): 848-851.

ZHEN Jiaqi, SI Xicai, WANG Tong, et al. Fast approach for direction estimation coherent signals based on virtual array[J]. Journal of Jilin University (Engineering and Technology Edition), 2010, 40(3): 848-851.

[8]張輝，葛臨東. 利用時間平滑技術(shù)估計(jì)相干信號二維波達(dá)方向[J]. 信息工程大學(xué)學(xué)報, 2007, 8(1): 68-71.

ZHANG Hui, GE Lindong. UCA based 2-D DOA estimation of multi-source multipath signals[J]. Journal of Information Engineering University, 2007, 8(1): 68-71.

[9]杜剛，張永順，姜新迎. 相干信號源的空時三維參數(shù)分離估計(jì)方法[J]. 現(xiàn)代雷達(dá)，2008, 30(3): 67-69, 73.

DU Gang, ZHANG Yongshun, JIANG Xinying. Space-time three-dimension parameter separable estimation method of coherent signals[J]. Modern Radar, 2008, 30(3): 67-69, 73.

[10]林亭，曾新吾，孫海洋. 離散傅里葉變換及變參考陣元特征法在相干信號估計(jì)中的應(yīng)用[J]. 信號處理, 2012, 28(3): 404-409.

LIN Ting, ZENG Xinwu, SUN Haiyang. DFT and shift reference sensor method applied in coherent signal parameter estimation[J]. Signal Processing, 2012, 28(3): 404-409.

[11]甄佳奇，丁群，趙冰. 虛擬陣列下的相干信號測向算法[J]. 系統(tǒng)工程與電子技術(shù), 2013, 35(10): 2032-2035.

ZHEN Jiaqi, DING Qun, ZHAO Bing. Direction finding algorithm for coherent signals by virtual array[J]. System Engineering and Electronics, 2013, 35(10): 2032-2035.

[12]刁鳴，安春蓮. 獨(dú)立信號與相干信號并存的二維DOA估計(jì)新方法[J]. 西安電子科技大學(xué)學(xué)報(自然科學(xué)版), 2013, 40(5): 66-71,98.DIAO Ming, AN Chunlian. 2-D DOA estimation of coexisting uncorrelated and coherent signals[J]. Journal of Xidian University(Natural Science), 2013, 40(5): 66-71,98.

石榮男，1974年生，博士，研究員。研究方向?yàn)殡娮訉?，雷達(dá)與通信系統(tǒng)等。

李瀟女，1993年生，碩士研究生。研究方向?yàn)殡娮訉埂?/p>

劉暢女，1988年生，碩士。研究方向?yàn)殡娮訉埂?/p>

Direction Finding Model for Coherent Multi-signals by Interferometer Based on Vectors Composition

SHI Rong，LI Xiao，LIU Chang

(Science and Technology on Electronic Information Control Laboratory, Chengdu 610036, China)

The conventional direction finding model by interferometer is based on the phase difference measurement among its channels. When it is used for the coherent multi-signals, there are huge errors and even failure sometimes. After the reason is given for this problem, the direction finding novel model for coherent multi-signals by interferometer is established for geometric configuration of its baselines, utilizing the method of vectors composition. Then the resolving method is discussed in detail. Not only the direction of coherent multi-signals can be found, and the number of signals is deduced by the result analysis, but also it is feasible for the single signal. Its compatibility downward is good for the conventional processing. The validity of this model is demonstrated by the simulations in the end. It is important for the development of direction finding theory and anti-jamming application of the interferometer.

interferometer；coherent multi-signals；vectors composition；deducing number of signals；direction finding model；compatibility downwards

10.16592/ j.cnki.1004-7859.2016.09.005

石榮Email：wyx1719@sina.com

2016-04-20

2016-06-22

TN971

A

1004-7859(2016)09-0023-05