郝才勇，劉恒

（國家無線電監(jiān)測中心深圳監(jiān)測站，廣東 深圳 518120）

采用空間稀疏性的單星無源定位方法

郝才勇，劉恒

（國家無線電監(jiān)測中心深圳監(jiān)測站，廣東 深圳 518120）

針對傳統(tǒng)單星定位中，參數(shù)估計和目標(biāo)定位被分為兩個獨立階段造成定位不優(yōu)化的問題，提出了一種采用空間稀疏性的單星無源定位方法。 該方法通過持續(xù)測量目標(biāo)信號的多普勒頻率，利用干擾源位置在地面分布的空間稀疏性，通過構(gòu)建干擾信號的多普勒頻率與地理位置滿足的凸優(yōu)化問題來實現(xiàn)干擾源定位。 該定位方法利用了干擾信號的累積時間和干擾源位置的空間分布信息，在一個階段中完成衛(wèi)星定位。 仿真證實該方法在單星定位中可得到較高的定位精度。

衛(wèi)星干擾源定位；多普勒頻率；空間稀疏性；凸優(yōu)化

1 引言

衛(wèi)星通信經(jīng)常受到地面輻射源的干擾，全球每年大約發(fā)生上百起較為嚴(yán)重的衛(wèi)星干擾事件，而且隨著衛(wèi)星通信業(yè)務(wù)的擴(kuò)展、用戶數(shù)量的增加以及無人值守站的增多，衛(wèi)星 干 擾 問 題 將 越 來 越 嚴(yán) 重［1］。目 前 解 決 衛(wèi) 星 通 信 干 擾 問 題時最普遍使用的技術(shù)是多星定位，如雙星時差／頻差定位、三 星 時 差 定 位 等［2-5］。多 星 定 位 的 優(yōu) 點 是 定 位 精 度 高 、定 位速度快，但要求定位系統(tǒng)必須找到滿足定位條件的鄰星?？捎绵徯潜仨殱M足非常嚴(yán)格的限制條件：軌道位置與受干擾衛(wèi)星（主星 ）接 近 （通 常 C 頻段間 隔 為 10°之 內(nèi) ，Ku 頻 段間隔為 7°之內(nèi)）；信號所在頻段、極化方式與主星一致；地面上行覆蓋范圍與主星一致；在干擾信號對應(yīng)的頻段上沒有在用信號；已知鄰星轉(zhuǎn)發(fā)器的本振頻率等。這導(dǎo)致了在某些干擾情況下很難找到甚至沒有可用鄰星，從而無法完成干擾源定位。

單星定位是解決多星定位系統(tǒng)中鄰星限制條件嚴(yán)格和降低系統(tǒng)復(fù)雜度的一種方法。傳統(tǒng)的單星定位方法主要利用了衛(wèi)星信號在時間狀態(tài)過程中的多普勒頻率變化特性，通過測量不同時刻接收信號的多個多普勒頻率，從而根據(jù)干擾源地理位置與多普勒頻移之間的對應(yīng)關(guān) 系 來 進(jìn) 行 定 位［6，7］。 傳 統(tǒng) 的 衛(wèi) 星 干 擾 源 定 位 方 法 （包 括單星定位和多星定位）在實現(xiàn)過程中一般分為兩個獨立的階段：信號參數(shù)估計階段和目標(biāo)干擾源定位階段。在信號參數(shù)估計階段，首先根據(jù)接收的衛(wèi)星下行信號，測量目標(biāo)信號的多普勒頻率，然后通過計算互相關(guān)函數(shù)來估計信號的到達(dá)頻率差（FDOA），該過程至少需要獲得兩組 FDOA；在目標(biāo)干擾源定位階段，首先利用上一階段獲 取 的 FDOA 和 已 知 的 衛(wèi) 星 星 歷 ，建 立 多 普 勒頻 差 與 干擾源在地面位置對應(yīng)關(guān)系的方程組，然后采用最小二乘算法在地面衛(wèi)星覆蓋區(qū)域中搜索，從而得到干擾源的最可 能 位 置［8，9］。

然而，傳統(tǒng)的定位方法并不優(yōu)化。因為在信號參數(shù)估計階段，假設(shè)前提是目標(biāo)干擾源可能存在于衛(wèi)星信號在地面覆蓋區(qū)域的任意位置，因此每次對信號進(jìn)行參數(shù)測量都是一個獨立過程。在目標(biāo)干擾源定位階段，對目標(biāo)干擾源的定位計算中僅使用了信號參數(shù)估計階段中傳遞來的多普勒頻率參數(shù)，而沒有對定位過程進(jìn)行系統(tǒng)性分析，沒有考慮到干擾源的地理位置在整個定位過程中具有唯一性。事實上，在每次參數(shù)測量和目標(biāo)定位中，目標(biāo)干擾源都位于地面同一位置，而且在多次測量過程中多普勒頻率的變化具有時間連續(xù)性，傳統(tǒng)的定位方法并沒有利用這些先驗信息。

為了優(yōu)化定位過程，提高單星干擾源定位的精度，提出了采用空間稀疏性的單星定位方法。該方法的優(yōu)勢在于，對于地面靜止干擾源，所有的信號測量針對的是相對處在同一位置的干擾源。網(wǎng)格化干擾源所在的地面區(qū)域，不單獨進(jìn)行參數(shù)估計，能夠?qū)⒉煌瑫r刻的信號從發(fā)射到接收作為系統(tǒng)整體來考慮。本文首先介紹了采用空間稀疏性單星無源定位的原理與數(shù)學(xué)模型，然后闡述了定位算法，通過持續(xù)測量目標(biāo)信號來得到多個 FDOA，并采用信號的空間稀疏性來建立信號模型，進(jìn)而構(gòu)造凸優(yōu)化問題以求解干擾源位置，最后通過系統(tǒng)仿真證實該方法在單星定位中可得到較高的定位精度。

2 定位原理與數(shù)學(xué)模型

該定位方法通過建立位于地面靜止位置的干擾輻射源與測量得到的衛(wèi)星信號在多個時間狀態(tài)下的多普勒頻率的關(guān)系，實現(xiàn)衛(wèi)星干擾無源定位。由于衛(wèi)星相對地面的運(yùn)動產(chǎn)生了多普勒頻移，在不同的觀測時刻，衛(wèi)星的運(yùn)動速度不同，從而導(dǎo)致衛(wèi)星接收信號的頻率各不相同，而干擾 源 的 地 理 位 置 影 響 著 多 普 勒 頻 率［10］。將 多 普 勒 頻 率 與 干擾源位置結(jié)合起來就可以得到二者之間滿足的約束條件。定位原理如圖1所示。

圖1 定位原理

對于位于地面的靜止干擾源，將其可能存在的區(qū)域網(wǎng)格化。選取 M個觀測時刻，即衛(wèi)星處于 M 個不同的位置，分別接收干擾源發(fā)射的信號。衛(wèi)星的位置和速度可從星歷中獲取，信號傳播速度已知 （為電磁波在空中的傳播速度），信號接收頻率可以通過地面衛(wèi)星監(jiān)測站測量，需要求解未知干擾源的位置（由于信號傳播過程中衛(wèi)星下行鏈路已知，因此只需考慮信號的上行鏈路傳播）。

該定位方法利用了衛(wèi)星接收信號頻率中的多普勒頻率參數(shù)，為了簡化模型，假設(shè)衛(wèi)星接收的信號僅為載波信號 ， 可 通 過 倍 頻 等 方 法 提 取 調(diào) 制 信 號 中 的 載 波 信 號［11］。 第m個觀測時刻的接收信號可表示為：

其 中 ，am是 衛(wèi) 星 接 收 信 號 的 傳 輸 衰 減 系 數(shù) ，s（t）是 干 擾源 發(fā) 射 的 載 波 信 號 復(fù) 包 絡(luò) ，Wm（t）是 零 均 值 高 斯 白 噪 聲 ，fm是 第 m 個 觀 測 時 刻 的 多 普 勒 頻 率 。信 號 的 多 普 勒 頻 率 fm與信號干擾源的地面位置相關(guān)，可表示為：

在每個觀測時刻，接收機(jī)截取信號并采樣，假設(shè)信號的采樣長度為 N，那么可以將衛(wèi)星接收的信號表示為：

其中，下標(biāo)m表示在第m個觀測時刻對應(yīng)的信號參數(shù) 。rm為 衛(wèi) 星 接 收 的 干 擾 信 號 ，采 樣 個 數(shù) 為 N；s 為 干 擾 源發(fā) 射 的 載 波 信 號 ；diag｛x1，… ，xn｝表 示 ｛x1，… ，xn｝的 主 對 角 矩陣 ；fm為 信 號 的 多 普 勒 頻 率 。

將干擾源在地面可能存在的區(qū)域（通常為衛(wèi)星在該頻段的覆蓋范圍）進(jìn)行網(wǎng)格劃分，為每個網(wǎng)格分配一個布爾值 zx，y，如 果 干 擾 源 位 于 該 網(wǎng) 格 內(nèi) ，設(shè) 置 zx，y為 1；在 其 他 情況 下 設(shè) 置 為 0，完 成 網(wǎng) 格 劃 分 。這 樣 ，zx，y可 以 表 示 為 ：

為了減少計算量，網(wǎng)格劃分過程采取遞進(jìn)細(xì)化的方法，步驟如下：

（1）將該頻段的衛(wèi)星覆蓋范圍劃分為粗粒度大網(wǎng)格，執(zhí)行初步定位，確定干擾源所在的粗略網(wǎng)格；

（2）將該粗略網(wǎng)格進(jìn)一步劃分為較小的子網(wǎng)格，再進(jìn)行定位；

（3）重復(fù)步驟（2），直到定位結(jié)果滿足精度要求。

進(jìn)行網(wǎng)格化處理之后，第 m 個觀測時刻的接收信號可表示為：

其 中 ，Lm，x，y表 示 干 擾 源 位 于 網(wǎng) 格 （x，y）的 多 普 勒 頻 率 參數(shù)。為了方便計算，將二維向量改寫為一維列向量，式（7）可表示為：

下面建立信號滿足的定位約束條件。假設(shè)載波信號 sk為確定性未知信號，可通過最大似然估計得到信號 sk的表達(dá)式［12］。從式（8）得到的信號觀測量為：

則相對于第m個觀測時刻的傳輸信號為：

考慮到所有 M 個觀測時刻（干擾源位于 網(wǎng)格 k），載波信號表示為：

假設(shè)所有的接收信號均來自于位于網(wǎng)格k的同一個干擾源，定義為相對所有M時刻的多普勒頻率算子：

考慮到所有的網(wǎng)格，定義矩陣Ψ包含所有的向量Ψk（k：1，…，K）為 ：

相對于不同觀測時刻的衛(wèi)星，干擾源始終位于地面上同一位置。因此，可以聯(lián)合所有觀測時刻的接收信號來估計出干擾源位置。由式（11）～式（14）可將接收信號表示如下：

其中，r是所有 M 個觀測 時刻接收信 號的向量，z是網(wǎng)格的空間稀疏向量，W是高斯白噪聲。

經(jīng)過這樣的轉(zhuǎn)換，就將衛(wèi)星干擾源定位問題轉(zhuǎn)化為了求解向量中的最稀疏非 0元素的位置。最后，通過建立一個凸優(yōu)化問題來求解：

式 （17）是 一 個 經(jīng) 典 的 凸 優(yōu) 化 問 題［13］，可 使 用 BPIC （basis pursuit with inequality constraint）［14］來 求 解 ，即 完 成 定 位 方 程計算。

3 定位仿真

為了評估提出的單星無源定位方法的可行性，使用蒙特 卡 洛 （Monte Carlo）仿 真 來 計 算 不 同 信 噪 比 下 的 定 位 誤差分布。如果定位誤差在可接受的范圍內(nèi)，說明該定位方法具有可行性。根據(jù)長期以來的衛(wèi)星干擾定位處理經(jīng)驗，衛(wèi) 星 干 擾 上 行 地 球 站 定 位 若 在 150 km 范 圍 內(nèi) ，通 常 可 作為可接受的定位結(jié)果。

假設(shè)干擾源的實際位置坐標(biāo)為 p0，第 n 次實驗計 算 的目 標(biāo) 位 置 為，蒙 特 卡 洛 實 驗 次 數(shù) 為 Nexp，計 算 定 位 點和 目 標(biāo) 位 置 的 均 方 根 誤 差 （root mean square error，RMSE）來 描 述 精 度 性 能［15］，RMSE 定 義 為 ：

定位實驗參數(shù)設(shè)置見表 1。

表1 定位參數(shù)

根據(jù)定位算法，使用表1中設(shè)置的實驗參數(shù)進(jìn)行仿真，計算得到定位結(jié)果及誤差分 布。圖 2和 圖 3表示信噪比 一 定 （分 別 為 7 dB 和 20 dB）時 定 位 仿 真 得 到 的 位 置 與干擾源實際位置距離的分布情況（X 和 Y 分別為地面二維平面的兩個坐標(biāo)軸，定位點分 布 在 XY 平 面 上 ），表 明 在 信噪比低時，定位點的收斂聚合度明顯比信噪比高時低。圖 4表示在不同信噪比下的定位誤差變化趨勢。從仿真結(jié)果可以 看 出 ，當(dāng) 干 擾 信 號 信 噪 比 接 近 6 dB 時 ，定 位 誤 差 迅 速 降低 ；信 噪 比 高 于 11 dB 時 ，定 位 誤 差 降 低 到 150 km 之 內(nèi) ，為可接受的定位結(jié)果。進(jìn)一步提高信噪比，定位精度可穩(wěn)定 在 80 km 左 右 ，且 定 位 點 的 收 斂 度 較 高 ，說 明 該 方 法 能實現(xiàn)精度相對較高的單星無源定位。

另 外 ，衛(wèi) 星 的 攝 動［16］會 造 成 衛(wèi) 星 相 對 速 度 較 小 ，因 而單星定位難以實現(xiàn)較高的定位精度。其定位精度依賴于衛(wèi)星速度、測頻精度、測頻次數(shù)（累積時間）、網(wǎng)格劃分粒度。

圖2 定位點的位置分布（SNR=7 dB）

圖3 定 位 點 的 位 置 分 布 （SNR=20 dB）

圖4 不同信噪比下定位均方根誤差距離

4 結(jié)束語

通過利用衛(wèi)星干擾源位置在地面分布的空間稀疏性對信號定位，能夠?qū)⒉煌瑫r間狀態(tài)下的觀測信號從發(fā)射到接收作為整體來考慮，從而避免了將定位過程分為參數(shù)估計和目標(biāo)定位兩個獨立的階段造成的定位不優(yōu)化問題。本文提出的單星定位方法通過持續(xù)采集衛(wèi)星接收的目標(biāo)信號，建立信號在不同觀測時刻的多普勒頻率和干擾源位置的關(guān)系，然后通過網(wǎng)格化初始區(qū)域形成稀疏向量。將干擾源定位問題轉(zhuǎn)化為求解該向量的最稀疏解的位置，這是一個經(jīng)典的凸優(yōu)化問題，可用很多成熟的算法來解算。由于該方法采用的是單星模式，很難實現(xiàn)高精度衛(wèi)星干擾源定位，但是可用于雙星定位、三星時差定位時無法找到滿足定位條件的鄰星的情況，或?qū)Χㄎ痪纫蟛皇呛芨叩那闆r。

［1］ HAWORTH D，SMITH N，BARDELLI R，et al.Interference localization for EUTELSAT satellites-the first European transmitter location system ［J］.International Journal of Satellite Communications，1997，15（3）：155-183.

［2］ HARWORTH D，SMITH N，BARDELLI R，et al.Interference localization for EUTELSAT satellites-the first European transmitter location system ［J］.International Journal of Satellite Communications，1997，15（3）：155-183.

［3］ 郝才勇，劉 恒，劉宏 立. 基 于 TDOA 的 衛(wèi) 星 干 擾 源定 位 方 法 的 研究［J］.重慶郵電大學(xué)學(xué)報（自然科學(xué)版），2012，24（4）：442-446. HAO C Y，LIU H，LIU H L.Research on satellite interference location based on TDOA ［J］.Journal of Chongqing University of Posts and Telecommunications（Natural Science Edition），2012，24（4）：442-446.

［4］ 李 文 華. 三 星 構(gòu) 型 設(shè) 計 與 時 差 定 位 精 度 研 究 ［J］. 宇 航 學(xué) 報 ，2010，31（3）：701-706. LI W H.Research on configuration of tri-satellites and location precision of TDOA ［J］.Journal of Astronautics，2010，31 （3）：701-706.

［5］ SMITH W W，STEFFES P G.Time delay techniques for satellite interference location system ［J］.IEEE Transactions on Aerospace and Electronic Systems，1989，25（2）：224-230.

［6］ 陸 安 南 ，孔 憲 正. 單 星 測 頻 無 源 定 位 法 ［J］. 通 信 學(xué) 報 ，2004，25（9）：160-168. LU A N，KONG X Z.Passive localization from frequency measurements by single satellite［J］.Journal on Communications，2004，25（9）：160-168.

［7］ 陸 安 南 ，楊 小 牛. 單 星 測 頻 測 相 位 差 無 源 定 位 ［J］. 系 統(tǒng) 工 程與電子技術(shù)，2010，32（2）：244-247. LU A N，YANG X N.Passive location from the combined set of frequency and phase difference measurements by single satellite［J］. Systems Engineering&Electronics，2010，32（2）：244-247.

［8］ BARDELLIR， HAWORTH D， SMITH N.Interference localization forthe EUTELSAT satellite system ［C］／IEEE Globecom’95，November 13-17，1995，Raffles，Singapore.New Jersey：IEEE Press，1995：1641-1651.

［9］ 夏 暢 雄 ，葉 尚 福. 衛(wèi) 星 干 擾 源 定 位 系 統(tǒng) 中 的 融 合 定 位 技 術(shù) ［J］.電波科學(xué)學(xué)報，2008，23（6）：1106-1110. XIA C X，YE S F.Transmitter location system using fusion technology ［J］.Chinese Journal of Radio and Science，2008，23（6）：1106-1110.

［10］瞿 文 中 ，葉 尚 福，孫 正 波. 衛(wèi) 星 干 擾 源 精 確 定 位 的 位 置 校 正算 法［J］. 電 波 科 學(xué) 學(xué) 報，2005，20（3）：342-346. QU W Z，YE S F，SUN Z B.Algorithm of position calibration for satellite interference location ［J］.Chinese Journal of Radio and Science，2005，20（3）：342-346.

［11］CHEN J.Carrier extraction ofQPSK baseband digital demodulation ［J］.Journal of Electronics & Information Technology，1991，13（4）：544-547.

［12］AMAR A，WEISSA J.Localizationofnarrowbandradio emitters based on doppler frequency shifts ［J］. IEEE Transactions on Signal Processing，2008，56（11）：5500-5506.

［13］BOYD S， VANDENBERGHE L.Convexoptimization ［M］. London：Cambridge University Press，2004.

［14］DUARTE M F，ELDAR Y C.Structured compressed sensing：from theory to applications ［J］.IEEE Transactions on Signal Process，2011，59（9）：4053-4085.

［15］YANG K，WANG G，LUO Z Q.Efficient convex relaxation methods for robust target localization by a sensor network using time differences of arrivals ［J］.IEEE Transactions on Signal Processing，2009，57（7）：2775-2784.

［16］肖 峰. 人 造 地 球 衛(wèi) 星 軌 道 攝 動 理 論 ［M］. 長 沙 ：國 防 科 技 大 學(xué)出版社，1997. XIAO F.Artificial satellite orbit perturbation theory［M］.Changsha：National University of Defense Technology Press，1997.

Method of single satellite passive geolocation using spatial sparsity

HAO Caiyong，LIU Heng
Shenzhen Station of State Radio Monitoring Center，Shenzhen 518120，China

Traditional single satellite geolocation is divided into two separate phases：parameter estimation and target positioning，which causes non-optimal localization.To address this issue，a method of single satellite passive geolocation using spatial sparsity was presented.After continuously measuring the Doppler frequency of the target signal and using the sparse distribution of the interference on the ground，a convex problem achieving geolocation through establishing the relation between Doppler frequency and target position of the interfering signal was introduced.The method leveraged cumulative time and spatial distribution of the interference and geolocation was completed in one phase.The simulations confirm that the proposed method can get high accuracy in single satellite geolocation.

satellite interference geolocation，Doppler frequency，spatial sparsity，convex optimization

TN927

：A

10.11959／j.issn.1000-0801.2016124

郝才勇（1985-），男，國家無線電監(jiān)測中心深圳監(jiān)測站工程師，英國倫敦帝國理工學(xué)院訪問學(xué)者，主要研究方向為衛(wèi)星監(jiān)測與無線電定位。

劉恒（1962-），男，國家無線電監(jiān)測中心深圳監(jiān)測站高級工程師、站長，主要研究方向為無線電監(jiān)測與信號分析。

2016-02-16；

2016-04-07