水下多目標跟蹤技術現狀與展望

張思宇, 何心怡, 張 馳, 祝 琳, 陳 雙

?

水下多目標跟蹤技術現狀與展望

張思宇, 何心怡, 張 馳, 祝 琳, 陳 雙

(海軍研究院, 北京, 100161)

水下多目標跟蹤技術在軍事和民用方面均有重要作用, 是軍民融合發展中的重點技術。針對水下水文條件復雜、作用距離相對較小等不利因素, 文中根據跟蹤濾波算法原理的不同, 分別論述了基于數據關聯的水下多目標跟蹤技術和基于隨機有限集的水下多目標跟蹤技術, 詳細闡述了其目標運動模型、跟蹤濾波方法和應用現狀, 梳理了在水下實現多目標跟蹤關鍵的數據關聯技術和隨機有限集技術的性能, 分析了由于漏報和虛警導致觀測信息的不確定、跟蹤過程中目標數量不確定和運動狀態以及跟蹤算法實時性差等3類技術瓶頸, 突顯了建立統一的隨機有限集框架描述跟蹤問題解決該類瓶頸的優勢。在此基礎上, 根據作戰使用和海洋開發的要求, 展望了水下多目標跟蹤技術發展方向, 供相關研究人員參考。

水下多目標跟蹤; 數據關聯; 隨機有限集

0 引言

水下目標跟蹤技術在海洋資源勘探、水下工程作業、海戰場監視及水中兵器精確自導等方面用途廣泛, 有著廣闊的軍民融合發展前景。該技術通過對量測的目標數量和狀態信息實時處理, 維持對目標數量和運動狀態準確估計的過程。根據探測器的數量和被跟蹤目標數量的不同可分單個探測器跟蹤單個目標(one-to-one, OTO)、單個探測器跟蹤多個目標(one-to-multiple, OTM)、多個探測器跟蹤單個目標(multiple-to-one, MTO)及多個探測器跟蹤多個目標(multiple-to-multiple, MTM)等4類情形, 其中OTM、MTO、MTM等3類情形要通過多個探測器之間的數據傳輸、數據融合和對多個目標的量測信息進行數據關聯, 技術實現較為復雜。與雷達系統目標跟蹤應用場景相比, 聲吶系統由于應用場景的水文條件復雜、作用距離相對較小、量測信息誤差大, 同時又存在和電磁波相比聲速很慢、系統數據刷新率低、能觀性差以及濾波收斂慢等問題, 使得跟蹤實現難度更大。

經過50多年的研究發展, 多目標跟蹤理論在實際工程應用中取得重大突破, 部分最新服役的魚雷等精確制導水中兵器具有多目標跟蹤能力[1]。如美國的MK54魚雷、法意聯合研制的MU90魚雷及意大利研制的“黑鯊”魚雷等均具備定位且同時對多個目標進行跟蹤的能力, 可自動判斷最具威脅的目標并辨識出真假目標。

多目標跟蹤問題最早由Wax[1]于1955年提出, 到上世紀70年代, 隨著計算機技術的發展, 逐步受到學者的廣泛重視和研究。Singer等[2]提出了最近鄰算法, 以目標統計距離最小值作為目標跟蹤門限, 實現多目標跟蹤, 但僅適用于高信噪比的跟蹤背景。Fortmann等[3]于1983年提出了聯合概率數據關聯算法, 通過計算每個量測信息和所有可能目標的相關概率, 篩選出最佳關聯概率, 進行多目標跟蹤, 但存在計算量巨大的缺陷。Reid[4]基于Singer和Fortmann等的貢獻提出了多假設跟蹤算法, 通過時延邏輯, 實現有效的多目標航跡保持, 該算法計算量小, 但存在對多目標的先驗信息依賴程度高的缺陷。2003年, Mahler[5]首次提出了基于隨機有限集的多目標跟蹤算法, 受到學界的廣泛關注。Vo等[6]對相關算法進行改進, 使基于隨機有限集的多目標跟蹤算法取得了長足的進步。

文中系統梳理了水下多目標跟蹤技術相關文獻, 總結評述了傳統基于數據關聯的多目標跟蹤技術和新興的基于隨機有限集的多目標跟蹤技術的目標運動模型、跟蹤算法、應用現狀及當前的技術瓶頸, 并展望了技術發展方向, 以期為相關研究提供參考。

1 基于數據關聯的多目標跟蹤

該跟蹤方法除了涵蓋單目標跟蹤理論的建立目標運動模型、選取合適的跟蹤坐標系、跟蹤濾波與狀態預測等關鍵環節外, 還產生了跟蹤門、數據關聯與維持、跟蹤的起始與終結、虛警與漏報等主要處理步驟, 其核心是確定探測器接收到量測信息和目標源對應關系的數據關聯[7]。多目標跟蹤利用貝葉斯方法將多目標跟蹤問題轉化為單目標跟蹤問題, 進而遞推處理, 其基本原理如圖1所示。

首先, 利用目標的數量和狀態的量測信息更新已有軌跡, 通過設定的跟蹤門判定量測信息和已有軌跡的配對情況, 進而用數據關聯確定量測數據與已有軌跡的最佳配對, 通過機動辨識和跟蹤濾波保證整個跟蹤過程的維持。來自新目標或雜波的量測信息和已有軌跡是互不相關的, 通過跟蹤起始的方法進行判別新目標的真假, 對真目標建立新的軌跡檔案。當目標離開跟蹤區域, 通常及時結束跟蹤, 并刪除目標軌跡檔案, 避免無效計算。在整個過程中, 跟蹤門的大小和中心位置要根據目標狀態估計實時調整, 確保跟蹤精度。以下對多目標跟蹤的關鍵環節進行逐一概述。

1.1 目標運動模型的建立

建立合適的目標運動模型要從兩方面考慮, 一是模型要準確地描述目標的運動狀態, 即盡可能涵蓋運動的全面特性, 如目標位置、速度、加速度和描述外界引起的不可預知現象的狀態噪聲等; 二是建立的模型數學處理簡單, 要求狀態向量的維數盡可能少, 對一些變量進行近似和簡化, 保證計算效率和跟蹤的實時性。因此, 如何建立合理的模型經歷了由簡單模糊到復雜精準的過程, 主要有勻速運動(constant velocity, CV)和勻加速運動(constant aceleration, CA)模型[8]、辛格(singer)模型[9]、半馬爾科夫過程模型[10-11]、目標運動“當前”統計模型和交互式多模型(interacting multiple model, IMM)[12]等描述目標運動的模型。

1) CV和CA模型

2) Singer模型

Singer模型采取有色的噪聲分布描述機動加速度, 更貼近目標的實際運動狀態, 故跟蹤效果較好, 應用廣泛; 但是潛艇等水下目標發現受到攻擊, 一般會采取大機動規避動作。此時, Singer模型與實際偏差較大, 很難準確跟蹤目標, 甚至會失跟。

3) 半馬爾科夫模型

為了克服Singer模型在做目標大機動規避時, 估計誤差較大的局限, Moose[11]提出了應用半馬爾可夫過程描述目標機動過程。

4) 目標運動“當前”統計模型

由于該模型的“當前”加速度是基于機動加速度統計量均值的預測值, 且服從修正瑞利分布的時間相關模型, 故對目標運動狀態描述更符合實際。

5) IMM模型

由于機動目標的運動狀態的不確定性, 單一目標運動模型難以準確描述目標實際運動, 故Blom等[12]提出用不同的機動目標運動模型來描述不同的運動狀態, 以馬爾可夫過程來描述不同機動模型之間的轉移概率, 用卡爾曼濾波實現機動模型和目標運動狀態估計。從模型理論上講, 如果構建所有運動模型則可實現對機動目標狀態的最佳估計。但是, 目標機動存在不確定性, 工程實踐中通常運用以上幾個典型的機動目標模型近似描述機動目標實際運動狀態。該模型不需要對機動目標進行檢測, 可以實現全面自適應選擇, 是一種能夠較好描述大機動目標的機動目標模型, 故可被廣泛應用[14]。

1.2 跟蹤濾波與狀態預測的確定

1.2.1 線性濾波

機動目標的運動過程一般可以用穩態過程和暫態過程來描述, 而上述濾波算法都有其對應的合適的運動過程, 需要依據實際情況選擇。另外, 根據濾波算法的實現精度和目標軌跡的實時性來鑒定系統整體性能, 從跟蹤精度來評價, 精度最高的是卡爾曼濾波, 但是計算時間最長; 從跟蹤實時性來評價, 兩點外推算法計算時間最短, 實時性強, 但是精度較低。隨著現代微型計算機性能大幅提升, 系統能夠滿足卡爾曼濾波的運算要求, 在線性高斯白噪聲模型中, 卡爾曼濾波是最佳選擇。

1.2.2 非線性濾波

非線性運動模型更符合目標的實際運動, 故非線性濾波方法有較高的精度。根據估計方法不同, 可將非線性濾波分為點估計和概率密度估計兩大類。其中, 點估計根據逼近方法不同可分為: 函數逼近類, 如擴展卡爾曼濾波; 統計逼近類, 如無跡濾波; 隨機模型逼近類, 如中心差分濾波。而概率密度估計應用最多的是粒子濾波。文中主要評述以下3種最常見的非線性濾波方法。

1) 擴展卡爾曼濾波(extended Kalman filter, EKF)算法是卡爾曼濾波在處理非線性跟蹤濾波問題的擴展, 一般用Taylor級數展開法對非線性的目標運動狀態方程或量測狀態方程線性近似[14]。該算法的優點在于Taylor級數展開將非線性問題線性近似, 算法較簡單, 能以較小的計算量取得較好的跟蹤濾波效果, 工程實現性好, 應用廣泛。其局限是由于該算法只是一種有限的線性化近似, 處理非線性度較高的系統時, 其模型的線性化誤差較大, 濾波精度低, 甚至可能出現濾波發散。

2) 無跡卡爾曼濾波(unscented Kalman filter, UKF)算法是利用機動目標模型初始分布的均值和協方差確定一系列sigma采樣點, 將采樣點通過非線性轉移, 進而估計出目標的均值和協方差[18]。該算法直接求出目標的均值和協方差, 不會產生非線性系統線性化的近似誤差, 故跟蹤精度要高于EKF。但是, 處理4維以上系統會出現非正定協方差, 導致濾波信息輸出不穩定, 甚至發生濾波發散, 所以難以應用到高維系統處理。

3) 粒子濾波(particle filter, PF)算法是利用序貫蒙特卡洛(sequential Monte Carlo, SMC)方法直接逼近非線性系統概率密度函數, 再加權求和估計出目標運動狀態的一種濾波算法[22]。該算法的優點是不受機動目標模型限制, 能夠處理在非線性非高斯動態過程中的機動目標跟蹤問題, 可減少計算量, 在未知目標數量的條件下, 實現有效估計。但該算法的粒子退化問題是序貫重要性采樣難以消除的, 在多次迭代計算后, 除少數粒子權值有效外, 絕大多數粒子權值趨零, 造成大量無效加權求和。且跟蹤精度和濾波收斂性只能依靠大量粒子, 計算量大, 難以保證跟蹤實時性。

綜合以上各種濾波算法的實質就是求取最大后驗概率的最優估計問題?？柭鼮V波是線性高斯白噪聲條件下的最優濾波器, EKF是非線性非高斯有色噪聲條件下的次優濾波器, 以上在未知目標數量的條件下, 濾波性能較差?；赟MC方法的PF能有效破解上述難點, 但目前計算量大和粒子衰退問題是需攻克的難點。從工程實際出發, 采用基于點估計的濾波算法, 可以以較高的精度估計目標狀態, 計算量小, 跟蹤實時性好, 系統成本較低。

1.3 多目標跟蹤的基本要素

1.3.1 跟蹤門的形成

跟蹤門作為跟蹤空間的一個子空間, 其大小取決于接收正確量測信息的概率大小, 以目標估計的運動狀態作為中心。跟蹤門是將觀測信息與已有目標軌跡或新建目標軌跡配對的一種預先粗檢手段, 如果探測器接收到的量測信息符合已設定的跟蹤門, 則將該量測信息作為已有目標的回波, 估計并更新目標狀態; 如果量測信息沒有落在跟蹤門內, 則將該信息視為虛警或者有新目標進入, 再進一步判斷采取摒棄該信息還是建立新的候選目標軌跡。

1.3.2 跟蹤起始與跟蹤終結

跟蹤起始是通過形成假設目標軌跡、初始化目標軌跡和確定目標軌跡等3個步驟, 對新進入監視區域的目標, 建立軌跡檔案的決策方法。跟蹤起始在量測信息與已有目標軌跡不能配對時發起, 通常采用次掃描中出現次相關的方法降低虛警概率, 保證跟蹤起始的有效性。

跟蹤終結是發現目標離開監視區域后, 刪除其軌跡檔案的決策方法, 可增加有效存儲, 提高計算速度, 增強跟蹤的實時性。

1.4 數據關聯算法

數據關聯是處理目標不確定性和量測不確定性的關鍵方法, 也是該跟蹤方法的最核心環節[23]。一般從面向目標、面向量測及面向航跡這3個方向設計數據關聯算法, 分別是基于觀測值是來自已知目標還是雜波、量測信息是來自已知目標還是新目標或雜波和目標軌跡是跟蹤起始、跟蹤維持還是跟蹤終結等3類情況考慮的。

1)“最近鄰”數據關聯算法

該算法是根據目標狀態估計值設計跟蹤門初步篩選所接收的量測信息作為初步候選回波, 進一步篩選出和目標狀態估計值的統計距離最小的量測信息作為候選回波, 然后將候選回波應用上述的單目標跟蹤濾波算法估計出目標運動狀態。該算法優點是能在高信噪比和目標稀疏的條件下, 以較小的運算量實現較高精度的跟蹤, 工程實現較簡單; 但是由于艦艇編隊密集或水文環境不佳時, 密集的量測信息回落于同一跟蹤門內, 則統計距離最小的候選回波不一定是目標回波, 容易出現虛警, 易產生誤跟或失跟的現象。

2) 聯合概率數據關聯(joint probability data association, JPDA)算法

由于利用概率進行關聯的算法只能用于多目標跟蹤門不重疊或者沒有多個目標回波落在同一跟蹤門的條件下, 通過對所有可能的量測信息加權平均估計出目標運動狀態的局限, 推廣出了JPDA算法。該算法通過計算每個量測信息和所有可能目標的相關概率, 篩選出最佳關聯概率, 進一步估計出目標運動狀態。該算法歷遍所有量測信息和所有可能目標, 在密集多目標和多雜波的情況下, 有較好的跟蹤效果; 但也存在目標數量較多時, 計算量太大導致組合爆炸的問題。

3) 軌跡分裂數據關聯算法

軌跡分裂數據關聯算法的基本原理是在接收到多余目標回波時, 按照現有新的目標軌跡數分裂原有軌跡與之對應, 進而求出每一條軌跡的似然函數, 僅保留高于門限的回波, 舍棄低于門限的回波。其優點是可以有效發現并舍棄虛警軌跡; 但是在密集多目標和多雜波情況下, 系統計算和存儲壓力很大, 有組合爆炸的可能。

在現有計算能力下, 由于目標不確定性和量測不確定性的限制, 數據關聯算法存在計算量小則跟蹤精度低, 跟蹤精度高則計算量大的矛盾。故可將模糊推理、模糊聚類、模糊數學結合神經網絡和目標屬性與狀態識別等新技術應用于數據關聯算法, 以期解決復雜水下環境背景的密集多目標跟蹤難題。

2 基于隨機有限集的多目標跟蹤

傳統基于數據關聯的多目標跟蹤算法, 實質仍是通過數據關聯將多目標跟蹤轉換化為單目標跟蹤, 難以實現多源多目標數據融合和突破計算量大的瓶頸。Mahler[24]于2007年首次提出了基于隨機有限集(random finite set, RFS)的多目標跟蹤算法, 將每一個目標狀態和量測信息視為隨機有限集結, 進而通過貝葉斯濾波算法直接估計出目標數量和運動狀態, 給出了相關近似算法, 有力地推動了該技術的發展。該算法能夠將5級的信息融合均統一于RFS的框架下, 較好實現數據關聯困難的復雜環境下密集多目標跟蹤問題, 受到業界的廣泛關注[15-28]。目前, 美國在該方法的研究取得了很大的理論進展和應用成果, 如洛馬公司的反潛多目標跟蹤識別系統[5, 24]、海軍研究辦公室的被動聲吶反潛信息融合和目標識別系統[6, 26]。國內也有許多研究機構開展了大量的探索工作, 取得了一些初步的研究成果, 但離實際應用仍有差距[29-36]。

該算法也是利用貝葉斯濾波遞推計算, 區別在于其是對隨機有限集進行積分, 如式(9)和式(10)所示[17]。

由上式可見, 對隨機有限集積分很難求解, 故Mahler等提出了許多近似的多目標濾波算法, 如: 概率假設密度(probability hypothesis density, PHD)濾波算法[26]、勢概率假設密度(cardinalized probability hypothesis density, CPHD)濾波算法[30]和多伯努利(multi-Bernoulli, MB)濾波算法[27]等。

2.1 基于RFS的跟蹤系統模型

一般情況下, 各集合之間相互獨立, 其模型如圖2所示。

2.2 基于RFS的跟蹤濾波算法

Mahler[24, 26-27, 30]首先給出了單傳感器探測多目標背景下的最優貝葉斯濾波器-多目標聯合跟蹤濾波器(joint target-detection and tracking, JoTT), 并針對目標數未知、最優貝葉斯濾波器無法實現的問題, 進一步推導出概率假設密度濾波器、勢概率假設密度濾波及多目標多伯努利濾波器。Vo[31]針對多目標多伯努利濾波器的目標數有偏估計問題, 推導了勢平衡多目標多伯努利濾波器, 提高了估計目標數的能力。

1) PHD濾波算法。該算法是將多目標狀態空間利用假設的后驗概率密度轉換到單目標狀態空間, 避免了計算復雜的數據關聯過程[15]。其核心是通過多目標狀態概率假設密度函數的1階矩求得后驗概率密度, 進而估計多目標狀態。對于任一監視區域, 其集值積分為

2) CPHD濾波算法。鑒于PHD算法的不足, Malher[39]提出了以遞推的形式同時傳播目標狀態概率假設密度函數和目標數的后驗概率分布, 且只需利用目標數分布函數即可估計出目標個數的算法, 即CPHD跟蹤濾波算法。該算法也是通過GM和SMC[41-42]方法來實現線性和非線性條件下多目標跟蹤, 并在具體問題有對應的改進算法, 如可通過EKF的思想, 將非線性函數近似線性化, 利用GM-CPHD濾波器在線性系統下處理; 也可運用粒子濾波的思想, 對SMC-CPHD濾波器進行優化, 提高濾波性能?？偟膩碚f, 運用CPHD濾波算法可以提高濾波和跟蹤的精度, 但隨著計算復雜程度和計算量的增加, 增大了工程的實現難度。

3) MB濾波算法[43-46]是通過傳播一個多伯努利隨機有限集參數來近似多目標后驗概率密度, 進而估計目標個數和狀態信息。但是Beard[43]和Schuhmacher[45]指出, MB濾波算法是一種有偏估計, 因而改進為勢平衡多伯努利濾波(cardinality balanced multi-target multi-Bernoulli, CBMeMBer)算法, 并提出了線性條件下GM實現和非線性非高斯條件下SMC實現方法。由于CBMeMBer算法提出時間不長, 仍有許多問題有待解決, 是當前的研究熱點之一。

綜上, 基于隨機有限集的多目標跟蹤濾波器, 根據目前公開發表的文獻來看, GM-CPHD在線性高斯條件下性能最好, 是最優貝葉斯濾波器; 在非線性非高斯程度很高的環境中, CBMeMBer濾波器能夠改進目標數的有偏估計, 且計算量較小, 跟蹤精度高和并可直接獲取目標狀態。故建議根據具體問題, 如在線性高斯條件下, 選擇基于GM-CPHD改進的濾波器; 在非線性非高斯條件下, 選擇基于SMC-CBMeMBer的次優貝葉斯濾波器。在水下多目標跟蹤實現上, Peter團隊[40]已有初步應用, 采用GM-PHD濾波器取得了較好的多傳感器的多目標跟蹤效果; 劉姝琴等[14]初步探索了如何利用水下無線傳感器網絡(underwater wireless sensor network, UWSN)實現多目標跟蹤。

2.3 基于RFS的跟蹤濾波性能評價指標

該跟蹤濾波方法沒有采用數據關聯過程, 故傳統將多目標跟蹤轉化為多個單目標跟蹤的評價方法不再適用?；诖? 有學者提出以合理量測2個隨機有限集之間的距離作為新的評價指標。主要有Hausdorff距離[45]、Wasserstein距離[46]和OPSA(optimal subpattern assignment)距離[47]等3種評價指標。

Schuhmacher等[45]提出了Hausdorff距離是度量RFS的性能的最初標準, 能夠表現出目標狀態和數目估計的局部性能。但是, 該度量方法沒有定義出現空集時Hausdorff距離的意義, 且對RFS之間的勢差異不敏感, 無法給出跟蹤性能的準確評價。Nagappa等[46]提出了Wasserstein距離是統計學中度量非空集合之間距離的方法, 其比Hausdorff距離對RFS的勢差更敏感, 當2個RFS的勢大小相同時, 是最優關聯下的距離; 但是也存在出現空集時, 距離無意義, 且對RFS的勢差敏感, 目標數量估計可能出現錯誤懲罰。Ristic等[47]提出了OPSA距離是對Wasserstein距離的一種改進度量方法, 其優點是在出現空集情況, OPSA距離都具有一定物理意義。Ristic等[47]和He等[48]在RFS的勢誤差和狀態估計誤差的基礎上, 并入估計質量信息, 提出了更高估計精度的Q-OPSA距離, 提高了度量多目標跟蹤精度。

綜述以上的度量方法, 由于目標估計質量信息一般很難保證, 故目前OPSA距離是認可度較高的基于隨機有限集多目標濾波評價指標。

3 水下目標跟蹤技術展望

攻防裝備發展是一個此消彼長、螺旋上升的發展過程, 隨著裝備的隱蔽性和機動性能的提高, 海區監視、魚雷打擊目標和反魚雷作戰的難度也不斷增大, 從當前和未來一個時期作戰使用需求來看, 水下跟蹤技術應重點朝以下幾個方向發展。

1) 隨著敵方艦艇戰場感知能力增強和海戰場大范圍監控的任務需要, 對隱蔽、大范圍跟蹤監視系統的需求也日益加強。因此, 應大力發展多基地聲吶監視系統, 發揮其跟蹤監視范圍大、可靠性高和己方艦艇隱蔽性強的優勢, 提高我國重點海區的水下跟蹤監視能力, 保證國家海洋安全。

2) 鑒于魚雷攻擊距離的增大, 對自身隱蔽性要求的提高, 以及打擊艦艇編隊和具備水聲對抗能力的單艇目標, 需要具備真假目標識別和多目標跟蹤的能力。因此, 因重點提高魚雷自導系統的測量精度, 克服水下可觀測性差, 跟蹤收斂慢的困難。同時, 應在自導系統引入目標識別功能, 有效提高魚雷發現并攻擊真目標和高價值目標的能力。

3) 隨著海洋開發的日益深入和海軍走向遠海大洋, 存在著許多不適合布放聲吶陣列的海域海況環境。因而, 應重點關注和發展UWSN, 發揮其布放簡單、成本低廉、容錯和抗毀傷能力強的特點, 拓展水下跟蹤監視的范圍。

4) 隨著多源多目標跟蹤技術的不斷發展, 傳統基于數據關聯的濾波算法難以適應信息融合的需要。因此, 要重點研究基于RFS的跟蹤濾波算法, 突破計算復雜度高和多目標狀態提取困難的限制, 追趕國際先進水平, 加快其工程應用進度。

5) 高實時性多功能目標跟蹤系統。隨著各方面技術的發展和作戰需求牽引, 水下目標跟蹤系統呈現向集成目標探測識別、預警跟蹤和水下通信等綜合系統的方向發展。應加強DSP和現場可編程門陣列(field-programmable gate array, FPGA)等專用高速處理器的使用, 進行相應的高速計算語言和算法的關注, 以突破水下目標跟蹤數據量巨大, 跟蹤實時性差的瓶頸。

4 結束語

現代海戰是體系與體系的對抗, 準確感知戰場環境, 精準毀傷目標是決定戰爭勝敗的關鍵因素之一, 而水下目標跟蹤技術是其關鍵環節。文中根據跟蹤濾波算法原理的不同, 分別論述了基于數據關聯的水下多目標跟蹤技術和基于隨機有限集的水下多目標跟蹤技術, 詳細闡述了其目標運動模型、跟蹤濾波方法和應用現狀及其技術性能, 分析了存在的技術瓶頸。并根據作戰使用和海洋開發的要求, 展望了水下多目標跟蹤技術發展方向。同時可以看出, 基于隨機有限集的多目標跟蹤技術能夠較有效克服傳統多目標跟蹤算法的3類瓶頸問題, 仍處在理論向工程實踐轉化的關鍵階段, 應為技術人員著力重點突破的方向。

[1] Wax N.Signal-to-noise Improvement and the Statistics of Track Populations[J].Journal of Applied Physics, 1955, 26(5): 586-595.

[2] Singer R, Stein, J.An Optimal Tracking Filter for Proces- sing Sensor Data of Imprecisely Determined Origin in Surveillance Systems[C]//IEEE Conference on Decision and Control, Miami Beach: IEEE, 1971: 171-175.

[3] Fortmann T E, Bar-Shalom Y, Scheffe M.Sonar Tracking of Multiple Targets Using Joint Probabilistic Data Association[J].IEEE Journal of Oceanic Engineering, 1983, 8(3): 173-184.

[4] Reid D B.An Algorithm for Tracking Multiple Targets [J].IEEE Transactions on Automatic Control, 1979, 24(6): 843-854.

[5] Mahler R.Multi-target Bayes Filtering Via First-order Multi-target Moments[J].IEEE Transactions on Aeros-p- ace and Electronic Systems, 2003, 39(4): 1152-1178.

[6] Vo B-N, Singh S, Doucet A.Sequential Monte Carlo Methods for Multi-target Filtering with Random Finite sets[J].IEEE Transactions on Aerospace and Electronic Systems, 2005, 41(4): 1224-1245.

[7] 劉妹琴, 蘭劍.目標跟蹤前沿理論與應用[M].北京: 科學出版社, 2015.

[8] 周宏仁, 敬忠良, 王培德.機動目標跟蹤[M].北京: 國防工業出版社, 1991.

[9] Singer R A.Estimating Optimal Tracking Filter Performance for Manned Maneuvering Target[J].IEEE Transa-c-tions on Aerospace and Electronic Systems, 1970, 6(4): 473-483.

[10] 黨建武.水下多目標跟蹤理論[M].西安: 西北工業大學出版社, 2009.

[11] Moose R L.An adaptive state estimation solution to the maneuvering target problem .IEEE Transactions on Automatic Control, 1975, 20(6): 359-362.

[12] Blom H A P, Bar-Shalom Y.The IMM Algorithm for Syst- em with Markovian Switching Coefficients[J].IEEE Tr- ansactions on Automatic Control, 1988, 33(8): 780-783.

[13] 周宏仁.機動目標當前統計模型與自適應跟蹤算法[J].航空學報, 1983, 4(1): 73-86.Zhou Hong-ren.Current Statistical Model and Adaptive Tracking Algorithm for Maneuvering Targets[J].Acta Aeronautica et Astronautica Sinica, 1983, 4(1): 73-86.

[14] Song T L.Observability of Target Tracking with Bearing- only Measurements[J].IEEE Transactions On Aerospace & Electronis Systems, 1996, 32(4): 1468-1471.

[15] 嚴衛生, 高磊, 張礴, 等.水下目標跟蹤技術[J].火力與指揮控制, 2000, 25(2): 2-6.Yan Wei-sheng, Gao Lei, Zhang Bo, et al.Underwater Target Tracking Technology[J].Fire Control & Command Control, 2000, 25(2): 2-6.

[16] Bar-Shalom Y, Tse E.Tracking in a Cluttered Environment with Probabilistic Data Association[J].Automatica , 1975, 11(5): 451-460.

[17] Ristic B, Arulampalam S, Gordon N.Beyond the Kalman Filter: Particle Filters for Tracking Application[M].Bos ton London: Artech House, 2004.

[18] 高劍, 徐德民, 嚴衛生, 等.無跡卡爾曼濾波及其在三維水下目標跟蹤系統中的應用[J].船舶工程, 2005, 27 (3): 24-28.Gao Jian, Xu De-min, Yan Wei-sheng, et al.Unscented Kalman Filter and Its Application in 3D Underwater Target Tracking System[J].Marine Engineering, 2005, 27(3): 24- 28.

[19] 吳盤龍, 孔建壽.基于平方根UKF的水下純方位目標跟蹤[J].南京理大學學報(自然科學版), 2009, 33(6): 751-755.Wu Pan-long, Kong Jian-shou.Underwater Pure Azimuth Target Tracking Based on Square Toot UKF[J].Journal of Nanjing Institute of Science Natural Sciences, 2009, 33(6): 751-755.

[20] Popp R L, Pattipati K R, Bar-Shalom Y.M-Best S-D Assignment Algorithm with Application to Multitarget Trac- king[J].IEEE Transactions on Aerospace and Electronic Systems, 2001, 37(1): 22-39.

[21] Bar-Shalom Y.Update without of Sequence Measure-me- nts in Tracking: Exact Solution[J].IEEE Transactions on Aerospace and Electronic Systems, 2002, 38(3): 769-778.

[22] Li T, Bolic M, Djuric P D.Resampling Methods for Particle Filtering[J].IEEE Signal Processing Magazine, 2015, 32(3): 70-86.

[23] Li X R, Jilkov V P.Survey of Maneuvering Target Tracking.Part I: Dynamic Models[J].IEEE Transa-ctions on Ae- ro--space and Electronic Systems, 2003, 39(4): 1333-1364.

[24] Mahler R.Statistical Multi-source Multi-target Informat- ion Fusion[M].Norwood, MA: Artech House, 2007.

[25] Mahler R.PHD Filters of Higher Order in Target Number[J].IEEE Transactions on Aerospace and Electronic Systems, 2007, 43(4): 1523-1543.

[26] Vo B N, Ma W K.The Gaussian Mixture Probability Hypothesis Density Filter[J].IEEE Transactions on Signal Processing, 2006, 54(11): 4091-4104.

[27] Vo B, Cantoni A.The Cardinality Balanced Multi-Target Multi-Bernoulli Filter and Its Implementations[J].IEEE Trans.on Signal Processing, 2009, 57(2): 409-423.

[28] Mahler R.CPHD and PHD Filters for Unknown Backg-rounds, II: Multitarget Filtering in Dynamic Clutter[C]// Proceedings of SPIE Sensors and Systems for Space Applications III, Orlando, FL, USA: SPIE, 2009.

[29] Vo B, Ma W.The Gaussian Mixture Probability Hypo-thesis Density Filter[J].IEEE Transactions on Signal Pro- ce-ssing, 2006, 54(11): 4091-4104.

[30] Clark D, Vo B, Bell J.GM-PHD Filter Multi-target Trac- king in Sonar Images[C]//Processings of Signal Pro--ce---ss- ing, Sensor Fusion , and Target Recognition XV, Orlando (Kissimmee), FL, USA: SPIE, 2006.

[31] 陳里銘.基于隨機有限集的多目標跟蹤方法研究[D].大連: 大連理工大學, 2013.

[32] 連峰.基于隨機有限集的多目標跟蹤方法研究[D].西安: 西安交通大學, 2009.

[33] 秦永.基于隨機有限集目標跟蹤算法研究[D].武漢: 華中科技大學, 2013.

[34] 李翠蕓, 王榮, 姬紅兵.基于變分貝葉斯勢均衡多目標多伯努利濾波的多擴展目標跟蹤算法[J].控制理論與應用, 2015, 32(2): 187-195.Li Cui-yun, Wang Rong, Ji Hong-bing.Based on Variable Decibel Equilibrium Multi-objective Bernoulli Filter, Multi Extended Target Tracking Algorithm[J].Control Theory and Application, 2015, 32 (2): 187-195.

[35] Vo B T, Vo B N, Hoseinnezhad R, et al.Multi-Bernoulli Filtering with Unknown Clutter Intensity and Sensor Fi- eld-of-view[C]//Proceedings of IEEE 45th Annual Con-fe-rence on Information Sciences and Systems, Baltimore, America: IEEE, 2011.

[36] Vo B T, Vo B N, Hoseinnezhad R.Robust Multi-Bernoulli Filtering[J].IEEE Journal of Selected Topics in Signal Processing, 2013, 7(3): 399-409.

[37] Mahler R, Vo B T.An Improved CPHD Filter for Unknown Clutter Backgrounds[C]//Proceedings of SPIE Signal Proce- ssing, Sensor/Information Fusion, and Target Recognition XXIII.Baltimore, Maryland, USA: SPIE, 2014.

[38] Georgescu R.Willett P.The GM-CPHD Tracker Applied to Real and Realistic Multistatic Sonar Data Sets[J].IEEE Journal of Oceanic Engineering, 2012, 37(2): 220-235.

[39] Beard M, Vo B-T, Vo B-N, et al.Arulampalam, A Partially Uniform Target Birth Model for Gaussian Mixture PHD/ CPHD Filtering[J].IEEE Transactions on Aerospace and Electronic Systems, 2013, 49(4): 2835-2844.

[40] Mahler R.Proceedings CPHD Filters for Unknown Clutter and Target-birth Processes[C]//Signal Processing, Sensor/In- formation Fusion, and Target Recognition XXIII.Baltimore, Maryland, USA: SPIE, 2014.

[41] Vo B N, Vo B T.Random Finite Set Multi-target Trackers: Stochastic Geometry for Space Situational Awareness[C]// Proceedings of Signal Processing, Sensor/Information Fusion, and Target Recognition XXIV: SPIE, 2015.

[42] Beard M, Vo B T, Vo B N.Bayesian Measurements Using Labelled Random Finite Multi-target Tracking with Sets [J].IEEE Transactions on Merged Signal Processing, 2015, 63(6): 1433-1447.

[43] Beard M, Vo B-T, Vo B-N.Bayesian Measurements Using Labelled Random Finitemulti-target Tracking with Sets [J].IEEE Transactions on merged Signal Processing, 2015, 63(6): 1433-1447.

[44] Hoffman J, Mahler R.Multitarget Miss Distance Via Optimal Assignment IEEE Transactions on Systems[J].Man and Cybernetic, 2004, 34(3): 327-336.

[45] Schuhmacher D, Vo B T, Vo B N.A Consistent Me-tric for Performance Evaluation of Multi-object Filters[J].IEEE Transactions on Signal Processing, 2008, 56(8): 3447- 3457.

[46] Nagappa S, Clark D E, Mahler R.Incorporating Track Uncertainty into the OSPA Metric[C]//14th International Conference on Information Fusion, Chicago IL, USA: IEEE, 2011: 1-8.

[47] Ristic B, Vo B N, Clark D, et al.A Metric for Perfor-ma- nce Evaluation of Multi-target Tracking Algorithms[J].IEEE Transactions on Signal Processing, 2011, 59(7): 3452-3457.

[48] He X, Tharmarasa R, Kirubarajan, et al.A Track Quality Based Metric for Evaluating Performance of Multi-target Filters[J].IEEE Transactions on Aerospace and Electronic Systems, 2013, 49(1): 610-616.

Present Situation and Prospect of Underwater Multi-Target Tracking Technologies

ZHANG Si-yu, HE Xin-yi, ZHANG Chi, ZHU Lin, CHEN Shuang

(Naval Research Academy, Beijing 100161, China)

Aiming at the unfavorable factors, such as the complexity of underwater hydrological conditions and the relatively small action distance, the key underwater multi-target tracking technologies based on respective data association and random finite set are discussed according to different principles of tracking filtering algorithms.The target motion models, tracking filtering algorithms and its applications status are expounded.The performances of two multi-target tracking technologies are analyzed.Three kinds of technical bottlenecks are analyzed, including the uncertainty of observation information due to failing to report and false alarm, the uncertainty of target number and motion state change in tracking process, and the poor real-time performance of tracking algorithm.The advantages of establishing a unified random finite set framework to describe the tracking problem and solve these bottlenecks are emphasized.Further, according to the requirements of operational application and ocean development, the development direction of underwater multi-target tracking technology is prospected to provide a reference for relevant researchers.

underwater multi-target tracking; data association; random finite set

TJ630.34; TN953

R

2096-3920(2018)06-0511-10

10.11993/j.issn.2096-3920.2018.06.001

2018-05-28;

2018-12-07.

張思宇(1992-), 男, 在讀碩士, 主要研究方向為魚雷自導與水聲信號處理.

張思宇, 何心怡, 張馳, 等.水下多目標跟蹤技術現狀與展望[J].水下無人系統學報, 2018, 26(6): 511-520.

(責任編輯: 楊力軍)