鄭高哲,劉偉駿,江建興
(國家無線電監測中心福建監測站,福建 廈門 361004)
傳統的無線電監測定位都是基于無線電測向的交會定位,無線電測向依賴于多信道測向接收機和各種測向體制、不同孔徑的測向天線陣,受其技術特性、設備復雜度的影響,在監測超視距信號尤其是沿海岸線布局時,會產生無法形成交會定位結果等情況。在當前我國各地市級別的無線電管理機構中,針對超視距信號也無法具備足夠數量的測向站。
相比之下,基于成本較低的分布式射頻傳感器監測網絡采用時差定位方法,在同樣場景尺度下更具可行性。相對于固定監測站、移動監測車等傳統的監測測向設備,其具有設備簡單、成本低廉、攜帶架設方便等優點。近年來,射頻傳感器應用越來越廣泛,包括城市內的網格化監測、沿海的水上業務監測、重大活動的無線電安全保障等。
到達時間差(Time Difference of Arrival, TDOA)無需依靠定向天線或測向天線陣,只要有至少3個監測點能夠接收到目標信號,即可通過3組時間差轉換為3組雙曲線的交點確定為信號源的位置。TDOA定位在算法上主要有兩個環節,即雙曲線方程解算和時差估計,而其性能受監測點布局形式的影響。
雙曲線方程解算主要有兩種經典算法:Taylor算法和Chan算法。
Taylor算法是利用泰勒級數展開的基于加權最小二乘法的迭代算法。該算法在給出定位估計初始值的情況下,通過循環迭代以逐漸收斂至某設定閾值范圍內,得到測量誤差的局部最小二乘解。該算法運算量大,定位精度取決于設定閾值,若定位估計初始值有缺失將無法保障該算法的收斂性。Chan算法是基于兩步加權最小二乘法的具有解析解的算法。該算法能夠利用多余的監測點提升定位精度,計算復雜度較低,不需要定位初始估計值,且能夠給出一個閉式解析解,但容易出現模糊與無解的情況。

圖1 協同算法流程圖
在工程實踐中通常使用二者的混合及協同算法,利用Chan算法給出目標初始位置,再利用Taylor算法進行迭代可減少迭代次數,快速給出高精度的解算結果[1]。
時差估計最為經典的即互相關算法,求兩個接收信號的互相關函數,峰值即為信號傳播的時間差。標準互相關是直接將兩路接收信號進行互相關處理,對其互相關函數進行峰值檢測。實際環境中由于信號特性不理想造成無法達到理論條件的峰值。廣義互相關在其基礎上,對兩路接收信號先進行預濾波,在頻域上作加權處理,以抑制噪聲功率使相關峰更加突出。

圖2 廣義互相關原理
1.3.1 幾何布局形式GDOP
幾 何 精 度 因 子(Geometric Dilution of Precision, GDOP)用于描述定位誤差與幾何的關系。定位目標的位置應盡可能靠近監測點布局形狀的幾何中心,這樣布局對定位精度的影響最小,從而監測點的布局應盡量在以信號源區域為中心的周邊形成較為規則的閉合區域。以及應盡量拉長監測點之間的基線長度,定位精度隨基線長度的增加而增大。

圖3 幾何精度因子原理[2]
1.3.2 圓概率誤差CEP
圓概率誤差(Circular Error Probable, CEP)的定義是以估計均值為圓心,包含50%定位估計結果的圓的半徑。CEP用于二維定位誤差分析,可簡單準確地表征定位性能[3]。
CEP是在對實測信號、環境 和系統進行簡化假設下計算出的最佳情況概率,它不包括環境、天線或干擾信號的影響。CEP半徑越小,則表明包含50%定位數據的范圍越小,定位結果越集中;反之,則定位結果代表的定位數據越分散。
本次超短波遠程接收實驗的目標信號是海洋對岸的調頻廣播信號,懷疑是由海洋蒸發波導類型的大氣波導發生超視距傳播。如圖4所示,地圖左側陸地距離海洋對岸最近處約330 km,最遠處約450 km。

圖4 測試場景地圖和GDOP分析
為了獲得更好的定位精度,部署5個監測點盡可能在南北方向上拉長監測點之間的基線長度,各監測點天線架設高度距離地面約10~30 m,具有可靠的供電條件和4G網絡覆蓋。GDOP仿真結果表明:接近直線排列情況下已對其東部海域形成最大監測覆蓋,這已經是測試環境下能夠利用海岸線擴大基線長度且盡可能保障共同收到大氣波導傳播信號的最佳布局形式。

圖5 單個監測點系統框圖及典型架設
接收設備采用Keysight公司的射頻傳感器N6841A,具有2個天線接口,配備1副調頻廣播頻段(87-108 MHz)專用的定向八木天線和1副全向天線(20 MHz-6 GHz),1副GPS天線,連接1套4G通信模塊。
在連續多天的測試中,僅有3天較為明顯地接收到了相關語種的調頻廣播信號,其余幾天均未收到,對照“赫本對流層指數”(Hepburn Tropo Index)網站[4]每天發布的分時大氣波導預測圖與測試接收情況完全相符。同時,判斷廣播的語種信息也驗證了超視距傳播。

圖6 大氣波導覆蓋預測圖
選取接收到信號的3天內出現時間最長、定位次數最多的某廣播信號進行定位分析,經業余無線電網站調研該頻率共有4個廣播臺址見圖7。功率信息見表1。

表1 93.9MHz電臺方位及功率

圖7 目標方向的定位擬合
直觀的統計分析可看出:
(1)定位結果中大部分位于靠近沿海5個監測點的一定距離范圍內。約73.5%的定位數據位于距海岸線100 km范圍內,約62.8%的定位數據位于距海岸線50 km范圍內。
(2)TDOA定位可對超視距信號產生示向線,并分辨不同方向的同頻信號。兩條匯集定位結果的紅色延長線分別朝向陸地的西北部、西南部電臺方向,且概率橢圓長軸方向逐漸偏向真實電臺的位置。
(3)兩條延長線的末端在真實的4個電臺位置附近均有一定數量的定位結果。
做進一步數據處理篩選50%較優的CEP半徑 和互相關數據的定位結果,見圖8。

圖8 目標方向50%較優數據的定位擬合
可以看出,兩條定位結果匯集形成的示向線依舊清晰,且仍有定位結果位于真實電臺位置附近。若TDOA單次定位無法準確定位信號發射源,可通過長時間持續定位的統計分布對信號傳播方向進行描跡。模糊預測當大氣波導陷獲層覆蓋至電臺的發射天線和監測點的接收天線高度時,TDOA定位方法可定位信號發射源的真實位置。
本文一方面驗證了接收經由大氣波導傳播的超視距信號,另一方面也驗證了沿海岸線部署射頻傳感器對超視距信號采用時差定位方法的可行性。后續有待繼續挖掘造成超視距傳播的大氣波導傳播特性,探求時差定位與電波傳播模型的結合點,以及研究將時差定位與傳統測向定位相結合的協同定位,為超視距信號的無線電監測定位提供參考。