專利名稱：Cdma無線通信的空間特征提取及信號到達方向估計方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種移動通信技術(shù)，更確切地說是涉及一種碼分多址(CDMA)無線通信的空間特征提取及信號到達方向(DOA)估計的方法。該方法不受信號相關(guān)性和信號DOA間距的影響，適于在CDMA移動通信系統(tǒng)中估計多徑信號的空間特征和DOA。所提取的無線信道空間特征，可用于對接收或發(fā)射的信號進行空間濾波，估計的DOA可以用于對移動目標進行定位和波束賦形。
背景技術(shù)：
在移動通信系統(tǒng)中應(yīng)用天線陣列可以有效地抑制干擾，提高系統(tǒng)容量，通過在具有天線陣列的移動通信系統(tǒng)中設(shè)置空間濾波器和/或波束賦形器來有效抑制干擾?？臻g濾波器的設(shè)計方法，大致可以分為“盲”和“非盲”兩大類。其中“非盲”的方法是利用訓(xùn)練序列估計信道的空間特征(或通過估計DOA得到信道的空間特征)來完成的。DOA估計還可以用于對移動目標進行定位和跟蹤。
對于單一的(無擴散的)一條無線信道傳播路徑，當天線陣列結(jié)構(gòu)確定以后，信道的空間特征與信號的DOA是一一對應(yīng)的(映射關(guān)系)，因此可以將信道的空間特征估計和信號的DOA估計作為同一個問題來處理。
估計信號DOA的方法有很多，最典型的有波束掃描(Classic)方法(一種傳統(tǒng)的DOA估計方法)，Capon的最小方差無畸變響應(yīng)(MVDR)方法，和特征結(jié)構(gòu)法(如MUSIC，ESPRIT及其擴展方法)。傳統(tǒng)的波束掃描方法是假設(shè)空間濾波器具有陣列方向矢量的形式，在空間一定的觀察范圍內(nèi)以一定的角度間隔搜索空間濾波器輸出功率最大的那些峰值，那些峰值所對應(yīng)的角度就是要估計的DOA。Capon的MVDR方法的優(yōu)化準則是在保證主瓣能量的約束條件下，使空間濾波器輸出功率最小，從而達到抑制干擾的目的。因此在分辨DOA相近的多個信號時，比傳統(tǒng)的波束掃描方法具有更好的分辨性能。特征結(jié)構(gòu)法通過對以陣列方式輸出的空間協(xié)方差矩陣進行特征值分解，定義信號子空間和噪聲子空間，進而利用信號子空間和噪聲子空間的正交性(對MUSIC及其擴展)或利用信號子空間的旋轉(zhuǎn)不變性(對ESPRIT及其擴展)來獲得對信號DOA的估計。特征結(jié)構(gòu)法在估計多個獨立信號的DOA時，其角度分辨能力大大超過了傳統(tǒng)的Classic算法，因此又被稱為高分辨技術(shù)。
上述的DOA估計算法在分辨多個信號的DOA時要求信號源的數(shù)目小于天線陣列中的天線單元總數(shù)。這一要求使這些算法不能直接用于CDMA通信系統(tǒng)，因為在CDMA通信系統(tǒng)中，信號源的數(shù)目一般大于天線單元總數(shù)。因此現(xiàn)有的CDMA信號的DOA估計方法均基于解調(diào)后的CDMA信號進行。
由于在CDMA移動通信系統(tǒng)中，DOA估計方法基于解調(diào)后的CDMA信號進行。而某個用戶解調(diào)后的信號中必然含有該用戶的多徑信號分量和加性噪聲，且具有如下的特征(1)信號的多徑分量之間存在很強的相關(guān)性；(2)信號的多徑分量之間的DOA間隔(間距)很??；(3)由于CDMA系統(tǒng)中擴頻碼的非正交性，解調(diào)后的噪聲往往是有色噪聲。對于高度相關(guān)的信號，采用傳統(tǒng)Classic算法的DOA估計方法時，會使角分辨能力嚴重下降，采用高分辨方法雖然可以結(jié)合空間平滑技術(shù)使其分辨相關(guān)信號的性能改善，但有損失陣列孔徑和加大算法計算量的不足。在移動通信系統(tǒng)中，多徑分量間的角度往往相隔很近，而受系統(tǒng)成本的影響，陣列天線中的天線單元總數(shù)又不可能很多，因此對于既相關(guān)又角度相近的多徑分量，高分辨的算法也將無能為力。上述分析說明，現(xiàn)有的空間特征提取及DOA估計方法，分辨多徑分量的能力受信號相關(guān)性和DOA間距影響很大。這是因為這些技術(shù)基于噪聲污染的用戶多徑數(shù)據(jù)的空間協(xié)方差矩陣，如果多徑信號相關(guān)，這些相關(guān)能量在分辨多徑時就會影響分辨性能。另一方面，在DOA間距較小時，現(xiàn)有技術(shù)受天線數(shù)及技術(shù)本身性能的影響，使其無法分辨這些DOA間距較小的多徑。
一個解決分辨多徑分量的途徑是利用RAKE接收機結(jié)構(gòu)上的特點進行的，即用一條路徑解調(diào)后的信號估計該徑的DOA(逐徑進行)，但這種方法仍然要受加性有色噪聲的嚴重影響，這是因為現(xiàn)有的技術(shù)是基于加性白噪聲發(fā)展起來的。此外，在分辨多個信號源時，高分辨DOA估計技術(shù)的優(yōu)點主要體現(xiàn)在與傳統(tǒng)的波束掃描技術(shù)相比較，但在計算上要花更大的消費，因為它要進行復(fù)數(shù)矩陣的特征值分解。當用單條路徑解調(diào)后的信號進行DOA估計時，高分辨技術(shù)就失去了它的優(yōu)越性。另外，用RAKE接收機結(jié)構(gòu)的每條徑解調(diào)后的信號進行DOA估計還對接收機有如下的要求(1)接受機必須是RAKE接收形式；(2)有分離每條徑的解調(diào)處理能力。否則，只為DOA估計而實行多徑分量的分離跟蹤和解調(diào)會給接收機帶來額外的復(fù)雜性。

發(fā)明內(nèi)容
本發(fā)明的目的是設(shè)計CDMA無線通信的空間特征提取及信號到達方向估計方法，不受信號相關(guān)性和信號到達方向間距的影響，提取的信道空間特征可以用于對接收和發(fā)射的信號進行空間濾波，估計的到達方向(DOA)可以用于對目標進行定位和波束賦形。
實現(xiàn)本發(fā)明目的的技術(shù)方案是這樣的一種CDMA無線通信的空間特征提取及信號到達方向估計方法，其特征在于包括以下步驟A.利用天線陣列接收的用戶數(shù)據(jù)估計所關(guān)心用戶信道的空時沖激響應(yīng)；B.從估計的該用戶信道的空時沖激響應(yīng)結(jié)果中分離并選取出其中主要傳播路徑的空間沖激響應(yīng)；C.用提取出的主要傳播路徑的空間沖激響應(yīng)，估計每條主要傳播路徑的空間特征以及對應(yīng)的到達方向。
所述步驟B，分離并選取出其中主要傳播路徑的空間沖激響應(yīng)，進一步包括B1.根據(jù)該用戶信道的空時沖激響應(yīng)結(jié)果，獲得其中每條傳播路徑在每個天線單元、每個時延抽頭上的沖激響應(yīng)強度；B2.利用沖激響應(yīng)強度，按沖激響應(yīng)平均功率最大原則確定與主要傳播路徑對應(yīng)的時延抽頭；B3.提取與主要傳播路徑對應(yīng)的時延抽頭的信道空間沖激響應(yīng)。
所述的步驟B2進一步包括求出該用戶對應(yīng)每個天線單元、每個時延抽頭的沖激響應(yīng)平均功率；選取對應(yīng)最大q個平均功率的時延抽頭作為與主要傳播路徑對應(yīng)的時延抽頭。
所述的q是在不同信道環(huán)境下，所存在的滿足一定信噪比門限的不少于q的主要傳播路徑。
所述的步驟B2進一步包括設(shè)定一功率門限值；選取滿足某一條徑對應(yīng)的空間沖激響應(yīng)的功率平均值大于等于功率門限值的所有時延抽頭作為與主要傳播路徑對應(yīng)的時延抽頭。
所述的功率平均值基于天線陣列中的全部天線單元或部分天線單元。
所述的功率門限值基于當前時隙信道空時沖激響應(yīng)數(shù)據(jù)確定，或聯(lián)合考慮當前時隙和當前時隙前面時隙的信道空時沖激響應(yīng)數(shù)據(jù)確定。
所述的聯(lián)合考慮當前時隙和當前時隙前面時隙的信道空時沖激響應(yīng)數(shù)據(jù)，是從當前時隙開始按指數(shù)衰減的原則考慮前面時隙的信道空時沖激響應(yīng)數(shù)據(jù)。
所述的步驟C按θi，a(θi)=argmaxθ∈
]>計算每條主要傳播路徑i的空間特征以及對應(yīng)的到達方向，<，>表示內(nèi)積運算，為陣列方向矢量a(θ)和其中一條主要傳播路徑i對應(yīng)的空間沖激響應(yīng)h(k)(τi)的內(nèi)積，k表示第k個用戶，τi表示與主要傳播路徑i對應(yīng)的時延抽頭。
本發(fā)明的技術(shù)方案，基于用戶信道空時沖激響應(yīng)，估計出每個用戶的信道空時沖激響應(yīng)后，再分離出每條信號主要傳播路徑的信道空間沖激響應(yīng)(有用信息)，只用每個用戶信號主要傳播路徑的信道空間沖激響應(yīng)，而不是利用所有傳播路徑的空間沖激響應(yīng)數(shù)據(jù)，來估計相應(yīng)徑的DOA，不僅避免了白噪聲的影響，還避免了多徑信息的相互干擾。即使多徑信號的相關(guān)性很高，DOA完全相同，本發(fā)明方法也是有效的。而其它現(xiàn)有的方法既不能有效地分辨高度相關(guān)的多徑信號，也不能分辨DOA相近的多徑。
此外，本發(fā)明方法計算簡單，適合工程應(yīng)用尤其是在CDMA系統(tǒng)中應(yīng)用。


圖1是8天線單元的天線陣列和沖激響應(yīng)窗長或抽頭數(shù)為16時的沖激響應(yīng)強度示意圖；圖2是對應(yīng)圖1的沖激響應(yīng)的平均功率示意圖；圖3是應(yīng)用本發(fā)明方法的DOA估計，以1°為步長計算的兩個方向角效果示意。
具體實施例方式
本發(fā)明提出了一種簡單有效的無線信道空間特征提取及信號到達方向(DOA)的估計方法，是基于信道空時沖激響應(yīng)的DOA估計方法，該方法不受信號相關(guān)性和信號DOA間距的影響，適合用在移動通信系統(tǒng)中估計多徑分量的信道空間特征和DOA。提取的信道空間特征可以用于對接收和發(fā)射的信號進行空間濾波，估計的DOA可以用于對目標進行定位和波束賦形。
本發(fā)明主要針對多用戶移動通信系統(tǒng)，尤其是CDMA系統(tǒng)，提出的一種基于空時信道估計的DOA估計方法。以下結(jié)合附圖詳細說明其技術(shù)方案，其主要技術(shù)步驟包括估計用戶信道空時沖激響應(yīng)；分離用戶多徑信號并按沖激響應(yīng)平均功率最大原則選取抽頭；和進行空間特征提取及DOA估計。
結(jié)合參見圖1，對于用戶k，假設(shè)估計后的用戶信道空時沖激響應(yīng)為H(k)＝[h(k)(0)h(k)(1)…h(huán)(k)(W-1)](1)
=h1(k)(0)h1(k)(1)h1(k)(W-1)h2(k)(0)h2(k)(1)h2(k)(W-1)hM(k)(0)hM(k)(1)hM(k)(W-1)----(2)]]>(1)式為矩陣表達式，(2)式是列矢量表達式，式中M為天線數(shù)，W為空時沖激響應(yīng)窗長或抽頭數(shù)。本發(fā)明方法中，可以采用任何現(xiàn)有的估計技術(shù)來估計用戶信道空時沖激響應(yīng)，即對估計用戶信道的空時沖激響應(yīng)的方法沒有限制。沖激響應(yīng)可以是任何形式歸一化后的結(jié)果。
定義對應(yīng)第m個天線單元第i+1個抽頭的沖激響應(yīng)強度為 (強度是絕對值)，圖1示意給出了M＝8和W＝16時，用戶k對應(yīng)每個天線單元、每個抽頭上的沖激響應(yīng)強度，用圓圈示出。經(jīng)過運算，假設(shè)其中兩條多徑的信號h(k)(0)和h(k)(3)(公式2中的第1、第4兩列，即第1條、第2條徑)分別對應(yīng)抽頭τ1＝0和τ2＝3(圖中橫軸表示時延抽頭τ，縱軸表示沖激響應(yīng)強度。
很明顯，此時對第k個用戶，信道的特征可由h(k)(0)和h(k)(3)完全確定。事實上，在每條路徑角度擴散比較小的情況下，我們可以將h(k)(0)和h(k)(3)近似表示成沖激響應(yīng)函數(shù) 式中α0、α1表示信道的幅度衰落系數(shù)，0、1分別是第1條、第2條徑的相位，表示成ej0，ek1，a(θ0)，a(θ1)分別表示第1條、第2條徑的空間特征，δ(τ)和δ(τ-3)沒有運算上的意義。因此經(jīng)過上述運算過程，從H(k)中提取了h(k)(0)和h(k)(3)，就可以估計這兩條路徑的空間特征a(θ0)和a(θ1)以及它們各自對應(yīng)的到達方向角θ0和θ1。
但是，在實際中，我們并不知道信道有幾個主要的傳播路徑(如上述的兩條多徑信號h(k)(0)和h(k)(3))，也不知道這些路徑所對應(yīng)的時延抽頭(如上述的抽頭τ1＝0和τ2＝3)。因此，必須按照某種準則判定并選取出信號的主要傳播路徑(即如何獲得上述兩條多徑信號h(k)(0)和h(k)(3)及分別對應(yīng)時延抽頭τ1＝0和τ2＝3)。可以采用兩種選取方法第一種選取方法是固定的選取q個最大平均功率的時延抽頭，或者說基于選取平均功率最大的q個徑。
對第k個用戶，定義對應(yīng)第i個時延抽頭的沖激響應(yīng)的平均功率為P(k)(i)=||h(k)(i)||2/M,i=0,···,W-1----(4)]]>對應(yīng)圖1，可得到8個天線信道空間沖激響應(yīng)的平均功率，如圖2中圓圈所示。實施時，計算平均功率可基于部分天線單元或全體天線單元進行(即M的取值)。
通過對無線信道進行測量和分析，假設(shè)在不同的信道環(huán)境下總存在不少于q(例如q＝2)個主要傳播路徑滿足一定的信噪比門限，因此，我們在計算了所有抽頭對應(yīng)的沖激響應(yīng)的平均功率后，就可選取對應(yīng)最大q個平均功率的抽頭，提取這些抽頭對應(yīng)的空間沖激響應(yīng)(對應(yīng)時延抽頭τ1＝0和τ2＝3的位置處)。
第二種選取方法是先設(shè)定一個功率門限，選取平均功率超過該門限的所有抽頭，或者說基于設(shè)定的某個門限值，選取平均功率大于該門限的徑。
在移動通信環(huán)境中，環(huán)境的變化和信道的衰落都會帶來多徑數(shù)和多徑強度的變化。方法1選擇固定q條徑就可能出現(xiàn)過估計(q太大)或欠估計(q太小)的情況。一種自適應(yīng)選取q的方法是根據(jù)估計的信道空間沖激響應(yīng)設(shè)定某個門限值 ，然后選取滿足下式(5) 的所有抽頭做為主徑抽頭，即某一條徑對應(yīng)的信道空間沖激響應(yīng)的功率平均值大于等于 的徑，而把低于門限 的徑當噪聲處理。門限可以只用本時隙估計的本用戶的沖激響應(yīng)數(shù)據(jù)來確定，也可以聯(lián)合考慮本時隙和前面時隙的沖激響應(yīng)數(shù)據(jù)來確定，比如從本時隙開始按指數(shù)衰減的原則考慮前面時隙數(shù)據(jù)的作用。
最后作空間特征提取及DOA估計。當我們由H(k)選出了某個信號多徑的空間沖激響應(yīng)h(k)(τi)(如對應(yīng)圖2的h(k)(0)，h(k)(3))后，就可以按下面的規(guī)則估計對應(yīng)路徑的空間特征(以a(θ)為變量)以及對應(yīng)的DOA(以θ為變量)
θi,a(θi)=argmaxθ∈
]>其中<，>表示內(nèi)積運算，為陣列方向矢量a(θ)(陣列方向矢量可以是以任何形式歸一化后的結(jié)果)和其中一條徑對應(yīng)的空間沖激響應(yīng)h(k)(τi)的內(nèi)積。由上式可見，DOA的估計精度由兩個因素確定。其一是h(k)(τi)的估計精度，其二是計算消費函數(shù)fi(θ)＝|<a(θ)，h(k)(τi)>|時，角度搜索步長的選取。
作為驗證舉例，圖3給出了對應(yīng)圖2中的兩條路徑在假設(shè)[θ0，θ3]＝[20°，60°]時，在0到2π范圍內(nèi)，以1°為步長計算(掃描)的f0(θ)＝|<a(θ)，h(k)(0)>|(圖中上半部分)和f1(θ)＝|<a(θ)，h(k)(3)>|的值(圖中下半部分)。圖中顯示出在θ分別為20°和60°時，計算出的f0(θ)和f1(θ)最大，說明了本發(fā)明方法的準確性。
權(quán)利要求
1.一種CDMA無線通信的空間特征提取及信號到達方向估計方法，其特征在于包括以下步驟A.利用天線陣列接收的用戶數(shù)據(jù)估計所關(guān)心用戶信道的空時沖激響應(yīng)；B.從估計的該用戶信道的空時沖激響應(yīng)結(jié)果中分離并選取出其中主要傳播路徑的空間沖激響應(yīng)；C.用提取出的主要傳播路徑的空間沖激響應(yīng)，估計每條主要傳播路徑的空間特征以及對應(yīng)的到達方向。
2.根據(jù)權(quán)利要求1所述的一種CDMA無線通信的空間特征提取及信號到達方向估計方法，其特征在于所述步驟B，分離并選取出其中主要傳播路徑的空間沖激響應(yīng)，進一步包括B1.根據(jù)該用戶信道的空時沖激響應(yīng)結(jié)果，獲得其中每條傳播路徑在每個天線單元、每個時延抽頭上的沖激響應(yīng)強度；B2.利用沖激響應(yīng)強度，按沖激響應(yīng)平均功率最大原則確定與主要傳播路徑對應(yīng)的時延抽頭；B3.提取與主要傳播路徑對應(yīng)的時延抽頭的信道空間沖激響應(yīng)。
3.根據(jù)權(quán)利要求2所述的一種CDMA無線通信的空間特征提取及信號到達方向估計方法，其特征在于所述的步驟B2進一步包括求出該用戶對應(yīng)每個天線單元、每個時延抽頭的沖激響應(yīng)平均功率；選取對應(yīng)最大q個平均功率的時延抽頭作為與主要傳播路徑對應(yīng)的時延抽頭。
4.根據(jù)權(quán)利要求3所述的一種CDMA無線通信的空間特征提取及信號到達方向估計方法，其特征在于所述的q是在不同信道環(huán)境下，所存在的滿足一定信噪比門限的不少于q的主要傳播路徑。
5.根據(jù)權(quán)利要求2所述的一種CDMA無線通信的空間特征提取及信號到達方向估計方法，其特征在于所述的步驟B2進一步包括設(shè)定一功率門限值；選取滿足某一條徑對應(yīng)的空間沖激響應(yīng)的功率平均值大于等于功率門限值的所有時延抽頭作為與主要傳播路徑對應(yīng)的時延抽頭。
6.根據(jù)權(quán)利要求5所述的一種CDMA無線通信的空間特征提取及信號到達方向估計方法，其特征在于所述的功率門限值基于當前時隙信道空時沖激響應(yīng)數(shù)據(jù)確定，或聯(lián)合考慮當前時隙和當前時隙前面時隙的信道空時沖激響應(yīng)數(shù)據(jù)確定。
7.根據(jù)權(quán)利要求6所述的一種CDMA無線通信的空間特征提取及信號到達方向估計方法，其特征在于所述的聯(lián)合考慮當前時隙和當前時隙前面時隙的信道空時沖激響應(yīng)數(shù)據(jù)，是從當前時隙開始按指數(shù)衰減的原則考慮前面時隙的信道空時沖激響應(yīng)數(shù)據(jù)。
8.根據(jù)權(quán)利要求5所述的一種CDMA無線通信的空間特征提取及信號到達方向估計方法，其特征在于所述的功率平均值基于天線陣列中的全部天線單元或部分天線單元。
9.根據(jù)權(quán)利要求1所述的一種CDMA無線通信的空間特征提取及信號到達方向估計方法，其特征在于所述的步驟C按θi，a(θi)=argmaxθ∈
]>計算每條主要傳播路徑i的空間特征以及對應(yīng)的到達方向，<，>表示內(nèi)積運算，為陣列方向矢量a(θ)和其中一條主要傳播路徑i對應(yīng)的空間沖激響應(yīng)h(k)(τi)的內(nèi)積，k表示第k個用戶，τi表示與主要傳播路徑i對應(yīng)的時延抽頭。
全文摘要
本發(fā)明涉及一種CDMA無線通信的空間特征提取及信號到達方向估計方法，可不受信號相關(guān)性和信號到達方向間距的影響。包括以下步驟A.利用天線陣列接收的數(shù)據(jù)估計所關(guān)心用戶信道的空時沖激響應(yīng)；B.從估計的該用戶信道的空時沖激響應(yīng)結(jié)果中分離并選取出其中主要傳播路徑的空間沖激響應(yīng)；C.用提取出的主要傳播路徑的空間沖激響應(yīng)，估計每條主要傳播路徑的空間特征以及對應(yīng)的到達方向。由于只用用戶信號主要傳播路徑的信道空間沖激響應(yīng)，而不是利用所有傳播路徑的空間沖激響應(yīng)數(shù)據(jù)，來估計相應(yīng)徑的DOA，不僅避免了白噪聲的影響，還避免了多徑信息的相互干擾。即使多徑信號的相關(guān)性很高，DOA完全相同，本發(fā)明方法也是有效的。
文檔編號H04L25/02GK1585318SQ03153929
公開日2005年2月23日 申請日期2003年8月20日 優(yōu)先權(quán)日2003年8月20日
發(fā)明者李平安, 楊貴亮, 康紹莉, 胡兵 申請人:大唐移動通信設(shè)備有限公司