0 引言

    廣播式自動相關監視(Automatic Dependent Surveillance-Broadcast，ADS-B)作為國際民航組織(International Civil Aviation  Organization，ICAO)主推的新一代監視技術，結合衛星導航、通信、機載設備以及地面設備等先進技術，能有效地提高運行安全水平，擴充空域容量，提升運行效率以及加強航空公司的運行控制能力^[1-2]。但隨著近年來ADS-B技術的快速發展，監視者與飛機使用ADS-B進行通信時，經常發生信號之間的交織現象，交織后的信號無法被接收機正常譯碼，給通信帶來困難^[3]?；诠_的文獻，現階段國內外的接收機主要通過拋棄交織信號中的其中一條信號來處理交織問題，目前，國內外研究院所正在也在積極展開對交織信號處理的研究工作，因此，根據實際需求，設計能夠進行實時信號解交織處理的接收機是非常有意義的。針對ADS-B信號交織問題，在交織檢測方面，大多是求得一個交織信號奇異值分解后的特征值關于時間的函數，根據曲線來確定合適的閾值并估計交織時刻^[4]；在交織信號分離方面，有ADS-B強FRUIT干擾環境中的解交織方法^[5]、累加分類^[6]和空域濾波^[7]等解交織方法。但上述算法有的要求緩存一整條交織信號、有的涉及高階統計量，計算復雜，不適用于硬件實時系統。

    基于現有公開文獻，發現基于二陣元特征值的交織檢測方法^[7]和投影算法^[8-10]兩種算法的研究都是建立在陣列信號模型的基礎上，效果穩定且優化后能夠降低計算量，適于硬件實現。本文基于現場可編程門陣列(Field Programmable Gate Array，FPGA)器件對上述兩種算法進行優化設計，以降低計算復雜度，減少硬件資源的占用，完成解交織系統的設計實現并通過天線接收ADS-B信號進行處理驗證其性能。

1 信號模型

    采用5個陣元組成的十字陣列天線接收ADS-B信號源，其中陣元天線在時刻t接收到的信號表示為：

2 算法原理與優化

2.1 算法原理

    完成信號解交織系統，需要實現兩個算法，首先在交織檢測方面，利用基于二陣元特征值的交織檢測方法判斷信號是否出現交織以及如果信號發生交織，判斷出交織時刻；然后在交織信號分離方面，利用投影解交織算法對信號進行解交織。

2.1.1 基于二陣元特征值的交織檢測方法

    基于二陣元特征值的交織檢測方法利用兩路陣元信號計算其協方差矩陣并分解得到特征值，設計算協方差矩陣的快拍數為K，則小特征值服從自由度為2K的卡方分布^[11]，其概率密度函數為：

    利用假設檢驗^[12-13]實現對信號的交織檢測，原假設H0：數據段內的信號不多于1個，備擇假設H1：數據段內至少包含兩條信號。以兩路陣元信號協方差矩陣的小特征值作為檢驗統計量^[14]，置信水平為α時的判決域為：

    如果小特征值位于判決區域，則拒絕原假設，認為數據段至少包含兩條信號。

2.1.2 投影算法

    投影算法是基于兩個ADS-B信號在時間上存在部分重疊，如圖1所示。從圖中可以看到在t₁到t₂只存在第一條源信號，t₃到t₄只存在第二條源信號。如果在交織檢測模塊中判斷出信號在t₂到t₃處出現交織，則分別對t₁到t₂和t₃到t₄兩段單獨信號源進行奇異值分解得到其最大奇異值對應的特征向量m₁和m₂，并分別求其特征向量對應的MOORE-PENROSE廣義逆矩陣m₁⁺和m₂⁺，由此得到原信號的估計值為：

2.2 算法優化

2.2.1 基于二陣元特征值的交織檢測方法的優化

    在基于二陣元特征值的交織檢測方法中，由于需要把兩路陣元信號協方差矩陣分解的小特征值作為檢驗統計量^[10]，因此要求對兩路陣元信號做歸一化和標準化處理，也就是說在硬件實現中需要緩存整條交織信號然后根據信號信息對信號完成歸一化和標準化之后才能求解特征值來對比判決域判斷信號是否出現交織。顯然不能滿足實時系統的要求，并且會占用大量資源從而為硬件實現增加難度。因此根據基于二陣元特征值的交織檢測方法設計出一種計算交織檢測動態門限值的方法以適用于硬件實時系統。

    在實際的硬件處理中，采用每隔一定的快拍數就對到來的信號做一次協方差矩陣計算并直接分解其特征值，并從ADS-B信號開始時對每個小特征值進行標準化處理，根據得到的小特征值和卡方分布的原理計算出交織檢測的動態門限值。根據卡方分布的原理：

其中等式右邊th為在卡方分布中置信度為α時的判決門限值；等式左邊為對小特征值標準化的實時處理，eig表示小特征值，m₁表示小特征值的均值，m₂表示小特征值的平方的均值。根據式(10)可計算出動態門限值：

    如果小特征值大于門限值即判斷信號出現交織。

2.2.2 投影算法的優化

    在2.1.2小節中提到的投影算法由于特征向量和廣義逆矩陣在FPGA中需要復雜的計算，將其優化為對協方差矩陣求逆，而且原算法中在估計第一條源信號時，首先需要得到第二條源信號的信息，因而不能滿足硬件實時系統的要求，因此，在本次解交織的硬件實現中通過計算第一條源信號一定快拍數的協方差矩陣m₁，并且計算其協方差矩陣的逆m₁⁺，通過計算估計出源信號：

3 系統實現

    本文采用Verilog語言實現前文描述的解交織系統，使用的軟件為ISE 14.4, FPGA芯片為Xilinx的Virtex5系列。

3.1 系統硬件結構

    解交織系統的硬件實現結構如圖2所示。系統在工作時，射頻信號經過十字陣列天線進入下變頻器，將1 090 MHz 的ADS-B射頻信號轉換為10 MHz的中頻信號，再經過模數轉換模塊形成數字中頻信號進入FPGA1，在FPGA1中，對數字中頻信號數據作預處理和ADS-B信號解交織系統的具體實現，3.2節將具體介紹FPGA1的功能設計。接著，FPGA1將處理后的數據傳輸到FPGA3，FPGA3中有ADS-B接收機功能^[15]，當它在某一時刻檢測到一條ADS-B信號時，將向FPGA1發送一條幀有效信號，FPGA1以此來確定ADS-B信號開始時刻，并從此時開始做ADS-B信號的交織檢測和解交織處理。

3.2 FPGA1功能設計

    本文所實現的ADS-B交織檢測和解交織系統的硬件設計與實現均在FPGA1中完成，本節將介紹FPGA1的功能設計。 

    如圖3所示，FPGA1在工作時，首先對經過模數轉換的數字中頻信號進行信號預處理，包括對數據做希爾伯特變換、去直流分量和下采樣數據復用；之后，預處理后的數據進入協方差矩陣計算模塊得到5×5的協方差矩陣，同時在交織檢測模塊中利用其中2×2的協方差矩陣分解求出雙天線下的特征值，并根據2.2.1節優化的計算方式進一步求出判決門限值，對ADS-B信號是否出現交織進行檢測；如果信號未出現交織，則原信號輸出至FPGA3接收機，一旦發現信號出現交織則協方差矩陣累加模塊根據幀有效信號和信號出現交織的時刻對交織前的五陣元協方差矩陣進行累加、求逆，然后根據2.2.2節中解交織的優化方式進行信號解交織計算，并將解交織后的兩條信號按照先后順序替換在交織時刻后的原信號中輸出。

4 系統測試

    分別利用不同功率的交織信號對實現后的解交織系統進行測試，信號源為ADS-B發射系統通過兩根發射天線發射的交織信號，用上位機軟件控制其發射功率，信號源的產生方式如圖4所示，同時使用5陣元十字陣列天線接收交織信號。接收交織信號并成功解交織后，FPGA3接收機會返回幀有效信號，可以在ChipScope中觀察并作為信號是否解交織成功的依據。

    圖5所示為經過信號預處理后的發射交織時間為60 ?滋s的ADS-B交織信號，隨后進入交織檢測模塊輸出jz_detect判斷交織時刻然后進行信號解交織后的信號如圖6所示。其中jz_detect表示交織檢測標志，當信號出現交織時，標志置1，反之為0；bit_decoding表示從FPGA3接收機中返回的幀有效信號，代表成功解碼ADS-B信號。由圖6可以看到，將圖5中的交織信號分離為兩條ADS-B信號，并解碼成功。

    利用ADS-B發射機發射不同功率的交織信號對解交織系統進行測試，發現影響解交織系統分離信號成功率的因素有交織信號的功率差和交織信號的首條信號功率的大小。如圖7所示，固定交織信號的首條信號功率不變然后改變第二條信號的功率，發現隨著交織的兩信號功率差的增大，解交織成功率逐漸增大至穩定在0.85左右。之后取交織中兩信號功率差在3 dB不變，通過同時增加交織信號兩信號的功率，如圖8所示，發現解交織成功率隨著交織信號的首條信號功率的增大而增大。

5 結論

    針對ADS-B信號的交織問題，對基于特征值的交織檢測方法和解交織方法進行深入分析并根據FPGA實時系統的特點對算法進行優化改進，同時設計一種計算交織檢測門限值的方法，減少了其硬件實現資源消耗，并據此在FPGA中設計實現解交織系統，在實測中性能穩定，證明了本文中設計方案的有效性與可行性，解決了目前國內外市場上接收機無法對ADS-B信號實時解交織的問題，可為未來ADS-B解交織系統在民航領域的實際應用提供參考，具有一定的現實意義。

參考文獻

[1] 高峰，楊秀峰.淺析ADS-B技術原理和常見故障分析[J].科技視界，2018(11)：39-50.

[2] 高永剛.ADS-B監視技術功能探討及特點分析——基于1090ES數據鏈[J].現代商貿工業，2018，39(4)：188-190.

[3] 張青竹，張軍，劉偉，等.民航空管應用ADS-B的關鍵問題分析[J].電子技術應用，2007(9)：72-74.

[4] 張英.基于陣列天線模型的雷達信號分離算法研究[D].成都：電子科技大學，2009.

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

[6] 吳仁彪，吳琛琛，王文益.基于累加分類的ADS-B交織信號處理方法[J].信號處理，2017，33(4)：572-576.

[7] 吳琛琛.ADS-B系統解交織算法研究[D].天津：中國民航大學，2017.

[8] PETROCHILOS N，GALATI G，PIRACCI E.Separation of SSR signals by array processing in multilateration systems[J].IEEE Transactions on Aerospace & Electronic Systems，2009，45(3)：965-982.

[9] PETROCHILOS N，GALATI G，PIRACCI E.Projection techniques for separation of multiple secondary surveillance radar sources in a real environment[C].Fourth IEEE Workshop on Sensor Array and Multichannel Processing，2006.

[10] PETROCHILOS N，GALATI G，MENE L，et al.Separation of multiple secondary surveillance radar sources in a real environment by a novel projection algorithm[C].Proceedings of the Fifth IEEE International Symposium on Signal Processing and Information Technology，2005.

[11] 蔡則元.三大抽樣分布的理解與具體性質[J].數學學習與研究，2019(12)：143.

[12] 王志福，潘旭，金姝，等.假設檢驗的原理及其應用[J].渤海大學學報(自然科學版)，2013，34(2)：101-105.

[13] 馬鳳鳴，王忠禮.假設檢驗方法分析及應用[J].長春大學學報，2012，22(2)：188-192，196.

[14] 王晶，劉彭.參數假設檢驗中統計量的選取問題[J].高師理科學刊，2019，39(4)：55-58.

[15] 劉丹陽.實時ADS-B信號接收機的設計與實現[D].天津：中國民航大學，2017.

作者信息:

胡鐵喬，韓  斌

(中國民航大學 天津市智能信號與圖像處理重點實驗室，天津300300)