文獻標識碼： A
DOI：10.16157/j.issn.0258-7998.174212
中文引用格式： 楊銳，黃海生，李鑫，等. 適用于北斗GNSS-R接收機的反射信號捕獲算法[J].電子技術應用，2018，44(8)：118-121，125.
英文引用格式： Yang Rui，Huang Haisheng，Li Xin，et al. A reflected signal acquisition algorithm for Beidou GNSS-R receiver[J]. Application of Electronic Technique，2018，44(8)：118-121，125.
0 引言
全球導航衛(wèi)星系統(tǒng)(Global Navigation Satellite System，GNSS)不僅可以為用戶提供導航定位信息、授時等功能，其反射信號也可以被接收與處理。利用GNSS反射信號進行目標遙感測量的技術，被稱為全球導航衛(wèi)星系統(tǒng)反射信號遙感技術(Global Navigation Satellite System-Reflections，GNSS-R)[1]。通過處理GNSS衛(wèi)星的反射信號，可以進行海面、地表、森林等參數(shù)的遙感探測與反演[2-4]。中國北斗衛(wèi)星導航系統(tǒng)(BeiDou Navigation Satellite System，BDS)是中國自行研制的全球衛(wèi)星導航系統(tǒng)。它可以進行全方位、多信號源、廣覆蓋、高分辨率的地表觀測實驗，為海洋、陸地、冰面等反射面的監(jiān)測提供了更為可靠有效的方法，這也必將成為北斗系統(tǒng)應用的一個重要領域[5]。
反射信號的捕獲是GNSS-R技術中的關鍵。然而信號經(jīng)反射面反射后，信號功率大幅度衰減，加大了反射信號的捕獲難度。對于反射信號的捕獲，通過加長積分時間提高信噪比。目前，常見的反射信號捕獲算法主要有相干非相干算法和差分相干捕獲算法[6]。其中，相干積分對于信號增益的提升最為有效，但相干積分受限于導航數(shù)據(jù)跳變[7]。本文針對北斗GNSS-R接收機，提出一種可以克服導航數(shù)據(jù)跳變的反射信號捕獲算法。
1 北斗衛(wèi)星信號概述
北斗衛(wèi)星導航系統(tǒng)簡稱北斗系統(tǒng)，其空間星座由5顆地球靜止軌道(GEO)衛(wèi)星、27顆中圓地球軌道(MEO)衛(wèi)星和3顆傾斜地球同步軌道(IGSO)衛(wèi)星組成。北斗衛(wèi)星1～5號星屬于GEO衛(wèi)星。GEO衛(wèi)星播發(fā)D2導航電文，數(shù)據(jù)率為500 b/s。北斗衛(wèi)星6～37號星屬MEO/IGSO衛(wèi)星，播發(fā)D1導航電文，數(shù)據(jù)率為50 b/s，并調制有速率為1 Kb/s的二次編碼。D1導航電文上調制的二次編碼是指在速率為50 b/s的電文上調制一個NH碼。NH碼周期為20 ms，共20 bit信息位（0 0 0 0 0 1 0 0 1 1 0 1 0 1 0 0 1 1 1 0），碼速率為1 Kb/s，碼寬為1 ms，以模二加的形式與擴頻碼和導航信息碼同步調制[8]。圖1為北斗信號二次調制原理圖。
2 反射信號捕獲算法研究
2.1 信號捕獲原理
衛(wèi)星信號的捕獲是對多普勒頻率及C/A碼的二維搜索過程，利用C/A碼的強自相關性，在二維搜索過程中找出一個自相關峰值，從而確定C/A碼相位偏移以及多普勒頻率信息。常見的捕獲搜索算法為頻率串行碼相位并行搜索，其處理框圖如圖2所示。
2.2 傳統(tǒng)反射信號捕獲算法
2.2.1 相干-非相干積分算法
相干積分利用的是C/A碼的強自相關性，對N個C/A碼周期的相關值矩陣直接累積，數(shù)學表達式如下：
式中，yk(t，fd)表示第k個C/A碼相關值矩陣，t和fd分別代表碼相位點和多普勒頻率。
相干積分增益可以表示為[9]：
式中，Bb=1/T，T為相干積分時間。對于北斗信號，則BW=4.092 MHz。式(2)表明，增加相干積分時間可以提高信號增益，由于導航數(shù)據(jù)和跳轉邊沿未知，使得相干積分時間不能無限增大，此時可以通過采用非相干積分的方式來提高信噪比。
非相干積分是對N個相干積分結果進行模值平方后再求和[10]，數(shù)學表達式如下：
非相干積分長度不受導航數(shù)據(jù)和跳轉邊沿影響，但是存在噪聲平方損失。文獻[11]分析了其增益和損耗，指出非相干次數(shù)在少的時候是有效的，但隨著累積次數(shù)的增多，損耗也會增加，總的處理增益會下降。
2.2.2 差分相干算法
差分相干算法，其原理可概括為如下表達式：
文獻[12]指出差分相干積分能夠減小平方損耗，對信噪比的提升優(yōu)于非相干積分。
2.3 直射輔助式反射信號捕獲算法
2.3.1 直射信號與反射信號關系研究
GNSS-R原理如圖3所示，接收機在接收衛(wèi)星直射信號的同時也將接收反射信號。反射信號經(jīng)由反射面反射后到達接收機。相對于直射信號而言，反射信號的傳播路徑有一個路徑延遲。北斗衛(wèi)星直射信號與反射信號的數(shù)字中頻表達式如下：
直射信號數(shù)學表達式：
式中，上標i表示衛(wèi)星編號，A、x、D、fd、τ、f和φi分別表示衛(wèi)星信號幅值、C/A碼、導航數(shù)據(jù)碼、直射多普勒頻移、反射相對于直射的碼片延遲、反射相對于直射的多普勒偏移和載波初始相位。
2.3.2 基于直射輔助的反射信號捕獲算法
GNSS-R接收機處理框圖如圖4所示。接收機基帶處理模塊分直射通道與反射通道。直射通道完成直射信號的捕獲跟蹤、導航解算與用戶定位。通過導航解算可求解出直射信號中的導航數(shù)據(jù)D(k)，D(k)取值為1或-1，使用D(k)剝離掉反射信號中的導航數(shù)據(jù)，則式(6)可改寫為：
去掉導航跳變的反射信號就可以進行長時間的相干積分。在使用FFT做相關處理前，對已剝離掉導航跳變的數(shù)據(jù)進行周期性累加，再執(zhí)行FFT相關操作，這樣可以降低FFT運算點數(shù)，從而極大降低長時間相干積分的運算量，提高反射信號捕獲速度。基于直射輔助的反射信號快速捕獲過程如下：
（1）以直射信號多普勒頻率為中心搜索頻率，對輸入的中頻北斗反射信號乘以正交本地載波，進行去載波操作。
（2）利用直射信號中提取的導航數(shù)據(jù)對反射信號進行導航數(shù)據(jù)剝離，根據(jù)北斗衛(wèi)星調制方式不同，分兩種情況：
①若為GEO衛(wèi)星，則導航數(shù)據(jù)直接乘以反射信號。
②若為MEO/IGSO衛(wèi)星，導航數(shù)據(jù)先調制NH碼，再與反射信號相乘。
（3）對不含導航跳變的信號以一個C/A周期進行累加操作。
（4）累加后信號進行FFT變換，與本地C/A的FFT變化的復共軛相乘，做IFFT變換之后取模值。
3 算法實驗仿真
3.1 實驗條件
實驗選用矩陣電子的多星座導航信號模擬器（型號為GNS-8332）產(chǎn)生衛(wèi)星信號。通過配置通道參數(shù)生成不同功率的北斗B1頻點衛(wèi)星信號，仿真場景中可以設置給衛(wèi)星加入多徑信號用以模擬該顆星的反射信號。反射信號相對于直射信號的功率衰減可以進行設置。采用萊特信息科技的多天線衛(wèi)星中頻信號采樣器（型號為LT-C-002）進行數(shù)據(jù)采集；采樣頻率為20 MHz。信號中心頻率為2.902 MHz。
3.2 算法仿真
3.2.1 實驗一
模擬器產(chǎn)生北斗1號衛(wèi)星信號，直射信號功率設為-130 dBm，反射信號相對直射信號衰減10 dB，即功率為-140 dBm。反射信號相對于直射信號的碼片延遲為20個碼片。3種算法對北斗1號星的捕獲結果分別如圖5～圖7所示。
實驗一仿真結果：3種捕獲方法所處理的數(shù)據(jù)長度均為20 ms。碼片分辨率均為1/4碼片。圖5～圖7所用捕獲方法分別為相干非相干積分、差分相干積分和直射輔助式相干積分，捕獲性能對比如表1所示。
3.2.2 實驗二
模擬器產(chǎn)生北斗7號衛(wèi)星信號，直射信號功率設為-130 dBm,反射信號相對直射信號衰減20 dB，即功率為-150 dBm。反射信號相對于直射信號的碼片延遲為20個碼片。3種算法對北斗7號星的捕獲結果分別如圖8～圖10所示。
實驗二仿真結果：3種捕獲方法所處理的數(shù)據(jù)長度均為100 ms。碼片分辨率均為1/4碼片。圖8～圖10所用捕獲方法分別為相干非相干積分、差分相干積分和直射輔助式相干積分，捕獲性能對比如表2所示。
4 結束語
本文在研究傳統(tǒng)反射信號捕獲算法的基礎上，利用直射與反射信號的關系提出一種適用于北斗GNSS-R接收機中的反射信號的捕獲算法，并進行了仿真，同時與相干非相干累積積分、差分相干算法做了對比。仿真結果表明，在捕獲-140 dBm信號時，3種算法的積分長度均為20 ms，本文算法的捕獲性能優(yōu)于相干非相干算法13.98 dB，優(yōu)于差分相干算法7.78 dB。捕獲-150 dBm的弱信號時，積分時間均為100 ms，傳統(tǒng)捕獲算法未能檢測到反射信號，本文算法正確檢測到反射信號。若采用更長的積分時間，則可以捕獲到更弱的反射信號。在利用GNSS-R技術進行目標反演時，該算法具有良好的應用前景。
參考文獻
[1] 李黃，夏青，尹聰，等.我國GNSS-R遙感技術的研究現(xiàn)狀與未來發(fā)展趨勢[J].雷達學報，2013，2(4)：389-399.
[2] CARDELLACH E，RUFFINI G，PINO D，et al.Mediterranean balloon experiment：ocean wind speed sensing from the stratosphere，using GPS reflections[J].Remote Sensing of Environment，2003，88(3)：351-362.
[3] RODRIGUEZ-ALVAREZ N，CAMPS A，VALL-LLOSSERA M，et al.Land geophysical parameters retrieval using the interference pattern GNSS-R technique[J].IEEE Transactions on Geoscience & Remote Sensing，2010，49(1)：71-84.
[4] 周旋，葉小敏，于暘，等.基于GNSS-R的海面風速探測技術研究[J].電子與信息學報，2013，35(7)：1575-1580.
[5] 楊東凱，張其善.GNSS反射信號處理基礎與實踐[M].北京：電子工業(yè)出版社，2012.
[6] ELDERS-BOLL H，DETTMAR U.Efficient differentially coherent code/Doppler acquisition of weak GPS signals[C].IEEE Eighth International Symposium on Spread Spectrum Techniques and Applications.IEEE，2005：731-735.
[7] 陳學，徐建華，葉甜春，等.基于FFT的微弱GPS信號捕獲算法[J].微電子學與計算機，2010，27(3)：98-101.
[8] 中國衛(wèi)星導航系統(tǒng)管理辦公室.北斗衛(wèi)星導航系統(tǒng)信號接口控制文件公開服務信號（2.1版）[S].2013.
[9] 魯郁.北斗/GPS雙模軟件接收機原理與實現(xiàn)技術[M].北京：電子工業(yè)出版社，2016.
[10] 徐曉蘇，范金洋.GPS軟件接收機捕獲算法[J].中國慣性技術學報，2009，17(2)：165-169.
[11] LIN D M，TSUI J B Y，LIOU L L，et al.Sensitivity limit of a stand-alone GPS receiver and an acquisition method[C].Proceedings of International Technical Meeting of the Satellite Division of the Institute of Navigation，2002：1663-1667.
[12] SCHMID A，NEUBAUER A.Performance evaluation of differential correlation for single shot measurement positioning[C].Proc.ION GNSS 2004，Long Beach，CA，2004：1998-2009.
作者信息:
楊 銳1，黃海生1，李 鑫1，曹新亮2
（1.西安郵電大學 電子工程學院，陜西 西安710121；2.延安大學 物理學與電子信息學院，陜西 延安716000）