彭金龍，陳建春，郭亞萍

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

基于差分時間平滑的多徑信號DOA估計

彭金龍，陳建春，郭亞萍

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

針對移動通信環(huán)境中，多徑信號與獨立信號共存的混合信號模型DOA估計問題，依據(jù)差分方法與時間平滑原理，提出一種改進(jìn)的DOA算法。該算法增加了可分辨信源個數(shù)，同時提高了相干信源DOA估計的性能，仿真實驗證明，該算法在混合信號模型下解相干性能優(yōu)于傳統(tǒng)空間平滑法。

時間平滑;差分;多徑信號;移動通信

在現(xiàn)代移動通信系統(tǒng)中，用戶方位對某些應(yīng)用具有重要作用，例如基于用戶方位的智能天線可提高通信系統(tǒng)頻率和時間復(fù)用率，從而滿足大容量與高質(zhì)量的信號傳輸。移動通信系統(tǒng)中，用戶經(jīng)常處于多徑信號環(huán)境中，例如同一移動用戶的信號經(jīng)過建筑物反射形成多徑傳輸。由于多徑信號來自同一信號源，當(dāng)多徑信號的傳播時延差較小時，通常認(rèn)為這些多徑信號是相干的，從而導(dǎo)致傳統(tǒng)基于特征子空間的DOA估計失敗［1］。

對于相干信源，為了實現(xiàn)DOA估計，常見的一種處理方式是首先采用空間平滑算法解相干［1－3］，然后進(jìn)行DOA估計。但在實際移動通信環(huán)境中，同一基站范圍內(nèi)存在多個用戶通話，空間平滑算法由于陣列孔徑的損失，降低了空間譜分辨率可檢測的信號數(shù)目，并且在獨立信號與相干信號共存環(huán)境中性能惡化。文獻(xiàn)［1，3］分析表明:前向、前后向空間平滑算法分別最多只能估計M/2和2M/3個多徑信號DOA(M為陣元數(shù))，而實際中多徑數(shù)目通常遠(yuǎn)大于此。

時間平滑算法［4－7］利用移動通信信號的衰落特性，在不同的相干時間段內(nèi)取多個數(shù)據(jù)塊，通過對數(shù)據(jù)塊的協(xié)方差矩陣重新整合，使得到的新協(xié)方差矩陣達(dá)到滿秩。時間平滑算法沒有陣列孔徑損失，但是對于獨立信源與相干信源的混合信號模型，無法有效區(qū)分，同時導(dǎo)致解相干的信源數(shù)降低?；诘染嗑€陣的平滑差分方法［8－9］，首先利用傳統(tǒng)的基于特征空間的波達(dá)方向估計方法實現(xiàn)獨立信號的空間譜估計，再利用空間平滑差分方法實現(xiàn)相關(guān)信號的空間譜估計。本文將文獻(xiàn)［8～9］的差分思想應(yīng)用到時間平滑中，提出一種基于差分方法的時間平滑的DOA算法。該方法可以進(jìn)一步增大可估計信源數(shù)目，同時估計性能相對于傳統(tǒng)空間平滑算法有明顯改善，計算量較小。

1 多徑混合信號模型

對于一組多徑信號，在時刻的采樣數(shù)據(jù)為Sm(t)，以其中一個信號數(shù)據(jù)s(t)為參考，根據(jù)多徑信號的衰落因子，可建立一組多徑信號的數(shù)據(jù)模型

其中，C為多徑信號對應(yīng)衰落因子組成的P×1維矢量(假設(shè)P條多徑信號)

其中，cm(t)= αm(t)e－jφm(t)。

推廣到多個多徑信號組，同時包含獨立信源，建立獨立信號與多徑并存的陣列接收數(shù)據(jù)模型，假設(shè)N個窄帶遠(yuǎn)場信號入射到M元陣列。入射信號中獨立信源數(shù)為Nu，多徑信號數(shù)為Nc，其中N=Nu+Nc。其中多徑信號分成K組相互獨立的相干信源，每一組多徑信源相干信源數(shù)為N，N=∑KN。陣列輸出噪聲為高kck=1k斯白噪聲，且陣元間噪聲相互獨立，其均值與功率分別為0、σ2N，陣列在t時刻接收數(shù)據(jù)為

其中，A為M×N維陣列流型矩陣;S(t)為N維入射信號源;N(t)為M×1維陣元輸出高斯白噪聲;Au為獨立信號對應(yīng)的M×Nu維陣列流型矩陣;Ac為所有多徑信源對應(yīng)的M×Nc維陣列流型矩陣;Su(t)為獨立信號數(shù)據(jù)組成的Nu×1維向量;Sc(t)為K組相互獨立的相干信源組中K個相互獨立的信源數(shù)據(jù)組成的K×1維向量;B為多徑信號的衰落因子組成Nc×K維塊對角矩陣。數(shù)學(xué)模型為

2 基于差分的時間平滑解相干

在移動通信環(huán)境中，對于移動的終端，通常用相干時間tcoh=c(vfc)表示時間選擇性信道衰落變化的快慢。其中，c、v、fc表示電磁波速度、終端移動速度和載波頻率。多徑信號與直射信號變化相對較慢，但由于信號波長短，復(fù)衰落向量C變化較快，尤其是相位。利用時間平滑算法可將接收數(shù)據(jù)協(xié)方差矩陣恢復(fù)為滿秩。假設(shè)接收數(shù)據(jù)快拍數(shù)為L，將所有快拍數(shù)據(jù)分成W(W＞P)個數(shù)據(jù)塊，每個數(shù)據(jù)塊的數(shù)據(jù)長度為D=L/W。對于同一傳播路徑，當(dāng)每個數(shù)據(jù)塊的采樣時間大于信道相干時間tcoh時，這些數(shù)據(jù)塊多徑衰落不相關(guān)［7］。

對于所有傳播路徑及陣列接收數(shù)據(jù)第j個數(shù)據(jù)塊，對陣列接收數(shù)據(jù)按列分塊，可表示為

式(8)中X為陣列接收數(shù)據(jù)，則數(shù)據(jù)塊對應(yīng)的協(xié)方差矩陣為

對W個數(shù)據(jù)塊的協(xié)方差矩陣求平均得到時域平滑的協(xié)方差矩陣為

基于時域平滑的解相干算法不存在陣列孔徑損失，對陣列結(jié)構(gòu)沒有要求，但為了解相干，數(shù)據(jù)塊需要足夠的觀察時間使得相鄰數(shù)據(jù)塊采樣起始時間間隔tcoh，使得不同采樣數(shù)據(jù)塊的衰減因子C不相關(guān)［5］。采用數(shù)據(jù)塊時域平滑算法，只要滿足數(shù)據(jù)塊個數(shù)W＞P，可以得到滿秩的協(xié)方差矩陣 Rts［7］。

針對多徑混合信號，為提高算法性能，同時增大可估計信號數(shù)目，本文提出一種差分時間平滑算法。對于Toeplitz矩陣Rt滿足:JRTtJ=Rt，其中，J 為反對角單位矩陣。式(7)中Ru為對角陣，滿足Toeplitz特性，有JRTuJ=Ru，但Rc不滿足Toeplitz特性。所以可以構(gòu)造多徑相干信號的數(shù)據(jù)協(xié)方差矩陣為

求得協(xié)方差矩陣不含獨立信源分量，與時間平滑算法結(jié)合，可以獲得滿秩的協(xié)方差矩陣，利用傳統(tǒng)子空間算法例如MUSIC算法可對多徑相干信源進(jìn)行有效DOA估計。

根據(jù)以上討論，基于差分方法的時間平滑算法可歸納如下:

(1)由陣列接收數(shù)據(jù)X及快拍數(shù)L，計算陣列協(xié)方差矩陣估計值

(3)根據(jù)式(8)與式(9)求出數(shù)據(jù)塊的協(xié)方差矩陣R^(j)。

(4)利用式(11)求得每個數(shù)據(jù)塊只包含多徑信源協(xié)方差矩陣分量的差分矩陣。

3 仿真分析

仿真1 均勻線陣陣元數(shù)M=8，陣元間距等于半波長。兩個多徑相干信號組成相干信源組，多徑信號波達(dá)方向與衰落系數(shù)分別為{0°，10°}，{1，exp(i2π(2rand(1)－1))}。一個獨立信源入射到天線陣列上，波達(dá)方向為30°，信噪比SNR=10 dB。分別用空間平滑，差分時間平滑算法進(jìn)行DOA估計，時間平滑數(shù)據(jù)塊數(shù)。對多徑信號的估計結(jié)果如圖1所示?？臻g平滑算法在混合信號模型由于天線孔徑損失中性能惡化，可以看出差分方法時間平滑DOA算法有效去除了獨立信源信息，同時估計多徑信號性能良好，增大了陣列天線可估計相干信源數(shù)目。

仿真2 陣列與仿真1相同，在4個多徑信號入射到天線陣列上，波達(dá)方向分別為{－20°，0°，20°，30°}，信噪比為SNR=10 dB。分別用空間平滑，差分時間平滑算法進(jìn)行DOA估計，時間平滑數(shù)據(jù)塊數(shù)W=10。對多徑信號的估計結(jié)果如圖2所示。根據(jù)仿真結(jié)果，可說明該算法對全相干信源模型依然有效，由于差分方法損失了部分多徑相干信號的信息，相對于空間平滑算法，該算法性能有一定損失。

圖1 估計混合信號

圖2 估計相干信號

仿真3 陣列同仿真1，為簡化信號模型，如仿真1混合信號模型，信號入射方向與衰落系數(shù)同仿真1，結(jié)合RootMUSIC進(jìn)行DOA估計，并研究估計多徑信號的性能，研究估計均方誤差RMSE與信噪比SNR關(guān)系。獨立進(jìn)行200次蒙特卡洛試驗的仿真結(jié)果，如圖3所示。

圖3 估計均方誤差

仿真4 陣列、混合信號模型、信號入射方向與衰落系數(shù)均與仿真1相同，結(jié)合RootMUSIC進(jìn)行DOA估計，并研究估計多徑信號的性能，研究估計多徑信號成功率與信噪比的關(guān)系。仿真結(jié)果如圖4所示。

圖4 估計成功概率

4 結(jié)束語

本文介紹了一種移動通信多徑信號的DOA算法，基于多徑信號的衰落特性，利用時間平滑與差分方法，對多徑信號有效解相干。該算法采用時域平滑，不存在孔徑損失，使可估計信源數(shù)增加，分離了獨立信源的影響，提高相干信源DOA估計的性能。理論與仿真分析表明，該算法性能在混合信號模型中優(yōu)于經(jīng)典空間平滑算法，對入射信號全部為相干信源的解相干能力有待改進(jìn)。

［1］Shan T J，Wax M，Kailath T.On spatial smoothing for direction of arrival estimation of coherent signals［J］.IEEE Transactions on ASSP，1985，33(4):806 －811.

［2］Piltai S U，Kwon B H.Forward backward spatial smoothing techniques for the coherent signal identification［J］.IEEE Transactions on Acoustics，Speech，and Signal Processing，1989，37(1):8 －15.

［3］WiIIiams R T，Prasad S，MahaIanabis A K，et aI.An improved spatial smoothing technigue for bearing estimation in a multipath environment［J］.IEEE Transactions on ASSP，1988，36(4):425 －432.

［4］van Rheeden D R，Gupta S C.A temporal smoothing approach to direction of arrival estimation of coherent signals in fading channels［C］.Newyork:Wireless Communications and Networking Conference，IEEE，1999:286 －290.

［5］Gu Z H，Gunawan E.A performance analysis of multipath direction finding with temporal smoothing［J］.IEEE Transactions on Signal Processing Letters，2003，10(7):200 －203.

［6］何子述，韓春林，唐斌.一種有效的多徑信號DOA估計算法［J］.電子學(xué)報，2003，3l(1):139 －141.

［7］張森，陳永革，郭藝奪，等.基于時域數(shù)據(jù)塊平滑的多徑信號DOA估計算法［J］.雷達(dá)科學(xué)與技術(shù)，2010，8(5):429 －421.

［8］Zhang Yufeng，Ye Zhongfu.Efficient method of DOA estimation for uncorrelated and coherent signals［J］.IEEE Antennas and Wireless Propagation Letters，2008(7):799 －802.

［9］Guo Yiduo，Zhang Yongshun，Tong Ningning.A new DOA estimation method for uncorrelated and coherent sources under nonstationary noise fields［C］.Changsha:IEEE 8th International Conference on Application Specific Integrated Circuit，2009:987－990.

［10］王永良，陳輝，彭應(yīng)寧，等.空間譜估計理論與算法［M］.2版.北京:清華大學(xué)出版社，2005.

DOA Estimation of Multipath Signals Based on Differential Temporal Smoothing

PENG Jinlong，CHEN Jianchun，GUO Yaping
(School of Electronic Engineering，Xidian University，Xi'an 710071，China)

An improved DOA algorithm is proposed for the mixed signal model where independent signals and multipath signals coexist in the mobile communication environment based on the differential method and temporal smoothing principle.The algorithm increases the number of signal sources resolved and promotes the DOA estimation performance of coherent signals.Simulation results show that the algorithm has better performance in the mixed signal model than the traditional spatial smoothing method.

temporal smoothing;differential;multipath signals;mobile communication

TN929.5

A

1007－7820(2015)06－162－04

2014-10-29

彭金龍(1989—)，男，碩士研究生。研究方向:基于差分時間平滑的多徑信號DOA估計。E-mail:525545597@qq.com

10.16180/j.cnki.issn1007－7820.2015.06.044