衛星導航長碼信號波形監測方法

劉建成，王 宇，宮 磊，徐曉燕

衛星導航長碼信號波形監測方法

劉建成，王 宇，宮 磊，徐曉燕

（北京衛星導航中心，北京100094）

由于信號微弱，如何獲得在軌導航衛星的清晰信號波形是衛星導航信號質量監測研究中的難點之一，為此提出了一種在軌導航衛星的信號波形監測方法。該方法基于Vernier采樣原理，利用大口徑拋物面天線對在軌衛星進行信號采集，經過消除初相和殘余頻率、累加平均和數據組合等處理，獲得清晰的碼片波形。對于相同碼速率的民用信號和長碼信號，可確定民用信號和長碼信號的偽碼相位偏差。利用大口徑拋物面天線對北斗衛星進行跟蹤，獲得了多顆北斗衛星B1頻點民用信號和長碼信號的碼片波形。結果表明，民用信號和長碼信號的碼片波形的輪廓差異較小，但偽碼相位存在偏差。

信號質量監測；長碼信號；波形監測；北斗衛星導航系統

引用格式：劉建成，王 宇，宮 磊，等.衛星導航長碼信號波形監測方法［J］.無線電工程，2016，46（5）：45-48.

0 引言

在軌導航衛星的信號發射通道故障會導致信號波形畸變［1］，而信號波形的畸變將影響接收機對信號的捕獲和跟蹤，因此需要對在軌衛星進行信號質量監測［2］。衛星導航信號質量監測項目包括：功率譜、眼圖、星座圖、相關峰、碼片波形、發播信號功率［3］、偽碼序列、I/Q支路載波相位正交性、I/Q支路偽碼相位一致性［4］、I/Q支路偽碼測量一致性［5］以及偽碼測量與載波相位測量一致性等。按評估方式可分為基于專家判讀的定性監測和基于自動判決的定量監測［6］。碼片波形監測是對導航信號的I支路碼片波形或Q支路碼片波形的監測，從時域上檢查信號是否發生畸變。通過比較實際接收信號碼片波形和理想信號碼片波形的寬度、上升沿和下降沿的差異等來判斷信號是否發生畸變。文獻［7］研究了基于清晰波形的多個評估指標的測試方法。

由于地面接收的導航信號載噪比低，碼速率高，因此普通的采集設備不足以描述波形細節，需要在數據處理階段對多個偽碼周期的采集數據進行平均。由于 GPS衛星民用信號一個信息碼內包含20個重復偽碼，文獻［8-10］提出了基于Vernier采樣原理的GPS民用信號采集方法和處理方法。采用大口徑天線跟蹤GPS衛星，實現對GPS衛星L1頻點C/A信號的采集，采樣率大小考慮了徑向速度影響。經過消除初相、殘余頻率后，對一個信息碼內的20個偽碼周期采集數據進行數據組合，再經累加減小噪聲影響后，獲得了多顆GPS健康衛星清晰的時域波形。北斗衛星民用信號一個信息碼內也包含20個重復偽碼，但偽碼周期又經過了二次編碼，一個信息碼內的偽碼符號不同，所以不能直接對一個信息碼內的 20個偽碼進行數據組合。文獻［11-12］提出了先對符號相同的偽碼進行累積平均，然后再數據組合，因此這種數據處理方法數據率提高倍數不受限制，同時選擇在徑向速度為零時進行信號采集，減小了徑向運動對碼片寬度的影響。

但對于長碼信號而言，在信號采集時間內擴頻碼是非周期性的，因此無法直接采用民用信號的處理方法實現多周期累積和提高數據率。民用信號的周期性表現為偽碼周期，對于長碼信號而言，只能利用碼片的周期性。在基于Vernier采樣原理進行數據采集的基礎上，先對符號相同的碼片進行累積平均，然后再數據組合。這種方法可以提高長碼信號的載噪比和數據率，獲得清晰的長碼信號的碼片波形。由于該方法也適用于民用信號，因此如果相同頻率的民用信號和長碼信號的碼速率相同，則可同時獲得2個支路的碼片波形，并根據碼片波形的起點差異，獲得I/Q支路偽碼相位偏差。

1 基于Vernier采樣的信號采集方法

對于周期重復性數據，如果每個周期的第1個采樣點相對周期起點的位置錯開固定間隔，如圖1所示，則經數據組合處理，可以把n個周期的數據合成為一個周期的數據，提高了數據率。

圖1　Vernier采樣原理的數據組合

若預提高數據率倍數為n，根據Vernier采樣原理確定采樣率為［8］：

在進行數據組合前，一般要對采集數據先消除初相和頻率殘差，然后對n個周期的數據進行累加以提高載噪比。

2 長碼信號采集方法

對民用信號而言，圖1中的一個周期對應的是一個偽碼周期，但對于長碼信號而言，圖1中的一個周期對應的是一個碼片寬度。由于相鄰碼片的符號不一定是相同的，因此不能直接進行數據組合處理。為改變這種特性，以n個周期的數據為平均處理周期，先對符號相同的采集數據進行累積平均，使每個周期的符號均相同。比如，第n+1個碼片的符號如果與第1個碼片的符號相同，則該周期進行平均，若符號不同則不做處理。完成累加平均后再進行數據組合，既提高了載噪比也提高了數據率。這樣的處理方法避開了由于長碼信號偽碼的周期太長導致無法直接采用民用信號的處理方法實現多周期累積和提高數據率的問題。

一般說來，衛星導航信號在相同的頻率上，I支路調制民用信號，Q支路調制長碼信號。如果民用信號和長碼信號的碼速率是相同的，則在同時獲得清晰碼片波形的基礎上，可獲得民用信號和長碼信號的偽碼相位偏差。

因此，長碼信號采集和數據處理步驟包括：

①根據式（1）確定采樣率，利用大口徑拋物面天線完成信號采集。

信號采集需采用大口徑拋物面天線，主要目的是進行接收信號放大和空域選擇，以減小非信號方向干擾和多徑的影響，同時接收鏈路采用帶通濾波器，在頻域上進行選擇，且濾波器帶寬與衛星載荷發射信號帶寬相匹配，以減小噪聲影響。為減小徑向運動對碼片寬度的影響，在徑向多普勒速度為0時實現信號采集。

②消除初相和殘余頻率，獲得民用信號和長碼信號的數據序列。

為消除采集數據的初相和殘余頻率，不斷調整頻率和相位，對采集數據進行反相位調制處理。然后查看數據的實部和虛部是否均表現為碼片形狀，且碼片寬度約為由碼速率決定的碼片寬度的倍數。調整頻率和相位直至滿足上述要求，反相位調制后，處理后數據的實部就是民用信號的數據序列，虛部是長碼信號的數據序列。

由于消除初相時存在π相位模糊，因此存在2個初相，分別對應衛星發播的原信號初相及原信號初相的反相，而不同初相獲得的是2組數據結果。由于信號采集時長不足以獲得電文幀同步頭，因此無法解決π相位模糊問題。下面對這2種數據分別進行處理，處理過程相同。

③確定碼片起始采樣點和每個碼片的符號。

根據民用信號數據序列的正負號變化粗略確定一個碼片的零交叉點位置，并把零交叉點后的第1個采樣點確定為該碼片的起始采樣點。該起始點同時作為長碼信號數據序列的起始點。根據起始采樣點延后半個碼片的數據的符號，確定為該碼片符號。若不考慮信息碼和二次編碼符號影響，則確定的該碼片符號就是偽碼的符號。

④對民用信號和長碼信號根據碼片符號分別進行累加平均。

單獨正碼片的特性是當前碼片符號為正，而前一個碼片和后一個碼片的符號均為負。按照符合該特性的碼片進行累加平均。

⑤對民用信號和長碼信號分別進行數據組合處理，獲得了清晰的碼片波形。

根據Vernier采樣原理，對n個碼片且長度為n·M-1的民用信號和長碼信號分別進行數據組合，分別形成一個碼片周期的數據序列。新的數據序列長度不變，但數據間隔縮小了n倍，可表示為：

式中，Tr為碼片寬度。

⑥確定民用信號和長碼信號不同的碼片起點位置，則2個起點位置差即為偽碼相位偏差。

分別對民用信號和長碼信號正碼片的上升沿部分進行線性擬合，令線性擬合獲得的多項式方程為零，解該方程獲得碼片起點位置。民用信號碼片和長碼信號碼片的起點位置差為偽碼相位偏差。

⑦分別處理π相位模糊的數據序列，獲得偽碼相位偏差。

3 北斗衛星B1頻點信號波形特性分析

大口徑拋物面天線分別對北斗衛星導航系統的1號、3號、4號、6號、7號和9號衛星進行跟蹤，利用具有矢量信號分析功能的Agilent E4445A在中頻實現了B1頻點信號的采集。

選擇n=50，Agilent E4445A最大分析帶寬為80 MHz，北斗衛星B1頻點碼速率為2.046 Mcps，根據式（1）計算采樣率為102.259 08 MHz。采集信號時長0.5 s，VSA 89600軟件存儲數據類型選擇mat類型，文件大小約為400 MB。

下面首先比較了文獻［11］方法和本文方法對民用信號的波形監測結果，然后對同一顆衛星采用新方法獲得的民用信號和長碼信號的波形進行了比較，最后比較了不同衛星的波形差異。

3.1 種處理方法對民用信號的波形監測結果比較

采用文獻［11］的采集的IGSO-1衛星數據，采樣率為102.399 95 MHz。采用文獻［11］的處理方法獲得了民用信號一個偽碼周期的碼片波形，選取其中一個正碼片，如圖2（a）所示。采用本文的采集數據和處理方法，獲得了民用信號的一個碼片波形，如圖2（b）所示。

圖2　2種處理方法比較

比較圖2波形可以看出，2種處理方法獲得的民用信號碼片波形的輪廓相同，但由于后者的累積次數更多所以碼片波形更清晰，因此本文的處理方法更具有優勢。

3.2 民用信號和長碼信號波形比較及偽碼相位偏差

采用本文方法對北斗衛星導航系統的1號、3號、4號、6號、7號和9號衛星的采集信號進行處理，獲得了民用信號和長碼信號碼片波形，如圖3所示。

圖3　民用信號和長碼信號碼片波形比較

由圖3可以看出，民用信號和長碼信號的碼片波形的輪廓差異較小，但存在2種波形間上升沿的超前或滯后現象，即偽碼相位存在偏差。

3.3 不同衛星的碼片波形比較

為比較不同衛星的波形差異，對本文處理方法獲得的多顆衛星的碼片波形做了19點移動平均處理，以改善噪聲影響。民用信號處理結果如圖4（a）所示，長碼信號處理結果如圖4（b）所示。由圖4可以看出，無論是民用信號還是長碼信號，不同衛星的碼片波形均有稍微差異。波形不同可能會導致測距偏差不同，進而影響測距精度。

圖4　不同衛星碼片波形比較

4 結束語

本文提出了一種在軌導航衛星的信號波形監測方法。利用大口徑拋物面天線對北斗衛星進行跟蹤，完成了多顆衛星的信號采集和數據處理，獲得了多顆北斗衛星B1頻點民用信號和長碼信號的碼片波形。下一步研究工作是進行足夠長時間的信號采集以解π相位模糊，進而獲得明確的民用信號和長碼信號的偽碼相位偏差。

［1］ 盧曉春，周鴻偉.GNSS空間信號質量分析方法研究［J］.中國科學：物理學力學天文學，2010，40（5）：528-533.

［2］ 劉建成，桑懷勝，徐 赟，等.衛星導航信號異常對信噪比的影響［J］.電訊技術，2013，53（1）：28-32.

［3］ 王 斌，龐 巖，劉會杰.導航信號有害波形檢測技術研究［J］.電子與信息學報，2011，33（7）：1 713-1 717.

［4］ MACABIAU C，CHATRE E.Impact of Evil Waveform on GBASPerformance［C］∥ ProceedingsoftheIEEE PLANS 2000，SanDiego California：IEEE PLANS，2000：22-29.

［5］ PHELTS R E.Multicorrelator Techniques for Robust Mitigation of Threats to GPS Signal Quality［D］.California：Standford University，2001.

［6］ 劉建成，溫日紅，王宏兵，等.衛星導航信號質量監測［J］.測繪科學與工程，2014（3）：43-46.

［7］ 楊再秀，郭曉峰，楊麗云.GNSS信號質量關鍵指標測試方法研究［J］.無線電工程，2015，45（6）：59-62.

［8］ MITELMAN A M.Signal Quality Monitoring for GPS Augmentation Systems［D］.California：Standford University，2004.

［9］ WONG Gabriel，PHELTS R E，WALTER Todd，et al. Bounding Errors Caused by Nominal GNSS Signal Deformations［C］∥Proceedings of the ION GNSS，Manassas，VA：ION GNSS，2011：2 657-2 664.

［10］PHELTS R E，WALTER Todd，ENGE Per.Characterizing Nominal Analog Signal Deformation on GNSS Signals ［C］∥Proceedings of the ION GNSS，Savannah，GA：ION GNSS 2009：1 343-1 350.

［11］LIU Jian-cheng，YANG Hua-feng，FAN Jian-jun，et al.Deformation Characteristics of COMPASS Satellite Signals ［C］∥China Satellite Navigation Conference（CSNC），Lecture Notes in Electrical Engineering，2014：171-180.

［12］LIU Jian-cheng，FAN Jian-jun，FENG Xiao-chao，et al. Ranging Bias of COMPASS Satellite Signals［C］∥China Satellite Navigation Conference（CSNC），Lecture Notes in Electrical Engineering，2015：347-355.

Long-code Signal Waveform Monitoring Method for Navigation Satellites

LIU Jian-cheng，WANG Yu，GONG Lei，XU Xiao-yan
（Beijing Satellite Navigation Center，Beijing 100094，China）

Due to the weakness of signal，signal waveform monitoring for navigation satellites in orbit is one of the difficulties in satellite navigation signal quality monitoring research，so a signal waveform monitoring method for navigation satellites in orbit is proposed.Based on the Vernier sampling principle，a large-diameter parabolic antenna is used for in-orbit satellite signal collection.After initial phase and residual frequency elimination，accumulation and combination，a clear chip waveform is obtained.For civilian and longcode signals with the same code rate，the PN code phase bias can be determined.By using a large-diameter parabolic antenna for COMPASS satellite tracking，the civilian and long-code chip waveforms of several COMPASS satellites in B1 band are obtained，and the PN code phase bias of the satellite signals are got.The results show that there is little difference between the civilian signal waveform and long-code signal waveform，but there is a code phase bias between them.

signal quality monitoring；long-code signal；waveform monitoring；COMPASS

TN96

A

1003-3106（2016）05-0045-04

10.3969/j.issn.1003-3106.2016.05.12

2016-01-06

劉建成 男，（1976—），博士，工程師。主要研究方向：衛星導航信號質量監測、衛星導航測量數據質量分析。

王 宇 男，（1975—），高級工程師。主要研究方向：衛星導航工程。