基于前導(dǎo)脈沖輔助的相位調(diào)制1090ES信號載頻同步算法

孫清清，王　洪，宋　妍，姒　強

(電子科技大學(xué) 電子工程學(xué)院,　成都 611731)

?

·信號處理·

基于前導(dǎo)脈沖輔助的相位調(diào)制1090ES信號載頻同步算法

孫清清，王洪，宋妍，姒強

(電子科技大學(xué) 電子工程學(xué)院,成都 611731)

摘要：基于相位調(diào)制擴容1090ES信號為ADS-B帶來新的功能和應(yīng)用，載波頻率同步是擴容信號需解決的一項關(guān)鍵技術(shù)。針對擴容1090ES信號的突發(fā)載頻同步問題，提出了一種基于頻偏捕獲和相位誤差跟蹤環(huán)的快速同步方法，頻偏捕獲采用前導(dǎo)輔助脈沖相關(guān)性判決進行二叉樹頻率搜索，捕獲載頻偏差并校正得到頻偏較小的基帶信號，相位誤差跟蹤環(huán)實現(xiàn)對小頻偏信號進行頻率估計和補償。仿真結(jié)果表明：該方法適用于基于相位調(diào)制的擴容1090ES信號接收系統(tǒng)。

關(guān)鍵詞：1090ES；相位調(diào)制；載頻同步；二叉樹搜索；相位跟蹤環(huán)

0引言

廣播式自動相關(guān)監(jiān)視(ADS-B)是下一代民航監(jiān)視系統(tǒng)的一項重要技術(shù)，在國內(nèi)外得到了廣泛應(yīng)用。隨著監(jiān)視技術(shù)的不斷發(fā)展和進步，未來的ADS-B將不再是執(zhí)行簡單的監(jiān)視任務(wù)，而將同時具備導(dǎo)航、監(jiān)視和通信的一體化功能，ADS-B信息的傳輸需求將不斷提高。為此，國際民航組織(ICAO)和歐控等權(quán)威機構(gòu)提出一種適用于新一代ADS-B技術(shù)的擴容1090ES數(shù)據(jù)鏈，該數(shù)據(jù)鏈在原有的1090ES數(shù)據(jù)鏈格式基礎(chǔ)上加入了多重相位調(diào)制(MPSK)，可在不增加頻點使用和廣播次數(shù)的前提下提高ADS-B信息的傳輸容量[1-4]。另外，擴容1090ES數(shù)據(jù)鏈也可一定程度上緩解多類設(shè)備使用1 090 MHz頻點造成的頻譜擁塞問題[5]?，F(xiàn)階段，擴容1090ES數(shù)據(jù)鏈正處于驗證階段，各項技術(shù)如調(diào)制、編碼、同步、解碼等都有待進一步研究和完善[6-7]。

擴容1090ES數(shù)據(jù)鏈載頻同步是難點。一般的1090ES信號解碼主要關(guān)注接收信號的幅度和位置，因此對載頻和相位的同步?jīng)]有嚴格的要求。擴容信號格式與現(xiàn)有1090ES信號的格式相同，前導(dǎo)脈沖主要用于信號檢測以及定時同步，加入相位調(diào)制后，提高了信號載頻、相位以及時間同步的精確度要求。ICAO規(guī)定1090ES發(fā)射信號的載頻為1 090 MHz±1 MHz，且機載信號多普勒頻移都可能達到幾百甚至幾千Hz，接收端載頻的同步必然很大程度上影響相位信息的解碼性能。另外，ADS-B地面接收機作用面臨的ADS-B發(fā)射機可能達到幾十到幾百部，擴容1090ES信號是一種特殊的單發(fā)多收的突發(fā)通信，難以采用常規(guī)的鎖相環(huán)實現(xiàn)同步。

現(xiàn)有針對具有前導(dǎo)數(shù)據(jù)輔助的突發(fā)同步方法有很多研究成果，典型的如：基于輔助數(shù)據(jù)的頻偏估計Kay 算法[8]、突發(fā)模式下數(shù)據(jù)輔助載波頻偏恢復(fù)Fitz算法[9]，以及改進的L&R 頻偏估計算法[10]。這些經(jīng)典算法在已知訓(xùn)練序列的輔助下，均能取得良好的同步效果，但是無法直接用于1090ES擴容信號的同步，原因在于1090ES信號前導(dǎo)四個脈沖是非連續(xù)調(diào)制的脈沖，即每個碼元之間的時間間隔不相同，另外前導(dǎo)輔助數(shù)據(jù)太少也將造成頻偏估計的誤差變大。

本文針對1090ES擴容信號載頻同步的特殊性和難點，提出一種新型的適用于1090ES擴容信號的載頻同步方法，主要分為兩步：頻偏捕獲和相位誤差跟蹤判決。頻偏捕獲首先基于前導(dǎo)輔助脈沖對信號載頻偏差的進行快速估計并校正，剩余的載頻偏差通過相位誤差跟蹤環(huán)路進行補償，提高了同步性能和解碼的準確率。

1系統(tǒng)模型

擴容1090ES信號使用一種復(fù)合調(diào)制方法，在現(xiàn)有的脈沖位置調(diào)制(PPM)高電平子脈沖上進行相位調(diào)制。圖1為對1090ES進行MPSK脈內(nèi)調(diào)相后的信號波形，四個前導(dǎo)脈沖調(diào)制為0相位作為參考脈沖，用于接收信號的相位和載頻同步。

圖1　相位調(diào)制1090ES信號

基于頻偏捕獲及誤差跟蹤環(huán)的擴容1090ES信號載頻同步系統(tǒng)模型如圖2所示。

圖2　擴容1090ES信號載頻同步系統(tǒng)模型

擴容1090ES信號接收端經(jīng)過混頻后基帶信號為x(n)，如式(1)

(1)

式中：Δfb為大頻率偏差，包含擴容信號經(jīng)過1 090 MHz混頻后剩余的載頻偏差及多普勒頻移；A為信號幅度；g(n)為1090ES信號的基帶幅度包絡(luò)，如圖1所示；2π·(m-1)/M為8PSK調(diào)制相位，M=8，m=1,2,…,8；θ0為載波相位。

基帶信號采用PPM的解碼信息進行時間同步，即由PPM信息的解碼結(jié)果判決每個數(shù)據(jù)位高電平子脈沖的位置，從而確定相位信息判決點位置，經(jīng)過頻偏捕獲及校正模塊后輸出的基帶信號為y(n)，如式(2)

(2)

式中：Δfs為經(jīng)過頻偏捕獲和校正后的小頻率偏差，為提高后續(xù)反饋判決的解碼性能，本文將小頻偏控制在5 kHz范圍內(nèi)。

相位誤差跟蹤環(huán)路通過跟蹤相位判決誤差對小頻偏信號進行同步并解碼，輸出得到擴容1090ES信號的相位信息。

2頻偏捕獲

頻率捕獲是載波頻率粗估計的一種方式，目的是將載頻的估計誤差控制在較小的范圍，以提高后續(xù)相位跟蹤判決的準確性。頻率捕獲應(yīng)用于載頻同步是較常見的，如鎖相環(huán)路實現(xiàn)載波同步的系統(tǒng)中，在環(huán)路中增加掃描捕獲輔助單元使用順序搜索進行頻率捕獲，目的是提高環(huán)路信噪比和縮短環(huán)路捕獲時間[11]；在Viterbi等非線性頻率估計方法完成之后使用頻率捕獲，目的是減小非線性算法估計的頻率誤差，提高系統(tǒng)性能[12]。

本文針對1090ES擴容信號的同步，設(shè)計了一種基于二叉樹折半搜索的頻率捕獲方法，主要的思想是：使用前導(dǎo)脈沖與一定數(shù)量的數(shù)據(jù)脈沖作為輔助判決序列，采用前導(dǎo)輔助脈沖相關(guān)性作為判決準則，通過二叉樹折半搜索原理進行迭代搜索，直到頻率捕獲精度滿足設(shè)定的要求。頻率捕獲算法的核心是前導(dǎo)輔助脈沖的相關(guān)性判決，前導(dǎo)輔助脈沖由1090ES信號前導(dǎo)脈沖與一定數(shù)量的數(shù)據(jù)脈沖組成，如采用四個數(shù)據(jù)脈沖與前導(dǎo)脈沖構(gòu)成的輔助脈沖如圖3所示，時間長度為12 μs，具體判決方法是本地產(chǎn)生特定中心頻率的輔助脈沖序列與接收信號的前12 μs進行相關(guān)，通過相關(guān)性大小判決與接收信號更接近的頻率。輔助脈沖長度加入數(shù)據(jù)脈沖的數(shù)量越多判決的準確率越高，但也損失了部分的數(shù)據(jù)信息，具體使用多少位數(shù)據(jù)脈沖加入需要進一步探討和權(quán)衡，具體見后面的仿真分析。另外，頻率捕獲算法的迭代次數(shù)應(yīng)在合理范圍，避免過長的捕獲時間導(dǎo)致突發(fā)信號同步性能下降。

圖3　前導(dǎo)輔助脈沖

頻偏捕獲基本原理可描述為：

假設(shè)：系統(tǒng)的最大頻偏為Δfmax(對于1090ES信號取值為±1 MHz)，要求的頻率捕獲精度為δf，則需要的迭代次數(shù)為N≈lb(Δfmax/δf)；第一次迭代采用中心頻率為-0.5 MHz和+0.5 MHz本地輔助判決序列與接收信號序列進行相關(guān)判決，第一次判決將捕獲誤差控制在±0.5 MHz范圍；如第一次判決為+0.5 MHz，則第二次迭代采用中心頻率為+0.25 MHz和+0.75 MHz本地輔助判決序列與接收信號序列進行相關(guān)判決，反之，則采用中心頻率為-0.75 MHz和-0.25 MHz本地輔助判決序列與接收信號序列進行相關(guān)判決，第二次判決將捕獲誤差控制在±0.25 MHz范圍；以此類推，直到捕獲誤差滿足要求，如要求將捕獲誤差控制在±5 kHz范圍，需要8次迭代搜索。為方便計算令Δfmax=1.024 MHz，δf=4 kHz，如實際頻偏為Δf=1.018 MHz，使用頻率捕獲方法，需經(jīng)過8次迭代搜索，得到捕獲的載頻偏差為1.020 MHz，如圖4所示。

圖4　二叉樹折半搜索頻率捕獲

頻偏捕獲可將載頻偏差控制在一個較小的范圍內(nèi)，具體需要控制的頻率范圍根據(jù)相位誤差跟蹤要求的精度決定。要求精度越高搜索迭代的次數(shù)也越多，這可能增加頻率捕獲時間影響同步實時性。針對該問題一般可采取并行計算的方式，即多個中心頻率序列與接收序列同時進行相關(guān)計算后判決，盡可能縮短頻率捕獲時間。

3相位誤差跟蹤原理

相位誤差跟蹤環(huán)路在頻偏捕獲和校正之后進行，采用的是2階科斯塔斯環(huán)(Costas)的原理，通過比較接收信號相位與理想判決相位產(chǎn)生相位誤差，誤差信號經(jīng)過環(huán)路濾波后進行前饋補償，流程如圖5所示。

圖5　相位誤差跟蹤環(huán)

(3)

μ(n+1)=μ(n)-γ(1+ρ)e(n+1)+γe(n)

(4)

(5)

相位誤差跟蹤環(huán)輸出了相位信息的判決結(jié)果，前提是在PPM脈沖位置已知條件下，即PPM解碼已經(jīng)正確完成，另外，前導(dǎo)輔助序列的參考相位(調(diào)制相位已知)起到相位同步作用，保證了前期判決的準確性。后文的仿真結(jié)果說明，相位誤差反饋判決方法在小頻偏情況下可以獲得良好的解碼性能。

4仿真實例

4.1頻偏捕獲算法仿真實驗

數(shù)據(jù)源說明：1090ES擴容信號載頻在1 089 MHz～1 091 MHz之間隨機產(chǎn)生，前導(dǎo)脈沖8 μs，數(shù)據(jù)脈沖112 μs，基帶數(shù)據(jù)隨機產(chǎn)生，信號脈沖上升沿和下降沿等參數(shù)符合ICAO標(biāo)準規(guī)范；擴容信號在原有1090ES信號高電平子脈沖上進行8PSK相位調(diào)制，前導(dǎo)輔助脈沖由8 μs前導(dǎo)脈沖和若干個數(shù)據(jù)脈沖構(gòu)成，前導(dǎo)輔助脈沖的相位調(diào)制為0，作為參考相位，擴容信號經(jīng)過本地1 090 MHz混頻后抽取到10 MSPS。PPM信息解調(diào)采用文獻[13]提出的基于簡化表格的多點判決法，并以解碼后的PPM信息作為時間同步基準。加入不同數(shù)據(jù)位作為前導(dǎo)輔助脈沖，隨機產(chǎn)生100 000次信號，分別進行頻偏捕獲算法仿真，偏差捕獲精度要求5 kHz，即迭代次數(shù)為8，統(tǒng)計仿真結(jié)果如圖6所示。data 0表示前導(dǎo)輔助脈沖不加入數(shù)據(jù)位，data 1表示加入1個數(shù)據(jù)位，data 7表示加入7個數(shù)據(jù)位，其他表示以此類推。

圖6　頻偏捕獲算法性能

加入數(shù)據(jù)位越多頻偏捕獲準確率越高，但使用數(shù)據(jù)位作為輔助脈沖將損失一部分有用信息，圖6的結(jié)果表明使用4個數(shù)據(jù)位時頻偏捕獲準確率已經(jīng)很高，繼續(xù)增加數(shù)據(jù)位捕獲的準確率提高并不明顯，因此選擇4個數(shù)據(jù)位加入前導(dǎo)輔助脈沖是最佳選擇。

需要說明的是，圖6的仿真結(jié)果包含了離散化精度導(dǎo)致的相關(guān)估計誤差因素，在相關(guān)度計算之前需要對接收序列進行歸一化操作，盡可能減小離散化精度造成的估計誤差。

4.2解碼性能仿真實驗

數(shù)據(jù)源的產(chǎn)生與實驗(1)相同，前導(dǎo)輔助脈沖由8 μs前導(dǎo)脈沖和4 μs數(shù)據(jù)脈沖構(gòu)成。隨機產(chǎn)生100 000次信號，首先，進行頻偏捕獲，捕獲誤差分別設(shè)置2.5kHz、5 kHz和10 kHz；然后，進入相位誤差跟蹤環(huán)路并輸出解碼信息，仿真結(jié)果如圖7所示。

圖7　擴容1090ES信號解碼性能

圖7的結(jié)果顯示載頻偏差小于5 kHz范圍的相位信息判決誤碼率接近于理想8PSK調(diào)制下(載頻偏差為0 kHz)的誤碼率，而捕獲后的載頻偏差為10 kHz時解碼性能將下降。因此，算法實際應(yīng)用時將頻偏捕獲精度控制在5 kHz范圍，再使用相位誤差反饋式解碼。

5結(jié)束語

仿真結(jié)果證明了本文提出的同步算法具有頻偏捕獲準確率高、相位解碼性能好、算法復(fù)雜度低和易于工程實現(xiàn)等優(yōu)點，可以用于擴容1090ES信號的同步。本文的算法是在時間定時同步準確的條件下完成，定時誤差是影響1090ES擴容信號解碼性能的重要因素，后期的研究需進一步分析定時同步誤差的影響；另外，1 090 MHz頻點面臨同頻信號A/C和S模式等信號的干擾，在FRUIT和GARBLE干擾環(huán)境下同步算法的性能也有待進一步研究和驗證。

參 考 文 獻

[1]王洪, 孫清清, 李華瓊, 等. 擴充1090ES數(shù)據(jù)鏈容量的相位調(diào)制技術(shù)[J]. 電訊技術(shù), 2015, 55(4): 385-389.

WANG Hong, SUN Qingqing, LI Huaqiong, et al. Techniques to increase 1090ES capacity based on phase modulation[J]. Telecommunication Engineering, 2015, 55(4): 385-389.

[2]王洪. 1 030/1 090 MHz頻譜的共享與干擾問題綜述[J]. 電訊技術(shù)，2013, 53(1): 105-109.

WANG Hong. Overview of sharing and interference at1 030/1 090 MHz[J]. Telecommunication Engineering, 2013, 53(1): 105-109.

[3]顧春平. 空中交通管制監(jiān)視新技術(shù)簡介[J]. 現(xiàn)代雷達, 2010, 32(9): 1-5.

GU Chunping. Introduction of new ATC surveillance techniques[J]. Modern Radar, 2010, 32 (9): 1-5.

[4]PIRACCI E G, GAKATI G, PAGNINI M. ADS-B signals reception: a software defined radio approach[C]//Metrology for Aerospace (MetroAeroSpace). [S.l.]: IEEE Press, 2014: 543-548.

[5]STAYTON G. ASP TSG WP 5-18: ATC overlay data link[R/OL]. [2015-01-08]. http://adsb.tc.faa.gov/ICAO-TSG.

[6]STAYTON G T. Systems and methods for enhanced ATC overlay modulation [P]. United States Patent: US83031105B2, 2011.

[7]PARANZINO J. ASP TSG WP16-22: update on the development of phase modulation of the 1 090 extended squitter[R/OL]. [2014-01-20].http://adsb.tc.faa.gov/ICAO-TSG.

[8]Kay S. A fast and accurate single frequency estimator[J]. IEEE Transactions on Acoustics, Speech and Signal Processing, 1989, 37 (12): 1987-1989.

[9]FITZ M P. Further results in the fast estimation of a single frequency[J]. IEEE Transactions on Communications, 1994, 42(234): 862-864.

[10]LUISE M, REGGIANNINI R. Carrier frequency recovery in all-digital modems for burst-mode transmissions[J]. IEEE Transactions on Communications, 1995, 43(2): 1169-1178.

[11]姜波, 王世練, 溫東. 高速率 8PSK 信號的載波同步及實現(xiàn)[J]. 電路與系統(tǒng)學(xué)報, 2009, 14(4): 21-26.

JIANG Bo, WANG Shilian, WEN Dong. Carrier synchronization for high rate 8PSK signals and implementation[J]. Journal of Circuits and Systems, 2009, 14(4): 21-26.

[12]李晶, 朱江, 張爾揚, 等. 高速8PSK調(diào)制信號的頻率捕獲及跟蹤算法研究[J]. 信號處理, 2005, 12(1): 66-69.

LI Jing, ZHU Jiang, ZHANG Eryang, et al. Algorithm for carrier frequency acquisition and tracking of high-speed 8PSK modulated signals[J]. Signal Processing, 2005, 12(1): 66-69.

[13]王洪，劉昌忠，汪學(xué)剛，等. 強干擾背景下S模式解碼方法[J]. 電子與信息學(xué)報, 2009, 31(12): 2876-2880.

WANG Hong, LIU Changzhong, WANG Xuegang, et al. Mode S decoding methods in heavy jamming environment[J]. Journal of Electronics & Information Technology, 2009, 31(12): 2876-2880.

孫清清男，1990年生，碩士研究生。研究方向為雷達信號處理、高速實時信號處理。

王洪男，1974年生，博士，副教授。研究方向為雷達信號處理、ADS-B、多點定位、數(shù)字接收機和高速實時信號處理。

宋妍女，1991年生，碩士研究生。研究方向為雷達信號處理高速實時信號處理。

姒強男 ，1973年生，博士，副教授。研究方向為雷達信號處理、數(shù)字波形產(chǎn)生和高速實時信號處理等。

Algorithm of Carrier Frequency Synchronization of Phase Modulated 1090ES Signals Based on Preamble Pulses Aided

SUN Qingqing，WANG Hong，SONG Yan，SI Qiang

(School of Electronic Engineering, University of Electronic Science and Technology of China,Chengdu 611731, China)

Abstract：Encoding phase information on 1 090 MHz extended squitter (1090ES) signals is a promising technical method for next generation ADS-B systems，carrier frequency synchronization is the key technique for phase modulated 1090ES signals. An carrier frequency synchronization algorithm is proposed in this paper, based on frequency acquisition and phase error tracking loop, in which binary search tree method based on correlation decision between local signals and input reply signals is used to achieve frequency acquisition for large frequency offset signals and phase error tracking loop is used to achieve small frequency offset estimation and correction. The simulation results show that this algorithm can be applied to the receiver system of phase modulated 1090ES.

Key words：1090ES; phase modulated; carrier frequency synchronization; binary search tree; phase error tracking loop

DOI：10.16592/ j.cnki.1004-7859.2016.05.006

基金項目：國家自然科學(xué)基金重點項目(61139003)

通信作者：孫清清Email：qingqingsunmail@163.com

收稿日期：2015-12-23

修訂日期：2016-02-22

中圖分類號：TN957

文獻標(biāo)志碼：A

文章編號：1004-7859(2016)05-0022-04