聯(lián)合DOA估計與主波束形成的干擾抑制方法

劉海濤, 劉亞洲, 張學(xué)軍

(1. 中國民航大學(xué) 天津市智能信號與圖像處理重點試驗室, 天津 300300； 2. 北京航空航天大學(xué) 電子信息工程學(xué)院, 北京 100191)

?

聯(lián)合DOA估計與主波束形成的干擾抑制方法

劉海濤1, 劉亞洲1, 張學(xué)軍2

(1. 中國民航大學(xué) 天津市智能信號與圖像處理重點試驗室, 天津 300300； 2. 北京航空航天大學(xué) 電子信息工程學(xué)院, 北京 100191)

為克服測距儀(DME)發(fā)射的高強度脈沖信號干擾L頻段數(shù)字航空通信系統(tǒng)1(L-DACS1)OFDM接收機的問題，提出聯(lián)合DOA估計與主波束形成的干擾抑制方法. 首先接收機通過波達方向矩陣算法估計接收信號來向；然后基于信號來向通過線性約束最小方差波束形成算法提取各個方向信號；隨后通過頻域功率與時域功率比較算法分辨各來向信號的類型；最后輸出OFDM直射徑信號. 仿真結(jié)果表明：該方法可顯著克服DME與OFDM散射徑信號干擾，提高L-DACS1系統(tǒng)鏈路傳輸?shù)目煽啃?

正交頻分復(fù)用; 測距儀脈沖干擾; 均勻圓陣; 波達方向估計; 線性約束最小方差

L頻段數(shù)字航空通信系統(tǒng)1(L-DACS1)是民航未來沿陸地航路部署的空-地蜂窩移動通信系統(tǒng)，主要為陸地航路、終端區(qū)及機場的航空器提供空中交通管制、航空公司運營管理數(shù)據(jù)通信服務(wù)，是民航未來航空移動通信系統(tǒng)的主要技術(shù)手段之一[1]. 依據(jù)國際民航組織的規(guī)劃，L-DACS1系統(tǒng)將以內(nèi)嵌的方式部署在L頻段測距儀(DME)的波道間，占用傳輸帶寬500KHz[2-3]. 由于DME與OFDM信號頻譜存在部分交疊，且DME系統(tǒng)發(fā)射功率較高，不可避免產(chǎn)生DME信號干擾L-DACS1系統(tǒng)OFDM接收機的問題[4-5]，因此開展OFDM接收機DME干擾抑制的研究具有重要意義.

目前，L-DACS1系統(tǒng)DME干擾抑制的研究主要集中在單天線干擾消除及陣列天線干擾抑制兩個方面. 在單天線干擾消除方面，文獻[6-7]利用DME干擾信號在時域呈現(xiàn)脈沖干擾的特點，提出脈沖熄滅干擾抑制法，該方法主要缺點：在時變信道環(huán)境下，脈沖熄滅門限的精確設(shè)置非常困難. 為解決脈沖熄滅導(dǎo)致OFDM信號產(chǎn)生子載波間干擾(ICI)的問題，文獻[8-10]提出脈沖熄滅ICI干擾補償方法，該方法主要缺點：ICI干擾補償法要求精確知曉各個子信道的衰落信息，而在干擾環(huán)境下OFDM接收機各個子信道衰落信息的精確獲得本身就非常困難. 在陣列天線在航空移動通信應(yīng)用方面，文獻[11]提出在無人機平臺中搭載移動基站為偏遠地區(qū)提供移動通信服務(wù)的構(gòu)想，文獻[12-14]構(gòu)建了機載實驗系統(tǒng)驗證了在無人機平臺使用陣列天線進行大容量空-地數(shù)據(jù)通信的可行性. 另外，在基于陣列天線的L-DACS1系統(tǒng)OFDM接收機DME干擾抑制方面，文獻[15]提出利用正交投影算法消除高強度DME干擾，然后采用盲波束形成方法提取OFDM直射徑信號的方法，該方法存在不足：當(dāng)測距儀干擾較小時，正交投影干擾抑制性能較差. 在正交投影消除脈沖干擾方法基礎(chǔ)上,文獻[16]進一步提出利用CLEAN算法估計OFDM直射徑信號來向的方法，并通過常規(guī)波束形成算法提取OFDM直射徑信號，該方法主要不足：當(dāng)接收OFDM信號功率較低且數(shù)據(jù)快拍較少時，提出方法的DOA估計性能較差.

針對高強度測距儀信號干擾L-DACS1系統(tǒng)OFDM接收機的問題，本文提出聯(lián)合DOA估計與主波束形成的干擾抑制方法. 首先接收機利用DME與OFDM信號空域來向的差異，通過波達方向矩陣法獲取各來向信號的DOA；然后利用各個來向信號的DOA信息通過線性約束最小方差(LCMV)波束形成算法確定各來向信號波束形成的權(quán)值，并利用波束形成權(quán)值對接收信號進行波束形成以分離出各個來向信號；最后利用OFDM與DME信號時頻特征的差異，通過頻域功率與時域功率比較算法獲得OFDM直射徑信號. 仿真結(jié)果顯示：提出的方法可有效抑制DME及OFDM散射徑信號干擾，提高鏈路傳輸?shù)目煽啃?

1 系統(tǒng)模型

1.1 聯(lián)合DOA估計與主波束形成的OFDM接收機 圖1為聯(lián)合DOA估計與主波束形成的OFDM接收機原理框圖. 接收機均勻圓陣輸出信號通過A/D轉(zhuǎn)換為數(shù)字基帶信號，數(shù)字基帶信號通過波達方向估計器估計各個來向信號的DOA信息. 隨后陣列接收信號及DOA信息同時送入主波束形成器，主波束形成器通過線性約束最小方差算法(LCMV)計算各來向信號的波束形成權(quán)值，然后通過波束形成提取各個來向信號. 提取的各個來向信號同時送入信號分類器，在信號分類器中通過頻域功率與時域功率比較算法分辨DME及OFDM信號，并輸出OFDM直射徑信號. 信號分類器輸出的OFDM直射徑信號經(jīng)多普勒頻偏估計與補償后移除循環(huán)前綴(CP)，送入FFT轉(zhuǎn)換器變換為頻域信號. 頻域信號經(jīng)頻域下采樣后送入信道估計及均衡器完成信道均衡，均衡器輸出信號經(jīng)解調(diào)器、解交織器后送入信道譯碼器恢復(fù)出原始發(fā)送比特序列估值.

圖1 聯(lián)合DOA估計與主波束形成的OFDM接收機原理

1.2 DOA估計算法

為解決相干信號環(huán)境下均勻圓陣二維DOA估計問題，采用文獻[17]提出的波達方向矩陣法解決相干信源的二維DOA估計問題. 接收機使用上下平行雙均勻圓陣(如圖2所示)，上下間隔d=λ/2，上下子陣的陣元位置一一對應(yīng)，每個子陣的陣元數(shù)均為M，圓半徑均為R，信號波長為λ.

m=0,1,…,M-1.

圖2 雙均勻圓陣結(jié)構(gòu)圖

式中: si(t)代表陣列天線接收的第i個信號源，{θi,φi}代表第i個信號源的方位角和俯仰角，n1m(t)代表信道輸入的窄帶高斯白噪聲，I代表陣列接收信號源總數(shù). 類似的，子陣2上坐標(biāo)為(xm,ym,-d)的陣元輸出信號為

利用文獻[17]的DOA方法計算信號俯仰角和方向角得

式中：R代表按照文獻[17]中方法構(gòu)造的空域平滑波達方向矩陣，ηi和εi分別代表R的特征值和對應(yīng)特征矢量，εi(j)代表εi中的第j個元素，Arg(·)代表幅角主值運算.

1.3 主波束形成算法

圖3 主波束形成器

其中，I代表陣列天線接收OFDM信號及測距儀信號源總數(shù)，wi代表第i個來向信號波束形成權(quán)值[19]：

1.4 信號分類算法

圖4給出接收機信號分類器的原理框圖. 信號分類器由頻域功率比較器和時域功率比較器兩個單元組成. 其中，頻域功率比較器用于分辨DME與OFDM信號，時域功率比較器用于分辨OFDM直射徑與散射徑信號.

圖4 信號分類器

1.2.1 頻域功率比較器

根據(jù)OFDM及DME信號的頻域特征[4]，OFDM信號的能量主要集中于[-250 kHz～+250 kHz]，測距儀信號的能量主要位于±250 kHz處,因此可通過比較接收信號在頻域±250 kHz及直流子信道附近平均功率的方法來分辨接收信號的類型.

(1)

(2)

(3)

1.4.2 時域功率比較器

如果待檢測信號被判定為OFDM信號，則還需要進一步判定待檢測信號是OFDM直射徑還是散射徑信號. 由于航空移動信道呈現(xiàn)為賴斯信道[18]，信道的賴斯因子取值為15～20 dB，因此OFDM直射徑信號功率遠高于散射徑，可直接通過比較各個OFDM信號的功率來分辨OFDM直射徑.

2 數(shù)值仿真及分析

2.1 仿真環(huán)境設(shè)置

為驗證聯(lián)合DOA估計與主波束形成干擾抑制方法的正確性，按照L-DACS1系統(tǒng)技術(shù)規(guī)范[3]設(shè)計實現(xiàn)L-DACS1系統(tǒng)陣列天線干擾抑制仿真系統(tǒng)，仿真參數(shù)如表1所示.

表1 L-DACS1系統(tǒng)仿真環(huán)境

2.2 DOA估計性能

圖5 給出不同信源的DOA估計性能，橫坐標(biāo)代表方位角(度)，縱坐標(biāo)代表俯仰角(度). 仿真實驗中，OFDM直射徑來向(60°,50°)，散射徑來向(140°, 40°)，信噪比SNR=10 dB；DME信號1來向(100°, 10°)，信干比SIR=-3 dB，載波偏置+500 KHz；DME信號2來向(200°, 30°)，信干比SIR=-10 dB，載波偏置-500 KHz. 由圖5可觀測到：1) OFDM直射徑與DME信號來向的估計值與預(yù)設(shè)值完全一致；2) OFDM散射徑來向的估計值與預(yù)設(shè)值存在微小的偏差. 圖5 DOA估計性能(SNR=10 dB，SIR= -3/-10 dB，賴斯因子=10 dB，200次仿真實驗)

Fig.5 DOA estimation (SNR=10 db,SIR=-3/10 dB,rice factor= 10 dB,simulations for 200 times)

圖6給出不同信源DOA估計的根均方誤差曲線，橫坐標(biāo)代表信噪比(dB)，縱坐標(biāo)代表根均方誤差值(度). 根均方誤差值定義為

圖6 DOA估計根均方誤差曲線(賴斯因子=10 dB)

Fig.6 Root mean square error curve of DOA estimation(Rice factor=10 dB)

2.3 LCMV波束形成性能

圖7和圖8分別給出LCMV算法波束圖，圖7橫坐標(biāo)代表方位角，圖8橫坐標(biāo)代表俯仰角，兩圖縱坐標(biāo)均代表歸一化波束增益(dB). 由圖7與圖8可觀測到：陣列天線的波束在干擾信號1來向(100°, 10°)、干擾信號2來向(200°, 30°)及OFDM散射徑來向(140°, 40°)形成零陷，零陷衰減達-300 dB，表明了提出方法可有效抑制DME及OFDM散射徑信號的干擾.

圖7 LCMV算法在方位角上的波束圖(SNR=10 dB, SIR=-3 dB/-10 dB,200次蒙特卡羅實驗)

Fig.7 LCMV beamforming on azimuth angle(SNR=10 dB,SIR=-3 dB/-10 dB for Montacarlo test for 200 times)

圖8 LCMV算法在俯仰角上的波束圖(SNR=10 dB, SIR=-3 dB/-10 dB,200次蒙特卡羅實驗)

Fig.8 LCMV beamforming on pitch angle(SNR=10 dB,SIR=-3 dB/-10 dB.Montacarlo test for 200 times)

圖9給出波束形成性能隨DOA估計誤差變化曲線，橫坐標(biāo)代表DOA估計誤差(方向角誤差+俯仰角誤差)，縱坐標(biāo)代表波束形成輸出信號(OFDM直射徑信號)信噪比. 由圖9觀測到：1)輸出信噪比隨DOA估計誤差的增大可逐漸降低；2)理想DOA時，波束形成輸出信噪比為1 dB，對輸入信噪比提升約5 dB.

圖9 DOA估計偏差對波束形成性能影響曲線(輸入信噪比為-4 dB)

Fig.9 Beamforming performance based on the DOA estimation deviation (input SNR =-4 dB)

2.4 信號分類器差錯性能

圖10給出信號分類器判定差錯概率的性能曲線，橫坐標(biāo)代表信噪比，縱坐標(biāo)代表信號分類器判定OFDM直射徑、OFDM散射徑及DME信號出現(xiàn)錯誤的概率. 由圖10觀測到：判定差錯概率隨輸入信噪比的增加而快速降低，當(dāng)輸入信噪比為6 dB時，判定出現(xiàn)差錯的概率為1.0×10-4. 考慮到典型情況下OFDM接收機輸入信噪比高于6 dB，因此提出的信號分類器能夠準確分辨OFDM直射徑、OFDM散射徑及DME信號.

圖10 信號分類器差錯概率性能(OFDM直射徑信號、散射徑信號、單個DME信號、信干比=-10 dB)

Fig.10 Error probability performance of signal classification

2.5 系統(tǒng)比特差錯性能

圖11給出OFDM接收機比特差錯性能曲線，橫坐標(biāo)代表信噪比(dB)，縱坐標(biāo)代表比特差錯概率. 圖11共包含4曲線：標(biāo)有“◆”的曲線代表無干擾信號時系統(tǒng)比特差錯性能；標(biāo)有“■”的曲線代表存在一個干擾信號時系統(tǒng)比特差錯性能；“●”的曲線代表存在兩個干擾信號時系統(tǒng)比特差錯性能；標(biāo)有“▼”的曲線代表按照文獻[15]方法仿真得到的比特差錯性能. 曲線比較表明：所提方法可完全克服DME干擾信號對OFDM接收機鏈路傳輸?shù)挠绊?，且性能?yōu)于文獻[15]提出的方法.

圖11 L-DACS1系統(tǒng)比特差錯性能

3 結(jié) 語

針對高強度測距儀脈沖信號干擾L-DACS1系統(tǒng)OFDM接收機問題，提出聯(lián)合DOA估計與主波束形成的干擾抑制方法. 該方法首先利用DME與OFDM信號空域來向的差異，通過線性約束最小方差波束形成算法分離各個信號，隨后借助OFDM與DME信號的時頻特征的差異，分辨并輸出OFDM直射徑信號. 仿真研究表明：提出的方法可有效克服DME信號及OFDM散射徑的干擾，提高接收機鏈路傳輸?shù)目煽啃?

[1] SCHELL M, EPPLE U, SHUTIN D, et al. L-DACS: future aeronautical communications for air-traffic management[J]. IEEE Communications Magazine, 2014,25(5):104-110.

[2] NEJI N, LACERDA DE R, AZOULAY A, et al. Survey on the future aeronautical communication system and its development for continental communications[J]. IEEE Transactions Vehicular Technology,2013,62(1): 182-191.

[3] SAJATOVIC M, HAINDL B, EHAMMER M, et al. L-DACS1 system definition proposal: deliverable D2[S]. Version 1.0. Brussels: Eurocontrol, 2009:1-173.

[4] EPPLE U, SCHELL M. Overview of legacy systems in L-band and its influence on the future aeronautical communication system LDACS1[J]. IEEE Aerospace and Electronic Systems Magazine, 2014(29)2: 31-37.

[5] EPPLE U, SCHELL M. Overview of interference situation and mitigation techniques for LDACS1[C]//2011 IEEE/AIAA 30th.Digital Avionics Systems Conference (DASC).Seattle,WA: IEEE, 2011: 4C5-1-4C5-12.

[6] EPPLE U, HOFFMANN F, SCHELL M. Modeling DME interference impact on LDACS1[C]//IEEE.Integrated Communications, Navigation and Surveillance Conference (ICNS). Herndon,VA:IEEE, 2012: G7-1-G7-13.

[7] EPPLE U, BRANDES S, GLIGOREVIC S, et al. Receiver optimization for L-DACS1[C]// IEEE/AIAA 28thDigital Avionics Systems Conference. Orlando,FL:IEEE, 2009: 4. B. 1-1-4. B. 1-12.

[8] BRANDES S, EPPLE U, SCELL M. Compensation of the impact of interference mitigation by pulse blanking in OFDM systems[C]//IEEE Global Telecommunications Conference. Honolulu, HI: IEEE, 2009: 1-6.

[9] EPPLE U, SHUTIN D, SCHELL M. Mitigation of impulsive frequency-selective interference in OFDM based systems[J]. IEEE Wireless Communications Letters, 2012, 1(5): 484-487.

[10]LI Q, ZHANG J, XIE J,et al. Iterative interference mitigation and channel estimation for LDACS1[C]// 2014 IEEE/AIAA 33rd Digital Avionics Systems Conference (DASC 2014). Colorado:IEEE,2014: 3B2-1- 3B2-11.

[11]DOVIS F, SELLONE F. Smart antenna system design for airborne GSM base-stations [C]//2000 IEEE Proceedings of the Sensor Array and Multichannel Signal Processing Workshop.Cambridge,MA: IEEE,2000:pp. 429- 433.

[12]MONDIN M, DOVIS F, MULASSANO P. On the use of HALE platforms as GSM base stations[J]. IEEE Personal Communications, 2001, 8(2): 37-44.

[13]AVAGNINA D, DOUIS F, GHIG LIONE A, et al. Wireless networks based on high-altitude platforms for the provision of integrated navigation/communication services[J]. IEEE Communications Magazine, 2002, 40(2): 119-125.

[14]KARAPANTAZIS S, PAVLIDOU F. Broadband communications via high-altitude platforms: a survey[J]. IEEE Communications Surveys & Tutorials, 2005, 7(1): 2-31.

[15]劉海濤，劉亞洲，成瑋，等.聯(lián)合正交投影與盲波束形成的干擾抑制方法[J]. 系統(tǒng)工程與電子技術(shù), 2015,37(8):1180-1186.

LIU H T,LIU Y Z,CHENG W,et al.Interference mitigation method based on subspace projection and blind adaptive beamforming[J]. Systems Engineering and Electronics, 2015,37(8):1180-1186.

[16]劉海濤，劉亞洲，張學(xué)軍. 聯(lián)合正交投影與CLEAN的測距儀脈沖干擾抑制方法[J].信號處理,2015,31(5):536-543.

LIU H T,LIU Y Z,ZHANG X J. DME impulse interference mitigation method based on subspace projection and CLEAN algorithm[J].Journal of Signal Processing,2015,31(5):536-543.

[17]毛維平，李國林，謝鑫.均勻圓陣相干信源二維波達方向估計[J].系統(tǒng)工程與電子技術(shù), 2013,35(8):1596-1601.

MAO W P,LI G L,XIE X.2D-DOA estimation of coherent signals based on uniform circular array[J]. Systems Engineering and Electronics, 2013,35(8):1596-1601.

[18]HAAS E. Aeronautical channel modeling[J]. IEEE Transactions on Vehicular Technology, 2002,51(2):254-264.

(編輯 王小唯 苗秀芝)

Interference mitigation method based on joint DOA estimation and main beam forming

LIU Haitao1, LIU Yazhou1, ZHANG Xuejun2

(1.Tianjin Key Lab for Advanced Signal Processing, Civil Aviation University of China, Tianjin 300300, China;2. College of Electronic and Information Engineering, Beihang University, Beijing 100191, China)

To mitigate the deleterious influence of Distance Measure Equipment (DME) interference on Orthogonal Frequency Division Multiplexing (OFDM) receiver of L-band Digital Aeronautical Communications System1(L-DACS1) operating as an inlay system, a new interference suppression method is proposed based on joint DOA estimation and main beamforming. Firstly, DOA matrix algorithm is used to estimate the DOA of the received signals. Secondly, with the DOA information, Linearly Constrained Minimum Variance(LCMV) beamforming algorithm is utilized to extract the signals in all directions. Finally, frequency-domain power comparison and time-domain power comparison method is proposed to distinguish the received signals and output the direct line-of-sight(LOS) path OFDM signal. Computer simulation results indicate that the proposed method can effectively overcome DME impulse interference and OFDM scattering signal and improve the reliability of the L-DACS1 system.

orthogonal frequency division multiplexing; distance measure equipment pulse interference; uniform circular array; DOA estimation; linearly constrained minimum varianc

10.11918/j.issn.0367-6234.2016.11.016

2015-04-01

國家自然科學(xué)基金(U1233117, 61271404)

劉海濤(1966—)，男，教授，碩士生導(dǎo)師; 張學(xué)軍(1971—)，男，教授，博士生導(dǎo)師

劉海濤，htliu@cauc.edu.cn

TN929.5

A

0367-6234(2016)11-0103-06