王帥祥

VICKF-IMM算法在機動目標跟蹤中的應用

王帥祥

（中北大學機電工程學院，太原，030051）

針對智能體移動方式復雜，對其進行觀測的傳感器測量的信息存在噪聲以及目標運動軌跡發(fā)生突然的改變會導致目標觀測失真甚至錯誤的問題，提出了一種變積容積卡爾曼濾波交互多模型算法（VICKF-IMM）。該算法將容積卡爾曼濾波與交互多模型算法相結(jié)合，并對容積卡爾曼濾波（CKF）中球面積分進行變積分轉(zhuǎn)換處理。優(yōu)化了其積分求解的方式，提高了整體的穩(wěn)定性。Monte-Carlo仿真分析，與CKF-IMM和UKF-IMM算法相比，該算法的跟蹤精度有明顯的提高，并在目標運動發(fā)生突變時有更高的穩(wěn)定性。

SCKF-IMM；機動目標跟蹤；容積卡爾曼濾波；交互多模型算法

0 引 言

信息化時代針對智能體的跟蹤定位是當前研究的新目標，但智能體的移動往往不是單一模式下的運動，其運動模式復雜多變，且其運動模式隨時間變化而改變。并且在對目標的觀測過程中往往會受到量測相關(guān)的噪聲影響，造成定位精度的下降。

為了使算法能準確描述目標的狀態(tài)，1988年，Blom等[1]提出交互多模型算法。但交互多模型算法一般會配合濾波算法進行使用，傳統(tǒng)的卡爾曼濾波針對非線性系統(tǒng)的解決能力不足，所以Arasaratnam將貝葉斯估計和容積變換有機結(jié)合提出了一種非線性濾波算法[2]，即容積卡爾曼濾波算法（Cubature Kalman Filter，CKF），它能計算出一些具有相同權(quán)重的點，并可直接通過非線性系統(tǒng)方程對系統(tǒng)中的狀態(tài)和誤差進行迭代的估計，其精度一般可達到3階展開的程度[3]。

但容積卡爾曼濾波算法必須將非線性函數(shù)分解成球面積分和徑向積分的形式，求解其積分的過程比較復雜，傳統(tǒng)方式精度不高。

針對這種類型情況，本文提出變積容積卡爾曼濾波交互多模型算法（Variable Cubature Kalman Filter Interactive Multiple Model Algorithm，VICKF-IMM）。該算法以CKF為基礎(chǔ)，改變了其中徑向積分的求解方式，再融合自適應交互多模型算法，使其更加適應復雜多變的智能體目標跟蹤領(lǐng)域。

1 系統(tǒng)模型

當目標個體在觀測區(qū)域進行運動時，一般的觀測器如雷達、紅外等可以將其在圖像中投影位置獲得的偏置角度作為其相關(guān)變量。其觀測目標如圖1所示。

圖1 觀測器觀測目標

由圖1中角度相關(guān)關(guān)系可以推導出觀測器觀測出來的估計夾角。其推導公式如下：

笛卡爾坐標系下對目標運動進行建模。在對移動智能體目標進行運動學分析時，可以假設目標以某種運動模型進行運動，對目標建模一般會達到更好的跟蹤效果，在考慮到觀測器觀測時，要考慮到其噪聲的存在，故建立的狀態(tài)方程和觀測方程如下：

2 變積容積卡爾曼濾波算法

2.1 容積卡爾曼濾波算法

CKF算法的核心思想是采用球面-徑向容積規(guī)則，近似處理非線性傳遞的后驗估計值和相應的協(xié)方差。

CKF存在計算容積點，而后利用2個容積位置計算出來的點進行加權(quán)求和來近似高斯積分。其數(shù)學表示為

2.2 傳統(tǒng)求解方式

其球面積分可以用容積準則來表示：

2.3 改進積分求解方式

由于傳統(tǒng)的徑向基求解方式，較為復雜，這里對徑向積分進行變積分求解。引用定理1——相徑積分可通過r個點的高斯形式的求積準則，進行求解工作，其求解公式為

徑向積分可以通過高斯-拉爾蓋積分進行表示：

并且因為球面容積準則的對稱性，可以將式（12）的結(jié)論部分進行簡化：

使用這種方式，可簡化徑向積分的求解方式，且與交互多模型算法相結(jié)合，就構(gòu)成VICKF-IMM算法的全部過程。

3 數(shù)字實驗仿真與分析

3.1 仿真條件

假設目標運動模式的模型集中含有勻速直線運動（CV）、勻加速直線運動（CA）和協(xié)同的圓周轉(zhuǎn)彎運動（CT）等運動模式集合，并且目標在運動過程中其運動模式，會隨時發(fā)生改變。

假定當前選定的目標的起始位置為（30 000， 20 000），目標的初始速度在軸分量V=80，在軸分量V=50，且目標將在上一節(jié)介紹的3種運動形式中隨機產(chǎn)生運動變化，預設目標運行的總時間為350 s。

模型狀態(tài)轉(zhuǎn)移矩陣的初始化設置為針對CA，CV，CT模型的概率狀態(tài)初始值設置為[0.6；0.2；0.2]。相應的交互多模型中的馬爾科夫概率轉(zhuǎn)移矩陣的初始為

本文將仿真中的噪聲簡化為高斯白噪聲，其均值為0，方差為Q=0.13I，采用蒙特卡洛方法進行120次仿真實驗。本次仿真中將CKF-IMM、UKF-IMM和VICKF-IMM算法進行對比。實際仿真過程目標運動的位置如圖2所示。

從圖2可分析出，3種算法都能對包含CV、CA和CT模型集的目標機動進行跟蹤。下面將詳細分析其算法的優(yōu)劣。

算法性能評價指標在大多數(shù)標準中為均方根誤差：

目標運動實際參數(shù)軌跡中加速度預設情況見表1。

表1 探測中的目標運動加速度變化情況

Tab.1 Changes in target motion acceleration during detection

目標機動時刻sX方向加速度m·s-2Y方向加速度m·s-2Z方向加速度m·s-2 13613160 150-680 2639-90 289430 302-580

3.2 仿真分析

仿真結(jié)果中，本文算法VICKF-IMM、CKF-IMM和UKF-IMM都可以對目標運動進行有效的跟蹤，其位置均方根誤差對比如圖3所示。由圖3分析可知，VICKF-IMM算法的定位精度比UKF-IMM算法高61.1%，VICKF-IMM算法的定位精度比CKF-IMM算法提高51.6%，其位置均方根誤差平均下降60.1%。

圖3 不同算法位置均方根誤差對比

CKF-IMM和UKF-IMM算法中，其定位位置在40 s發(fā)生突變，而本文算法的定位精度改變較小，所以VICKF-IMM算法的穩(wěn)定性比傳統(tǒng)的CKF-IMM算法和UKF-IMM算法的高，魯棒性表現(xiàn)更好。

圖4為VICKF-IMM、CKF-IMM和UKF-IMM算法的速度均方根誤差。圖4表明其速度均方根誤差下降32%。

圖4 不同算法速度均方根誤差對比

4 結(jié) 論

綜上所述，本文提出的方法優(yōu)點可總結(jié)為以下幾方面：

a）通過改變并簡化求其球面積分方式，簡化了計算難度，并提高了相關(guān)跟蹤精度，在CA、CV、CT傳統(tǒng)的模型級突變算例中，其跟蹤精度在位置屬性上提高約60%，速度屬性的精度提高了32%；

b）優(yōu)化了傳統(tǒng)算法的處理步驟，避免了傳統(tǒng)CKF算法中計算數(shù)值時誤差累計造成的過程誤差分解的錯誤，VICKF-IMM算法穩(wěn)定性（魯棒性）相比傳統(tǒng)型算法有所提高。

[1] 張冉, 杜清, 萬岳, 等. 一種基于交互多模型卡爾曼濾波的被動協(xié)同定位濾波算法研究[J]. 雷達與對抗, 2023, 43(1): 9-12.

ZHANG Ran, DU Qing, WAN Yue, et al. Research on passive collaborative positioning filtering algorithm based on interactive multi-model Kalman filter[J]. Radar and Countermeasures, 2023, 43(1): 9-12.

[2] 封普文, 黃長強, 曹林平, 等. 基于自適應馬爾可夫參數(shù)交互多模型算法的彈道導彈跟蹤研究[J]. 兵工學報, 2014, 35(12): 2041-2049.

FENG Puwen, HUANG Changqiang, CAO Linping, et al. Research on ballistic missile tracking based on adaptive markov parameter interactive multi-model algorithm[J]. Acta Armamentarii, 2014, 35(12): 2041-2049.

[3] 蘭義華, 任浩征, 張勇, 等. 一種基于“當前”模型的改進卡爾曼濾波算法[J]. 山東大學學報（工學版）, 2012, 42(5): 12-17.

LAN Yihua, REN Haozheng, ZHANG Yong, et al. An improved Kalman filter algorithm based on "current" model[J]. Journal of Shandong University (Engineering Science), 2012, 42(5): 12-17.

[4] Aryal A, Ghahramani A, Becerik Gerber B. Monitoring fatigue in construction workers using physiological measurements[J]. Automation in Construction, 2020, 82(10): 2999-3007.

Application of VICKF-IMM Algorithm in Maneuvering Target Tracking

WANG Shuaixiang

(College of Mechanical and Electrical Engineering, North University of China, Taiyuan, 030051)

In the process of tracking targets with complex and variable motion patterns during movement, there is noise in the information measured by the sensors that observe them, and sudden changes in the target's motion trajectory can lead to distortion or even errors in the target observation. A variable cubature Kalman filter interacting multiple model algorithm (VICKF-IMM) is proposed. This algorithm combines volumetric Kalman filtering with interactive multiple model algorithm, and performs variable integral transformation on the spherical surface integral in cubature Kalman filtering. The optimization of its integral solution method has improved overall stability. Monte Carlo simulation analysis shows that compared with CKF-IMM and UKF-IMM algorithms, the tracking accuracy of this algorithm is significantly improved, and it has higher stability in case of sudden changes in target motion.

SCKF-IMM; maneuvering target tracking; cubature Kalman filter; interactive multiple model algorithm

2097-1974(2023)02-0016-04

10.7654/j.issn.2097-1974.20230204

TN953

A

2023-04-01；

2023-04-03

王帥祥（1996-），男，碩士研究生，主要研究方向為探測與預測。