米波雷達(dá)低仰角目標(biāo)多徑模型及其反演方法研究

鄭軼松 陳伯孝

?

米波雷達(dá)低仰角目標(biāo)多徑模型及其反演方法研究

鄭軼松*陳伯孝

(西安電子科技大學(xué)雷達(dá)信號(hào)處理國(guó)家重點(diǎn)實(shí)驗(yàn)室 西安 710017)

現(xiàn)有的低仰角測(cè)高方法多采用鏡面反射模型，將直達(dá)波與多徑簡(jiǎn)化為兩遠(yuǎn)場(chǎng)點(diǎn)源；然而實(shí)際中不規(guī)則反射面使多徑回波波前畸變，遠(yuǎn)場(chǎng)點(diǎn)源模型難以完全描述多徑信號(hào)。針對(duì)此模型失配問(wèn)題，該文重點(diǎn)研究低仰角目標(biāo)多徑模型，首先分析經(jīng)典多徑模型，對(duì)反射系數(shù)和反射面高度進(jìn)行參數(shù)反演；然后提出一種擾動(dòng)多徑模型，將反射面對(duì)多徑回波的影響建模為擾動(dòng)反射系數(shù)，并利用最大似然算法反演擾動(dòng)反射系數(shù)。計(jì)算機(jī)仿真結(jié)果驗(yàn)證了參數(shù)反演方法的有效性；實(shí)測(cè)數(shù)據(jù)驗(yàn)證了復(fù)雜情形下所建模型的合理性和反演方法的有效性，提高了低仰角測(cè)高算法在實(shí)際陣地的適用性。

米波雷達(dá)；多徑模型；低仰角測(cè)高；擾動(dòng)反射系數(shù)

1 引言

低仰角測(cè)高問(wèn)題一直是米波雷達(dá)技術(shù)亟待解決的重要問(wèn)題之一。由于米波雷達(dá)波長(zhǎng)長(zhǎng)，主瓣寬，往往很難從空域、時(shí)域或頻域分辨直達(dá)信號(hào)和多徑信號(hào)。低仰角測(cè)高實(shí)質(zhì)是在低信噪比情況下，分辨和測(cè)量?jī)蓚€(gè)夾角很小的相干源來(lái)波方向的問(wèn)題。針對(duì)該問(wèn)題人們已經(jīng)進(jìn)行了大量的理論與實(shí)驗(yàn)研究。近年來(lái)隨著壓縮感知、稀疏貝葉斯等技術(shù)的興起，稀疏恢復(fù)類(lèi)算法也被用來(lái)提高低信噪比下的角度估計(jì)性能。但現(xiàn)有算法多采用鏡面反射模型，將直達(dá)波與反射波簡(jiǎn)化為兩點(diǎn)源模型，研究多集中于如何解相干和提高測(cè)角精度。然而實(shí)際中的復(fù)雜陣地(如崎嶇的山地、高海情的海面)并不滿足鏡面反射模型，復(fù)雜反射面會(huì)造成多徑信號(hào)的波前畸變，導(dǎo)致各陣元反射點(diǎn)不再處于同一反射面，或各陣元接收多徑信號(hào)的反射系數(shù)不再相同，即多徑信號(hào)導(dǎo)向矢量不再滿足遠(yuǎn)場(chǎng)點(diǎn)源模型。由于多徑信號(hào)導(dǎo)向矢量無(wú)法精確已知，即使直達(dá)波信號(hào)滿足遠(yuǎn)場(chǎng)點(diǎn)源模型，也會(huì)被淹沒(méi)在相干性很強(qiáng)的多徑信號(hào)中而難以分辨和測(cè)量。這也是低仰角測(cè)高問(wèn)題一直難以解決的原因之一，即多徑信號(hào)的模型失配問(wèn)題。眾多學(xué)者針對(duì)此問(wèn)題對(duì)低仰角測(cè)高算法進(jìn)行了一系列改進(jìn)。文獻(xiàn)[2,7]利用雷達(dá)高度、目標(biāo)斜距、地球曲率等先驗(yàn)信息，提出基于精確反射多徑信號(hào)模型的最大似然算法，但沒(méi)有考慮反射面高度對(duì)算法性能的影響，只能適用于非常平坦的陣地。文獻(xiàn)[11]闡述了反射面高度對(duì)算法的影響，提出基于地形匹配的合成導(dǎo)向矢量算法，但是該文沒(méi)有提出反演反射面高度的方法，也沒(méi)有考慮復(fù)雜反射面對(duì)多徑信號(hào)的影響。文獻(xiàn)[12]考慮漫反射對(duì)測(cè)高算法的影響，將漫反射占主導(dǎo)地位時(shí)的多徑信號(hào)建模為分布源，但是其分布模型難以確定，一般假設(shè)分布源服從高斯分布，因此其適用性不高。因此本文從多徑模型及反演方法入手，針對(duì)經(jīng)典多徑模型中反射面高度、反射系數(shù)提出基于最大似然算法的參數(shù)反演方法；并利用擾動(dòng)反射系數(shù)建模復(fù)雜反射面對(duì)多徑信號(hào)的擾動(dòng)，提出擾動(dòng)多徑模型；最后通過(guò)仿真數(shù)據(jù)分析和實(shí)測(cè)數(shù)據(jù)驗(yàn)證所提模型和反演方法的有效性。

2 經(jīng)典多徑模型及其反演

經(jīng)典多徑模型示意圖如圖1所示，假設(shè)接收陣列為線性均勻線陣，包含個(gè)陣元，快拍數(shù)為。陣列接收信號(hào)為

圖1 經(jīng)典多徑模型示意圖

多徑模型一般可分為不考慮地球曲率的平面反射模型和考慮地球曲率的曲面反射模型。這兩種模型之間具有轉(zhuǎn)化關(guān)系[11]，因此為了計(jì)算簡(jiǎn)便，本文采用平面反射模型。對(duì)于遠(yuǎn)場(chǎng)目標(biāo)，可將直達(dá)波導(dǎo)向矢量表示為

如圖1所示，隨著目標(biāo)的運(yùn)動(dòng)，不同點(diǎn)跡對(duì)應(yīng)不同的反射點(diǎn)，而實(shí)際陣地并不平坦，因此不同點(diǎn)跡其反射面高度也不相同，設(shè)反射面關(guān)于參考平面的高度為，文中以雷達(dá)架設(shè)陣地所在平面為參考平面，如圖1中虛線所示，則直達(dá)波和反射波波程分別為

其中，0為目標(biāo)與雷達(dá)的水平距離，h為目標(biāo)高度，為各陣元關(guān)于參考平面的高度，，其中，為最下邊陣元的高度。由文獻(xiàn)[13]可知，當(dāng)h,時(shí)，式(4)和式(5)可近似表示為

此時(shí)波程差可表示為

由式(8)可知，若h,和0都精確已知，則地面反射面高度h決定波程差。

將陣列接收信號(hào)重新表示為式(9)的形式：

陣列接收信號(hào)協(xié)方差矩陣的最大似然估計(jì)為

對(duì)于合成導(dǎo)向矢量最大似然算法[2]，首先構(gòu)造合成導(dǎo)向矢量矩陣及其投影矩陣：

通過(guò)DOA估計(jì)得到目標(biāo)仰角和斜距后，對(duì)于距離較近的目標(biāo)，目標(biāo)高度可由式(14)計(jì)算得到：

對(duì)于距離較遠(yuǎn)的目標(biāo)，必須考慮地球曲率和空氣折射的影響，修正后得到式(15)：

合成導(dǎo)向矢量最大似然算法實(shí)質(zhì)是一種模型匹配算法，對(duì)似然函數(shù)影響較大的參數(shù)包括直達(dá)波仰角，反射系數(shù)和反射面高度。現(xiàn)有算法多假設(shè)所有點(diǎn)跡的反射系數(shù)和反射面高度不變，而在實(shí)際陣地，由于不同點(diǎn)跡對(duì)應(yīng)不同的反射區(qū)，地面反射系數(shù)和反射面海拔高度時(shí)變而未知。因此需要根據(jù)二次雷達(dá)數(shù)據(jù)，在已知目標(biāo)仰角時(shí)，反演反射系數(shù)和反射面高度。當(dāng)直達(dá)波仰角確定時(shí)，只需構(gòu)造似然函數(shù)，對(duì)反射系數(shù)和反射面高度進(jìn)行2維聯(lián)合估計(jì)即可。反射系數(shù)的搜索范圍一般可取，而反射面高度的搜索范圍需由以下準(zhǔn)則確定：

進(jìn)而得到反射系數(shù)和反射面高度的最大似然估計(jì)：

3 擾動(dòng)多徑模型及其反演

當(dāng)雷達(dá)陣地周?chē)匦螐?fù)雜時(shí)，由于地面有效反射區(qū)(第1菲涅爾區(qū))內(nèi)反射面高低不一、反射介質(zhì)各不相同，各陣元接收的多徑信號(hào)的幅度不再一致，相位也不滿足線性變化，此時(shí)，反射系數(shù)、目標(biāo)仰角和反射面高度等參數(shù)都不足以完全描述多徑回波。另外對(duì)于同一點(diǎn)跡，若反射介質(zhì)均勻而各陣元地面反射區(qū)高度不同，其反射系數(shù)可以認(rèn)為近似相同，但是反射區(qū)高度不一會(huì)造成多徑信號(hào)的相位擾動(dòng)，此時(shí)傳統(tǒng)遠(yuǎn)場(chǎng)點(diǎn)源模型也難以匹配真實(shí)多徑模型。這也是傳統(tǒng)測(cè)高算法雖然具有較好角度分辨率且能解相干信號(hào)，卻難以在實(shí)測(cè)數(shù)據(jù)處理中取得較高測(cè)角測(cè)高精度的原因。因此本文針對(duì)此多徑信號(hào)導(dǎo)向矢量的模型失配問(wèn)題，提出擾動(dòng)多徑模型，如圖2所示，不同陣元接收的多徑信號(hào)受復(fù)雜反射面擾動(dòng)影響，反射系數(shù)各不相同，定義第個(gè)陣元接收信號(hào)的擾動(dòng)反射系數(shù)為,為復(fù)數(shù)，其幅度的均值為經(jīng)典多徑模型中的反射系數(shù)。擾動(dòng)多徑模型下陣列接收信號(hào)為

圖2 擾動(dòng)多徑模型示意圖

通過(guò)式(20)得到反射系數(shù)和反射面高度的最大似然估計(jì)后，為了得到更精確的擾動(dòng)多徑模型，需要進(jìn)一步對(duì)擾動(dòng)反射系數(shù)進(jìn)行估計(jì)。由式(21)可知，式中未知參數(shù)只有信號(hào)，噪聲和擾動(dòng)反射系數(shù)。若已知信號(hào)，可通過(guò)最大似然估計(jì)得到擾動(dòng)反射系數(shù)。由文獻(xiàn)[14]可知，信號(hào)的最大似然估計(jì)為

下面給出擾動(dòng)多徑模型參數(shù)反演的具體步驟：

(1)利用二次雷達(dá)關(guān)聯(lián)一次雷達(dá)數(shù)據(jù)，作為當(dāng)前點(diǎn)跡目標(biāo)仰角的真值，進(jìn)而得到目標(biāo)與雷達(dá)的水平距離和目標(biāo)高度，聯(lián)合各陣元關(guān)于參考平面的高度等先驗(yàn)信息通過(guò)式(19)計(jì)算反射面高度搜索范圍；同時(shí)確定反射系數(shù)搜索范圍，如；

(3)根據(jù)式(10)、式(11)和式(12)分別計(jì)算接收信號(hào)的協(xié)方差矩陣，合成導(dǎo)向矢量矩陣及其投影矩陣；

由以上步驟可知，經(jīng)典多徑模型反演是擾動(dòng)多徑模型反演的前提，只有通過(guò)經(jīng)典多徑模型反演得到正確的反射面高度和反射系數(shù)后，才能進(jìn)一步得到擾動(dòng)反射系數(shù)的估計(jì)。在已知擾動(dòng)反射系數(shù)和反射面高度后，可構(gòu)造復(fù)合導(dǎo)向矢量矩陣，結(jié)合最大似然、MUSIC等算法構(gòu)造空間譜進(jìn)行DOA估計(jì)。此時(shí)信號(hào)模型與真實(shí)模型匹配，因此能得到更高的測(cè)角測(cè)高精度。

4 仿真數(shù)據(jù)分析

仿真1 驗(yàn)證本文參數(shù)反演方法的正確性 需要說(shuō)明的是，由于經(jīng)典多徑模型的參數(shù)反演方法包含在擾動(dòng)多徑模型反演的步驟中，因此在仿真1中只針對(duì)擾動(dòng)多徑模型進(jìn)行參數(shù)反演，若擾動(dòng)多徑模型反演正確，也就說(shuō)明了經(jīng)典多徑模型反演方法的正確性。仿真參數(shù)如下：假設(shè)高斯白噪聲背景下，陣元的均勻線陣，波長(zhǎng)，陣元間隔半波長(zhǎng)，陣列斜向上傾斜，信噪比，快拍數(shù)為，陣列最下邊陣元高度為，反射面高度為，直達(dá)波入射角度為，擾動(dòng)反射系數(shù)幅度服從(0.6, 1.0)的均勻分布，相位服從的均勻分布。由式(19)計(jì)算得反射面高度搜索范圍為m，取反射系數(shù)搜索范圍為[－1.0,0]，仿真計(jì)算得到反射系數(shù)為，反射面高度為，與仿真所設(shè)參數(shù)相符，說(shuō)明該方法能精確反射系數(shù)和反射面高度。為了觀察似然函數(shù)關(guān)于反射面高度的周期性，取反射面高度搜索范圍為，繪制似然函數(shù)關(guān)于反射系數(shù)和反射面高度的聯(lián)合2維估計(jì)圖。如圖3(a)所示，在的范圍中似然函數(shù)均勻分布4個(gè)譜峰，可見(jiàn)似然函數(shù)是關(guān)于反射面高度的周期函數(shù)。取反射系數(shù)為，繪制似然函數(shù)關(guān)于反射面高度的剖面圖，如圖3(b)所示，圖中虛線為各譜峰所在位置，可見(jiàn)似然函數(shù)譜峰分別位于和處，周期為19.1 m，與理論分析相符。在反演得到反射面高度和反射系數(shù)估計(jì)結(jié)果后，代入擾動(dòng)多徑模型中反演擾動(dòng)反射系數(shù)，估計(jì)結(jié)果如圖4所示，可見(jiàn)擾動(dòng)反射系數(shù)估計(jì)值與真實(shí)值非常接近，說(shuō)明本文算法可有效反演擾動(dòng)反射系數(shù)。綜合仿真1的結(jié)果可知本文反演方法可有效反演擾動(dòng)多徑模型。

仿真2 考察擾動(dòng)大小及信噪比與擾動(dòng)反射系數(shù)估計(jì)均方根誤差的關(guān)系 信噪比由~變化，考察5組擾動(dòng)反射系數(shù)，其擾動(dòng)大小滿足表1中參數(shù)分布，表1中參數(shù)均為均勻分布，擾動(dòng)大

圖3 似然函數(shù)周期性仿真圖

表1 擾動(dòng)反射系數(shù)分布

小由序號(hào)1至序號(hào)5擾動(dòng)依次增大，蒙特卡羅次數(shù)=200次，其余仿真參數(shù)與仿真1一致。定義擾動(dòng)反射系數(shù)的均方根誤差RMSE(Root Mean SquareError)為，上標(biāo)表示蒙特卡羅次數(shù)，表示范數(shù)。RMSE統(tǒng)計(jì)結(jié)果如圖5所示，可見(jiàn)低信噪比時(shí)各組測(cè)試數(shù)據(jù)都無(wú)法有效估計(jì)擾動(dòng)反射系數(shù)，即便對(duì)于無(wú)擾動(dòng)的數(shù)據(jù)組，由于噪聲對(duì)信號(hào)影響較大，導(dǎo)致擾動(dòng)反射系數(shù)估計(jì)結(jié)果精度很差，因此本算法在低信噪比(小于0 dB)時(shí)無(wú)法適用；信噪比較高(大于4 dB)時(shí)，各數(shù)據(jù)組RMSE均小于0.2，估計(jì)結(jié)果較好，且序號(hào)較小的數(shù)據(jù)組的RMSE都小于序號(hào)較大的數(shù)據(jù)組，即反射系數(shù)擾動(dòng)越小可得到越好的估計(jì)結(jié)果。

5 實(shí)測(cè)數(shù)據(jù)驗(yàn)證

為了驗(yàn)證本文所建模型及反演方法在實(shí)測(cè)數(shù)據(jù)中的有效性，對(duì)某米波雷達(dá)實(shí)測(cè)數(shù)據(jù)進(jìn)行處理。某航班航跡圖如圖6所示。該航班自北向南平穩(wěn)飛行，高度約11 km，目標(biāo)距離由170~300 km變化，仰角在至之間，雷達(dá)3 dB波束寬度約為，因此該航班位于低仰角區(qū)。根據(jù)一次雷達(dá)和二次雷達(dá)關(guān)聯(lián)結(jié)果，將二次雷達(dá)的仰角信息作為真值反演多徑模型，反射面高度和反射系數(shù)幅度反演結(jié)果如圖7所示。可見(jiàn)該陣地較為平坦，反射面高度變化不大，反射系數(shù)在~之間。進(jìn)而對(duì)擾動(dòng)反射系數(shù)進(jìn)行參數(shù)反演，分別取點(diǎn)跡號(hào)30~40和80~90的點(diǎn)跡作擾動(dòng)反射系數(shù)估計(jì)，結(jié)果如圖8，圖9所示，圖中各曲線表現(xiàn)了不同點(diǎn)跡的擾動(dòng)反射系數(shù)估計(jì)結(jié)果?？梢?jiàn)相鄰點(diǎn)跡的擾動(dòng)反射系數(shù)相關(guān)程度較高，圖中表現(xiàn)為相鄰點(diǎn)跡曲線間聚合程度較高。

對(duì)于間隔較遠(yuǎn)的點(diǎn)跡(對(duì)比圖8和圖9)，隨著目標(biāo)運(yùn)動(dòng)，地面反射區(qū)域不再相同，擾動(dòng)反射系數(shù)也隨之變化，圖中反映為不同點(diǎn)跡間擾動(dòng)反射系數(shù)曲線的漸變特征。值得注意的是，圖7~圖9中出現(xiàn)了反射波強(qiáng)度大于直達(dá)波的現(xiàn)象，分析原因如下：如果反射面為完全光滑的水平面，一般假設(shè)反射波經(jīng)由地面反射點(diǎn)產(chǎn)生，此時(shí)反射波能量與直達(dá)波能量大小一致，不會(huì)大于直達(dá)波能量。而實(shí)際中地面反射并非只來(lái)自地面等效反射點(diǎn)一點(diǎn)，而是來(lái)自地面菲涅爾反射區(qū)。若地面菲涅爾反射區(qū)高低不一、或具有一定弧度，可能會(huì)對(duì)反射波造成聚焦效應(yīng)，類(lèi)似于凹面鏡對(duì)光源的聚焦，這種情形下反射波經(jīng)地面多點(diǎn)、多條路徑反射至雷達(dá)(這也是本文所分析地面擾動(dòng)的來(lái)源之一)，因此造成某些點(diǎn)跡反射波強(qiáng)度大于直達(dá)波的現(xiàn)象。且該現(xiàn)象應(yīng)該具有連續(xù)性，與反射面的地形有關(guān)。

對(duì)該航班進(jìn)行仰角估計(jì)，不同算法的DOA估計(jì)結(jié)果如圖10所示，其中“”表示前后向空間平滑算法[15]，“”表示數(shù)字波束形成，真值為二次雷達(dá)應(yīng)答的仰角信息，“反演前”表示沒(méi)有進(jìn)行參數(shù)反演的合成導(dǎo)向矢量最大似然算法，其反射面高度和反射系數(shù)取固定值，分別取0 m和，“反演后”表示通過(guò)前文所述的參數(shù)反演方法，獲得反射面高度和擾動(dòng)反射系數(shù)等先驗(yàn)信息的合成導(dǎo)向矢量最大似然算法。由圖10可見(jiàn)由于角度分辨率不夠和多徑信號(hào)的模型失配等問(wèn)題，傳統(tǒng)算法(DBF和SSMUSIC)無(wú)法在低仰角區(qū)得到正確結(jié)果。而SVML類(lèi)算法測(cè)角效果較好，由于SVML類(lèi)算法結(jié)果與真值接近，特繪

圖4 擾動(dòng)反射系數(shù)估計(jì)結(jié)果 圖5 擾動(dòng)反射系數(shù)均方根誤差隨信噪比變化圖 圖6 某航班航跡圖

制其與真值的仰角估計(jì)誤差圖以方便對(duì)照。如圖11所示。由圖可見(jiàn)參數(shù)反演后SVML算法的仰角估計(jì)精度較參數(shù)反演之前有明顯提升，參數(shù)反演前SVML算法的仰角均方根誤差為，而參數(shù)反演后SVML算法的多徑模型與真實(shí)多徑信號(hào)匹配，其仰角估計(jì)均方根誤差僅為，滿足實(shí)際需求。

為進(jìn)一步驗(yàn)證本文多徑模型與參數(shù)反演方法的有效性，對(duì)兩組同航線航班數(shù)據(jù)進(jìn)行處理。航跡圖如圖12所示，由圖可知航班1與航班2在航跡圖中高度重合。對(duì)航線所在區(qū)域進(jìn)行放大，繪制各航班方位角與距離，如圖13所示，不同航班距離相同的點(diǎn)跡方位角相差最大不超過(guò)，因此可認(rèn)為航班1與航班2為同一航線，即相同仰角區(qū)域?qū)?yīng)相類(lèi)似的地形。對(duì)航班1與航班2進(jìn)行參數(shù)反演，主要參數(shù)如圖14所示，可見(jiàn)在相同的仰角區(qū)域其模型參數(shù)變化趨勢(shì)一致，這也是后續(xù)利用反演參數(shù)提高算法性能的前提。將航班1反演的參數(shù)代入航班2中進(jìn)行仰角估計(jì)，結(jié)果如圖15所示，各算法標(biāo)識(shí)如前文所述。由圖15可知參數(shù)反演的優(yōu)勢(shì)在于低仰角，此時(shí)由于孔徑限制，大部分算法由于無(wú)法分辨直達(dá)波和反射波而算法失效，而SVML類(lèi)算法也需要更準(zhǔn)確的模型參數(shù)才能對(duì)目標(biāo)仰角進(jìn)行跟蹤與估計(jì)，反應(yīng)了參數(shù)反演的必要性。

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

本文研究了低仰角目標(biāo)多徑模型及其反演方法，實(shí)現(xiàn)了經(jīng)典多徑模型下的反射系數(shù)和反射面高

圖10 不同算法仰角估計(jì)結(jié)果對(duì)比圖　　　 圖11 仰角估計(jì)誤差結(jié)果對(duì)比圖　　　　 圖12 各航班航跡圖

圖13 各航班方位角與距離　 圖14 參數(shù)反演結(jié)果 圖15 航班2仰角估計(jì)結(jié)果

度的參數(shù)反演；同時(shí)為解決實(shí)際陣地中多徑信號(hào)的模型失配問(wèn)題，提出干擾多徑模型，將復(fù)雜反射面對(duì)多徑信號(hào)的影響建模為干擾反射系數(shù)，合理解釋了復(fù)雜陣地下傳統(tǒng)算法失效的原因；最后反演擾動(dòng)反射系數(shù)，得到其最大似然估計(jì)。仿真數(shù)據(jù)分析表明本文方法在較高信噪比時(shí)能獲得良好的參數(shù)反演效果；實(shí)測(cè)數(shù)據(jù)驗(yàn)證了擾動(dòng)多徑模型的合理性，并通過(guò)參數(shù)反演方法反演擾動(dòng)反射系數(shù)，反映出復(fù)雜情形下陣地反射特性的擾動(dòng)和漸變特征。本文研究成果表明，若需進(jìn)一步提高測(cè)高算法在實(shí)際陣地中的精度，必須考慮復(fù)雜反射面的對(duì)多徑的擾動(dòng)影響。

[1] BARTON D K. Low-angle radar tracking[J]., 1974, 62(6): 587-704. doi: 10.1109/PROC.1974. 9509.

[2] LO T and LITVA J. Use of a highly deterministic multipath signal model in low-angle tracking[J].,,1991, 138(2): 163-171. doi: 10.1049/ip-f-2.1991.0022.

[3] 陳伯孝, 胡鐵軍, 鄭自良, 等. 基于波瓣分裂的米波雷達(dá)低仰角測(cè)高方法及其應(yīng)用[J]. 電子學(xué)報(bào), 2007, 35(6): 1021-1025. doi: 10.3321/j.issn:0372-2112.2007.06.003.

CHEN Baixiao, HU Tiejun, ZHENG Ziliang,. Method of altitude measurement based on beam splitin VHF radar and its application[J]., 2007, 35(6): 1021-1025. doi: 10.3321/j.issn:0372-2112.2007.06.003.

[4] 洪升, 萬(wàn)顯榮, 柯亨玉. 空間色噪聲背景下雙基地多輸入多輸出雷達(dá)低仰角估計(jì)方法[J]. 電子與信息學(xué)報(bào), 2015, 37(1): 15-21. doi: 10.11999/JEIT140290.

HONG Sheng, WAN Xianrong, and KE Hengyu. Low- elevation estimation for bistatic MIMO radar in spatially colored noise[J].&, 2015, 37(1): 15-21. doi: 10.11999/JEIT140290.

[5] 徐振海, 黃坦, 熊子源, 等. 基于頻率分集的陣列雷達(dá)低角跟蹤算法[J]. 國(guó)防科技大學(xué)學(xué)報(bào), 2014(2): 93-98. doi: 10.11887 /j.cn.201402016.

XU Zhenhai, HUANG Tan, XIONG Ziyuan,. Low angle tracking algorithm using frequency diversity for array radar[J]., 2014(2): 93-98. doi: 10.11887/j.cn.201402016.

[6] PARK D, YANG E, AHN S,. Adaptive beamforming for low-angle target tracking under multipath interference[J]., 2014, 50(4): 2564-2577. doi: 10.1109/TAES.2014.130185.

[7] 胡曉琴, 陳建文, 王永良. 米波雷達(dá)測(cè)高多徑模型研究[J]. 電波科學(xué)學(xué)報(bào), 2008, 23(4): 651-657. doi: 10.3969/j.issn. 1005-0388.2008.04.011.

HU Xiaoqin, CHEN Jianwen, and WANG Yongliang. Research on meter-wave radar height-finding multipath model[J]., 2008, 23(4): 651-657. doi: 10.3969/j.issn.1005-0388.2008.04.011.

[8] ZHU Wei and CHEN Baixiao. Novel methods of DOA estimation based on compressed sensing[J]., 2013, 26(1): 113-123. doi: 10.1007/s11045-013-0239-2.

[9] 王園園, 劉崢, 曹運(yùn)合. 基于壓縮感知的米波雷達(dá)低空測(cè)角算法[J]. 系統(tǒng)工程與電子技術(shù), 2014(4): 667-671. doi: 10.3969/j.issn.1001-506X.2014.04.10.

WANG Yuanyuan, LIU Zheng, and CAO Yunhe. Low-angle estimation method based on compressed-sensing for meter-wave radar[J]., 2014(4): 667-671. doi:10.3969/j.issn.1001-506X.2014.04.10.

[10] YANG Zai, XIE Lihua, and ZHANG Cishen. A discretization-free sparse and parametric approach for linear array signal processing[J]., 2014, 62(19): 4959-4973. doi: 10.1109/TSP. 2014.2339792.

[11] ZHU Wei and CHEN Baixiao. Altitude measurement based on terrain matching in VHF array radar[J]., 2013, 32(2): 647-662. doi: 10.1007 /s00034-012-9472-4.

[12] HAN Yinghua, WANG Jinkuan, ZAHO Qiang,. A blind central DOA estimation algorithm in radar low-angle tracking environment[C]. 2009 IET International Radar Conference, Guilin, 2009: 20-22.

[13] BOMAN K and STOICA P. Low angle estimation: models, methods, and bounds[J]., 2001, 11(1): 35-79. doi: 10.1006/dspr.2000.0373.

[14] STOICA P and NEHORAI A. MUSIC, maximum likelihood, and Cramer-Rao bound[J].,, 1989, 37(5): 720-741. doi: 10.1109/29.17564.

[15] SHAN Tiejun, WAX M, and KAILATH T. On spatial smoothing for direction-of-arrival estimation of coherent signals[J].,, 1985, 33(4): 806-811. doi: 10.1109/TASSP. 1985.1164649.

鄭軼松： 男，1990年生，博士生，研究方向?yàn)殛嚵行盘?hào)處理、米波雷達(dá)低仰角測(cè)高方法.

陳伯孝： 男，1966年生，教授，博士生導(dǎo)師，研究方向包括新體制雷達(dá)系統(tǒng)設(shè)計(jì)及其實(shí)現(xiàn)、雷達(dá)信號(hào)處理、目標(biāo)精確制導(dǎo)與跟蹤等.

Multipath Model and Inversion Method for Low-angle Target in Very High Frequency Radar

ZHENG Yisong CHEN Baixiao

(,,710071,)

The existing methods of altitude measurement for low-angle targets adopt the specular reflection surface model, and the direct and multipath signals are considered as two correlated far-field point sources. However, in reality, the wavefront of multipath signal is distorted by irregular reflection surface, and the far-field point source model is not enough to describe the multipath signal. To deal with this model mismatch problem, the low-angle multipath model is mainly studied. This paper begins with a discussion of classical multipath model and is followed by the inversion method of reflection coefficient and the height of reflection surface. Then the perturbation of the multipath signal caused by irregular reflection surface is modeled as perturbational reflection coefficient and a perturbational multipath model is developed with a maximum likelihood method to invert the proposed parameter. Simulation data processing results validate the effectiveness of the inversion method. The effectiveness of the proposed model and inversion method are validated by measured data processing results. These research results can provide valuable information for enhancing the applicability of the low-angle altitude measurement method in practical situations.

High frequency radar; Multipath model; Altitude measurement of low-angle target; Perturbational reflection coefficient

The National Natural Science Foundation of China (61571344)

TN958

A

1009-5896(2016)06-1468-07

10.11999/JEIT151013

2015?09?09；改回日期：2016-01-22；網(wǎng)絡(luò)出版：2016-03-29

鄭軼松 zhengys90@163.com

國(guó)家自然科學(xué)基金(61571344)