代培培，王騰飛，姚 錚，邢建平

GNSS接收機晶振頻率穩(wěn)定度仿真分析

代培培1，王騰飛2,3，姚 錚2,3，邢建平1

（1. 山東大學 微電子學院，濟南 250100；2. 清華大學 北京信息科學與技術(shù)國家研究中心，北京 100084；3. 清華大學 電子工程系，北京 100084）

針對全球衛(wèi)星導航系統(tǒng)（GNSS）提供的高精度定位、導航及授時（PNT）服務以精密測量的傳輸時間為基礎，然而GNSS接收機中普遍采用的頻率源為穩(wěn)定度較差的晶振的問題，提出一種對GNSS接收機晶振頻率穩(wěn)定度進行分析和建模的方法：根據(jù)GNSS接收機晶振頻率穩(wěn)定度受隨機頻率誤差及老化率的影響的原理，給出晶振隨機頻率誤差及老化率的數(shù)學模型；并根據(jù)不同隨機頻率誤差成分仿真相應的各項晶振隨機頻率誤差序列；同時基于LAG1自相關法辨識不同時期主導的隨機頻率誤差類型，定性分析其對晶振頻率穩(wěn)定度的影響。實驗結(jié)果表明，提出的方法可為GNSS接收機晶振設計及頻率穩(wěn)定度能評估提供參考。

全球衛(wèi)星導航系統(tǒng)（GNSS）接收機；晶振；隨機頻率誤差；老化率；阿倫（Allan）方差；LAG1自相關

0 引言

在高精度導航定位授時中，精確的時間是十分重要的。如今在全球衛(wèi)星導航系統(tǒng)中，為了能夠給用戶提供高精度的定位、導航及授時（positioning，navigation and timing，PNT）服務，衛(wèi)星端均搭載了穩(wěn)定度更高的高精度原子鐘，比如銣鐘、銫鐘或氫鐘[1]。然而，用戶端的導航定位授時性能依賴于全球衛(wèi)星導航系統(tǒng)（global navigation satellite system，GNSS）接收機對衛(wèi)星信號的高精度觀測[2-3]，因此，接收機的時鐘穩(wěn)定度有著至關重要的作用。GNSS接收機時鐘在滿足精度要求的同時，還須兼顧體積、功耗和成本。由于高精度原子鐘的價格較高，并且對運行環(huán)境要求極高，因此具備低功耗、低價格及體積更小等優(yōu)勢的晶振獲得了廣泛的關注[4]。目前大多數(shù)GNSS接收機的時鐘為晶振，常用于校準和頻率傳遞；然而晶振頻率穩(wěn)定度變化會導致接收機對衛(wèi)星信號的測量產(chǎn)生誤差[5]，從而疊加在偽距和載波相位的測量值上。其中，隨機頻率誤差[6-7]及老化率[8-9]是引起晶振頻率穩(wěn)定度變化的重要因素。隨機頻率誤差主要影響晶振的短期頻率穩(wěn)定度，老化率主要影響晶振的長期頻率穩(wěn)定度。因其誤差特性十分復雜，很難進行實時精確建模來消除對GNSS接收機服務精度的影響。因此，研究基于GNSS接收機晶振的誤差特性及評估分析其對晶振頻率穩(wěn)定度的影響具有十分重要的意義。

本文通過不同隨機頻率誤差成分仿真相應的各項晶振隨機頻率誤差序列，同時基于LAG1自相關法辨識不同時期主導的隨機頻率誤差類型，從而定性分析其對晶振頻率穩(wěn)定度的影響，為GNSS接收機性能評估提供參考。首先闡述晶振頻率穩(wěn)定度的分析方法及隨機頻率誤差特性；然后針對晶振的隨機頻率誤差，即5種獨立的冪率噪聲建立數(shù)學模型，并采用LAG1自相關法對冪率噪聲進行識別，仿真分析其對晶振頻率穩(wěn)定度的影響；最后針對晶振的老化率，采用可更準確擬合晶振老化率的經(jīng)典老化率對數(shù)模型，定性分析老化率對晶振頻率穩(wěn)定度的影響。

1 晶振頻率穩(wěn)定度分析方法

1.1 晶振信號的數(shù)學模型

晶振的輸出信號一般為正弦波，可表示為[10-11]

為方便分析晶振頻率穩(wěn)定度，建立晶振模型的輸出頻率數(shù)據(jù)隨機數(shù)，采用隨機生成的方法仿真晶振輸出的頻率數(shù)據(jù)，此數(shù)據(jù)未添加任何影響[12]。如圖1所示的仿真晶振輸出頻率數(shù)據(jù)的時間-頻率圖可以看出晶振頻率在標稱頻率附近波動，因此能夠仿真包含噪聲的晶振輸出頻率。

圖1 仿真晶振輸出頻率數(shù)據(jù)的時間-頻率

1.2 穩(wěn)定度分析方法

2 晶振隨機頻率誤差數(shù)學模型及分析方法

2.1 晶振隨機頻率誤差數(shù)學模型

晶振的隨機頻率誤差為冪率噪聲，主要由白調(diào)相噪聲、白調(diào)頻噪聲、閃爍調(diào)相噪聲、閃爍調(diào)頻噪聲以及隨機游走調(diào)頻噪聲5種獨立的噪聲組成。美國學者大衛(wèi)·萊森（David Lesson）給出了經(jīng)典的冪率噪聲模型[14]，可表示為

因此，晶振噪聲頻譜分量與Allan方差之間的轉(zhuǎn)換可表示為[13]

結(jié)合式（5）和式（6），可以得到噪聲頻譜密度的Allan方差的時域，表示為

根據(jù)前面所提到的5種噪聲之間是相互獨立的隨機過程，晶振噪聲的Allan方差可表示為[15]

2.2 晶振隨機頻率誤差仿真方法

本文根據(jù)典型時間間隔的Allan方差值，仿真晶振5種冪率噪聲頻率序列。其中，白調(diào)相噪聲、白調(diào)頻噪聲、隨機游走噪聲可采用蒙特卡洛仿真方法，閃爍噪聲為1/噪聲[16]，其功率譜為冪函數(shù)，可通過利用頻域波形，將白噪聲整形成為1/噪聲，逆變換生成時域1/噪聲序列。

由式（10）可得頻率數(shù)據(jù)序列為

由式（10）可得相位數(shù)據(jù)序列為

2.3 晶振噪聲類型識別方法

3 實驗與結(jié)果分析

3.1 晶振隨機頻率誤差仿真分析

本文仿真的對象為10 MHz TCXO晶振，根據(jù)文獻[15]提供的該晶振不同時間間隔的經(jīng)典Allan方差平方根值進行仿真分析，如表1所示。

表1 典型Allan方差平方根值

圖3 各項噪聲的頻率穩(wěn)定度比較

在實際應用中，更關注的是不同時期噪聲的識別；因此，本文采用LAG1自相關法通過估計其冪率指數(shù)對上述仿真噪聲進行分類識別。圖4給出了LAG1自相關法識別各項噪聲自相關函數(shù)圖以及相位數(shù)據(jù)的滯后散點圖，同時給出了估計的冪律指數(shù)值，可通過冪律噪聲系數(shù)的識別區(qū)間得到相應時期的主導噪聲類型。其中，相位數(shù)據(jù)的滯后散點圖為表示數(shù)據(jù)相關程度的另一種方式，散點的線性擬合斜率和LAG1密切相關。滯后散點圖是白噪聲的一組不相關點，對于相關性更強的閃爍噪聲和隨機游動噪聲表現(xiàn)出線性特征。當滯后間隔等于或接近周期的1/4時，具有正弦分量的噪聲呈現(xiàn)出圓形分布特征。

3.2 老化率對晶振頻率穩(wěn)定度的影響分析

在中短期影響晶振頻率穩(wěn)定度的主要因素是冪率噪聲；對于晶振長期頻率穩(wěn)定度的影響占主導地位的是晶振的老化率，其會造成晶振輸出頻率產(chǎn)生緩慢漂移。Filter.R.L[9]通過大量的實驗研究，發(fā)現(xiàn)晶振的老化率隨時間的增加而減??；因此，線性模型不能準確表征晶振的老化率。同時，晶振的老化率以對數(shù)曲線形式呈現(xiàn)，給出了對數(shù)函數(shù)模型，能夠更準確地擬合晶振老化率的實際情況。Viq.J.R[8]給出了經(jīng)典老化率數(shù)學模型，如圖5所示。

圖5 經(jīng)典的晶振老化率數(shù)學模型

采用隨機生成的方法仿真受老化率影響的晶振輸出頻率數(shù)據(jù)。由如圖6（a）所示的老化率影響的晶振輸出時間-頻率圖可以看出，晶振頻率前期在標稱頻率附近波動，隨著時間的增加，老化率影響逐漸顯現(xiàn)，呈現(xiàn)出不太標準的正態(tài)分布。由如圖6（b）所示的老化率影響的晶振輸出頻率直方圖可知，該方法能夠仿真受老化率影響的晶振輸出頻率數(shù)據(jù)。

基于上述晶振輸出頻率數(shù)據(jù)，分析其頻率穩(wěn)定度，如圖7所示給出了老化率影響下晶振頻率穩(wěn)定度結(jié)果。從圖中可以看出，前期晶振老化率隨著時間的增加而減小，在后期又逐漸增大并趨于平穩(wěn)，符合實際晶振老化率的情況。

圖6 仿真老化率影響的晶振輸出頻率數(shù)據(jù)的時間-頻率及頻率直方圖

圖7 老化率影響下晶振頻率穩(wěn)定度結(jié)果

4 結(jié)束語

本文針對隨機頻率誤差及老化率對GNSS接收機晶振頻率穩(wěn)定度的影響進行仿真分析。首先推導了晶振的隨機頻率誤差數(shù)學模型，通過對晶振隨機頻率誤差的仿真分析，得出晶振頻率偏移程度隨著采樣間隔的增大而減小，并且隨著采樣時間的增加，在不同時期主導噪聲影響不同的結(jié)論：在中期，隨機游走調(diào)頻噪聲和閃爍調(diào)頻噪聲為主導噪聲；在短期，隨機游走噪聲和白調(diào)相噪聲為主導噪聲。其次，本文采用LAG1自相關法通過估計其冪率指數(shù)對上述仿真噪聲進行分類識別，驗證了晶振隨機頻率誤差數(shù)學模型的正確性。最后，采用經(jīng)典老化率對數(shù)數(shù)學模型，通過仿真分析老化率對晶振頻率穩(wěn)定度的影響，得出晶振的老化率隨時間的增加而減小最后趨于平緩，符合實際情況的結(jié)果。因此，通過仿真分離老化率及隨機頻率誤差對晶振頻率穩(wěn)定度的影響，可以為GNSS接收機晶振設計及頻率穩(wěn)定度評估提供參考。

[1] 郭海榮. 導航衛(wèi)星原子鐘時頻特性分析理論與方法研究[D]. 鄭州: 信息工程大學, 2006.

[2] 黃觀文. GNSS星載原子鐘質(zhì)量評價及精密鐘差算法研究[D]. 西安: 長安大學, 2012.

[3] 王嘉琛, 樊多盛, 肖波, 等. 基于晶振的遠程時間同步參考源性能分析[J]. 時間頻率學報, 2020, 43(4): 270-278.

[4] 彭慧麗, 陳金和, 韓艷菊, 等. 高頻高穩(wěn)恒溫晶體振蕩器設計[J]. 宇航計測技術(shù), 2020, 40(1): 46-50．

[5] Allan D W, Levine J. A historical perspective on the development of the Allan variances and their strengths and weaknesses[J]. IEEE Transactions on Ultrasonics Ferroelectrics＆ Frequency Control, 2016, 63 (4) : 513-519.

[6] Rutman J. Characterization of phase and frequency instabilities in precision frequency sources: Fifteen years of progress[J]. Proceedings of the IEEE, 1978, 66(9): 1048-1075.

[7] Vig J R. Quartz crystal resonators and oscillators for frequency control and timing applications: A tutorial[J]. Nasa Sti/recon Technical Report N, 1994, 95: 19519.

[8] Vig J R. The aging of bulk acoustic wave resonators, filters and oscillators[C]// Proceedings of the 45th Annual Symposium on Frequency Control 1991. [S.l.]: IEEE, 1991: 77-101.

[9] Filter R L, Vig J R. Long-Term aging of oscillators[J]. IEEE Transactions on Ultrasonics, Ferroelectrics and Control, 1993, 40(4): 387-394.

[10] 孫江濤. 基于GPS的1PPS的二級頻標馴服技術(shù)[D]. 西安:西安電子科技大學, 2010.

[11] 魯剛. 基于GNSS的晶振馴服與保持技術(shù)研究[D]. 北京: 中國科學院大學(中國科學院國家授時中心), 2021.

[12] 賈麗娜, 彭慧麗, 于德江, 等. 高穩(wěn)晶振短期頻率穩(wěn)定度的仿真分析[J]. 宇航計測技術(shù), 2021, 41(3): 27-33.

[13] Riley W J. Handbook of frequency stability analysis[Z]. National Institute of Standards and Technology: Gaithersburg, MD, USA, 2008.

[14] Guo W. Estimations of frequency and its drift rate[J]. IEEE Transactions on Instrumentation and Measurement, 1997, 46(1): 79-82.

[15] 寇艷紅, 張其善. GPS接收機中晶振誤差的模擬方法[J]. 電子與信息學報, 2004, 26(8): 1319-1324.

[16] 馮愛明, 林敏. 利用幾何插值法合成1/f噪聲[J]. 中國計量學院學報, 2002, 13(2): 109-112.

[17] 孫宏偉. 原子鐘噪聲的蒙特卡絡模擬方法[J]. 時間頻率學報, 1996, 19(0): 33-37.

[18] Riley W J. User manual: Stable32 frequency stability analysis[Z]. Hamilton Technical Services, Beaufort, South Carolina, 2007.

Simulation analysis on frequency stability of crystal oscillator in GNSS receivers

DAI Peipei1, WANG Tengfei2,3, YAO Zheng2,3, XING Jianping1

（1. School of Microelectronics, Shandong University, Jinan 250100, China;2. Beijing National Research Center for Information Science and Technology, Tsinghua University, Beijing 100084, China;3. Department of Electronic Engineering, Tsinghua University, Beijing 100084, China）

Aiming at the problem that the high-precise positioning, navigation and timing (PNT) services provided by the global navigation satellite system (GNSS) are based on precisely measured transmission time, but the frequency source commonly used in the GNSS receivers is crystal oscillators with poor stability, the paper proposed an analysis and modelling method of the frequency stability of crystal oscillatorin GNSS receivers: according to the principle that the crystal frequency stability of GNSS receivers is affected by random frequency error and aging rate, the mathematical model of the random frequency error and the aging rate was given; and the corresponding random frequency error sequences of crystal oscillator were simulated according to different random frequency error components; meanwhile, the dominant random frequency error types in different periods were identified based on the LAG1 autocorrelation method, to qualitatively analyze their influence on the frequency stability of crystal oscillator. Experimental result showed that the proposed method could provide a reference for the design of crystal oscillator in GNSS receivers and the performance evaluation of the frequency stability.

global navigation satellite system (GNSS) receiver; crystal oscillator; random frequency error; aging rate; Allan variance; LAG1 autocorrelation

代培培, 王騰飛, 姚錚, 等. GNSS接收機晶振頻率穩(wěn)定度仿真分析[J]. 導航定位學報, 2023, 11(4): 90-96.（DAI Peipei, WANG Tengfei, YAO Zheng, et al. Simulation analysis on frequency stability of crystal oscillator in GNSS receivers[J]. Journal of Navigation and Positioning, 2023, 11(4): 90-96.）DOI:10.16547/j.cnki.10-1096.20230413.

P228

A

2095-4999(2023)04-0090-07

2022-10-24

國家重點研發(fā)計劃項目（2021YFA0716603）；國家自然科學基金面上項目（42274018）；第六屆中國科協(xié)青年人才托舉工程支持項目（YESS20200226）。

代培培（1993—），女，山東聊城人，博士研究生，研究方向為GNSS定位授時算法。