趙滿川，王玉文，尹　婷，陳　捷

(電子科技大學(xué) 航空航天學(xué)院，四川 成都 611731)

?

基于FFT和ESPRIT融合的單星測頻定位技術(shù)

趙滿川，王玉文，尹婷，陳捷

(電子科技大學(xué) 航空航天學(xué)院，四川 成都 611731)

針對我國應(yīng)急通信服務(wù)主要依靠國外衛(wèi)星通信系統(tǒng)的重大隱患，提出了基于FFT和ESPRIT融合的單星測頻定位技術(shù)。通過FFT和ESPRIT算法詳細論述了頻移估計過程，分析了2種算法各自的優(yōu)缺點，提出了基于FFT和ESPRIT融合的新的頻移估計算法。用蒙特卡羅方法模擬仿真了在不同信噪比下各算法的測頻精度，比較了各個算法的性能，并對頻移的誤差進行了分析。基于應(yīng)急通信與導(dǎo)航定位的需求，提出了單星測頻定位體制，用提出的融合算法所測得的多普勒頻移實現(xiàn)單星測頻定位，對定位結(jié)果進行了仿真，結(jié)果表明，提出的融合算法有助于提高單星測頻定位的精度，能滿足應(yīng)急條件下的導(dǎo)航定位需求。

單星；多普勒頻移；融合；導(dǎo)航定位；ESPRIT

0　引言

近些年，各種自然災(zāi)害和人為事故頻發(fā)，實施搶險救災(zāi)的首要任務(wù)就是為救援人員和部門等提供準確的導(dǎo)航定位服務(wù)。為此，研究具有高精度的應(yīng)急響應(yīng)導(dǎo)航定位系統(tǒng)對挽救人民群眾的生命財產(chǎn)安全、維護國家的國土安全具有十分重要的現(xiàn)實意義[1]。目前，我國的單星導(dǎo)航定位技術(shù)發(fā)展還處于初級研制階段，由于技術(shù)不夠成熟，導(dǎo)致功能與費用上都存在問題，不能在災(zāi)難面前用作一種應(yīng)急導(dǎo)航定位手段[2]。為了彌補我國在單星導(dǎo)航定位技術(shù)領(lǐng)域的缺欠，解決我國應(yīng)急響應(yīng)的需求，本文提出的單星測頻定位關(guān)鍵技術(shù)研究，使得終端設(shè)備能夠配合應(yīng)急響應(yīng)衛(wèi)星實現(xiàn)單星定位的功能。

多普勒頻率是任何具有相對運動物體之間電磁信號的頻率變化特性，對于運動目標來說，可以利用運動目標的運動所引起的頻率變化來確定目標的位置和運動特性，對于靜止的目標，可以人為的產(chǎn)生相對運動來確定目標的位置[3]。由于衛(wèi)星和目標輻射源之間始終存在著相對運動，因此衛(wèi)星上接收到的信號是受多普勒效應(yīng)影響的。顯然，信號的瞬時頻率和脈沖信號的脈沖重復(fù)間隔(PRI)會發(fā)生變化，且信號頻率的改變又產(chǎn)生了頻率變化率等參數(shù)。這些參數(shù)反映著多普勒效應(yīng)的強度，實際上就是反映了當(dāng)前衛(wèi)星與目標輻射源之間的相對位置關(guān)系和相對運動關(guān)系[4]。因此，衛(wèi)星的多普勒頻移測量精度直接決定著單星測頻定位的精度。針對單星快速定位精度不高的原因，本文提出了基于快速傅里葉變換(Fast Fourier Transformation，F(xiàn)FT)和旋轉(zhuǎn)不變子空間(Estimating Signal Parameters via Rotational Invariance Techniques，ESPRIT)相融合的算法，以提高測量的多普勒頻移精度，改善單星快速定位的精度。

1　多普勒頻移估計算法

1.1基于FFT算法的頻移估計

FFT屬于經(jīng)典的頻域分析方法，該技術(shù)具有相對較小的計算量，并且在一定的噪聲環(huán)境下仍能穩(wěn)定地工作。將FFT應(yīng)用在多普勒頻移估計中要先對接收端信號進行下變頻，分離出采樣數(shù)據(jù)中的多普勒頻率部分，然后直接估計多普勒頻率值[5]。通過計算得出FFT算法求出的頻率與真實頻率存在一定的誤差fe：

(1)

式中，fs=1/Ts表示采樣率；nmax為高斯白噪聲；N為采樣點。

同時傅氏變換也存在許多問題，比如柵欄效應(yīng)和頻譜泄漏。這些問題來自于對信號周期性的采樣上硬性問題，包括了對無限長信號的有限采樣和對有限信號數(shù)據(jù)準同步采樣[6]。

1.2基于ESPRIT算法的頻移估計

ESPRIT算法最基本的是假設(shè)2個具有相同結(jié)構(gòu)和陣元數(shù)(M)且間距為已知的子陣[7]。當(dāng)接收到一個入射信號時，因為陣元對的一致性，使得陣元輸出的數(shù)據(jù)間只相差一個相位[8]。這里用X1和X2來表示陣元對的輸出，于是有

(2)

式中，S(t)為入射信號源；N1(t)和N2(t)分別為子陣1和子陣2的白噪聲矢量，且N1(t)和N2(t)的均值為0，方差為σ2。子陣1的陣列流型矩陣A1=[a(θ1),…,a(θk)]，子陣2的陣列流型矩陣A2=A1φ，φ=diag[e-jφ1,…,e-jφk]。

由上式可知，只要得到2個子陣間的旋轉(zhuǎn)不變關(guān)系φ就能得到需要的信號參數(shù)。將2個子陣模型合并得到[9]：

(3)

計算協(xié)方差并進行特征分解得到：

(4)

(5)

φ可以由φ的特征值來構(gòu)成，同時φ對應(yīng)的特征向量可以構(gòu)成T，進而可以得到旋轉(zhuǎn)不變關(guān)系矩陣φ，得到需要的信號參數(shù)。

1.3基于FFT和ESPRIT融合算法的頻移估計

經(jīng)典譜算法本質(zhì)上都相當(dāng)于使用帶通信號與源信號進行時域乘積，所以會產(chǎn)生頻譜能量泄露的效應(yīng)，ESPRIT估計算法在運行中需要分離出噪聲空間。文獻[10-11]中提出一種FFT、MUSIC聯(lián)合測頻算法，在該算法中利用FFT算法的旋轉(zhuǎn)不變性，采用經(jīng)過FFT變換后的譜峰系數(shù)xk構(gòu)造數(shù)據(jù)矩陣，利用FFT對噪聲的不敏感性提高MUSIC算法的性能。算法雖然提高了頻率估計精度，但是還有提升的空間。為此利用FFT算法的旋轉(zhuǎn)不變性，得到大致的頻譜范圍，然后利用粗估計的頻率進行ESPRIT算法，實現(xiàn)精估計。

為了能夠和噪聲信號更好地分離開來，這里考慮加入自相關(guān)技術(shù)。因為信號是完全獨立于噪聲的，所以通過自相關(guān)運算可以達到去除噪聲的目的。經(jīng)過FFT變換的正弦信號，保留譜峰范圍為(k0-Δ,k0+Δ)的數(shù)據(jù)，可以得：

(6)

(7)

從式(7)可以看出，自相關(guān)運算和FFT變換在信號的形式上具有一定的相似性，所以可以直接在FFT變換過程中使用譜線進行自相關(guān)運算，即在頻域進行自相關(guān)，然后再進行重構(gòu)信號。針對譜峰附近的數(shù)據(jù)進行重新挑選，可以將白噪聲的干擾減小，提高了信噪比，達到提高ESPIRT算法輸入數(shù)據(jù)的信噪比的目的，使得ESPRIT算法在中低信噪比時具有更好的抗噪特性。

① 初始化

② 迭代，求出ES1=W(n)。

這里將新的算法命名為Alg1。

1.4測頻算法仿真

對Alg1的算法性能進行仿真實驗，設(shè)置頻偏900 kHz，采樣點數(shù)為256，采樣率設(shè)置為2 MHz。在(-15 dB，20 dB)的范圍進行仿真，對每個SNR點做2 000次Mento Carlo實驗，結(jié)果如圖1所示。從圖1可以看出，新的算法整體性能比ESPRIT算法更好，在中高信噪比下更是超越了TLS-ESPRIT算法。算法Alg1對多普勒頻移的估計更能逼近真實值，頻移估計的誤差也更小。

圖1　算法測頻誤差與信噪比關(guān)系

多普勒頻率是任何具有相對運動物體之間電磁信號的頻率變化特性，對于運動目標來說，可以利用運動目標的運動所引起的頻率變化來確定目標的位置和運動特性，對于靜止的目標，可以人為地產(chǎn)生相對運動來確定目標的位置。由于衛(wèi)星和目標輻射源之間始終存在著相對運動，因此衛(wèi)星上接收到的信號是受多普勒效應(yīng)影響的。當(dāng)衛(wèi)星運行在高軌道的遠地點附近，需要用多普勒頻移對目標進行定位解算時，多普勒頻移估計越精確，定位誤差越小。為此用本文提到的多普勒頻移估計技術(shù)對20 000 km以上的單星測頻定位進行仿真，驗證新的算法對單星測頻定位性能的改善情況。

2　基于FFT和ESPRIT的單星測頻定位

仿真

單星測頻定位，即只利用一顆衛(wèi)星收集的衛(wèi)星數(shù)據(jù)和多普勒頻移數(shù)據(jù)實現(xiàn)對目標位置信息的解算，如圖2所示，文獻[12]詳細介紹了單星測頻定位的原理，其實質(zhì)上就是利用時間信息換取衛(wèi)星在空間位置信息的差異。當(dāng)衛(wèi)星處于遠地點軌道時，衛(wèi)星的相對速度較小，能在較長時間內(nèi)對地面目標區(qū)域進行覆蓋，有利于在特殊條件下的應(yīng)急通信。但是此時在短時間內(nèi)衛(wèi)星的位置信息基本保持不變，要實現(xiàn)對目標區(qū)域的位置信息解算就需要精度較高的多普勒頻移測量算法。

圖2　單星測頻定位模型

根據(jù)多普勒效應(yīng)，某一時刻ti衛(wèi)星接收到的信號頻率fi與輻射源的發(fā)射頻率fT有如下關(guān)系：

(8)

式中，fi為衛(wèi)星在ti時刻接收到的頻率；Vi為ti時刻衛(wèi)星的速度；θi為衛(wèi)星和輻射源連線與衛(wèi)星速度方向的夾角；c為光速；fT為輻射源發(fā)射的頻率。

根據(jù)3個不同時刻的測量值可以得到3個不同的方程，聯(lián)立這3個方程組成一個非線性方程組，利用解析的方法即可求解這個方程組即可得到輻射源的坐標(x,y,z)。

利用STK軟件構(gòu)造高軌道衛(wèi)星運行軌跡，獲取遠地點附近的衛(wèi)星軌道數(shù)據(jù)和多普勒頻移數(shù)據(jù)，并把此頻移作為基準頻移。用各個測頻算法算出多普勒頻偏并與基準頻移進行比較，結(jié)果如圖3所示?；诟鳒y頻算法的單星測頻定位結(jié)果如圖4所示(頻移測量時間為60 s)。

圖3　各測頻算法的性能比較

圖4　各測頻算法下的單星定位均方誤差

從圖3可以看出，用本文提出的測頻算法得出的頻移能更好地擬合基準頻移，相比TLS-ESPRIT和ESPRIT測頻算法，新算法平穩(wěn)性更好。圖4仿真了在遠地點用單星測頻定位得到的均方誤差，可以看出基于新算法的測頻定位精度有一定的提高。但是由于衛(wèi)星運行在20 000 km左右的遠地點附近，且只有一顆衛(wèi)星，衛(wèi)星數(shù)據(jù)之間的相關(guān)性較強，得到的測頻定位精度還有很大的提高空間。

3　結(jié)束語

本文分別使用經(jīng)典的FFT算法和子空間類型的ESPRIT算法對多普勒頻移進行估計，提出了基于FFT和ESPRIT融合的多普勒頻移估計算法，通過仿真驗證了新算法有助于提高ESPRIT類型算法在低信噪比環(huán)境下的頻移估計性能和穩(wěn)定性。在單星測頻定位系統(tǒng)中，特別是當(dāng)衛(wèi)星運行在遠地點軌道時，在短時間內(nèi)衛(wèi)星的位置信息差異很小，此時要利用多普勒頻移定位，對多普勒頻移的測量精度有較高的要求。為此把本文提出的新算法得到的多普勒頻移用到單星測頻定位中，通過仿真驗證了新的測頻算法提高了單星測頻定位的精度。同時也看到，當(dāng)單星衛(wèi)星運行在遠地點附近軌道時，雖然能在長時間內(nèi)覆蓋目標區(qū)域，有利于在應(yīng)急條件下與地面建立通信服務(wù)，但是要在短時間內(nèi)實現(xiàn)對地面目標位置信息的高精度解算還有進一步提高的空間。

[1]王玉峰,孟蘇英,劉英.應(yīng)急通信保障情況綜合處理平臺研究[J].無線電通信技術(shù),2009,35(3):55-58.

[2]劉振嶺,張水興,李駿.短波通信網(wǎng)在應(yīng)急通信領(lǐng)域的應(yīng)用[J].無線電通信技術(shù),2014,40(2):82-85.

[3]張騫丹,田紅心.GPS系統(tǒng)多普勒頻移估算的研究[J].無線電工程,2007,37(4):21-23.

[4]劉麗哲,馮東升.多普勒頻移測量方法研究[J].無線電工程,2009,39(5):51-53.

[5]于飛,吳文全.基于FFT頻率測量誤差分析[J].電子測量技術(shù),2004(6):3-5.

[6]齊國清.幾種基于FFT的頻率估計方法精度分析[J].振動工程學(xué)報,2006,19(1):86-92.

[7]司偉建,林晴晴.基于延時相關(guān)處理的ESPRIT算法[J].系統(tǒng)工程與電子技術(shù),2012,34(3):439-443.

[8]SUN W T,JUN L V,TONG L I.An Algorithm of Frequency Estimation for ESPRIT Frequency-Hopping Signal Based on FFT Rotational Invariance[J].Measurement & Control Technology,2014，23(6):347-353.

[9]XU B,SUN L,XU L,et al.Improvement of the Hilbert Method via ESPRIT for Detecting Rotor Fault in Induction Motors at Low Slip[J].IEEE Transactions on Energy Conversion,2013,28(28):225-233.

[10]ZHANG S,ZHU Y,KUANG G.Imaging of Downward-looking Linear Array Three-dimensional SAR Based on FFT-MUSIC[J].IEEE Geoscience & Remote Sensing Letters,2015,12(4):885-889.

[11]蔣毅.基于非參量法的諧波頻率估計研究[D].成都：電子科技大學(xué),2008.

[12]陸安南,孔憲正.單星測頻無源定位法[J].通信學(xué)報,2004,25(9):160-168.[13]李獻斌.單星測頻無源定位技術(shù)研究[D].長沙：國防科學(xué)技術(shù)大學(xué),2009.

[14]李楊,周彬.基于概率的單星測頻無源定位算法[J].電子信息對抗技術(shù),2007,22(3):7-10.

[15]李璟璟,張迎春,鄭靖,等.基于信息融合的大橢圓軌道衛(wèi)星組合導(dǎo)航方法[J].宇航學(xué)報,2012,33(9):1 233-1 240.

趙滿川男，(1990—)，碩士研究生。主要研究方向：單星導(dǎo)航定位、空天地一體化測控通信。

王玉文男，(1962—)，副教授,高級工程師。主要研究方向：空天地一體化測控通信與導(dǎo)航定位技術(shù)、空間系統(tǒng)運行控制與交通管理技術(shù)。

Single-satellite Frequency Measurement Positioning Technology Based on FFT and ESPRIT Fusion

ZHAO Man-chuan,WANG Yu-wen,YIN Ting,CHEN Jie

(SchoolofAstronauticsandAeronautics，UniversityofElectronicScienceandTechnologyofChina,ChengduSichuan611731,China)

To address the issue that the emergency communication services in China mainly rely on overseas satellite communication system,this paper presents the single-satellite frequency measurement positioning technology based on the fusion algorithm of FFT and ESPRIT.The frequency shift estimation process is discussed in detail through FFT and ESPRIT algorithms,then the advantages and disadvantages of the two algorithms are analyzed,and a new algorithm based on FFT and ESPRIT fusion is proposed.The precision of frequency measurement is simulated by Monte-Carlo method under different SNRs,the performances of the algorithms are compared,and the errors of Doppler frequency shift are analyzed.Based on the demands of emergency communication,positioning and navigation,the paper proposes a single-satellite positioning system,and uses the Doppler frequency shift measured by the fusion algorithm proposed to achieve the single-satellite frequency measurement positioning.Finally,the simulated results show that the proposed algorithm helps to improve the precision of single-satellite frequency measurement positioning,which meets the needs of navigation and positioning under emergency conditions.

single-satellite;Doppler frequency shift;fusion;navigation and positioning;ESPRIT

10.3969/j.issn.1003-3106.2016.10.11

2016-07-06

TN96.1

A

1003-3106(2016)10-0047-04

引用格式：趙滿川，王玉文，尹婷，等.基于FFT和ESPRIT融合的單星測頻定位技術(shù)[J].無線電工程,2016,46(10):47-50，68.