專利名稱：一種移動(dòng)通訊系統(tǒng)來波方向的高分辨率估計(jì)方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種TD-SCDMA系統(tǒng)超多DOA的高分辨率估計(jì)的方法或裝置，本方法同時(shí)適用于其它CDMA通信領(lǐng)域。
背景技術(shù)：
在無線通訊系統(tǒng)中，往往需要估計(jì)移動(dòng)用戶的來波方向(DOA)，以滿足各種不同的應(yīng)用需要，如移動(dòng)用戶定位，波束賦形，等等，DOA估計(jì)一般是通過對(duì)多陣元的天線陣列的接收信號(hào)進(jìn)行聯(lián)合處理來實(shí)現(xiàn)的。傳統(tǒng)的DOA估計(jì)技術(shù)，如MUSIC算法，ESPRIT算法和以MUSIC算法為基礎(chǔ)改進(jìn)的Root-MUSIC等，所能估計(jì)的DOA數(shù)受天線陣列的陣元數(shù)的限制，M根陣元只能估計(jì)M-1個(gè)DOA。而在實(shí)際中，考慮成本與復(fù)雜度，基站所采用的天線陣列一般為4-8根陣元，也就是說，在實(shí)際的8陣元的天線陣列系統(tǒng)中，傳統(tǒng)的DOA估計(jì)方法只能估計(jì)最多7個(gè)DOA。這是不能滿足需要的，因?yàn)樵跓o線通信環(huán)境中，一個(gè)基站往往需要同時(shí)與十幾到幾十個(gè)用戶同時(shí)通信，且這些用戶信號(hào)的傳播環(huán)境還存在多徑，需要估計(jì)的DOA個(gè)數(shù)常常多于7個(gè)?；贛USIC算法改進(jìn)的循環(huán)MUSIC算法根據(jù)譜相關(guān)性，能估計(jì)的DOA數(shù)不受陣元數(shù)的限制，但要求的數(shù)據(jù)量較大，計(jì)算復(fù)雜。
因此，現(xiàn)有技術(shù)存有缺陷，而有待于改進(jìn)和發(fā)展。

發(fā)明內(nèi)容
本發(fā)明的目的在于提供一種移動(dòng)通訊系統(tǒng)來波方向的高分辨率估計(jì)方法，能夠高分辨率地估計(jì)出各用戶的多徑來波方向，且最多能估計(jì)的DOA數(shù)數(shù)倍于天線陣列的陣元數(shù)，為(M-1)K個(gè)，其中M是天線陣列的陣元數(shù)，K是用戶個(gè)數(shù)，而且本發(fā)明提出的方法計(jì)算量少，易于工程實(shí)現(xiàn)。
本發(fā)明的技術(shù)方案為一種移動(dòng)通訊系統(tǒng)來波方向的高分辨率估計(jì)方法，其包括以下步驟估計(jì)出各天線陣元上每個(gè)用戶的信道信息，然后利用各天線陣元上每個(gè)用戶的信道信息構(gòu)成各個(gè)用戶的信號(hào)協(xié)方差矩陣，再對(duì)每個(gè)用戶的信號(hào)協(xié)方差矩陣進(jìn)行特征分解，分離信號(hào)子空間和噪聲子空間，然后在噪聲子空間搜索來波方向，將每個(gè)用戶的多徑方向準(zhǔn)確搜索出來。
所述的方法，其中，所述方法還包括以下步驟a)估計(jì)用戶信道響應(yīng)值利用每個(gè)時(shí)隙中已知的Midamble碼構(gòu)成G矩陣，估計(jì)所有用戶在M元陣元上的信道響應(yīng)值，第m根陣元上的信道響應(yīng)值為h^(m)=(GHG)-1GH·r(m)----(4);]]>b)根據(jù)用戶k的midamble碼移位取得該相應(yīng)位置上的信道矩陣；c)求取用戶k的信道矩陣的相關(guān)矩陣Rhh(k)=H(k)·H(k)H=h(k,1)h(k,2)···h(k,M)h(k,1)*h(k,2)*···h(k,M)*;]]>=h(k,1)h(k,1)*h(k,1)h(k,2)*···h(k,1)h(k,M)*h(k,2)h(k,1)*h(k,2)h(k,2)*···h(k,2)h(k,M)*············h(k,M)h(k,1)*h(k,M)h(k,2)*···h(k,M)h(k,M)*----(7)]]>d)對(duì)該相關(guān)矩陣進(jìn)行特征分解；e)根據(jù)特征值的大小排序，取最小的M-Lk個(gè)特征值對(duì)應(yīng)的特征向量組成噪聲空間Vn(k)=g1(k)g2(k)···gM-Lk(k)----(9);]]>f)用戶k的第l徑的方向向量a(θl(k))，l＝1∶Lk，可表示為a(θl(k))=ej2πrλcosθl(k)ej2πrλcos(θl(k)-2π(ka-1)Ka)···ej2πrλcos(θl(k)-2π(Ka-1)Ka)T---(10)]]>由于用戶k的第l徑屬于信號(hào)空間，因此有a(θl(k))HVnVnHa(θl(k))=0]]>于是，用戶k的Lk個(gè)入射信號(hào)的來波方向可通過確定MUSIC空間譜的峰值而作出估計(jì)，具體如下令P(θl(k))=1al(k)HVn(k)Vn(k)Hal(k)----(11)]]>或P(θl(k))=al(k)Hal(k)al(k)HVn(k)Vn(k)Hal(k)---(12)]]>令θ從0到360度取值，計(jì)算P(θ)值，取P(θ)的Lk個(gè)峰值對(duì)應(yīng)的Lk個(gè)θ值為θl(k)，l＝1…Lk；g)重復(fù)步驟b)-f)，估計(jì)出所有用戶(k＝1…K)的來波方向。
本發(fā)明所提供的一種移動(dòng)通訊系統(tǒng)來波方向的高分辨率估計(jì)方法，結(jié)合傳統(tǒng)MUSIC的思想，利用各天線陣元上每個(gè)用戶的信道信息構(gòu)成各個(gè)用戶的信號(hào)協(xié)方差矩陣，再對(duì)每個(gè)用戶的信號(hào)協(xié)方差矩陣進(jìn)行特征分解，分離信號(hào)子空間和噪聲子空間，然后在噪聲子空間搜索來波方向，估計(jì)的DOA數(shù)遠(yuǎn)遠(yuǎn)突破天線陣列陣元數(shù)的限制，并且計(jì)算量少，易于工程實(shí)現(xiàn)。


圖1、圖2是本發(fā)明的移動(dòng)通訊系統(tǒng)來波方向的高分辨率估計(jì)方法的效果示意圖。
具體實(shí)施例方式
下面以實(shí)施例子具體說明本發(fā)明方法。
本發(fā)明方法的核心思想為，首先估計(jì)出各天線陣元上每個(gè)用戶的信道信息，然后結(jié)合傳統(tǒng)MUSIC算法的思想，利用各天線陣元上每個(gè)用戶的信道信息構(gòu)成各個(gè)用戶的信號(hào)協(xié)方差矩陣，再對(duì)每個(gè)用戶的信號(hào)協(xié)方差矩陣進(jìn)行特征分解，分離信號(hào)子空間和噪聲子空間，然后在噪聲子空間搜索來波方向，將每個(gè)用戶的多徑方向準(zhǔn)確搜索出來。
下面是對(duì)本發(fā)明方法做進(jìn)一步詳細(xì)描述。
對(duì)一個(gè)離散輸入輸出的CDMA系統(tǒng)，設(shè)系統(tǒng)共有K個(gè)用戶接入，第k個(gè)用戶所發(fā)送的信號(hào)可表示為s(k)(t)=Σn=1Ndn(k)·Σq=1Qcq(k)·gc{t-[(n-1)Q+q]Tc}----(1)]]>其中，dn(K)表示用戶發(fā)送的數(shù)據(jù)符號(hào)；N為用戶發(fā)送的符號(hào)數(shù)，Q為擴(kuò)頻系數(shù)，Tc為碼片持續(xù)時(shí)間，gc(τ)為碼片沖激響應(yīng)，cq(k)為第k個(gè)用戶的擴(kuò)頻碼片值。
設(shè)基站使用M元的天線陣列(均勻線陣、圓陣或其它流形)，則第m個(gè)陣元上的接收數(shù)據(jù)為x(m)(t)=Σk=1KΣl=1Lkal(k,m)sk(t)·hl(k,m)(τ,t)+n(t)----(2)]]>其中，Lk是第k個(gè)用戶信號(hào)到達(dá)接收端的多徑數(shù)；hl(k，m)(τ，t)是第k個(gè)用戶的第l徑到達(dá)天線第m陣元的等效信道沖激響應(yīng)，它包含了空中衰落信道和接收端信道濾波器的影響，含有幅度的瑞利衰落、多徑時(shí)延、因移動(dòng)臺(tái)高速移動(dòng)引起的多譜勒頻移、角度擴(kuò)展等信息，是時(shí)變信道；al(k，m)是第k個(gè)用戶的第l徑到達(dá)天線第m陣元的陣列因子，它主要與陣列間距和信號(hào)的波達(dá)角有關(guān)；n(t)為加性噪聲。
則M根天線的接收信號(hào)合成一個(gè)矩陣為X＝[x(1)(t) x(2)(t) ... x(M)(t)]T(3)其中，T為轉(zhuǎn)置。
根據(jù)TD-SCDMA物理層協(xié)議，一個(gè)時(shí)隙中的K個(gè)用戶的Midamble碼由一基本的Midamble碼移位形成，移位的窗長(zhǎng)度W＝[128/K]，這個(gè)W值對(duì)于基站和用戶都是事先已知的。根據(jù)Steiner信道估計(jì)方法，利用每個(gè)時(shí)隙中已知的Midamble碼構(gòu)成G矩陣，估計(jì)所有用戶在第m根陣元上的信道響應(yīng)值h^(m)=(GHG)-1GH·r(m)----(4)]]>其中r(m)是接收信號(hào)x(m)中相應(yīng)于Midamble碼的部分。
聯(lián)合M根陣元上的信道響應(yīng) 每個(gè)用戶根據(jù)自己的Midamble碼的移位，在 中有相應(yīng)的M×W的信道值，如式(5)的標(biāo)記所示。第k個(gè)用戶的信道矩陣H(k)=h1(k,1)···hW(k,1)h1(k,2)···hW(k,2)·········h1(k,M)···hW(k,M)----(6)]]>由式(4)估計(jì)的信道響應(yīng)值，不僅包含時(shí)延和衰落信息，還包含角度信息，是式(2)中的hl(k，m)(τ，t)和al(k，m)兩項(xiàng)的乘積。求用戶k的信道相關(guān)矩陣Rhh(k)=H(k)·H(k)H=h(k,1)h(k,2)···h(k,M)h(k,1)*h(k,2)*···h(k,M)*]]>=h(k,1)h(k,1)*h(k,1)h(k,2)*···h(k,1)h(k,M)*h(k.2)h(k,1)*h(k,2)h(k,2)*···h(k,2)h(k,M)*············h(k,M)h(k,1)*h(k,M)h(k,2)*···h(k,M)h(k,M)*----(7)]]>容易驗(yàn)證，協(xié)方差矩陣Rhh(k)為Hermitian矩陣，它的特征值為非負(fù)值。令特征值λ1≥λ2≥...≥λM≥0，其特征分解值可寫成Rhh(k)=UΣUH=Σi=1MλiuiuiH.]]>式中，U＝[u1，u2...，uM]為特征向量組成的酉矩陣(即每個(gè)向量的模為1，向量之間相互正交，且UUH＝I，是復(fù)數(shù)的正交陣)?！疲絛iag[λ1，λ2，...，λM]為特征值構(gòu)成的對(duì)角矩陣。若將M個(gè)特征值按大小依次排列，則前Lk(Lk≤W)個(gè)與用戶k的多徑功率有關(guān)。從第Lk+1開始的小特征值完全確定于噪聲，其值約等于σ2，即λLk+1≈λLk+2≈...≈λM≈σ2.]]>假設(shè)各用戶有效的多徑數(shù)Lk小于陣元數(shù)M(這是非常合理的假設(shè)，一般多徑數(shù)在2-4間，陣元數(shù)在4-8間)，Rhh(k)有Lk個(gè)大特征值對(duì)應(yīng)多徑信號(hào)數(shù)，有M-Lk個(gè)最小特征值，對(duì)應(yīng)于噪聲功率。將與Lk個(gè)大特征值對(duì)應(yīng)的特征向量ul(k)，l＝1，...，Lk，它們構(gòu)成多徑信號(hào)空間，與M-Lk個(gè)最小特征值對(duì)應(yīng)的特征向量gj(k)，j＝1∶M-Lk，它們是噪聲矢量，構(gòu)成噪聲空間。多徑信號(hào)空間與噪聲空間正交，也就是說u1(k)u2(k)···uLk(k)⊥g1(k)g2(k)···gM-Lk(k)----(8)]]>構(gòu)造噪聲矩陣Vn(k)=g1(k)g2(k)···gM-Lk(k)----(9)]]>根據(jù)(8)式，我們可得ul(k)HVnVnHul(k)=0.]]>用戶k的第l徑的方向向量a(θl(k))，l＝1∶Lk可表示為
a(θl(k))=ej2πrλcosθl(k)ej2πrλcos(θl(k)-2π(ka-1)Ka)···ej2πrλcos(θl(k)-2π(Ka-1)Ka)T----(10)]]>由于用戶k的第l徑屬于信號(hào)空間，因此有a(θl(k))HVnVnHa(θl(k))=0]]>于是，用戶k的Lk個(gè)入射信號(hào)的DOA可通過確定MUSIC空間譜的峰值而作出估計(jì)。
具體如下令P(θl(k))=1al(k)HVn(k)Vn(k)Hal(k)---(11)]]>或P(θl(k))=al(k)Hal(k)al(k)HVn(k)Vn(k)Hal(k)---(12)]]>令θ從0到360度取值，計(jì)算P(θ)值，取P(θ)的Lk個(gè)峰值對(duì)應(yīng)的Lk個(gè)θ值為θl(k)，l＝1∶Lk。
一般來說，與Lk個(gè)大特征值對(duì)應(yīng)的特征向量al(k)=a(θl(k)),]]>l＝1∶Lk，為導(dǎo)向矢量，可以將它們直接作為下行的加權(quán)矢量用于波束形成。但在仿真中發(fā)現(xiàn)，此方法形成的波束有一個(gè)缺點(diǎn)，就是求得的每個(gè)權(quán)矢量并沒有嚴(yán)格地與來波方向一一對(duì)應(yīng)，它們?cè)诔俗约悍较虻钠渌较蛞灿幸欢ǚ鹊仨憫?yīng)，這樣使得形成的波束用于發(fā)射時(shí)，則把目標(biāo)信號(hào)同時(shí)發(fā)射到其它徑向上和其它用戶方向上，形成干擾。另外，在宏小區(qū)的環(huán)境中，角度擴(kuò)展一般為45°，使得同一用戶的兩個(gè)多徑方向可能很相近，此時(shí)形成的波束則不能兩徑信號(hào)分開。如果用本發(fā)明的方法高精度地估計(jì)出各徑的來波方向DOA，采用此方向的陣列響應(yīng)矢量作為加權(quán)矢量，則會(huì)很大程度地避免上述問題。
本發(fā)明方法的具體操作步驟如下1.用Steiner方法估計(jì)用戶信道響應(yīng)值，見(4)式；2.根據(jù)用戶k的midamble碼移位取得相應(yīng)位置上的信道矩陣；3.求取用戶k信道矩陣的相關(guān)矩陣，見(7)式；4.對(duì)相關(guān)矩陣特征分解；5.根據(jù)特征值的大小排序，取最小的M-Lk個(gè)特征值對(duì)應(yīng)的特征向量組成噪聲空間(9)式；6.取θ等于0到360度，計(jì)算a(θ)和P(θ)，見(10)式和(11)式，取P(θ)的Lk個(gè)峰值對(duì)應(yīng)的Lk個(gè)θ值為用戶k的Lk個(gè)入射信號(hào)的DOA估計(jì)，記為θl(k)，l＝1∶Lk；重復(fù)步驟2到6，得出所有用戶的DOA估計(jì)。
假設(shè)TD-SCDMA系統(tǒng)同時(shí)通信的用戶數(shù)為5，每用戶有4個(gè)多徑，各徑能量分配為
dB，這樣總的來波數(shù)＝5*4＝20個(gè)，8陣元圓陣天線接收，陣列半徑0.5λ，用戶來波方向隨機(jī)產(chǎn)生，角度擴(kuò)展范圍45，用戶移動(dòng)速度0km/h，輸入信噪比Eb/NO為10dB。
則用戶DOA的估計(jì)過程如下1分別用8陣元數(shù)據(jù)估計(jì)信道沖激響應(yīng)；2根據(jù)第1用戶的midamble碼位移對(duì)應(yīng)的位置，取得相應(yīng)的信道矩陣。
3求信道矩陣的相關(guān)矩陣。
4對(duì)相關(guān)矩陣特征分解。
5根據(jù)特征值的大小排序，將特征向量分成信號(hào)空間和噪聲空間。
6利用噪聲空間矩陣在360度空間搜索DOA，得到第1用戶多徑來波方向，如圖1所示。
7第2用戶至第5用戶的DOA估計(jì)重復(fù)第2步到第6步。
如圖1和圖2所示，分別顯示了其中第1用戶和第2用戶的所有4徑的來波方向，圖上注明了實(shí)際的DOA以便與結(jié)果對(duì)比。其它3個(gè)用戶的結(jié)果相似，就不再贅述。從仿真結(jié)果可以看出本發(fā)明能有效估計(jì)全部的5個(gè)用戶的20個(gè)DOA。
應(yīng)當(dāng)理解的是，上述針對(duì)具體實(shí)施例的描述過于具體，不能因此而限制專利請(qǐng)求保護(hù)范圍，專利請(qǐng)求保護(hù)范圍應(yīng)以所附權(quán)利要求為準(zhǔn)。
權(quán)利要求
1.一種移動(dòng)通訊系統(tǒng)來波方向的高分辨率估計(jì)方法，其包括以下步驟估計(jì)出各天線陣元上每個(gè)用戶的信道信息，然后利用各天線陣元上每個(gè)用戶的信道信息構(gòu)成各個(gè)用戶的信號(hào)協(xié)方差矩陣，再對(duì)每個(gè)用戶的信號(hào)協(xié)方差矩陣進(jìn)行特征分解，分離信號(hào)子空間和噪聲子空間，然后在噪聲子空間搜索來波方向，將每個(gè)用戶的多徑方向準(zhǔn)確搜索出來。
2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述方法還包括以下步驟a)估計(jì)用戶信道響應(yīng)值利用每個(gè)時(shí)隙中已知的Midamble碼構(gòu)成G矩陣，估計(jì)所有用戶在M元陣元上的信道響應(yīng)值，第m根陣元上的信道響應(yīng)值為h^(m)=(GHG)-1GH·r(m)---(4);]]>b)根據(jù)用戶k的midamble碼移位取得該相應(yīng)位置上的信道矩陣；c)求取用戶k的信道矩陣的相關(guān)矩陣Rhh(k)=H(k)·H(k)H=h(k,1)h(k,2)···h(k,M)[h(k,1)*h(k,2)*···h(k,M)*]]]>=h(k,1)h(k,1)*h(k,1)h(k,2)*···h(k,1)h(k,M)*h(k,2)h(k,1)*h(k,2)h(k,2)*···h(k,2)h(k,M)*············h(k,M)h(k,1)*h(k,M)h(k,2)*···h(k,M)h(k,M)*---(7)]]>d)對(duì)該相關(guān)矩陣進(jìn)行特征分解；e)根據(jù)特征值的大小排序，取最小的M-Lk個(gè)特征值對(duì)應(yīng)的特征向量組成噪聲空間vn(k)=g1(k)g2(k)···gM-Lk(k)---(9);]]>f)用戶k的第1徑的方向向量a(θl(k)),l=1:Lk,]]>可表示為a(θl(k))=ej2πrλcosθl(k)ej2πrλcos(θl(k)-2π(ka-1)Ka)···ej2πrλcos(θl(k)-2π(Ka-1)Ka)T---(10)]]>由于用戶k的第1徑屬于信號(hào)空間，因此有a(θl(k))HVnVnHa(θl(k))=0]]>于是，用戶k的Lk個(gè)入射信號(hào)的來波方向可通過確定MUSIC空間譜的峰值而作出估計(jì)，具體如下令P(θl(k))=1al(k)HVn(k)Vn(k)Hal(k)---(11)]]>或P(θl(k))=al(k)Hal(k)al(k)HVn(k)Vn(k)Hal(k)---(12)]]>令θ從0到360度取值，計(jì)算P(θ)值，取P(θ)的Lk個(gè)峰值對(duì)應(yīng)的Lk個(gè)θ值為θl(k)，l＝1...Lk；g)重復(fù)步驟b)-f)，估計(jì)出所有用戶(k＝1...K)的來波方向。
全文摘要
本發(fā)明的一種移動(dòng)通訊系統(tǒng)來波方向的高分辨率估計(jì)方法，其包括以下步驟估計(jì)出各天線陣元上每個(gè)用戶的信道信息，然后利用各天線陣元上每個(gè)用戶的信道信息構(gòu)成各個(gè)用戶的信號(hào)協(xié)方差矩陣，再對(duì)每個(gè)用戶的信號(hào)協(xié)方差矩陣進(jìn)行特征分解，分離信號(hào)子空間和噪聲子空間，然后在噪聲子空間搜索來波方向，將每個(gè)用戶的多徑方向準(zhǔn)確搜索出來。本發(fā)明方法估計(jì)的DOA數(shù)遠(yuǎn)遠(yuǎn)突破天線陣列陣元數(shù)的限制，并且計(jì)算量少，易于工程實(shí)現(xiàn)。
文檔編號(hào)H04J13/00GK1595190SQ20041002799
公開日2005年3月16日 申請(qǐng)日期2004年7月6日 優(yōu)先權(quán)日2004年7月6日
發(fā)明者周用芳 申請(qǐng)人:中興通訊股份有限公司