基于FFT/IFFT的弱信號捕獲方法研究

袁先舉，盧艷娥

（1.貴州大學 大數(shù)據(jù)與信息工程學院，貴州 貴陽　550025；2.中國科學院 重慶綠色智能技術(shù)研究院北斗導航中心，重慶 400714）

基于FFT/IFFT的弱信號捕獲方法研究

袁先舉1，2，盧艷娥2

（1.貴州大學 大數(shù)據(jù)與信息工程學院，貴州 貴陽550025；2.中國科學院 重慶綠色智能技術(shù)研究院北斗導航中心，重慶 400714）

為了縮短捕獲時間，提高對微弱信號的檢測能力，利用FFT/IFFT方法代替?zhèn)鹘y(tǒng)時域串行相干累積，相干累積之后的處理，提出采用各個相干累積值相乘取模的方法提升檢測量，最后進行門限比較。該方法在提高信噪比的同時且計算量更少，平方損耗小。仿真實驗表明，檢測量提升比非相干方法快，檢測概率比非相干大，門限設(shè)置范圍更廣，可用于捕獲GPS和北斗等GNSS衛(wèi)星微弱信號。

微弱信號捕獲；提升檢測量；相干值相乘；FFT/IFFT

隨著我國北斗系統(tǒng)、美國GPS、歐洲Galileo系統(tǒng)、俄羅斯GLONASS系統(tǒng)的發(fā)展，國內(nèi)外許多專家學者都在進行GNSS接收機的研制和改善。其中接收機研制的一個重要方向就是對微弱信號的捕獲［1］。目前接收機的弱信號捕獲方法主要是向輔助GNSS［2］和高靈敏度GNSS［3］兩方向發(fā)展。于高靈敏度接GNSS收機而言，通常通過加長相干積分時間和非相干積分時間提高信噪比來捕獲微弱信號。而導航數(shù)據(jù)比特的跳變，會影響到相干累積結(jié)果，為消除導航數(shù)據(jù)比特跳變帶來的影響，估計比特位翻轉(zhuǎn)的全比特法和半比特法得到應用，從而得到更長時間的相干累積結(jié)果，提高微弱信號的信噪比。由于FFT方法運算的有效快速性，許多學者都基于FFT方法進行弱信號捕獲的研究［4-5］，部分匹配濾波和快速傅里葉變換（PMF-FFT）結(jié)合的方法［6-7］與FFT/IFFT方法得到廣泛應用。采用接收到的數(shù)據(jù)進行每10 ms數(shù)據(jù)分組，然后數(shù)據(jù)進行奇偶分組，加長了相干積分時間，消除了導航比特跳變的影響。相干積分之后的處理，采用提出的各個相干值相乘法，迅速提升檢測量。

1　基于FFT/IFFT的相干值相乘捕獲方法

1.1基于FFT/IFFT的相干值相乘法原理

如果接收到長度為N點的衛(wèi)星信號偽碼序列為a（n），長為N點本地復制碼序列為b（n），則a（n）與b（n）的相關(guān)值z（n）可表示為

由于在做長時間的相干累積時，計算量較大。利用FFT/ IFFT方法代替?zhèn)鹘y(tǒng)相干累積操作，若對（1）式進行離散傅里葉變換，可得到

式（2）中，A（K）與B（K）分別表示a（n）與b（n）傅里葉變換，“*”表示取復共軛；根據(jù)傳統(tǒng)半比特非相干捕獲［8］的思想，相干后續(xù)處理用相干值相乘，取相乘后最大峰值與次最大峰值的比值（以下稱峰-峰比值）作為門限判決量，可得到如圖1所示的捕獲框圖。若接收數(shù)據(jù)奇數(shù)分組和偶數(shù)分組分別為m組，第m個相干累積值為Y（m），則其中一組相干相乘輸出為

采用FFT/IFFT方法，N ms的 數(shù)據(jù)，每1 ms數(shù)據(jù)j個采樣點經(jīng)過FFT/IFF運算后，第i ms得到每個碼相位處的相干累積值數(shù)組Ai（1，j）（i=1，2，3…，N），最后得到的相干相乘輸出值為

與傳統(tǒng)幾種微弱信號捕獲方法相比，相干值相乘法計算量更小。表1給出了幾種方法運算量的比較，把1 ms GPS數(shù)據(jù)采樣成2 048個點，設(shè)置多普勒頻移搜索單元為20，PMFFFT以32點FFT，64點相干累加為例。表2給出了非相干操作與相干值相乘法的運算量情況，以M次相干值為例。結(jié)合表1和表2數(shù)據(jù)可以分析出，相同數(shù)據(jù)長度，F(xiàn)FT/IFFT相干值相乘法運算次數(shù)明顯低于其他幾種算法，當相干積分次數(shù)越多，相干值相乘法計算量比起非相干法乘積計算量減少越明顯。

1.2檢測量提升分析

經(jīng)過FFT/IFFT后，在每個碼相位處輸出相干積分值I（n）+jQ（n），其中：

式中，P為接收信號功率，R（·）為偽碼的自相關(guān)函數(shù)，τn為相位估計誤差，Tc為相干積分時間，且有TC=NcTs，Nc為Tc內(nèi)采樣的點數(shù)。Δfd為多普勒頻移估計誤差，Δφ為載波相位殘差，滿足Δφ=2πΔfd［（k-1）Tc+Tc/2］+θe；以FFT/IFFT輸出相乘值作為檢測量，相干后的I/Q支路含有均值為零的高斯噪聲，I（n）+jQ（n）表示為（I（n）+nI）+j（Q（n）+nQ），其中，nI，nQ～N（0，σ2），有：

式中，nI+jQ=nI+jnQ～N（0，2σ2），R（n）=I（n）+jQ（n）；

若每次FFT/IFFT后得到的各個碼相位處的相干幅值都相同，N次相干值相乘法得到下式：

表1　幾種方法計算量比較Tab.1 Comparison of several methods

表2　非相干與相干值相乘運算量Tab.2 Incoherent and coherent value multiplied

圖1　相干值相乘法捕獲框圖Fig.1　Block diagram of Coherent multiplication value

式中第二個等號右邊第三項是均值為零的高斯分布，在累加過程中可被濾除；等號右邊第二項為（nI+jQ）N～（0，σ2N），此項噪聲均值為零，可被濾除，無平方損耗。

根據(jù)以上分析，可以得到M個相干值相乘值為：

以上相干值相乘取模得：

非相干積分中：

式（10）中nNC是均值不為零的高斯噪聲。這種噪聲不能被積分器濾除，這是導致非相干積分平方損耗的原因。非相干積分在計算相關(guān)信號功率的時候，同一噪聲也被平方了，平方后的噪聲均值不為零，這導致了噪聲不能被積分器濾除［9］。

假設(shè)每次FFT/IFFT過后正確碼相位處的相干幅值A(chǔ)和其他碼相位處的最大幅值B（A≥B）不變，實際情況下，每一個相干積分后A、B值不是固定不變的，B值變動會抑制其增長速度，峰-峰比值會得到迅速增長。如圖2得到不同M值和不同A/B值的情況下，峰峰比值提升的變化情況。

圖2　不同M和A/B值的峰峰比值Fig.2　Peak ratios of different M and A/B values

2　實驗仿真與驗證

2.1用模擬數(shù)據(jù)對算法仿真

利用MATLAB軟件平臺，模擬產(chǎn)生GPS 1號衛(wèi)星數(shù)據(jù)，加入適當高斯白噪聲，對200 ms模擬GPS數(shù)據(jù)進行捕獲，設(shè)定接收碼相位延遲721個碼片，仿真相干相乘捕獲結(jié)果與傳統(tǒng)半比特捕獲法做比較，得到碼相位捕獲圖如圖3所示。從圖3中可以看出，相干值相乘法的碼相位捕獲圖互相關(guān)干擾峰值更小，效果更佳，且峰峰比值遠大于傳統(tǒng)半比特非相干捕獲法的比值，從碼相位捕獲結(jié)果來看，在721碼相位處取得最大值，可以在正確碼相位處成功捕獲信號。

圖3　半比特交替與新算法碼相位捕獲Fig.3　Half bit and the new algorithm acquisition results

2.2實際接收的衛(wèi)星數(shù)據(jù)進行算法的驗證

設(shè)置射頻前端下變頻后的中頻頻率為11.029 6 MHz，采樣頻率為49.999 06 MHz，C/A碼碼速率1.023 MHz，則一個周期C/A碼碼片長為1 023碼片。基于Darius Plausinaitis和Dennis M.Akos的軟件接收機，結(jié)合提出的相干值相乘法捕獲衛(wèi)星信號，通過檢測峰-峰比值來判斷是否成功捕獲到衛(wèi)星信號。如圖4和圖5分別給出了相干值相乘法以及傳統(tǒng)半比特算法捕獲衛(wèi)星的結(jié)果，為證明相干值相乘法成功捕獲，根據(jù)捕獲到的衛(wèi)星，解調(diào)出的各顆衛(wèi)星的載波頻率值，碼相位補償值，以及當前衛(wèi)星的多普勒頻移值，在表3給出。結(jié)合圖6捕獲結(jié)果和表3數(shù)據(jù)可以得出，相干值相乘法可成功捕獲衛(wèi)星信號，且較非相干方法門限范圍設(shè)置更寬，在相同捕獲門限下，檢測概率更大，捕獲結(jié)果峰值更高。

圖4　相干值相乘法捕獲結(jié)果Fig.4　The results of coherent multiplication value

圖5　半比特非相干法捕獲結(jié)果Fig.5　The results of incoherent method

3　結(jié)束語

提出的相干值相乘法比起非相干積分，計算量更小，同時對噪聲的提高較小，平方損耗低于非相干積分。從模擬數(shù)據(jù)仿真結(jié)果看出，獲得的碼相位自相關(guān)幅值大于傳統(tǒng)半比特非相干捕獲法，檢測量提高快。實際接收的數(shù)據(jù)捕獲驗證結(jié)果，處理200 ms數(shù)據(jù)可成功捕獲八顆GPS衛(wèi)星，捕獲載噪比可低至 23 dBHz。相干值相乘法也可以用來捕獲北斗等GNSS衛(wèi)星信號，能快速提升正確碼相位處的相干幅值，捕獲結(jié)果明顯優(yōu)于半比特非相干捕獲法，在對微弱信號的捕獲方面有一定意義。

表3　相干值相乘法捕獲解算結(jié)果Tab.3 Coherent multiplication value capture solver results

［1］莫凡，王新龍.高靈敏度GPS捕獲技術(shù)研究綜述［J］.航空兵器，2014（2）:19-26.

［2］Wang X，Ji X，F(xiàn)eng S.A scheme for weak GPS signal acquisition aided by SINS information［J］.GPS Solutions，2013，18 （2）:243-52.

［3］Shanmugam S.New enhanced sensitivity detection techniques for GPS L1 C/A and modernized signal acquisition［D］. Canada，University of Calgary（Canada），2008.

［4］任記達.基于FFT/iFFT的GPS信號捕獲算法研究 ［J］.現(xiàn)代導航，2013（6）:407-410.

［5］姜士奇.一種改進的基于FFT的擴頻碼快速捕獲方法［J］.黑龍江科技信息，2013（17）:134.

［6］齊華，施小茜，冀樂，等.PMF-FFT的PN碼捕獲方法分析及仿真［J］.西安工業(yè)大學學報，2010（1）:57-61.

［7］方科.基于PMF-FFT的擴頻信號多通道并行捕獲［J］.通信技術(shù)，2013（7）:8-16.

［8］Elders B H，Derrmar U.Efficient Differentially Coherent Code/ Doooler Acquisition of Weak GPS［C］//Spread Spectrum Techniques and Applications，2004 IEEE Eighth International asaymposium on.IEEE，2004:731-735.

［9］謝鋼.GPS原理與接收機設(shè)計［M］.北京:電子工業(yè)出版社，2009.

The research of weak signal acquisition algorithm based on FFT/IFFT

YUAN Xian-ju1，2，LU Yan-e2
（1.College of Big Data and Information Engineering，Guizhou University，Guiyang 550025，China；2.Beidou Navigation Center of Chongqing Institute of Green and Intelligent Technology，Chinese Academy of Sciences，Chongqing 400714，China）

In order to shorten acquisition time and improve the ability to detect weak signals using FFT/IFFT method instead of the traditional time-domain serial coherent accumulation，After coherent accumulated processing is proposed using various coherence cumulative value multiplied modulo method to enhance the detection volume，F(xiàn)inally，to compare with the threshold.which at the same time to improve the signal to noise ratio and the calculation of the amount of less-squares loss. Simulation results show that，using FFT/IFFT method can make detection of increased faster，detection probability larger and threshold setting range wider compared with non-coherent detection method.It can be used to capture GPS and compass GNSS satellite weak signal.

weak signal acquisition；enhance detectable quantity；coherent value multiplied；FFT/IFFT

TN966

A

1674－6236（2016）05-0126-04

2015-04-17稿件編號：201504184

袁先舉（1990—），男，貴州遵義人，碩士研究生。研究方向：電路系統(tǒng)設(shè)計自動化。