石 榮，胡來(lái)招

(電子信息控制國(guó)家重點(diǎn)實(shí)驗(yàn)室，成都 610036)

0 引 言

干涉儀測(cè)向是常用的角度測(cè)量方法，工程應(yīng)用十分廣泛，一般情況下都使用多基線干涉儀，用長(zhǎng)基線獲得高精度的相位差測(cè)量值，用短基線來(lái)解相位差模糊，并且形成了一套完整的多基線干涉儀的無(wú)模糊測(cè)向傳統(tǒng)理論，以此來(lái)確定干涉儀中各條基線的長(zhǎng)度［1～4］。如果在應(yīng)用中僅采用單基線干涉儀，為了避免相位模糊通常要求基線長(zhǎng)度不超過(guò)電磁波的半個(gè)波長(zhǎng)，即d≤λ/2，但此時(shí)測(cè)向精度又無(wú)法保證;如果d?λ/2，雖然能獲得高的測(cè)向精度，但又會(huì)產(chǎn)生多值模糊。顯然在傳統(tǒng)理論指導(dǎo)下單基線干涉儀的兩種性能無(wú)法兼顧。

針對(duì)上述問(wèn)題，分析了傳統(tǒng)理論條件下單基線干涉儀產(chǎn)生測(cè)向模糊的原因，并在作者已經(jīng)研究的干涉儀時(shí)差分析理論基礎(chǔ)上［5］，提出了單基線干涉儀無(wú)模糊測(cè)向新方法，從而使得在單基線長(zhǎng)度d?λ/2 的條件下，巧妙地避免了測(cè)向多值性的發(fā)生，并對(duì)測(cè)向過(guò)程中的整個(gè)信號(hào)處理流程進(jìn)行了闡述，最后以雷達(dá)脈沖信號(hào)和數(shù)字通信相位調(diào)制信號(hào)兩類典型的工程實(shí)際信號(hào)的測(cè)向應(yīng)用為例，通過(guò)仿真驗(yàn)證了單基線無(wú)模糊測(cè)向理論的有效性與實(shí)用性。這一研究結(jié)果將為干涉儀測(cè)向更加廣泛的應(yīng)用提供新的理論指導(dǎo)。

1 單基線干涉儀測(cè)向模糊原因分析

單基線相位干涉儀測(cè)向原理如圖1 所示。圖中有一信號(hào)波長(zhǎng)為λ 的平面電磁波從與天線視軸夾角為θ 的方向達(dá)到測(cè)向天線A 和B。

A 與B 兩天線接收到信號(hào)的相位差為

圖1 單基線相位干涉儀測(cè)向原理圖

如果兩通道的相位響應(yīng)完全一致，則由接收機(jī)輸出信號(hào)的相位差不變，經(jīng)過(guò)鑒相器之后取出相位差信息為

仔細(xì)分析上述單基線干涉儀的傳統(tǒng)測(cè)向理論模型可以發(fā)現(xiàn)，該模型利用了一個(gè)假設(shè)條件，即“來(lái)波信號(hào)為單頻平面電磁波”，這意味著可以利用的唯一信息只有電磁波的載波頻率(對(duì)應(yīng)信號(hào)波長(zhǎng))，除此之外，其他信息均不能利用，于是干涉儀輸出的信號(hào)中只含有與信號(hào)載波相位相關(guān)的信息，而載波相位是以2π 為周期的，因此便形成了單基線測(cè)向求解的多值性。這就是單基線干涉儀測(cè)向模糊的本質(zhì)原因。

2 單基線干涉儀無(wú)模糊測(cè)向理論模型

如前所述，傳統(tǒng)的單基線干涉儀分析模型假設(shè)波長(zhǎng)為λ 的電磁波是一個(gè)單頻信號(hào)，其上沒(méi)有攜帶信息。但是對(duì)于一個(gè)實(shí)際信號(hào)來(lái)說(shuō)，信號(hào)上都承載有信息［6］，即使是單頻雷達(dá)脈沖信號(hào)，脈沖的幅度包絡(luò)也是一個(gè)調(diào)制信息，所以，如果從另外一個(gè)信息提取角度來(lái)分析，利用信號(hào)上調(diào)制的其他信息，就可以從單基線干涉儀的2 個(gè)通道輸出的信號(hào)中提取出更多的信息，從而消除測(cè)向模糊，具體分析如下。

2.1 干涉儀時(shí)差信息與信號(hào)來(lái)波方向的關(guān)系

除了從傳統(tǒng)的相位差處理角度來(lái)分析干涉儀測(cè)向應(yīng)用外，還可以從雙站時(shí)差測(cè)量的角度來(lái)分析單基線干涉儀測(cè)向過(guò)程，這就是干涉儀的時(shí)差分析理論［5］。如圖1 所示，當(dāng)電磁波到達(dá)干涉儀A，B 兩天線時(shí)會(huì)產(chǎn)生一個(gè)時(shí)間差Δt，且有

式中，c 為電磁波傳播速度。如果Δt 已知，顯然可以求得到電磁波的來(lái)波方向θ 為

由式(4)可見(jiàn):只要將干涉儀的兩通道中信號(hào)時(shí)間差測(cè)量值Δt 得到，就可以實(shí)現(xiàn)電磁信號(hào)來(lái)波方向的測(cè)量。在這一理論模型中，有如下兩點(diǎn)值得特別關(guān)注。

(1)理論模型中信號(hào)的來(lái)波方向與測(cè)量的時(shí)間差之間有一一映射關(guān)系。因?yàn)槭?3)中Δt·c≤d是由實(shí)際單基線干涉儀2 個(gè)天線布置位置條件完全決定的，所以測(cè)角過(guò)程不存在模糊解。

(2)新的理論模型沒(méi)有對(duì)信號(hào)的實(shí)際調(diào)制形式做出任何限定。無(wú)論是對(duì)于雷達(dá)信號(hào)，通信信號(hào)，還是其他類型的調(diào)制信號(hào)，均可以應(yīng)用。

2.2 從調(diào)制信號(hào)的頻域相位提取時(shí)差信息

設(shè)承載有信息的時(shí)域調(diào)制信號(hào) f (t)的傅里葉變換為F(ω)，根據(jù)傅里葉變換的位移特性，對(duì)于任何實(shí)因子a，f (t－a)的傅里葉變換為e－jaωF(ω)，其物理意義是:時(shí)間位移對(duì)應(yīng)為頻率調(diào)制。如果將上面的實(shí)因子a 理解為干涉儀兩個(gè)通道接收到的信號(hào)的時(shí)間差，于是就可以在頻域內(nèi)來(lái)求解時(shí)間差a。

設(shè)干擾儀兩個(gè)通道所接收到的信號(hào)分別為:f (t)和f (t－Δt)，將上述兩路信號(hào)從時(shí)域變換到頻域，進(jìn)行共軛相乘和幅度歸一化處理

由式(5)可知，經(jīng)過(guò)信號(hào)處理之后，將在被測(cè)信號(hào)所在頻段(ωd＜ω ＜ωu，其中ωd是信號(hào)的低端頻率，ωu是信號(hào)的高端頻率)的相位函數(shù) φ (ω) 上可得到一條相位直線方程，為

該直線的斜率直接對(duì)應(yīng)了2 個(gè)通道中信號(hào)之間的時(shí)間差Δt。將Δt 代入式(4)，便可以求得到調(diào)制信號(hào)的來(lái)波方向θ。

2.3 新理論對(duì)傳統(tǒng)理論的兼容性說(shuō)明

前面在單基線干涉儀無(wú)模糊測(cè)向理論模型分析過(guò)程中，針對(duì)的都是實(shí)際承載有信息的調(diào)制信號(hào)。對(duì)于沒(méi)有調(diào)制任何信息的單頻電磁波，載波頻率為ω0，新理論模型同樣適用。單頻電磁波是在頻域上表現(xiàn)為一個(gè)沖擊脈沖，信號(hào)帶寬無(wú)限窄，即ωd=ω0=ωu，這樣式(6)就不再表現(xiàn)為一段直線，而退化為1 個(gè)點(diǎn)，于是有

由式(5)所求得的相位函數(shù)φ (ω0)的值域是被限制在［－π，π］范圍內(nèi)，在單頻電磁波條件下同樣存在相位模糊的問(wèn)題。所以對(duì)于單頻電磁波的測(cè)向，新理論模型與傳統(tǒng)的單基線干涉儀相位差分析模型是完全兼容的。但是一旦信號(hào)被調(diào)制有信息之后，就可以利用新理論來(lái)解決測(cè)向模糊問(wèn)題。

綜上所述，在實(shí)際應(yīng)用中信號(hào)都是具有調(diào)制特性的，所以新提出的理論模型正是利用了這一點(diǎn)來(lái)實(shí)現(xiàn)了單基線干涉儀的無(wú)模糊測(cè)向。

3 無(wú)模糊測(cè)向信息處理流程

根據(jù)前面的單基線干涉儀無(wú)模糊測(cè)向理論模型，按照如下的處理流程來(lái)完成整個(gè)測(cè)向過(guò)程。

(1)利用信號(hào)檢測(cè)結(jié)果，確定實(shí)施測(cè)向的電磁目標(biāo)對(duì)象，同時(shí)采集單基線干涉儀A 與B 兩天線接收到的信號(hào)樣本;

(2)對(duì)采集后的時(shí)域信號(hào)樣本進(jìn)行補(bǔ)零后，變換到頻域，以利于后續(xù)的求解過(guò)程;

(3)利用式(5)求得到頻域相位函數(shù) φ (ω)，此時(shí)φ (ω)的值域被限制在［－π，π］范圍內(nèi)，雖然φ (ω)是一條直線方程，但是其縱坐標(biāo)也會(huì)被強(qiáng)制折疊在［－π，π］范圍內(nèi)，所以需要對(duì) φ (ω)進(jìn)行相位解擾處理，這一步驟可以直接 用MATLAB 的unwrap 函數(shù)來(lái)實(shí)現(xiàn)。

(4)根據(jù)信號(hào)所占帶寬:ωd＜ω ＜ωu，在該帶寬范圍內(nèi)對(duì)相位函數(shù) φ (ω)進(jìn)行直線擬合，所擬合直線的斜率即對(duì)應(yīng)時(shí)差數(shù)值Δt;

(5)將Δt 代入式(4)，即可求得到調(diào)制信號(hào)的來(lái)波方向θ。

關(guān)于上述單基線干涉儀無(wú)模擬測(cè)向信息處理流程的細(xì)節(jié)將在后續(xù)仿真驗(yàn)證中進(jìn)一步說(shuō)明。

4 仿真驗(yàn)證

本節(jié)將以兩類典型的實(shí)際信號(hào):雷達(dá)脈沖信號(hào)和數(shù)字通信相位調(diào)制信號(hào)的測(cè)向應(yīng)用為例，來(lái)對(duì)單基線無(wú)模糊測(cè)向理論進(jìn)行驗(yàn)證。

4.1 對(duì)單頻雷達(dá)脈沖信號(hào)的無(wú)模糊測(cè)向仿真

仿真條件:S 波段單頻雷達(dá)脈沖信號(hào)的載波頻率為3 GHz，對(duì)應(yīng)的信號(hào)波長(zhǎng)為0.1 m，矩形脈沖寬度為0.9 μs，單基線干涉儀的長(zhǎng)度為2 m，顯然這一長(zhǎng)度遠(yuǎn)遠(yuǎn)超過(guò)了信號(hào)的半波長(zhǎng)，信號(hào)來(lái)波方向θ=30°，SNR≥20 dB。采樣頻率取為10 GHz，采樣時(shí)間長(zhǎng)度取1 μs，則樣本點(diǎn)數(shù)為1e4，干涉儀A、B 兩天線接收到的信號(hào)時(shí)域波形，如圖2 所示。

圖2 單基線干涉儀兩天線接收到的信號(hào)時(shí)域波形

按照傳統(tǒng)的干涉儀測(cè)量理論，A 與B 兩天線接收到信號(hào)的理論相位差φ =0°，根據(jù)基線長(zhǎng)度與波長(zhǎng)關(guān)系，一共有39 個(gè)來(lái)波方向都會(huì)產(chǎn)生φ =0°的結(jié)果，所以傳統(tǒng)理論將無(wú)法獲得唯一的測(cè)向解。

按照本文提出的理論模型，對(duì)采集到的1 μs 長(zhǎng)信號(hào)進(jìn)行補(bǔ)零，將樣本點(diǎn)數(shù)擴(kuò)展到1e5，并按照式(5)求解頻域相位函數(shù) φ (ω) ，根據(jù)仿真條件，被測(cè)信號(hào)所在頻段的高/低端頻率分別為:ωu=3000.5 MHz，ωd=2999.5 MHz。在此頻段區(qū)域內(nèi)對(duì)應(yīng)的信號(hào)的頻域幅度曲線與相位函數(shù) φ (ω) 的曲線如圖3 所示(圖中實(shí)線)。為了進(jìn)行對(duì)比，在相位函數(shù)中，還將臨近的在傳統(tǒng)干涉儀測(cè)向理論中會(huì)產(chǎn)生模糊的2 個(gè)來(lái)波方向，分別是26.74°與33.37°，對(duì)應(yīng)的相位曲線圖中虛線所示。

圖3 信號(hào)的局部頻域幅度曲線與相位函數(shù) φ(ω) 的曲線

如圖3 所示，相位函數(shù)曲線在對(duì)應(yīng)區(qū)域表現(xiàn)為一條直線，按照式(6)進(jìn)行直線擬合，可求得直線的斜率Δt=3.294 ns，然后利用式(4)可求得信號(hào)的來(lái)波方向θ=29.61°。在傳統(tǒng)干涉儀測(cè)向理論中會(huì)產(chǎn)生模糊的2 個(gè)來(lái)波方向的相位曲線函數(shù)所擬合的2條直線的斜率分別為:2.892 ns 與3.503 ns，這一數(shù)值是可以做到多條不同斜率直線的區(qū)分。由此可見(jiàn)，按照新的理論模型可以在d?λ/2 時(shí)實(shí)現(xiàn)單基線干涉儀對(duì)雷達(dá)脈沖信號(hào)的無(wú)模糊測(cè)向。

4.2 對(duì)BPSK 通信信號(hào)的無(wú)模糊測(cè)向仿真

仿真條件:VHF 波段通信信號(hào)的載波頻率為250 MHz，對(duì)應(yīng)的信號(hào)波長(zhǎng)為1.2 m，BPSK 調(diào)制的符號(hào)速率為20 MHz，單基線干涉儀的長(zhǎng)度為20 m，顯然這一長(zhǎng)度遠(yuǎn)遠(yuǎn)超過(guò)了半波長(zhǎng)，信號(hào)來(lái)波方向θ =13°，SNR≥20 dB。帶通采樣頻率取為200 MHz，采樣時(shí)間長(zhǎng)度取250 μs，則樣本點(diǎn)數(shù)為5e4，干涉儀A、B 兩天線接收到的信號(hào)時(shí)域波形如圖4 所示。

圖4 單基線干涉儀兩天線接收到的信號(hào)時(shí)域波形

按照傳統(tǒng)的干涉儀測(cè)量理論，A 與B 兩天線接收到信號(hào)的理論相位差φ= －90.3°，根據(jù)基線長(zhǎng)度與波長(zhǎng)關(guān)系，一共有33 個(gè)不同的來(lái)波方向都會(huì)產(chǎn)生同樣的結(jié)果，所以按照傳統(tǒng)理論將無(wú)法獲得唯一的測(cè)向解。

按照本文提出的理論模型，對(duì)采集到的250 μs長(zhǎng)信號(hào)進(jìn)行補(bǔ)零，將樣本點(diǎn)數(shù)擴(kuò)展到1e5，并按照式(5)求解頻域相位函數(shù) φ (ω) ，根據(jù)前面的信號(hào)載波頻率，被測(cè)信號(hào)所在頻段的高低端頻率(按帶通采樣后的值來(lái)計(jì)算)分別為:ωu= 60 MHz，ωd=40 MHz。整個(gè)信號(hào)所在頻段的頻域幅度曲線與局部相位函數(shù)曲線如圖5 所示。

圖5 信號(hào)的頻域幅度曲線與局部相位函數(shù) φ(ω) 的曲線

如圖5 所示，相位函數(shù)曲線在對(duì)應(yīng)區(qū)域表現(xiàn)為一條直線，按照式(6)進(jìn)行直線擬合，可求得直線的斜率Δt=15.024 ns，然后利用式(4)可求得信號(hào)的來(lái)波方向θ =13.02°。由此可見(jiàn)，按照新的理論模型，可以在d?λ/2 時(shí)實(shí)現(xiàn)單基線干涉儀對(duì)通信信號(hào)的無(wú)模糊測(cè)向。

5 結(jié) 語(yǔ)

本文主要針對(duì)單基線干涉儀無(wú)模糊測(cè)向問(wèn)題進(jìn)行了理論研究與仿真驗(yàn)證，分析了傳統(tǒng)的單基線干涉儀相位差測(cè)量理論產(chǎn)生測(cè)向多值性模糊的本質(zhì)原因，利用干涉儀時(shí)差分析理論，揭示了干涉儀時(shí)差信息與信號(hào)來(lái)波方向之間的對(duì)應(yīng)關(guān)系，并從調(diào)制信號(hào)的頻域相位來(lái)提取時(shí)差信息，建立了單基線干涉儀無(wú)模糊測(cè)向的理論模型，分析了新理論對(duì)傳統(tǒng)理論的兼容性，闡述了單基線干涉儀無(wú)模糊測(cè)向的完整信息處理流程。并通過(guò)對(duì)兩類典型信號(hào)的仿真，驗(yàn)證了新理論的有效性與實(shí)用性。這一研究結(jié)果將為干涉儀測(cè)向更加廣泛的應(yīng)用提供新的理論指導(dǎo)。下一步，還會(huì)根據(jù)所提出的新理論，進(jìn)一步開(kāi)展單基線干涉儀無(wú)模糊測(cè)向的應(yīng)用性研究工作。

［1］DAVID L ADAMY. EW101:A First Course in Electronic Warfare［M］.Beijing:Publishing House of Electronic Industry，2009.

［2］ DAVID L ADAMY. EW102:A Second Course in Electronic Warfare［M］. Beijing:Publishing House of Electronic Industry，2009.

［3］DAVID L ADAMY. EW103:Tactical Battlefield Communication Electronic Warfare［M］.Beijing:Publishing House of Electronic Industry，2010.

［4］栗平，趙國(guó)慶. 信息對(duì)抗技術(shù)［M］. 北京:清華大學(xué)出版社，2008.

［5］石榮，胡來(lái)招. 干涉儀的時(shí)差分析理論及應(yīng)用［J］. 航天電子對(duì)抗，2012 年1 月錄用。

［6］胡來(lái)招.信號(hào)與信息［M］.北京:電子工業(yè)出版社，2010.