基于反演策略的星基ADS-B信號譯碼方法

馮 濤，梁 俊

(1.空軍工程大學信息與導航學院 西安 710000；2.四川九洲空管科技有限責任公司 四川 綿陽 621000)

廣播式自動相關監視(ADS-B)是航空監視體系中對空管監視雷達的替代和擴展技術[1-3]，是國際民航組織(International Civil Aviation Organization,ICAO)主推的新一代監視技術。星基ADS-B覆蓋范圍遠、無視地形等障礙，可真正實現全球實時連續無縫覆蓋。世界各航空大國均投入大量資源開展相關研究，美國、歐洲等已經成功發射搭載ADSB系統的通信衛星。

星基ADS-B系統接收信號能量弱、信噪比低，目前沒有一種專門適用于星基ADS-B信號的信號檢測與譯碼方法，這是后續解交織[4]、解碼[5-6]等信號處理的基礎。最常用的方法是單脈沖檢測技術[7-8]，通過檢測報文中每個碼脈沖信號上升沿、下降沿、平坦度，再根據一定的門限準則判斷該碼脈沖是否有效。對于星基ADS-B信號，其主要特征是能量弱、信噪比低，適用于陸基ADS-B信號的檢測技術對星基ADS-B信號并不適用。常規單脈沖檢測譯碼技術效率低的原因之一是，該方法只利用了接收端信號而丟失了發射端ADS-B信號中的可用信息，如ADS-B報頭具有穩定的結構，充分利用該信息可以提高報文的檢測概率。反演技術正是基于這種思路的一種檢測技術，在工業CT檢測[9]、醫學成像[10]、地質勘探[11]等領域均有大量成功應用案例。

本文針對低信噪比ADS-B信號譯碼誤碼率大的難題，提出了一種基于反演策略的星基ADSB譯碼方法。該方法假設傳輸一幀ADS-B信號的空間在傳輸信號的時間段內是一個穩定的線性系統，接收端接收信號是發射端信號與該線性系統的卷積和。首先通過滑動窗檢測方法完成低信噪比的ADS-B報頭檢測，利用檢測到的報頭數據與理論報頭數據進行反演計算，估計空間系統參數；然后利用該幀ADS-B接收信號和空間系統參數進行反演計算，得到反演信號；最后將反演信號進行波形校正，實現了ADS-B信號譯碼。通過一組仿真實驗展示了本文譯碼方法的完整操作流程，仿真實驗與常規單脈沖檢測方法對比，結果表明反演譯碼方法能夠實現對ADS-B信號中每個碼脈沖準確定位，譯碼正確率遠高于常規單脈沖檢測方法。

1 ADS-B系統

1.1 ADS-B信號格式

ADS-B信號有3種，分別是1090ES數據鏈、UAT數據鏈和VDL-4數據鏈，其中1090ES數據鏈是全球運輸航空標準，我國運輸航空和通用航空ADS-B系統均采用1090ES數據鏈。

1090ES 模式 ADS-B 信號調制頻率 1 090 MHz，載波頻率允許最大偏差 1 MHz，幀長 120 μs，其中，報頭(前導信號)8 μs，數據域 112 μs。一個完整的1090ES模式ADS-B信號基本結構如圖1所示。報頭信號由4個脈沖組成，脈沖前沿時刻分別為 0、1、3.5 和 4.5 μs，每個脈沖持續時間為 0.5 μs。數據從第8 μs開始，120 μs處結束，采用脈沖位置調制進行信息調制。每一個比特位的碼元有兩個碼片，比特1的碼元由一個前置脈沖和低電平碼片組成，比特0由一個低電平和一個后置脈沖碼片組成，每個碼片持續時間為0.5 μs。數據域中第33～88比特為信息位，包含了ADS-B信息的主要信息，第89～112比特為CRC校驗位。

1.2 星基ADS-B技術挑戰

根據某衛星系統運行軌道高度和星座布局，所需探測距離為3 000 km，信號在傳輸過程中存在路徑損耗，根據式(1)估算最大損耗約為162.7 dB。

當前機載ADS-B信號的發射功率為250 W，該信號經過最大損耗后到達星載ADS-B系統的功率只有?109 dBm，有用信號受到噪聲干擾嚴重。如何在低信噪比下完成對ADS-B信號脈沖碼識別與準確譯碼是星基ADS-B技術發展的重要基礎。

2 反演譯碼技術

2.1 ADS-B信號報頭檢測

每一個1090ES模式ADS-B都有固定的信號報頭格式，如4個相鄰脈沖間隔分別為1、3.5和4.5 μs，則可以判斷該信號為一個有效的ADS-B信號報頭。根據ICAO規定，1090ES模式有效脈沖判決的依據是，脈內平坦度在1～2 dB以內。但是在低信噪比下，信號脈沖受到噪聲干擾嚴重，低電平位置采樣點與高電平采樣點幅值相差不大。當參考電平選取不當時，1～2 dB的平坦度可能無法區分信號中的高低電平，但平坦度過小則會導致檢測容錯力低而丟失真實脈沖，因此，采用單脈沖檢測技術可能無法檢測出完整的報頭。本文在單脈沖檢測的基礎上，利用ADS-B信號報頭結構特征進行綜合匹配判斷，實現低信噪比中報頭檢測。

ADS-B報頭檢測包括單脈沖檢測和報頭結構滑動匹配判斷兩個部分。在單脈沖檢測中，假設ADS-B接收系統的采樣頻率為FS，在接收機收到一幀信號后，選取一個參考脈沖電平HR，以及有效脈沖平坦度范圍 ΔH，從接收信號第1個采樣點開始逐點判斷其電平幅值H的有效性，如果時刻ts處采樣點幅值H＞HR?ΔH，則認為該采樣點幅值有效；若在近鄰NS=FS×0.5μs+ΔN采樣范圍內，有效采樣點個數N的占有率p=p0滿足N/NS≥p，則認為該NS個采樣點組成了一個有效的ADS-B碼片脈沖 P ，脈沖前沿時間為t0=ts。參數HR可根據實際接收信號實時計算，參數 ΔH、ΔN、p可根據實際應用環境人為設置，通常HR由第1個采樣點確定，ΔH≤3dB ，ΔN≤FS×0.1μs ，p≥0.75。

在報頭結構滑動匹配判斷中，首先將有效碼片脈沖 P看作ADS-B報頭第1個脈沖，依次從時刻t0+1μs±ΔT、t0+3.5μs±ΔT和t0+4.5μs±ΔT處作為起始時刻計算時間窗NS內 有效采樣點個數p1、p2和p3，若p1、p2和p3中至少有兩個滿足pi≥p，i=1,2,3，則P1=p1+p2+p3，否則P1=0；再將有效碼片脈沖 P看作ADS-B報頭第2個脈沖，依次從時刻t0?1μs±ΔT、t0+2.5μs±ΔT和t0+3.5μs±ΔT處作為起始時刻計算時間窗NS內 有效采樣點個數p1、p2和p3，若p2和p3中 滿 足pi≥p，i=2,3，則P2=p1+p2+p3，否則P2=0。

在ADS-B信號報頭最終判斷中，若P1=P2=0，則判決脈沖 P 為無效報頭信號；若P1≥P2，且P1＞0，則判決脈沖 P為有效ADS-B報頭信號的第1個脈沖，ADS-B信號報頭起始時刻為t0=ts；若P2≥P1，且P2＞0，則判決脈沖 P為有效ADS-B報頭信號的第2個脈沖，ADS-B信號報頭起始時刻為t0=ts?0.1μs。

2.2 反演譯碼原理

在本文星基ADS-B信號譯碼應用中，假設ADS-B信號傳輸空間為一個線性系統，在一幀ADS-B信號傳輸過程中，該空間系統參數恒定，則從機載ADS-B發射機的無噪信號x(t)、空間系統參數s(t)和星基ADS-B接收系統收到的信號r(t)之間滿足如下關系：

式中，*表示卷積和計算；噪聲ne(t)滿足高斯白噪聲分布。

因此，在接收到信號r(t)后，若能夠估計出s(t)和ne(t)，則使用反卷積算法即可完成原始無噪信號x(t)的估計。為便于計算，將卷積計算公式轉換成矩陣形式：

式中，R、X和n分 別為r(t)、x(t)和ne(t)的離散序列表示；若將s(t)離散化表示為一維序列S，則S為S形成的二維計算矩陣。

根據卷積定義，若x(t)和s(t)的采樣點個數分別為M和L，則r(t)的采樣點個數應該為M+L，并且離散序列滿足：

估計誤差為：

根據最小均方誤差準則，最佳的參數Xopt應該使上述誤差的均方誤差最小，即：

取得最小值或者極小值，式(6)稱為均方誤差性能函數。

2.3 反演譯碼流程

反演譯碼的前提是能夠具備系統參數S的估計值。ADS-B信號報頭不包含任何有用信息，結構穩定，因此，在完成ADS-B信號報頭檢測后，根據反演計算方法[12]可以利用從檢測到的報頭數據和標準的報頭數據得到系統參數S的估計值。

假設檢測出ADS-B信號的報頭部分為RP，標準ADS-B信號報頭數據為XP，將式(6)中的待求參數換成S，則：

令：

直接計算可得：

理論上講，標準報頭數據的幅值不影響最終的參數估計，因為工程應用感興趣的是參數的波形而不是數值。但是根據式(8)進行反演計算時涉及到矩陣求逆或者迭代求逆計算，檢測數據與標準數據之間的幅值相差太大會影響計算結果，因此實際計算中應根據式(1)的損耗計算方法對標準信號進行修正，使其幅度值與實際檢測信號幅值相近。

一幀ADS-B報文在傳輸過程中經歷的信道環境完全一樣，因此，ADS-B信息報頭和數據域兩個部分受到空間系統的影響相同。假設該幀ADSB信號數據部分的一個碼片采樣數據為RD，期望得到的標準ADS-B信號報頭數據碼片采樣數據為XD，則兩者與系統參數S，滿足卷積計算式(3)，其中S由S形成。由于RD和XD長度相同，RD實際選取XD與S卷積和的中間部分，相對應的，矩陣S行也應進行截斷處理。即從原始矩陣S的第L/2行開始，保留包括該行在內的后續M行數據，使其矩陣行數與RD長度相匹配。在反演譯碼計算中，直接將式(9)中的XP替換為S，RP替換為RD，采用矩陣求逆方法或者線性迭代方法。一幀ADS-B信號的數據域包括了112個碼元，反演譯碼分為224個窗口進行。

最后對反演得到的整段ADS-B報文數據進行碼元脈沖波形修正，使各碼元具有標準的方波脈沖特征。修正方法仍采用單脈沖檢測技術，經過波形修正后的波形即為完整的ADS-B報文脈沖序列，計算機讀取該標準方波序列即可實現對該ADS-B報文中有用信息讀取。

3 仿真實驗

3.1 實驗場景設置

為了驗證本文譯碼方法的有效性，將該方法應用到一幀1090ES模式ADS-B信號的前40 μs部分譯碼中，并將該譯碼方法與常規單脈沖方法比較，本文譯碼方法過程包括報頭檢測、空間系統參數估計、數據域譯碼3個步驟，各步驟中的參數如表1所示，信道干擾為加性隨機高斯白噪聲，常規單脈沖檢測方法相關參數與表1中相同。

表1 仿真參數

仿真實驗中模擬的ADS-B信號數據域第1～5位設置為10001，第6～8位為010，第9～32位地址位設置為0。為數據觀測方便，本文仿真實驗只對上述部分的ADS-B報文做譯碼測試，測試結果并不影響整體性能分析。另外，本文實驗中使用的模擬ADS-B報文碼元組成僅為了仿真實驗觀測方便，并不代表真實ADS-B報文結構。

3.2 仿真結果分析

3.2.1 報頭檢測

仿真ADS-B信號按表1中信號生成參數產生，報頭信號幅度值為1，信噪比為?1.58 dB，仿真實驗中數據采樣率為40 MHz，ADS-B信號報頭到達時間為5 μs，分別采用常規單脈沖檢測技術和本文滑動窗匹配檢測技術提取該信號報頭，兩方法中的參考脈沖電平HR為含噪接收信號0 μs時刻開始前20個波峰平均值，兩方法其他參數與表1中報頭檢測部分設置參數一致。

實際報頭檢測結果如圖2所示，由圖可見，常規單脈沖檢測技術將含噪接收信號中疑似碼片的信號段提取出來，并標記為有效碼片脈沖，但是缺乏報頭整體判斷，該方法在第4～12 μs時間段內從該段信號中共檢測出7個碼片脈沖，實際只有3個與理論ADS-B報頭的第2、3和4碼片脈沖位置吻合。滑動窗匹配方法在單脈沖提取結果的基礎上識別出ADS-B的第1個碼片脈沖位置，其他無效位置上提取的疑似碼片脈沖全部被抑制，實測第1個檢測碼脈沖前沿時刻為5 μs，第4個碼脈沖結束時刻為10 μs，滑動窗匹配檢測方法提取出的信號報頭與理論信號報頭碼片脈沖位置完全一致。

3.2.2 空間系統參數估計

按照ADS-B信號報頭檢測結果，選取圖2中4.75～10.25 μs段含噪接收數據為RP，標準ADS-B信號報頭部分為XP，信號傳輸空域的狀態參數為S，根據反演計算式(9)可得持續時間0.5 μs的空間系統狀態參數分布如圖3所示。

空間系統參數并沒有固定的數值結構，只與信號傳輸當時的空間狀態有關，在不同時段和不同空間位置傳輸的ADS-B信號所經歷的空間系統參數不能通用。

3.2.3 數據域反演譯碼

假設同一幀ADS-B信號在空間傳輸過程中空間狀態穩定，該幀廣播信號各段數據經歷的空間系統參數完全相同，構造反演系統參數矩陣為：

在對接收信號數據域反演計算中，一個碼元接收數據RD和待求碼元數據XD長度均為L，則三者的卷積計算可表示為：

則第i個碼元處的反演計算式為：

式中，RDi為第i個碼元接收數據；XDi為第i個碼元反演數據。

從接收到的ADS_B信號第t=t0+8處開始依次截選時間窗為0.5 μs的數據作為數據域一個碼元長度數據RD，窗口內數據采樣點數為L=Fs×Δt=40×0.5=20。按照式(12)通過滑動窗口分步完成對數據域的反演計算。

為了保持反演計算結果數據的完整性，再對報頭部分的數據進行反演計算，即根據接收到的ADS-B信號第0～t=t0+8 μs處數據RP、系統參數矩陣S，按照式(12)反演計算報頭數據XP。

對于一幀完整的ADS-B報文信號，反演計算得到的ADS-B信號格式為對反演得到的ADS-B信號進行波形修正得到最終的報文脈沖序列。

按照本文提供的譯碼操作流程，對該段模擬ADS-B信號進行譯碼，得到的計算結果如圖4所示，圖4a中虛線為接收到的含噪ADS-B信號，實線為該段ADS-B信號的理論波形，即無噪聲ADS-B信號；圖4b中虛線為通過反演計算式得到的ADS-B信號，實線為對反演結果進行波形修正后的ADS-B報文脈沖序列。由圖可見，反演計算結果能將ADS-B報文中低電平處的噪聲很好地壓制，提升報文中高電平處的脈沖特征，通過對反演結果進行簡單的波形修正即可得到與理論ADS-B報文完全一致的脈沖序列。

把常規單脈沖檢測技術應用到模擬ADS-B信號的譯碼中，將檢測結果與本文的反演譯碼結果進行對比。圖5為兩種方法得到的譯碼結果，其中，圖5a為反演譯碼結果及其理論ADS-B報文參照，圖5b為常規單脈沖檢測技術譯碼結果及其理論ADS-B報文參照。由圖可見，反演譯碼結果的每個碼元脈沖位置均與理論ADS-B報文脈沖序列完全吻合，單脈沖檢測技術得到的譯碼脈沖則很少能與理論ADS-B報文脈沖位置重合，并且該方法得到了大量的無效脈沖。單脈沖檢測技術在本文實驗中譯碼效率低下的主要原因是模擬信號的信噪比低，噪聲高能量和高平穩特性掩蓋了原始報文的高低電平特征，限制了單脈沖檢測技術的性能。通過該組對比實驗進一步驗證了本文提出的反演譯碼方法的有效性。

4 結 束 語

低信噪比下的星基ADS-B信號檢測與譯碼是星基ADS-B中的難點。本文提出了一種基于反演策略的譯碼技術，假設傳輸一幀ADS-B信號的空間在傳輸信號的時間段內是一個穩定的線性系統，接收端接收信號是發射端信號與該現行系統的卷積和。利用1090ES模式ADS-B信號報頭特點，本文設計了一種反演譯碼算法流程，并通過仿真試驗進行了驗證，結果表明本文方法能在信噪比為?1.58 dB的情況下實現對ADS-B信號各碼脈沖準確譯碼，正確率遠高于常規單脈沖檢測技術。

雖然將反演技術應用于星基ADS-B信號譯碼是本文首創，也僅是對該方法進行了理論仿真驗證。但反演技術是在多個領域經過實際工程驗證的可靠技術，因此，基于反演技術的ADS-B譯碼方法具有理論可行性。在下一步的研究中，可將該方法應用于實際星基ADS-B信號的譯碼中，研究適用于實際應用環境特征(如動態性、非平穩性、非高斯性等)的反演模型、反演計算方法、算法適用邊界、算法實施流程、不同幀信號空間參數區分等，提高該方法的適用性和穩定性。