陶瓊+孟李林+黃海生+李鑫+魏婷+楊銳+邢永強

摘 要： 針對目前導(dǎo)航衛(wèi)星應(yīng)用領(lǐng)域中一個新的熱點研究內(nèi)容，利用GNSS衛(wèi)星反射信號作為無源雷達進行目標探測和定位，基于GNSS?R技術(shù)進行車輛探測定位實驗研究，并給出反射信號相對于直射信號的延遲時間求解方法。實驗對接收到的直射信號進行處理，提取直射信號信息作為先驗條件，進行反射信號的細捕獲，縮小反射信號捕獲的搜索范圍，采用基于FFT的并行碼相位捕獲算法將碼相位延遲時間精確到碼片。將實驗數(shù)據(jù)與理論值進行對比分析，驗證了該方法的正確性。

關(guān)鍵詞： GNSS； 車輛探測； 反射信號； 直射信號； 相位捕獲算法； 相位延遲

中圖分類號： TN967?34； P228.4 文獻標識碼： A 文章編號： 1004?373X（2018）05?0106?04

Abstract： The reflected signal of global navigation satellite system （GNSS） is used to perform the target detection and location for the passive radar. The vehicle detection and localization experiment based on GNSS?R technology is studied. The delay time solving method of reflected signal relative to direct signal is given. The received direct signal is processed in the experiment. The direct signal information is extracted as the prior condition to carry out the fine capture of the reflected signal， so as to narrow the search range of the reflected signal capture. The FFT?based parallel code phase capture algorithm can correct the code phase delay time to chip. The experimental data is compared with the theoretical values. The correctness of the method was verified.

Keywords： global navigation satellite system； vehicle detection； reflected signal； direct signal； phase capture algorithm； phase delay

0 引 言

基于GNSS?R（Global Navigation Satellite System Reflections）目標探測技術(shù)是一種新興的遙感技術(shù)，其本質(zhì)是利用傳統(tǒng)導(dǎo)航定位系統(tǒng)中認為是“有害的”多徑信號，遙感地表參數(shù)、探測外來目標等[1]。四大衛(wèi)星導(dǎo)航系統(tǒng)在L波段的GPS L1，GPS L2，GPS L5，BDS B1，BDS B2，BDS B3，GLONASS L1以及伽利略L1等多個頻點上發(fā)射導(dǎo)航服務(wù)信號，可作為輻射源。該技術(shù)具有成本低、覆蓋寬、全天候等眾多優(yōu)勢，在海面測高[2]、海面測風[3]、積雪厚度探測、陸地濕度[4?5]、海洋和陸地表面成像、空間飛行器以及海洋或陸地移動目標探測[6?7]等方面都展現(xiàn)出應(yīng)用前景。

利用目標反射的GNSS信號，探測目標所在的空間位置和運動變化軌跡的雷達稱為非合作式“雙基地”（或“多基地”）無源雷達。該技術(shù)可廣泛應(yīng)用于探測戰(zhàn)斗機、導(dǎo)彈、無人機等移動目標的出現(xiàn)和運動變化現(xiàn)象中，以及對陸地和海洋上移動目標的探測和識別[8]。GNSS?R探測儀同時接收來自衛(wèi)星發(fā)射的直射信號和從目標物體上反射的衛(wèi)星信號，并進行信號處理，通過對比直射信號與反射信號二維相關(guān)峰值的時間延遲，推算目標物體與接收機之間的距離，從而確定被測物體的空間位置。計算反射信號相對于直射信號的時間延遲是利用GNSS?R探測目標的核心技術(shù)。

本文研究GNSS的目標物體回波特性，利用回波進行目標定位。具體實驗是用固定車輛的GPS反射信號計算汽車的空間位置。

1 GNSS?R目標探測原理

在地球上任何地方任何時刻都能同時觀測到多顆導(dǎo)航衛(wèi)星，因此，可利用一臺接收機接收多顆衛(wèi)星的反射信號，就可以實現(xiàn)對反射目標的定位?；贕NSS反射信號的目標探測系統(tǒng)原理如圖1所示。接收機同時接收直射和反射兩路信號，在圖1中分別用虛線和實線標出。

通過對直射信號和反射信號的分析處理，獲得反射信號與直射信號的傳播路徑差，結(jié)合衛(wèi)星及接收機的位置，計算反射目標的位置坐標。具體做法如下：

在地心直角坐標系下，假設(shè)滿足條件的導(dǎo)航衛(wèi)星的坐標分別為；目標位置坐標為；接收機的位置坐標為。在圖1中，對于1號導(dǎo)航衛(wèi)星，其反射信號相對于直射信號到達接收機的路徑差為：

式中：導(dǎo)航衛(wèi)星與接收機位置可通過接收機直射信號的定位解算模塊求解得到，c為真空中光速；為反射信號相對于直射信號到達接收機的時間延遲，可將其轉(zhuǎn)化為接收機直射信號與反射信號碼相位延遲時間（以下簡稱碼延遲時間）。同理，可得到滿足條件的其他3顆衛(wèi)星的路徑差方程，綜合上述方程建立方程組，便可解算出探測目標的位置坐標。

2 碼延遲時間的求解

由于直射信號與反射信號之間存在波程差，二者之間必然存在碼相位延遲，則兩路信號的二維相關(guān)功率如圖2所示。設(shè)導(dǎo)航信號碼周期為對于GPS信號而言，一個碼周期有1 023個碼片，則有如下關(guān)系式：endprint

式中：為反射信號相對于直射信號的延遲碼片數(shù)（單位：個）；與式（1）中的意義相同，為碼延遲時間。求碼延遲時間的問題可轉(zhuǎn)換成信號碼相位的求解問題。

2.1 碼延遲求解過程

GPS?R測量碼延遲時間的電路系統(tǒng)結(jié)構(gòu)圖如圖3所示。硬件部分由天線和中頻信號采集器組成，直射天線為右旋圓極化（RHCP）天線，為普通的商業(yè)GPS接收機天線，用來接收直射信號。反射天線為左旋圓極化（LHCP）天線，具有10 dB增益，指向反射目標，接收經(jīng)目標物體反射的反射信號。GN2012為GPS中頻信號采集器，對天線接收到的模擬信號進行降頻并采樣，得到數(shù)字中頻信號。將直射信號存儲在通道A，反射信號存儲在通道B，用圖3中所示數(shù)據(jù)格式交替存儲。

軟件部分為信號處理單元，完成直射信號和反射信號的數(shù)據(jù)處理，求出碼延遲時間。通過設(shè)置選取A路或B路信號進行分析，直射信號的處理與常規(guī)GPS接收機相同，包括捕獲、跟蹤和定位解算。將獲得的衛(wèi)星號、碼相位和多普勒頻率作為先驗條件，送入反射信號估算模塊，對反射信號的碼相位及多普勒頻移進行估算，設(shè)置搜索范圍和搜索步長，進入反射信號細測量模塊，進行精度小于一個碼片的精細捕獲，求得反射信號碼相位。最后將二者碼相位進行比較，得到延遲碼片數(shù)由公式（2）求得

2.2 信息估算

信息估算模塊主要包括載波多普勒估算與碼相位估算。對于靜止接收機且反射面為陸地上的車輛而言，其運動速度相比于導(dǎo)航衛(wèi)星運動速度來說很小，其反射信號與直射信號的多普勒頻移在同一數(shù)量級上，因此直接用直射信號載波NCO的輸出作為反射信號的多普勒頻移[9]。反射信號本地C/A碼相位的估算為：

式中：為先驗信息中直射信號碼相位；為接收機高度；為衛(wèi)星高度角；為一個C/A碼片的持續(xù)時間。

2.3 反射信號細測量

對于反射信號的處理如圖4所示，在基于FFT的并行碼相位捕獲算法[9]的基礎(chǔ)上進行改進，將估算模塊獲得的碼相位和載波相位分別送入可調(diào)本地載波和本地C/A碼產(chǎn)生模塊，輸入數(shù)據(jù)與可調(diào)本地載波混頻，通過一個由16點累加器實現(xiàn)的低通濾波器，濾除和頻分量，實現(xiàn)載波的剝離；而與C/A碼的相關(guān)運算，采用快速傅里葉變換實現(xiàn)，將剝離載波后的信號做FFT并取共軛，與本地C/A碼做FFT后的結(jié)果相乘，再取IFFT得到反射信號的相關(guān)值波形；記錄相關(guān)峰值及相關(guān)峰所在的碼片位置。由于對碼相位分辨率的要求，此處將本地C/A碼按位擴展16倍，以達到碼延遲時間精確到碼片的目的。

3 車輛探測實驗

探測對象為校園空曠處停放的車輛，數(shù)據(jù)采集環(huán)境如圖5所示，主要包括左旋天線、右旋天線、GN2012中頻信號采樣器及筆記本電腦。圖中直射天線為低增益RHCP天線，指向天頂方向；反射天線為四陣列高增益LHCP天線，指向車輛方向；GN2012中頻信號采集器輸出數(shù)據(jù)是采樣率為16.369 MHz，信號中心頻率為4.123 MHz，2 bit量化的數(shù)字中頻信號；筆記本電腦裝有中頻信號采集軟件，進行數(shù)字中頻信號的傳送和保存。實驗中，天線在車輛正南面，與接收機直線距離為35 m，將實驗計算結(jié)果與此對比，驗證該方法的正確性。

GPS信號中C/A碼長度為1 023個碼片，周期為1 ms，則一個C/A碼周期所對應(yīng)的距離約為：

一個碼片對應(yīng)的距離為 m。在2.3節(jié)的反射信號細捕獲算法中，將每一個碼片擴展了16倍，即采樣頻率為16.368 MHz，由此推算出理論測距精度約為 m。

先選4 ms直射信號數(shù)據(jù)進行處理，如圖6所示，結(jié)果給出數(shù)據(jù)采集時刻接收機捕獲到的直射信號的衛(wèi)星號、碼相位及多普勒頻率的粗略值。圖7為捕獲的衛(wèi)星星空視圖，同心圓為等仰角線，從外到內(nèi)分別代表0°，45°，90°。根據(jù)實驗場景設(shè)置，低仰角衛(wèi)星的反射信號會更強[10]，天線朝向偏西北方向，因此，應(yīng)選取仰角為10.30°的6號衛(wèi)星進行碼延遲時間的求解。

由圖6可知，6號衛(wèi)星多普勒頻移為2 600 Hz，碼相位為第32個碼片，直射信號跟蹤結(jié)果顯示載波多普勒頻移為2 597 Hz。由于反射信號強度較弱，這里采用增加相干積分時間的方法，而GPS導(dǎo)航電文數(shù)據(jù)比特寬度為20 ms，因此做反射信號細捕獲時，截取10 ms數(shù)據(jù)，以保證在連續(xù)兩段數(shù)據(jù)中至少有一段數(shù)據(jù)無導(dǎo)航數(shù)據(jù)跳[11]。根據(jù)估算，碼延遲時間在一個碼片范圍內(nèi)，因此以第32號碼片為起點，搜索步進為碼片，向后1個碼片，按2.3節(jié)所述進行搜索。將反射信號相關(guān)峰值與直射信號進行對比，結(jié)果如圖8所示。

從中選取某一頻率點如2 600 Hz，對直射和反射信號碼相位相關(guān)峰值做歸一化處理并進行對比，結(jié)果如圖9所示。

從圖9中可以看出6號衛(wèi)星反射信號相比于直射信號而言，有較大衰減，且相關(guān)峰延遲個碼片，換算成延后距離約為55 m，理論路徑延遲為m，誤差13.87 m，在測量精度范圍之內(nèi)。

4 結(jié) 語

為探究用GNSS衛(wèi)星反射信號作為無源雷達進行目標探測和定位的可行性，本文以GPS信號為例，進行基于GNSS?R技術(shù)車輛探測定位實驗。實驗給出反射信號相對于直射信號延遲時間的求解方法，用直射信號信息作為先驗條件，進行反射信號細捕獲，縮小了反射信號捕獲的搜索范圍，減少了捕獲所需的時間。采用基于FFT的并行碼相位捕獲算法，將碼相位延遲時間精確到碼片，即18.33 m。但實驗場地仍存在地面等其他反射面對實驗結(jié)果的影響，本次實驗并未考慮這些因素，實驗場地其他多徑信號對探測目標反射信號的影響將是今后的研究方向。

參考文獻

[1] 楊東凱，張其善.GNSS反射信號處理基礎(chǔ)與實踐[M].北京：電子工業(yè)出版社，2012.

YANG Dongkai， ZHANG Qishan. GNSS reflected signal proces?sing： fundamentals and applications [M]. Beijing： Publishing House of Electronics Industry， 2012.endprint

[2] 邵連軍，張訓械，劉經(jīng)南，等.GNSS?R海面測高算法[J].海洋測繪，2010（2）：1?3.

SHAO Lianjun， ZHANG Xunxie， LIU Jingnan， et al. The algorithm of sea surface altimetry on GNSS?R [J]. Hydroaphic surveying and charting， 2010（2）： 1?3.

[3] LOFGREN J S， HASSR， JOHANSSON J M. Monitoring coastal sea level using reflected GNSS signals [J]. Advances in space research， 2011（47）： 213?220.

[4] CARDELLACH E， FABRA F， NODUES?CORREIG O， et al. GNSS?R ground?based and airborne campaigns for ocean， land， ice and snow techniques： application to the GOLD?RTR datasets [J]. Radio science， 2011， 46（6）： 3604?3629.

[5] CHEW C， SMALL E E， LARSONK M. An algorithm for soil moisture estimation using GPS?interferometric reflectometry for bare and vegetated soil [J]. GPS solutions， 2016， 20（3）： 525?537.

[6] SUBERVIOLA I， MAYORDOMO I， MENDIZABAL J. Experimental results of air target detection with a GPS forward?scattering radar [J]. IEEE geoscience and remote sensing letters， 2012， 9（9）： 47?51.

[7] 裴振寧，葉偉，曲衛(wèi)，等.基于GPS信號的無源雷達性能分析[J].雷達與對抗，2012（2）：9?12.

PEI Zhenning， YE Wei， QU Wei， et al. Performance analysis of passive radar based on GPS signal [J]. Radar & ecm.， 2012（2）： 9?12.

[8] CLARIZIA M P， BRACA P， RUF C S， et al. Target detection using GPS signals of opportunity [C]// 18th International Confe?rence on Information Fusion. Washington， DC： IEEE， 2015： 1429?1436.

[9] 騰小霽，孫越強，吳迪，等.GNSS?R接收機信號捕獲算法設(shè)計[J].科學技術(shù)與工程，2010（10）：2278?2283.

TENG Xiaoji， SUN Yueqiang， WU Di， et al. Signal acquisition algorithm design for GNSS?R receiver [J]. Science technology and engineering， 2010（10）： 2278?2283.

[10] 馬小東.GNSS?R反射信號特征分析及仿真[D].北京：北京化工大學，2013.

MA Xiaodong. GNSS?R reflected signal feature analysis and application simulation [D]. Beijing： Beijing University of Chemical Technology， 2013.

[11] 黃海生，張偉，趙煥煥.GPS并行頻率域捕獲算法分析[J].西安郵電大學學報，2016（3）：38?41.

HUANG Haisheng， ZHANG Wei， ZHAO Huanhuan. Analysis on GPS parallel frequency acquisition algorithm [J]. Journal of Xian University of Posts and Telecommunications， 2016（3）： 38?41.

[12] 謝鋼.GPS原理與接收機設(shè)計[M].北京：電子工業(yè)出版社，2009.

XIE Gang. Principles of GPS and receiver design [M]. Beijing： Publishing House of Electronics Industry， 2009.endprint