摘 要:目前TD-SCDMA系統(tǒng)多采用智能天線和上行同步技術(shù)對(duì)移動(dòng)用戶的位置信息進(jìn)行聯(lián)合估計(jì)?，F(xiàn)有的智能天線技術(shù)在測(cè)角過程中，空間譜搜索計(jì)算量大且用戶定位容量受到天線陣元數(shù)目的限制。針對(duì)上述缺點(diǎn)，提出一種用戶數(shù)目不受限制的快速定位算法。新算法主要利用Midamble碼的非相關(guān)特性對(duì)角度估計(jì)算法進(jìn)行了改進(jìn)。仿真實(shí)驗(yàn)證明了該算法的有效性。

關(guān)鍵詞:時(shí)分同步碼分多址; 定位; 波達(dá)方向; 智能天線

中圖分類號(hào):TN914 文獻(xiàn)標(biāo)識(shí)碼:A

文章編號(hào):1004-373X(2010)07-0001-05

Research on Rapid Non-limited Wireless Location Algorithm in TD-SCDMA

JING Yang1，2， DING Qian-jun2， CHEN Hui2， WANG Fei-xue1

(1. Satellite Navigation RD Center National University of Defense Technology， Changsha 410073， China;

2. Communication Institute， Air Force Armament Research Academy， Beijing 100085， China)

Abstract: The estimation of the users′ location is based on smart antenna and uplink synchronization in TD-SCDMA. Considering the large computation and users number limitation of spatial spectrum searching， the existing smart antenna technology does not adapt to the requirement. Therefore， a new non-limited algorithm with less computational complexity is proposed. The direction vector of the users can be estimated by decorrelation of the midambles of users， and then the Direction of Arrival (DOA) of the users can be computed owing to the receiving array manifold. The excellent performance of the new technique is verified through experimental simulation.

Keywords: TD-SCDMA; location; direction of arrival; intelligent antenna

0 引 言

TD-SCDMA是我國(guó)自主研發(fā)的第三代移動(dòng)通信標(biāo)準(zhǔn)之一，其本身有自己的特點(diǎn)和優(yōu)勢(shì)。移動(dòng)用戶的定位作為當(dāng)前TD-SCDMA數(shù)據(jù)服務(wù)中的關(guān)鍵性業(yè)務(wù)，將成為未來3G市場(chǎng)的新利潤(rùn)增長(zhǎng)點(diǎn)。移動(dòng)定位業(yè)務(wù)是指運(yùn)營(yíng)商向用戶提供滿足一定服務(wù)質(zhì)量、一定私密性的移動(dòng)用戶位置信息。運(yùn)營(yíng)商、服務(wù)提供商或安全部門可利用位置信息向公眾提供各種增值服務(wù)和安全保障。常用的定位業(yè)務(wù)有導(dǎo)航、交通調(diào)度、網(wǎng)絡(luò)選擇增強(qiáng)、緊急救援等。3GPP也正加速TD-SCDMA的移動(dòng)業(yè)務(wù)標(biāo)準(zhǔn)化進(jìn)程[1-3]。目前，TD-SCDMA定位技術(shù)根據(jù)測(cè)量實(shí)體的不同，可分為基于網(wǎng)絡(luò)定位、基于移動(dòng)用戶定位、TD-SCDMA特有定位方式[4，5]。

(1) 基于網(wǎng)絡(luò)的定位技術(shù) (Observed Time Difference of Arrival，OTDOA)是一種典型的基于網(wǎng)絡(luò)定位方式?；驹砣缦?移動(dòng)用戶不斷測(cè)量多個(gè)基站(三個(gè)以上)的導(dǎo)頻信號(hào)，得到不同基站的到達(dá)時(shí)間差;根據(jù)測(cè)量結(jié)果及基站的位置坐標(biāo)就可以解算出用戶的位置信息。該算法雖不需要增加定位設(shè)備，但定位精度較差，很難滿足3G網(wǎng)絡(luò)中用戶的高精度定位需求。

(2) 基于移動(dòng)用戶定位，主要是利用手機(jī)上的GPS測(cè)量信息而實(shí)現(xiàn)的?；驹硎?移動(dòng)用戶集成GPS接收芯片，直接接收衛(wèi)星導(dǎo)航信號(hào)，解算出用戶位置信息。此外可以將內(nèi)置GPS的部分接收功能轉(zhuǎn)移到小區(qū)基站上進(jìn)行，例如，基站根據(jù)自身位置獲取當(dāng)前可見衛(wèi)星清單、衛(wèi)星多普勒參數(shù)等輔助信息，通過網(wǎng)絡(luò)發(fā)送給小區(qū)內(nèi)的移動(dòng)用戶。這些數(shù)據(jù)縮短用戶終端內(nèi)置GPS的定位時(shí)間，移動(dòng)用戶將自身位置發(fā)送給中心網(wǎng)絡(luò)。該方法定位精度高，但也有易受遮擋干擾，耗電量較大等明顯的缺點(diǎn)，同時(shí)還需要更新硬件設(shè)備。

(3) TD-SCDMA特有定位方式[6]，主要是利用智能天線、上行同步等一系列關(guān)鍵技術(shù)，對(duì)移動(dòng)用戶的方位角和距離進(jìn)行聯(lián)合估計(jì)，最終解算出位置信息。智能天線技術(shù)多采用陣列信號(hào)處理中的空間譜估計(jì)方法來估計(jì)用戶到基站的方位角;上行同步技術(shù)通過往返定時(shí)同步的方法得到用戶到基站的距離值，如圖1所示。

圖1 上行同步過程

上行同步過程:假設(shè)移動(dòng)用戶在延期了ts時(shí)間后收到基站發(fā)送的子幀，并從這個(gè)子幀的Dw-PTS知道了Up-PTS的位置。為了確保各個(gè)移動(dòng)用戶同時(shí)到達(dá)基站進(jìn)行處理，Up-PTS必須提前發(fā)出，提前量由上次測(cè)得的ts時(shí)延決定，由此可知，Up-PTS上行發(fā)送的時(shí)間為t1-2ts。顯而易見，在估計(jì)出基站到移動(dòng)用戶的時(shí)間為ts之后，由光速在空氣中傳播的速度即可得出基站到移動(dòng)用戶的距離。

TD-SCDMA特有定位方法原理簡(jiǎn)明，定位精度高，可以克服OTDOA多基站在偏僻區(qū)域很難得到滿足的困難，也克服了移動(dòng)用戶定位方法需增添GPS接收設(shè)備的困難。當(dāng)然，TD-SCDMA傳統(tǒng)的特有定位方式也具有一定的缺點(diǎn)，如:智能天線中的空間譜搜索過程計(jì)算量大，而解算的工作基本是由小區(qū)基站完成的，這就加大了基站處理的負(fù)荷;還有傳統(tǒng)測(cè)角技術(shù)使得智能天線中的陣元數(shù)目限制了待估計(jì)用戶數(shù)目。

針對(duì)上述不足，本文在保留原有測(cè)距算法的基礎(chǔ)上，對(duì)傳統(tǒng)的測(cè)角算法進(jìn)行了改進(jìn)，提出一種非受限單基站快速定位算法。實(shí)驗(yàn)表明，該算法在與傳統(tǒng)測(cè)向算法的測(cè)向精度保持一致的同時(shí)，大大減少了計(jì)算量，并且估計(jì)的用戶數(shù)目不受智能天線信號(hào)源數(shù)目的限制，因此，新算法是一種較理想的定位方式。

1 傳統(tǒng)的角度估計(jì)算法[1]

在智能天線處理中，設(shè)N個(gè)移動(dòng)信號(hào)源(波長(zhǎng)為λ)入射到空間M元測(cè)向均勻圓陣上，每個(gè)陣元天線形成一個(gè)接收通道，則陣列接收數(shù)據(jù)模型的向量表達(dá)為:

X(t)=A*S(t)+N(t)(1)

A=[a(θ1，φ1)，a(θ2，φ2)，…，a(θi，φi)，…，a(θN，φN)](2)

式中:X(t)為M×1維快拍數(shù)據(jù)矢量;N(t)為M×1維噪聲數(shù)據(jù)矢量;S(t)為空間信號(hào)快拍數(shù)據(jù)矢量;a為空間陣列的M×1維導(dǎo)向矢量;第i個(gè)移動(dòng)信號(hào)源的入射方位角為θi，俯仰角為φi。

如圖2所示，數(shù)字代表陣元天線的編號(hào)，這里用k來表示，k=1，2，…，M。對(duì)于均勻圓陣，設(shè)圓陣半徑為r，以圓心為參考點(diǎn)，則第i個(gè)移動(dòng)用戶的導(dǎo)向矢量為:

a(θi，φi)=[ejβ1(θi，φi)，ejβ2(θi，φi)，…，ejβM(θi，φi)]T(3)

βk(θi，φi)=2πrλcos2π(k-1)M-θicos φi(4)

式中:βk(θi，φi)為第i個(gè)用戶入射到第k個(gè)陣元天線的相位值。對(duì)于Music算法[8，9]，需要構(gòu)造陣列協(xié)方差矩陣R:



R=E[XXH]=AE[SSH]AH+σ2I

=UsΣsUHs+UnΣnUHn(5)

圖2 智能天線的天線接收模型

噪聲和信號(hào)是相互獨(dú)立的，對(duì)協(xié)方差矩陣進(jìn)行分解。其中，AE[SSH]AH是信號(hào)部分，σ2I是噪聲部分。因此可對(duì)R進(jìn)行特征分解，如式(5)所示。式(5)中，Us是由大特征值對(duì)應(yīng)特征矢量組成的信號(hào)子空間，而Un是由小特征值對(duì)應(yīng)特征矢量組成的噪聲子空間。理想條件下信號(hào)子空間與噪聲子空間是相互正交的。實(shí)際中，協(xié)方差矩陣R為復(fù)Hermitian矩陣，一般通過Jacobi算法或QR算法迭代對(duì)其進(jìn)行特征分解，得到特征值和特征向量，再通過AIC準(zhǔn)則、MDL準(zhǔn)則、CCT等方法來估計(jì)信號(hào)源數(shù)目，進(jìn)而確定信號(hào)子空間和噪聲子空間。Music算法就是搜索式(6)的極大值來實(shí)現(xiàn)的對(duì)DOA的估計(jì):

PMusic=1aH(θ，φ)UnUHna(θ，φ)(6)

Music算法原理說明較大特征值對(duì)應(yīng)的向量為信號(hào)特征向量，因此估計(jì)信源數(shù)目肯定不超過陣元數(shù)目，如五元圓陣估計(jì)的移動(dòng)用戶數(shù)目不能同時(shí)超過五個(gè)，這對(duì)于3G中移動(dòng)定位業(yè)務(wù)的拓展是個(gè)極大的限制。其次，譜搜索是方位角從0°~360°的搜索過程，較大的計(jì)算量增加了基站處理的負(fù)荷，也影響了定位業(yè)務(wù)的實(shí)時(shí)性。綜上考慮，急需一種更大容量、更快速度的測(cè)角處理算法。

2 改進(jìn)的角度估計(jì)算法

TD-SCDMA采用時(shí)分的制式，上下行空中數(shù)據(jù)在同一個(gè)工作頻段。在一幀(10 ms)長(zhǎng)度中可認(rèn)為信道是不變的，因此由上行信道的接收數(shù)據(jù)估計(jì)出的信道參數(shù)可用于下行發(fā)射信號(hào)的波束成形，使下行信號(hào)對(duì)準(zhǔn)移動(dòng)用戶，這就是智能天線的測(cè)向技術(shù)。另外，TD-SCDMA對(duì)上行同步有著嚴(yán)格的要求(同步精度可達(dá)1/8 Chip，對(duì)應(yīng)時(shí)間為97 ns)。不同用戶的數(shù)據(jù)都要以基站時(shí)間為基準(zhǔn)，在預(yù)定時(shí)刻到達(dá)Node B。在上行同步的過程中，基站也同時(shí)估計(jì)出信號(hào)從移動(dòng)臺(tái)到基站的傳輸時(shí)間。在獲取移動(dòng)用戶相對(duì)基站的方位和距離的條件下，根據(jù)極坐標(biāo)原理，即可求出移動(dòng)用戶的相對(duì)位置。最終由基站的位置得到移動(dòng)用戶的絕對(duì)位置。

TD-SCDMA在時(shí)隙中設(shè)計(jì)一個(gè)midamble訓(xùn)練序列[10]。序列在系統(tǒng)中起著十分重要的作用，除了估計(jì)信道之外，還用于功率控制測(cè)量、上行同步保持、波束賦形、頻率校正等。這就為聯(lián)合檢測(cè)和方位角估計(jì)提供了物理基礎(chǔ)。圖3為midamble所在時(shí)隙結(jié)構(gòu)的位置。

圖3 TD-SCDMA系統(tǒng)突發(fā)結(jié)構(gòu)

TD-SCDMA系統(tǒng)一共有128個(gè)長(zhǎng)度為144位的midamble序列，midamble在每個(gè)常規(guī)時(shí)隙中發(fā)送，它的發(fā)送功率與數(shù)據(jù)符號(hào)相同，但不進(jìn)行擴(kuò)頻和加擾。設(shè)互不相關(guān)的midamble碼分別為m(1)，m(2)，…，m(k)，(k=1，2，…，128)，其中m(k)=[m(k)1，m(k)2，…，m(k)144]。

假設(shè)同時(shí)有N個(gè)用戶接入到小區(qū)網(wǎng)絡(luò)中，即智能天線同時(shí)接收到N個(gè)不同的midamble碼(N可以大于智能天線陣元數(shù)目M)時(shí)，第i個(gè)天線上(i=1，2，…，M)的接收數(shù)據(jù)為:



yi=∑Nl=1m(l)1ejβi(θl，φl)，∑Nl=1m(l)2ejβi(θl，φl)，…，

∑Nl=1m(l)144ejβi(θl，φl)+[n1，n2，…，n144](7)

[n1，n2，…，n144]為對(duì)應(yīng)的midamble序列碼片處的噪聲(設(shè)為高斯白噪聲)。令表示為第i個(gè)天線上的接收數(shù)據(jù)與本地m(k)序列的相關(guān)輸出:



=∑Nl=1∑144p=1m(l)pm(k)pejβi(θl，φl)+

=∑Nl=1ejβi(θl，φl)+

=∑Nl=1clkejβi(θl，φl)+(8)



那么各個(gè)天線上的接收信號(hào)與N個(gè)本地midamble序列的相關(guān)輸出為:



………

=ejβ1(θ1，φ1)ejβ1(θ2，φ2)…ejβ1(θN，φN)ejβ2(θ1，φ1)ejβ2(θ2，φ2)…ejβ2(θN，φN)ejβM(θ1，φ1)ejβM(θ2，φ2)…ejβM(θN，φN)#8226;

c11c12…c1Nc21c22…c2NcN1cN2cNN+H

=[a(θ1，φ1)，a(θ2，φ2)，…，a(θN，φN)]*C

+(9)

式中:C為midamble碼互相關(guān)矩陣;yi為各個(gè)陣元的接收數(shù)據(jù);仍為噪聲。因此，對(duì)于第k個(gè)移動(dòng)用戶方向矢量的最小二乘估計(jì)[11]為:

(θk，φk)=C-1(10)

由于用戶的攜帶midamble是已知的，C矩陣的逆可預(yù)先求出并存儲(chǔ)到DSP的數(shù)據(jù)空間中，以減少計(jì)算量。從式(9)可以看出，對(duì)用戶方向矢量的估計(jì)不受接收陣元數(shù)目的限制，只取決于midamble序列數(shù)目，最大可達(dá)到128個(gè)，遠(yuǎn)可以滿足用戶密集區(qū)域移動(dòng)定位業(yè)務(wù)的容量需求。

抽取單個(gè)用戶k的估計(jì)方向矢量進(jìn)行分析，由上文可知，利用估計(jì)的方向矢量可求出陣元間的相對(duì)相位差。用ψij表示第i個(gè)陣元與第j個(gè)陣元的相位差，則有:



ψij=βi(θk，φk)-βj(θk，φk)

=2πrλcos φkcos2πM(i-1)-θk-

cos2πM(j-1)-θk(11)

對(duì)于基站上的智能天線，陣元的結(jié)構(gòu)和數(shù)量是固定的。根據(jù)已知陣元的相對(duì)布局，利用最小二乘擬合方法，可以推導(dǎo)出波達(dá)方向和各個(gè)相位差的換算關(guān)系。下面以五元圓陣為例，闡述推導(dǎo)過程。

設(shè)4πrλsinπMcos φk=V，cos θk=c，sin θk=s，則對(duì)于各相鄰天線相位差化簡(jiǎn)得:

ψ21=V(s#8226;cosπ5-c#8226;sinπ5)

ψ32=V(s#8226;sin3π5-c#8226;cos3π5)

ψ43=Vs

ψ54=V(s#8226;cos3π5+c#8226;sin3π5)

ψ15=V(s#8226;cosπ5+c#8226;sinπ5)(12)

由導(dǎo)向矢量推出的實(shí)測(cè)相位差為:21，32，43，54，…，1M。采用擬合處理，取實(shí)測(cè)值與理想值之差的平方和為目標(biāo)函數(shù)，則有:

ε(θ)=(21-ψ21)2+(32-ψ32)2+…+(1M-ψ1M)2(13)



將式(13)按照M=5展開:

ε(θ)=2V-cosπ521-cos3π532+43-

cos3π554-cosπ515#8226;s+

2Vsinπ521+sin3π532+0-

sin3π554-sinπ515#8226;c+K(14)



式中:K是與s，c無關(guān)的常數(shù)。再設(shè)s前系數(shù)為a;c前系數(shù)為b，則ε(θ)可簡(jiǎn)寫成:

ε(θ)=as+bc+K

=K-a2+b2cos(θk-k)(15)



式中:k滿足:

cos k=-ba2+b2

sin k=-aa2+b2(16)

可推出cos(θ-)=1時(shí)，ε(θ)取得最小值。此時(shí)，=θ，即所估計(jì)出的方位角k為用戶k入射的方位角θk。將方位角代入任意一個(gè)相位差表達(dá)式即可求出入射俯仰角。

該算法是利用陣元間的相位差關(guān)系直接解算出用戶的方位角信息，省略了Music算法的譜搜索過程，極大地減少了計(jì)算量。算法雖然以五元均勻圓陣為例，但這種推導(dǎo)過程適用于任何均勻圓陣，可以滿足TD-SCDMA智能天線的要求，有很強(qiáng)的拓展性。

3 實(shí)驗(yàn)仿真

針對(duì)改進(jìn)的測(cè)角算法，通過實(shí)驗(yàn)來驗(yàn)證其有效性和準(zhǔn)確性，并從測(cè)向精度、測(cè)向耗時(shí)等方面與空間譜估計(jì)的Music算法進(jìn)行對(duì)比。

3.1 算法精度

仿真條件:假設(shè)智能接收天線為五元圓陣;單個(gè)移動(dòng)用戶的入射俯仰角為10°，方位角在0°~360°之間變化;快拍數(shù)等于midamble碼片數(shù)(144個(gè));信噪比變化范圍為0~10 dB;統(tǒng)計(jì)結(jié)果為100次獨(dú)立實(shí)驗(yàn)的平均。

對(duì)于單個(gè)移動(dòng)用戶源的入射，分別采用傳統(tǒng)測(cè)角算法和改進(jìn)測(cè)角算法進(jìn)行方位角測(cè)量。圖4~圖6分別為入射方位角為60°，130°，300°時(shí)兩種算法的估計(jì)均方根誤差對(duì)比結(jié)果。由三幅圖可知，采用改進(jìn)的角度估計(jì)算法可以達(dá)到空間譜估計(jì)的測(cè)角精度，即兩種測(cè)角算法的精度基本一致。

圖4 方位角為60°時(shí)兩種測(cè)向方法對(duì)比

圖5 方位角為130°時(shí)兩種測(cè)向方法對(duì)比

圖6 方位角為300°時(shí)兩種測(cè)向方法對(duì)比

考慮待估計(jì)用戶數(shù)量超過智能天線陣元數(shù)目時(shí)基站對(duì)移動(dòng)用戶的測(cè)向情況。仿真條件:假設(shè)待定位的移動(dòng)用戶數(shù)目為六個(gè)，入射方位角分別為20°，75°，90°，120°，130°，140°;智能天線為五元均勻圓陣，對(duì)這六個(gè)用戶同時(shí)進(jìn)行方位角的估計(jì);表1中數(shù)據(jù)為100次獨(dú)立實(shí)驗(yàn)的平均統(tǒng)計(jì)結(jié)果;信噪比分別為-5 dB，0 dB，5 dB，10 dB。

表1為五元圓陣下采用改進(jìn)的測(cè)角算法估計(jì)六個(gè)用戶的方位角測(cè)量精度。由于原理上的限制，此時(shí)空間譜估計(jì)測(cè)角算法將失效。

表1 五元圓陣估計(jì)六用戶的方位角精度

SNR /dB20°75°90°120°130°160°

-522.7478.3389.11118.09132.39158.50

020.7875.6289.80118.94129.64158.94

520.2274.9089.95119.46130.01159.86

1019.9874.9690.02120.00129.86159.98

從表 1可以看出，此時(shí)改進(jìn)的測(cè)角算法是有效的，并且保持很高的測(cè)量精度。

3.2 算法耗時(shí)

測(cè)向時(shí)間的比較分別在Matlab虛擬環(huán)境和DSP實(shí)際平臺(tái)上進(jìn)行。

Matlab虛擬環(huán)境:采用tic-toc命令來分別記錄空間譜估計(jì)測(cè)向和干涉儀測(cè)向的算法耗用平均時(shí)間。仿真時(shí)使用的Matlab版本為7.1;電腦配置為Pentium 4 2.0 GHz，384 MB內(nèi)存。DSP實(shí)際環(huán)境:利用CCS軟件，通過對(duì)處理器的cycle計(jì)數(shù)來估計(jì)兩種測(cè)向算法所耗的時(shí)間。仿真時(shí)使用的是TI公司的6416系列處理芯片，代碼使用了-o3的優(yōu)化功能。在兩種仿真環(huán)境中對(duì)單個(gè)移動(dòng)用戶的入射角進(jìn)行估計(jì)，記錄估計(jì)所耗時(shí)間如表2所示。

表2 兩種測(cè)向算法所耗時(shí)間比較

傳統(tǒng)測(cè)角算法改進(jìn)測(cè)角算法

Matlab記錄時(shí)間 /s1.277 80.020 3

DSP記錄cycle數(shù) /個(gè)621 03816 796

表2說明，無論在虛擬環(huán)境還是在實(shí)際環(huán)境中，改進(jìn)的測(cè)角算法較傳統(tǒng)的測(cè)角算法(Music)都大幅度地縮短了運(yùn)行時(shí)間。

4 結(jié) 語

主要針對(duì)原有定位方法中的測(cè)角算法進(jìn)行了改進(jìn)。實(shí)驗(yàn)結(jié)果表明:

精度方面:改進(jìn)的測(cè)角算法并沒有對(duì)相位差信息造成損耗，因此與原有的Music算法保持相同的測(cè)角精度。

定位用戶數(shù)目方面:原有的測(cè)角算法中采用空間譜估計(jì)技術(shù)，用戶數(shù)目不允許超過智能天線陣元數(shù)目;而

改進(jìn)的算法中，用戶的數(shù)目只受到midamble碼數(shù)目的限制，TD-SCDMA系統(tǒng)共有128個(gè)midamble碼，遠(yuǎn)可以滿足小區(qū)定位數(shù)目的需求。

計(jì)算量方面:新算法不需要Music中的譜搜索過程，而是由方向矢量利用陣元結(jié)構(gòu)的關(guān)系直接解算出方位角和俯仰角，因此計(jì)算量成數(shù)量級(jí)的減少。

本文依據(jù)TD-SCDMA系統(tǒng)特有的關(guān)鍵技術(shù)，提出了一種用戶數(shù)目不受限的快速定位算法，具有廣闊的應(yīng)用前景。

參 考 文 獻(xiàn)

[1]3GPP 25.302. Services provided by the physical layer[S].

[2]3GPP.TS 22.071. Location services(LCS) service description stage 1[S].

[3]3GPP2 PN4747.Location services enhancements[S].

[4]宋宇恒，林金朝，謝顯忠.TD-SCDMA系統(tǒng)定位算法研究與仿真[J].重慶郵電學(xué)院學(xué)報(bào)，2003，15(1):1-4.

[5]崔專城，孫肖子，徐莉，等.3G定位業(yè)務(wù)研究及在TD-SCDMA系統(tǒng)中的初步實(shí)現(xiàn)[J].數(shù)據(jù)通信，2004(4):39-41.

[6]李校林，趙勇，高飛.基于聯(lián)合參數(shù)估計(jì)法的TD-SCDMA單基站定位技術(shù)的研究[J].通信技術(shù)，2009，42(4):172-174.

[7]王永良，陳輝，彭應(yīng)寧，等.空間譜估計(jì)理論與算法[M].北京:清華大學(xué)出版社，2004.

[8]SCHMIDT R O.Multiple emitter location and signal parameter estimation [J].IEEE Trans.，1986，34(3):276-280.

[9]BARABELL A J.Improving the resolution performance of eigenstructure-based direction-finding algorithms[C].[S.l.]: ICASSP，1983.

[10]李小文.TD-SCDMA第三代移動(dòng)通信系統(tǒng)、信令及實(shí)現(xiàn)[M].北京:人民郵電出版社，2003.

[11]郭梯云，鄔國(guó)揚(yáng)，李建東.移動(dòng)通信[M].西安:西安電子科技大學(xué)出版社，2000.