基于二次相差法的多頻連續(xù)波測(cè)距雷達(dá)仿真與研究

柏業(yè)鑫

(西安電子科技大學(xué)電子工程學(xué)院，陜西 西安 710071)

多頻連續(xù)波雷達(dá)具有較高的測(cè)距精度，還可提供速度信息，并具有結(jié)構(gòu)簡(jiǎn)單、造價(jià)低、便于應(yīng)用FFT等現(xiàn)代數(shù)字信號(hào)處理方法等優(yōu)點(diǎn)，主要應(yīng)用于靶場(chǎng)測(cè)量，因其優(yōu)良性能逐漸應(yīng)用于各種場(chǎng)合。

多頻連續(xù)波是基于雙頻比相測(cè)距［1］，是對(duì)雙頻比相測(cè)距的改進(jìn)和提高。雙頻連續(xù)波雷達(dá)測(cè)距的理論精度［2］為 σR=c/(4πΔf)，其中，Δf為雙頻測(cè)距時(shí)所用的頻差;SNR為FFT輸出信噪比;c為光速。同時(shí)在該雙頻頻差下的最大不模糊距離［3］為 Runamb=c/(2Δf)。可見雙頻連續(xù)波雷達(dá)存在著測(cè)距精度和最大不模糊距離不可調(diào)和的矛盾。而為解決這一矛盾，提出了多頻連續(xù)波測(cè)距法，主要有二次相差法和參差多頻法。傳統(tǒng)多頻連續(xù)波雷達(dá)基于雙頻比相測(cè)距，都是先由雙頻比相分離目標(biāo)，并解出相位差。

二次相差法相對(duì)于參差多頻法來(lái)說(shuō)，距離的解模糊算法較為簡(jiǎn)單，運(yùn)算量相對(duì)較小，易于工程實(shí)現(xiàn)。且只要頻點(diǎn)選擇適當(dāng)，目標(biāo)就一直處在最大不模糊距離之內(nèi)。

1 二次相差法基本原理

首先發(fā)射雙頻信號(hào)f0和fn+1，然后發(fā)射下一組雙頻信號(hào)f0和fn，…，最后發(fā)射雙頻信號(hào)f0和f1。多頻信號(hào) fn+1、fn、…、f1的值依次增大。

由頻率f0和fn+1可得到他們的相位差Δφ0，n+1為

其中，Δfn+1=fn+1－f0，同理可以通過(guò)頻率f0和fn得到它們的相位差 Δφ0，n為

由上面兩式可以得到頻率fn和頻率fn+1的二次相差Δ2φn，n+1為

若令補(bǔ)償相位Δφ補(bǔ)償=2πΔfn+1代入式(3)得

則對(duì)應(yīng)的最大不模糊距離為

其中，vn為目標(biāo)在頻率fn所對(duì)應(yīng)測(cè)量周期的徑向速度，它可以通過(guò)FFT后譜峰搜索的方法得到，T表示相鄰兩次測(cè)量之間的時(shí)間間隔。由此，得到對(duì)應(yīng)于頻率fn的距離Rn。在上式的推導(dǎo)過(guò)程中，由于頻率fn、fn+1發(fā)射時(shí)間間隔較短，可以認(rèn)為目標(biāo)是勻速運(yùn)動(dòng)的，Δφ補(bǔ)償就是由于目標(biāo)的運(yùn)動(dòng)所引起的。

頻點(diǎn)選擇時(shí)，二次頻差按二次頻差系數(shù)N遞增，如式(6)所示

式中，N為二次頻差系數(shù)。最小二次頻差Δ2fn，n+1=Δfn－Δfn+1的值需要根據(jù)本次測(cè)量所要求的測(cè)距范圍來(lái)確定，如Runamb=30 km，則由上式可以得到Δ2fn，n+1=Δfn－Δfn+1=c/2Runamb=5 kHz。最大的雙頻頻差Δf1由本次測(cè)量的精度確定。

2 仿真實(shí)驗(yàn)

取基準(zhǔn)頻率f0=1.5 GHz，發(fā)射7組點(diǎn)頻信號(hào)，N=5，因此可知

因此，f0=1.5 GHz、f1=1.515 GHz、f2=1.512 GHz、f3=1.5114 GHz、f4=1.51128 GHz、f5=1.511256 GHz、f6=1.5112512 GHz。最小的二次頻差為Δf4－Δf5=4800 Hz，其對(duì)應(yīng)的最大不模糊距離為 Runamb=31.25 km。最大的雙頻頻差為f1－f0=15 MHz，其對(duì)應(yīng)的最大不模糊距離為10 m。

首先通過(guò)仿真來(lái)分析二次相差法多頻連續(xù)波測(cè)距體制的測(cè)距均方根誤差與目標(biāo)距離之間的關(guān)系。這里假設(shè)目標(biāo)徑向速度vr=100 m/s不變，但目標(biāo)距離是隨著每次仿真實(shí)驗(yàn)而變化的，從0.1～30 km，每次遞增100 m。對(duì)于每個(gè)目標(biāo)距離的仿真實(shí)驗(yàn)，均進(jìn)行1000次蒙特卡羅仿真實(shí)驗(yàn)，圖1給出了二次相差法多頻連續(xù)波測(cè)距體制下，1000次蒙特卡羅仿真實(shí)驗(yàn)統(tǒng)計(jì)的測(cè)距均方根誤差與距離曲線。但這里忽略了信噪比隨距離增大而減小的影響，僅假設(shè)信噪比SNR始終不變。可以看出，測(cè)量均方根誤差較小，這是因?yàn)樽詈笠粋€(gè)最大不模糊距離為10 m，根據(jù)雙頻比相測(cè)距法的理論公式可以得出測(cè)量誤差會(huì)非常小。仿真信噪比為30 dB。二次相差法多頻連續(xù)波雷達(dá)可以發(fā)射很少的點(diǎn)頻連續(xù)波信號(hào)，以達(dá)到最大不模糊距離和高測(cè)距精度。

其次通過(guò)仿真分析了二次相差法多頻連續(xù)波測(cè)距體制的測(cè)距均方根誤差與FFT輸出信噪比之間的關(guān)系。這里依然假設(shè)目標(biāo)徑向速度vr=100 m/s不變，目標(biāo)距離為30 km，但SNR從－5～35 dB依次增大，每次增加0.1 dB。對(duì)于每個(gè)SNR與測(cè)距均方根誤差的仿真實(shí)驗(yàn)，均進(jìn)行1000次蒙特卡羅仿真實(shí)驗(yàn)。圖2給出了二次相差法多頻連續(xù)波測(cè)距體制下，1000次蒙特卡羅仿真實(shí)驗(yàn)統(tǒng)計(jì)的測(cè)距均方根誤差與輸出信噪比曲線。

由圖2可以看出，信噪比對(duì)二次相差法測(cè)距的影響較大，測(cè)距均方根誤差隨FFT輸出信噪比增大而迅速降低，當(dāng)信噪比＜0 dB時(shí)，測(cè)量誤差急劇增大。但二次相差法對(duì)信噪比的要求不高，10 dB以上的信噪比都能達(dá)到測(cè)距精度要求。

3 二次相差法多頻連續(xù)波雷達(dá)多目標(biāo)測(cè)距

由于二次相差法多頻連續(xù)波測(cè)距雷達(dá)是基于雙頻比相測(cè)距體制的，下面研究下雙頻比相測(cè)距體制下，多目標(biāo)分離的問(wèn)題。

對(duì)于發(fā)射信號(hào)頻率為f0和f1的回波差拍采樣信號(hào)分別為

其中，ωdA0=2πfdA0Ts、ωdA1=2πfdA1Ts、ωdB0=2πfdB0Ts、ωdB1=2πfdB1Ts。EA、EB分別為 A、B 目標(biāo)的回波幅度，與目標(biāo)距離和截面積等有關(guān)系。

對(duì)上式做FFT變換得

(1)當(dāng)兩個(gè)目標(biāo)A、B具有不同的速度時(shí)，即ωdA0=ωdB0且 ωdA1=ωdB1。

此時(shí)X0(k)和X1(k)的譜峰都會(huì)有兩個(gè)極大值kA0、kA1和 kB0、kB1，分別對(duì)應(yīng)為 A、B 目標(biāo)在 f0和 f1上的多普勒頻率。此時(shí)A、B兩個(gè)目標(biāo)就可以通過(guò)FFT譜峰搜索分離，根據(jù)f0和f1的FFT譜峰處的相位值就可以分別求出A、B的距離引起的相位差信息，然后就可以通過(guò)二次相差法解模糊求出A、B兩個(gè)目標(biāo)的精確距離值。

(2)當(dāng)兩個(gè)目標(biāo)A、B具有相同的速度時(shí)，即ωdA0≠ωdB0且 ωdA1≠ωdB1。

此時(shí)X0(k)和X1(k)都會(huì)只有一個(gè)FFT譜峰k0、k1，A、B兩個(gè)目標(biāo)就無(wú)法通過(guò)FFT后譜峰搜索分離，也就是把A、B兩個(gè)目標(biāo)看成了一個(gè)目標(biāo)。此時(shí)由于EA、EB分別為A、B目標(biāo)的回波幅度，與目標(biāo)距離和截面積等有關(guān)，而兩個(gè)目標(biāo)A、B的距離和截面積等都不相等，則EA、EB，求解出的FFT譜峰相位值是兩個(gè)目標(biāo)距離信息引起的相位項(xiàng)的一種非線性和，從而通過(guò)相位差求解的距離是RA和RB的某種非線性和。

圖3和圖4說(shuō)明，不同速度的兩個(gè)目標(biāo)的FFT譜峰是分開的，而同速度的兩個(gè)目標(biāo)的FFT譜峰是重合的。其實(shí)，當(dāng)多個(gè)速度相當(dāng)接近的目標(biāo)同在波束內(nèi)時(shí)，基于雙頻比相的多頻連續(xù)波雷達(dá)是無(wú)法分離的。

總之，雙頻比相測(cè)距體制只能分離不同速度的目標(biāo)，無(wú)法分離同一速度但距離不同的目標(biāo)。也就是說(shuō)基于雙頻比相測(cè)距的二次相差法多篇連續(xù)波雷達(dá)，只能求解速度不同的目標(biāo)的距離值，無(wú)法求解同速度不同距離目標(biāo)的距離值。

4 結(jié)束語(yǔ)

二次相差法的多頻連續(xù)波測(cè)距體制的算法簡(jiǎn)單，實(shí)現(xiàn)容易，對(duì)信噪比要求較低。若參數(shù)選用得當(dāng)，二次相差法完全能滿足雷達(dá)測(cè)距戰(zhàn)術(shù)指標(biāo)要求，具有廣泛的應(yīng)用前景。但是傳統(tǒng)的多頻連續(xù)波雷達(dá)都是基于雙頻比相測(cè)距體制的，而雙頻比相測(cè)距有著明顯的缺陷，就是無(wú)法分離速度十分接近的目標(biāo)，需要其他的方法進(jìn)行分離目標(biāo)，比如加入線性調(diào)頻信號(hào)進(jìn)行運(yùn)動(dòng)目標(biāo)剝離［4］，或者多頻連續(xù)波測(cè)距體制棄用雙頻比相法，改用其他點(diǎn)頻連續(xù)波測(cè)距方案。

［1］BOYER W D.A diplex doppler phase comparison radar［J］.IEEE Trans Actions on Aerospace and Navigational Electronics，1963，ANE －10(1):27 －33.

［2］SKOLNIK M I.雷達(dá)系統(tǒng)導(dǎo)論［M］.3版.左群聲，徐國(guó)良，馬林，等，譯.北京:國(guó)防工業(yè)出版社，1992.

［3］袁俊泉，龔享銥，皇甫堪.基于二次差頻的多頻連續(xù)波測(cè)距方法研究［J］.電子學(xué)報(bào)，2004，32(12):2056 －2058.

［4］戴奉周，馮維婷，沈福民.數(shù)字多頻連續(xù)波雷達(dá)信號(hào)處理中的關(guān)鍵算法［J］.雷達(dá)科學(xué)與技術(shù)，2005，3(3):177－184.