李立峰,馮曉東

(1.裝備工程技術(shù)研究實(shí)驗(yàn)室,河北 石家莊 050081；2.中國電子科技集團(tuán)公司第五十四研究所,河北 石家莊 050081；3.裝備發(fā)展部軍事代表局駐石家莊地區(qū)軍事代表室，河北 石家莊050081)

低信噪比小樣本下的短波測向技術(shù)

李立峰1,2,馮曉東3

(1.裝備工程技術(shù)研究實(shí)驗(yàn)室,河北 石家莊 050081；2.中國電子科技集團(tuán)公司第五十四研究所,河北 石家莊 050081；3.裝備發(fā)展部軍事代表局駐石家莊地區(qū)軍事代表室，河北 石家莊050081)

針對短波信道衰落導(dǎo)致接收信噪比低引起測向性能下降的問題，提出了譜線增強(qiáng)技術(shù)，利用窗函數(shù)對陣元接收信號在測向前進(jìn)行復(fù)指數(shù)加權(quán)處理，達(dá)到提高信噪比的目的。針對傳統(tǒng)空間譜測向算法在小樣本數(shù)或信號源相近等條件下分辨率嚴(yán)重下降的問題，提出了基于譜函數(shù)二階偏導(dǎo)數(shù)的處理方法，達(dá)到提高測向分辨率的目的。最后對測向性能進(jìn)行了分析驗(yàn)證，證實(shí)了該方法的有效性和適用性。

短波信道；低信噪比；譜線增強(qiáng)；二階偏導(dǎo)數(shù)

0 引言

短波測向[1]存在信號密集度高、接收場強(qiáng)差異大以及示向度游動大等特點(diǎn)，由于短波信道[2]存在著嚴(yán)重的衰落現(xiàn)象，有可能采集到的信號正好達(dá)到衰落最低點(diǎn)，達(dá)不到測向所需的信噪比要求，導(dǎo)致測向結(jié)果出現(xiàn)較大的偏差[3]。傳統(tǒng)空間譜測向算法[4]為了能對樣本空間更加準(zhǔn)確地估計(jì)，一般采用比較大的樣本數(shù)，但是樣本數(shù)的過分增多會帶來計(jì)算量大的問題，陣列分辨空間信號的能力與陣列孔徑有關(guān)，陣列孔徑越大，陣列的分辨能力越強(qiáng)。所以，要提高算法的分辨率，就必須增大天線的孔徑，但是在實(shí)際測向系統(tǒng)中，陣列孔徑往往是有限制的，那么如何在陣列孔徑一定且信號環(huán)境非理想的條件下獲得較高的分辨率，是值得研究的問題。

1 提高接收信號信噪比的處理方法

在低信噪比環(huán)境下，空間譜測向性能下降[5]。為了提高測向信噪比[6]，可以在測向前對采樣信號進(jìn)行譜線增強(qiáng)[7]的處理，如圖1所示。

譜線增強(qiáng)處理就是對每個(gè)接收通道的采樣信號進(jìn)行加權(quán)的過程，處理過程如下：

① 根據(jù)采樣信號的信噪比，選擇適當(dāng)?shù)臅r(shí)間窗長度L；

② 產(chǎn)生復(fù)指數(shù)時(shí)間窗函數(shù)w：

w=[1,ej2πf0Ts,ej2πf02Ts,……,ej2πf0(L-1)Ts]，

式中，Ts為采樣信號的采樣周期；

③ 運(yùn)用窗函數(shù)w對采樣信號進(jìn)行復(fù)指數(shù)加權(quán)處理。

對譜線增強(qiáng)后的數(shù)據(jù)進(jìn)行如文獻(xiàn)[8-11]所述的非均勻取樣，進(jìn)一步減弱噪聲信號間的相關(guān)性，從而進(jìn)一步提高對微弱信號[12]的估計(jì)性能。

圖1 接收信號譜線增強(qiáng)過程

2 基于二階偏導(dǎo)數(shù)的高分辨處理方法

為了提高在低信噪比和小樣本數(shù)下二維DOA的分辨率和估計(jì)性能，提出了一種基于二階偏導(dǎo)數(shù)[13]的二維DOA估計(jì)新方法。

假設(shè)方位角θ范圍為Rθ，搜索間隔為Δθ，仰角φ范圍為Rφ，搜索間隔為Δφ，令N=Rθ/Δθ，K=Rφ/Δφ，空間譜函數(shù)P可表示[14]為：

由上式可以看出，對于每個(gè)給定的坐標(biāo)(i,j),i=1、2、…、N,j=1、2、…、K，都對應(yīng)著一個(gè)確定的實(shí)數(shù)pij，設(shè)自變量

θ∈D1={θ1,θ2,…,θN},θi+1-θi=Δθ，

φ∈D2={φ1,φ2,…,φK},φj+1-φj=Δφ，

這里當(dāng)i≠j時(shí)，θi與θj互不相同，φi與φj互不相同，則：

pij=f(θi,φj)，

式中，f為二元離散函數(shù)。由文獻(xiàn)[15]可知，任意二元有限離散函數(shù)在其定義域內(nèi)都存在偏導(dǎo)，且其性質(zhì)與連續(xù)函數(shù)的常規(guī)偏導(dǎo)數(shù)非常相似。因此，在點(diǎn)(θk,φl)處，可得二元有限離散函數(shù)p對自變量θ的一階偏導(dǎo)數(shù)為：

利用MUSIC譜函數(shù)極大點(diǎn)對方位角和俯仰角的二階偏導(dǎo)數(shù)都小于零的特性，對空間譜進(jìn)行負(fù)向譜峰搜索可得到目標(biāo)的DOA。

3 性能分析

① 譜線增強(qiáng)處理對測向性能的影響

8單元均勻圓陣，真實(shí)來波方向65°，信號頻率10.83 MHz，接收信號信噪比低于5 dB，對譜線增強(qiáng)處理前后的測向性能進(jìn)行對比。

圖2中，橫坐標(biāo)表示方向，單位為(°)，縱坐標(biāo)表示合成增益，單位dB。黑實(shí)線表示譜線增強(qiáng)后的波束圖，波束指向64°，黑虛線表示譜線增強(qiáng)前的波束圖，波束指向60°。從試驗(yàn)結(jié)果可以看出，譜線增強(qiáng)后空間譜圖變得很尖銳，測向分辨率得到明顯提高，測向結(jié)果更加接近真實(shí)方位。

圖2 譜線增強(qiáng)處理前后的測向?qū)Ρ?/p>

② 譜函數(shù)二階導(dǎo)數(shù)處理對測向性能的影響

8單元均勻圓陣，真實(shí)來波方向65°，信號頻率10.83 MHz，樣本點(diǎn)數(shù)100，對譜函數(shù)二階導(dǎo)數(shù)處理前后的測向性能進(jìn)行對比。

圖3中，橫坐標(biāo)表示測向次數(shù)，縱坐標(biāo)表示測向結(jié)果，單位為(°)。黑Ο線表示譜二階導(dǎo)數(shù)處理后的測向結(jié)果，測向結(jié)果都在65°附近，黑*線表示譜二階導(dǎo)數(shù)處理前的測向結(jié)果，相對于真值浮動較大。從試驗(yàn)結(jié)果可以看出，經(jīng)過對譜函數(shù)作二階偏導(dǎo)數(shù)處理后，出現(xiàn)較大偏差的概率大大減小，測向精度、測向分辨率有所提高。

圖3 譜函數(shù)二階導(dǎo)數(shù)處理前后的測向結(jié)果對比

③ 聯(lián)合處理對測向性能的影響

綜合利用以上2種處理方法，對處理前后的測向性能進(jìn)行對比。

圖4中，橫坐標(biāo)表示方向，單位為(°)，縱坐標(biāo)表示合成增益，單位dB。黑實(shí)線表示聯(lián)合處理后的波束圖，波束指向64°，黑虛線表示譜線增強(qiáng)后的波束圖，波束指向64°，黑點(diǎn)劃線表示譜線增強(qiáng)前的波束圖，波束指向68°。從試驗(yàn)結(jié)果可以看出，通過譜線增強(qiáng)和對譜函數(shù)二階導(dǎo)數(shù)處理后，可以在波達(dá)方向附近形成更加尖銳的譜峰，測向精度和分辨率具有明顯的提高。

圖4 聯(lián)合處理前后的測向?qū)Ρ?/p>

4 結(jié)束語

本文通過譜線增強(qiáng)技術(shù)有效地提高了接收信號的信噪比，解決了短波信道的衰落導(dǎo)致接收信號太弱而不能正確測向的技術(shù)難題，明顯提高了測向精度和穩(wěn)定性。通過對譜函數(shù)求二階偏導(dǎo)數(shù)處理，有效解決了小樣本條件下測向分辨率低的問題，可以明顯提高測向分辨率和測向精度。綜合2種處理方法，可以同時(shí)提高低信噪比、小樣本數(shù)下的測向精度和分辨率。新方法在短波偵測系統(tǒng)中得以應(yīng)用，在天波試驗(yàn)中得到了充分的驗(yàn)證。

[1] Rubsaman M,Gershman A B.Direction-of-arrival Estimation for Non Uniform Sensor Arrays：from Manifold Separation to Fourier Domain MUSIC Methods[J].IEEE Trans on Signal Processing,2009,57(2)：588-599.

[2] Watterson C C,Juroshek J,Bensema W.Experim Ental Confirmation of an HF Channel Model[J].IEEE Trans.Commun Technol,1970,18(6)：792-803.

[3] 涂旭東.寬帶短波信道特性和建模研究[D].成都：電子科技大學(xué)，2004.

[4] Schmidt R O.Multiple Emitter Location and Signal Parameter Estimation[J].IEEE Trans.on Antennas and Propagation,1986,34(3)：276-280.

[5] 蔡臻祥，馬義德，李紅領(lǐng).多徑效應(yīng)對測向精度試驗(yàn)的影響[J].無線電工程，2006，36(2)：37-39.

[6] 潘玉平，齊青茂.波束形成技術(shù)在短波偵測中的應(yīng)用[J].無線電工程，2007，37(5)：33-34.

[7] 劉輝濤，叢衛(wèi)華，潘 翔.窄帶弱信號的線譜檢測[J].浙江大學(xué)學(xué)報(bào)，2007，41(12)：2048-2051.

[8] 李 安，王 殊，屈曉旭，等.基于非均勻采樣的空時(shí)級聯(lián)頻率和到達(dá)角估計(jì)[J].系統(tǒng)工程與電子技術(shù)，2010，32(3)：478-482.

[9]LeeJY，GreengardL.TheType3NonuniformFFTandItsApplications[J].JournalofComputationalPhysics，2005，206(1)：1-5.

[10]KubatJA，ChouJJ，RovnyakD.NonuniformSamplingandMaximumEntropyReconstructionAppliedtotheAccurateMeasurementofResidualDipolarCouplings[J].JournalofMagneticResonance，2007，186(2)：201-211.

[11]LiuYJ，MengHD，WangDS，etal.AdaptiveStaggeringTimeEstimationforTargetTrackinginPeriodicNonuniformSamplingSystem[J].ElectronicsLetters，2007，43(24)：1385-1387.

[12] 張家平，劉 青.低信噪比中MUSIC算法研究[J].通信技術(shù)，2009，42(1)：87-89.

[13] 司偉建，藍(lán)曉宇.提高空間譜算法分辨率的DOA估計(jì)[J].應(yīng)用科學(xué)學(xué)報(bào)，2013，31(4)：382-386.

[14] 劉 劍，于紅旗.基于二階預(yù)處理的共軛擴(kuò)展MUSIC算法[J].系統(tǒng)工程與電子技術(shù)，2008，30(1)：57-60.

[15] 王文麗.有限離散函數(shù)的導(dǎo)數(shù)和性質(zhì)[J].大學(xué)數(shù)學(xué)，2006，21(3)：110 -113.

HF DOA Technique under Low SNR and Small Number of Snapshots

LI li-feng1,2，F(xiàn)ENG Xiao-dong3

(1.Equipment Engineering Technology Research Laboratory,Shijiazhuang Hebei 050081,China;2.The 54th Research Institute of CETC,Shijiazhuang Hebei 050081,China;3.Military Representative Office of Equipment Development Department in Shijiazhuang,Shijiazhuang Hebei 050081,China)

The performance of DOA decreases under low signal noise ratio(SNR) because of shortwave channel fading.In view of this problem,a method named spectral-line enhancement algorithm is proposed.In order to improve SNR,the array arriving data is processed by complex-exponential weighting using window function before DOA.The resolution decreases in direction of arrival(DOA) estimation under severe environments such as small number of snapshots or incident waves from close angles.In view of this problem,a high-resolution estimation method is proposed based on calculation of the second derivative of spectrum.Finally,the test results show that this method is available and applicable.

shortwave channel;low SNR;spectral-line enhancement;second derivative

10.3969/j.issn.1003-3114.2017.02.11

李立峰,馮曉東.低信噪比小樣本下的短波測向技術(shù)[J].無線電通信技術(shù),2017,43(2)：44-46,54.

2016-11-02

國家部委基金資助項(xiàng)目

李立峰(1977—)，男，高級工程師，主要研究方向：通信對抗、陣列信號處理。馮曉東(1977—)，男，工程師，主要研究方向：通信工程、通信信號處理。

TN911

A

1003-3114(2017)02-44-3