CN103487800B - 基于残差反馈的多模型高速高机动目标跟踪方法 - Google Patents
基于残差反馈的多模型高速高机动目标跟踪方法 Download PDFInfo
- Publication number
- CN103487800B CN103487800B CN201310404989.8A CN201310404989A CN103487800B CN 103487800 B CN103487800 B CN 103487800B CN 201310404989 A CN201310404989 A CN 201310404989A CN 103487800 B CN103487800 B CN 103487800B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- matrix
- target
- sampling point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 239000011159 matrix material Substances 0.000 claims abstract description 84
- 239000013598 vector Substances 0.000 claims abstract description 35
- 238000005070 sampling Methods 0.000 claims description 57
- 238000004422 calculation algorithm Methods 0.000 claims description 25
- 230000007704 transition Effects 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 description 8
- 230000001133 acceleration Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000013179 statistical model Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
- G01S13/522—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves
- G01S13/524—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi
- G01S13/526—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi performing filtering on the whole spectrum without loss of range information, e.g. using delay line cancellers or comb filters
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种基于残差反馈的多模型高速高机动目标跟踪方法,主要解决现有技术计算量大的问题。其步骤是:1)初始化雷达观测值,得到初始状态向量和初始协方差矩阵;2)选择若干运动模型,确定运动模型的状态转移矩阵;3)由运动模型的个数假定初始权值;4)利用状态向量和状态转移矩阵作一步预测;5)利用初始权值对一步预测值加权求和,得到状态估计向量和估计协方差矩阵;6)根据观测值和状态估计向量得到预测误差;7)由预测误差和估计协方差矩阵计算增益矩阵;8)利用增益矩阵和状态估计向量更新状态向量和状态协方差矩阵;9)利用预测误差和状态估计向量更新权值。本发明跟踪性能稳定、计算量小,可用于对临近空间目标的跟踪。
Description
技术领域
本发明属于雷达技术领域,具体涉及一种对强机动目标的跟踪方法,可用于对临近空间目标的跟踪。
背景技术
由于目标跟踪在国防和民用各领域的重要价值,一直以来很多学者致力于研究目标跟踪算法。随着现代航空的迅速发展,各种飞行器的机动性和运动速度越来越高,但是可靠而精确的跟踪目标始终是目标跟踪***设计的主要目的与难点。目标跟踪实际上就是对目标状态的跟踪滤波问题,而实现目标跟踪首先要使所建立的目标运动模型与实际的目标运动模型匹配,建立目标的运动模型和自适应滤波是机动目标的两个关键部分。
在当前目标机动能力日益增强的情况下,单模型算法很难满足跟踪速度和精度的需求,而多模型算法可以避免采用单模型时由于目标机动而造成模型的不准确,提高机动目标跟踪性能,从而实现对强机动目标的精确跟踪。经过国内外学者几十年来的不断研究,提出了很多目标模型与算法。目前广泛采用的算法是1984-1989年Blom和Bar-Shalom提出的交互式多模型IMM算法,该算法具有Markov转移概率,且算法中有多个模型并行工作,多个滤波器交互作用得到目标状态估计的结果,由于该模型有较好的自适应能力,因而比较理想。但是IMM算法中由于其模型转移概率是人为设置的,因此引入了人为误差,影响了跟踪性能。除此之外,因为IMM算法需要对每一个模型都进行滤波,计算量过大,对目标跟踪会有延迟。
发明内容
本发明的目的在于针对上述已有技术的不足,提出一种基于残差反馈的多模型高速高机动目标跟踪方法,以降低计算复杂度,提高雷达在目标强机动下的跟踪精度。
实现本发明的技术思路是:利用LMS算法,通过信号的预测残差调整滤波器的权值,从而使估计信号逼近期望信号,以提高跟踪性能,其实现步骤包括如下:
(1)雷达对接受到的机动目标的位置检测值进行N次采样,得到一个长度为N的观测值序列{Z(k)},k=1,2,…,N;
(2)从观测值序列{Z(k)}中取出前三个采样点的观测值Z(1),Z(2),Z(3),并利用差分法初始化机动目标的运动状态,得到机动目标在第3个采样点的状态向量和初始协方差矩阵P(3);
(3)根据机动目标的运动特性选择M个目标运动模型,得到M个目标运动模型的状态转移矩阵Fi,i=1,2,…M,M≥2;
(4)根据选择的目标运动模型个数M,确定目标运动模型在第3采样点的权值:
W(3)=[1/M 1/M … 1/M]1×M,
其中,[·]1×M表示该矩阵为1行M列的矩阵;
(5)利用步骤(2)中得到第k采样点的状态向量和步骤(3)中得到的目标运动模型的状态转移矩阵Fi进行一步预测,得到对应模型的一步预测值
其中,表示第i个目标运动模型在第k采样点预测的机动目标在第k+1采样点的状态向量;
(6)利用步骤(4)中得到的权值对步骤(5)中得到的各个目标运动模型的一步预测值和协方差矩阵P(k)进行加权求和,得到基于LMS算法的状态估计向量和估计协方差矩阵P(k+1/k):
其中,wi(k)表示在第k采样点权值W(k)的第i个元素;[·]T表示矩阵的转置;P(k)表示机动目标在第k采样点的协方差矩阵;Qi(k)表示第i个目标运动模型在第k采样点的过程噪声协方差;
(7)根据步骤(1)中得到的雷达对机动目标位置的观测值Z(k+1)和步骤(6)中得到的基于LMS算法的状态估计向量通过下式求解机动目标在第k+1采样点的预测误差
其中,其表示第k采样点预测机动目标在第k+1采样点的观测值;Z(k+1)表示雷达在第k+1采样点对机动目标位置的观测值;H表示机动目标的量测矩阵;
(8)由步骤(7)中得到的预测误差和步骤(6)中得到的估计协方差矩阵P(k+1/k)确定预测误差的协方差矩阵S(k+1),继而得到增益矩阵G(k+1):
S(k+1)=HP(k+1/k)HT+R(k+1)
G(k+1)=P(k+1/k)HTS-1(k+1),
其中,S(k+1)是预测误差在第k个采样点的协方差矩阵;R(k+1)是观测噪声在第k+1采样点的协方差矩阵;G(k+1)表示第k+1采样点的滤波增益矩阵;[·]-1表示矩阵的逆;
(9)利用步骤(8)中得到的增益矩阵G(k+1)和步骤(6)中得到的状态估计向量完成对状态向量和协方差矩阵P(k)的更新,即:
P(k+1)=P(k+1/k)-G(k+1)S(k+1)GT(k+1)
其中,表示机动目标更新后在第k+1采样点的状态向量;P(k+1)表示机动目标在第k+1采样点的协方差矩阵;[·]T表示矩阵的转置;
(10)通过步骤(7)中得到的预测误差和步骤(5)中得到的目标运动模型的一步预测值对权值W(k)进行更新,得到更新后的权值:
其中,W(k+1)表示第k+1采样点目标运动模型的权值,μ表示步长因子,M为采用的目标运动模型的个数;
(11)判断跟踪是否完成,即判断k是否满足k≤N-1,若条件满足,则k递增,并继续执行步骤(5);若条件不满足,则跟踪过程结束。
本发明与现有技术相比的优点在于:
1.跟踪性能高
本发明中引入了类似IMM算法中转移概率的权值,但与传统IMM算法中转移概率固定不同的是,本发明中滤波器可以自适应调节权值,从而提高了跟踪精度,减小了误差,同时提高了对机动目标的跟踪性能。
2.计算简单
现有IMM算法中每次滤波时对每个目标运动模型都要进行卡尔曼滤波,即每次滤波过程需要进行M次卡尔曼滤波;本发明中每个目标运动模型只需要进行一步预测,每次滤波过程都只需要进行一次卡尔曼滤波,大幅度降低了计算复杂度。
附图说明
图1为本发明的流程图;
图2为本发明对机动目标跟踪的效果示意图;
图3为本发明与IMM方法对机动目标跟踪三维位置误差局部比较图;
具体实施方式
参照图1,本发明实现步骤如下:
步骤1,对运动状态和模型参数进行初始化。
1a)雷达对其接收到的机动目标的位置检测值进行N次采样,得到一个长度为N的观测值序列{Z(k)},k=1,2,…,N;
1b)从观测值序列中取出前三个采样点的观测值Z(1),Z(2),Z(3),其中:
Z(1)=[xx(1) yy(1) zz(1)]T
Z(2)=[xx(2) yy(2) zz(2)]T
Z(3)=[xx(3) yy(3) zz(3)]T
式中xx(·),yy(·),zz(·)分别表示机动目标在x轴、y轴、z轴的位置,[·]T表示矩阵的转置;
1c)利用差分法对Z(1),Z(2),Z(3)进行运算,获得初始状态向量如下:
其中,(xx(3)-xx(2))/T、(yy(3)-yy(2))/T、(zz(3)-zz(2))/T分别表示机动目标在第3采样点处的x轴、y轴、z轴方向速度;((xx(3)-xx(2))/T-(xx(2)-xx(1))/T)/T、((yy(3)-yy(2))/T-(yy(2)-yy(1))/T)/T、((zz(3)-zz(2)/T)-(zz(3)-zz(2)/T))/T分别表示机动目标在第3采样点的x轴、y轴、z轴方向加速度;T为采样间隔;
1d)根据初始状态向量获得初始协方差矩阵P(3):
其中,E[·]表示数学期望;[·]T表示矩阵的转置。
步骤2,选取目标运动模型,确定状态转移矩阵。
常用的目标运动模型包括二阶匀速CV模型,三阶匀加速CA模型,协调转弯CT模型,singer模型和“当前”统计模型,根据机动目标的运动特性,在常用的目标运动模型中选取M个目标运动模型,得到M个目标运动模型的状态转移矩阵Fi,i=1,2,…M,M≥2。
步骤3,由目标运动模型的个数,确定目标运动模型在第3采样点的权值W(3):
W(3)=[1/M 1/M … 1/M]1×M,
其中,[·]1×M表示该矩阵为1行M列的矩阵。
步骤4,利用状态向量和状态转移矩阵进行一步预测,得到对应模型的一步预测值
其中,表示第i个目标运动模型在第k个采样点预测的机动目标在第k+1采样点的状态向量;表示机动目标在第k个采样点的状态向量;Fi表示第i个目标运动模型的状态转移矩阵。
步骤5,计算基于LMS算法的最优状态估计向量和估计协方差矩阵。
5a)对各个目标运动模型的一步预测值进行加权求和,得到基于LMS算的最优状态估计向量
其中,wi(k)表示在第k采样点权值W(k)的第i个元素;
5b)对协方差矩阵P(k)进行加权求和,得到基于LMS算法的估计协方差矩阵P(k+1/k):
其中,P(k)表示机动目标在第k采样点的协方差矩阵;Qi(k)表示第i个目标运动模型在第k采样点的过程噪声协方差:
Qi(k)=Ci(k)qCi(k)T,
其中,Ci(k)表示第i个目标运动模型在第k采样点对应的控制矩阵,由选取的运动模型不同而取不同的矩阵;q为过程噪声方差,由机动目标位置的数量级确定。
步骤6,计算机动目标观测值的新息。
根据雷达对机动目标位置的观测值Z(k+1)和基于LMS算法的最优状态估计向量用下式求解机动目标预测值在第k+1采样点的预测误差
其中,其表示第k采样点预测机动目标在第k+1采样点的观测值;Z(k+1)表示雷达在第k+1采样点对机动目标位置的观测值;H表示机动目标的量测矩阵:
步骤7,计算增益矩阵G(k+1)。
7a)由估计协方差矩阵P(k+1/k)确定预测误差的协方差矩阵S(k+1)
S(k+1)=HP(k+1/k)HT+R(k+1),
其中,S(k+1)是预测误差在第k个采样点的协方差矩阵;R(k+1)是观测噪声在第k+1采样点的协方差矩阵;
7b)由预测误差的协方差矩阵S(k+1)和估计协方差矩阵P(k+1/k),计算增益矩阵G(k+1):
G(k+1)=P(k+1/k)HTS-1(k+1),
其中,G(k+1)表示第k+1采样点的增益矩阵;[·]-1表示矩阵的逆。
步骤8,对机动目标状态向量和状态协方差矩阵的进行更新。
8a)利用增益矩阵G(k+1)、状态估计向量和预测误差完成对状态向量的更新,即:
其中,表示机动目标更新后在第k+1采样点的状态向量;
8b)利用增益矩阵G(k+1)、估计协方差矩阵P(k+1/k)和预测误差的协方差矩阵S(k+1)完成对协方差矩阵P(k)的更新,即:
P(k+1)=P(k+1/k)-G(k+1)S(k+1)GT(k+1)
其中,P(k+1)表示机动目标更新后在第k+1采样点的协方差矩阵。
步骤9,通过预测误差和目标运动模型的一步预测值完成对权值W(k)的更新,即:
其中,W(k+1)表示更新后得到的第k+1采样点机动目标运动模型的权值,μ表示步长因子,保证算法收敛的μ的取值范围为:0<μ<1/λmax,λmax为矩阵xin的自相关矩阵的最大特征值,
步骤10,判断跟踪是否完成,即判断k是否满足k≤N-1,若条件满足,则k递增,并继续执行步骤(4);反之,则跟踪过程结束。
本发明的效果可以通过以下仿真进一步说明:
1.仿真条件:
假设目标的起始位置为[50 125 75]km,起始速度为[0.7 0.65 0.6]km/s,x轴方向的加速度a1为0.15sin(0.04πt)km/s2,y轴方向的加速度a2为0.1cos(0.04πt)km/s2,z轴方向的加速度a3为0.12sin(0.02πt)km/s2。距离误差的标准差为方位角误差的标准差为仰角测量误差的标准差为对目标进行30s的观测,采样间隔T为0.25s,Monte Carlo仿真100次。
选择的运动模型是singer模型、“当前”统计模型和匀加速CA模型。IMM算法中的概率转移矩阵p为:
2.仿真内容
仿真1,采用本发明方法对仿真条件中所述的目标进行跟踪仿真,结果如图2。从图2可以看出,本发明方法的跟踪轨迹很接近机动目标运动的真实轨迹。
仿真2,对用本发明方法与传统IMM方法的跟踪误差进行比较,结果如图3。其中图3(a)表示本发明方法和传统IMM算法在x轴上的跟踪误差比较;图3(b)表示本发明方法和传统IMM算法在y轴上的跟踪误差比较;图3(c)表示本发明方法和传统IMM算法在z轴上的跟踪误差比较。
图3表明本发明方法的跟踪误差在整体趋势上小于IMM算法的跟踪误差,且跟踪稳定,更适合对高速高机动目标的跟踪。
Claims (1)
1.一种基于最小均方误差的强机动目标跟踪方法,包括如下步骤:
(1)雷达对接收到的机动目标的位置检测值进行N次采样,得到一个长度为N的观测值序列{Z(k)},k=1,2,",N;
(2)从观测值序列{Z(k)}中取出前三个采样点的观测值Z(1),Z(2),Z(3),并利用差分法初始化机动目标的运动状态,得到机动目标在第3个采样点的状态向量和初始协方差矩阵P(3);
(3)根据机动目标的运动特性选择M个目标运动模型,得到M个目标运动模型的状态转移矩阵Fi,i=1,2,…M,M≥2;
(4)根据选择的目标运动模型个数M,确定目标运动模型在第3采样点的权值:
W(3)=[1/M 1/M … 1/M]1×M,
其中,[·]1×M表示该矩阵为1行M列的矩阵;
(5)利用第k采样点的状态向量和步骤(3)中得到的目标运动模型的状态转移矩阵Fi进行一步预测,得到对应模型的一步预测值
其中,表示第i个目标运动模型在第k采样点预测的机动目标在第k+1采样点的状态向量;
(6)利用权值对步骤(5)中得到的各个目标运动模型的一步预测值和协方差矩阵P(k)进行加权求和,得到基于LMS算法的状态估计向量和估计协方差矩阵P(k+1/k):
其中,wi(k)表示在第k采样点权值W(k)的第i个元素;[·]T表示矩阵的转置;P(k)表示机动目标在第k采样点的协方差矩阵;Qi(k)表示第i个目标运动模型在第k采样点的过程噪声协方差;
(7)根据步骤(1)中得到的雷达对机动目标位置的观测值Z(k+1)和步骤(6)中得到的基于LMS算法的状态估计向量通过下式求解机动目标在第k+1采样点的预测误差
其中,其表示第k采样点预测机动目标在第k+1采样点的观测值;Z(k+1)表示雷达在第k+1采样点对机动目标位置的观测值;H表示机动目标的量测矩阵;
(8)由步骤(7)中得到的预测误差和步骤(6)中得到的估计协方差矩阵P(k+1/k)确定预测误差的协方差矩阵S(k+1),继而得到增益矩阵G(k+1):
S(k+1)=HP(k+1/k)HT+R(k+1)
G(k+1)=P(k+1/k)HTS-1(k+1),
其中,S(k+1)是预测误差在第k个采样点的协方差矩阵;R(k+1)是观测噪声在第k+1采样点的协方差矩阵;G(k+1)表示第k+1采样点的滤波增益矩阵;[·]-1表示矩阵的逆;
(9)利用步骤(8)中得到的增益矩阵G(k+1)和步骤(6)中得到的状态估计向量完成对状态向量和协方差矩阵P(k)的更新,即:
P(k+1)=P(k+1/k)-G(k+1)S(k+1)GT(k+1)
其中,表示机动目标更新后在第k+1采样点的状态向量;P(k+1)表示机动目标在第k+1采样点的协方差矩阵;[·]T表示矩阵的转置;
(10)通过步骤(7)中得到的预测误差和步骤(5)中得到的目标运动模型的一步预测值对权值W(k)进行更新,得到更新后的权值:
其中,W(k+1)表示第k+1采样点目标运动模型的权值,μ表示步长因子,M为采用的目标运动模型的个数;
(11)判断跟踪是否完成,即判断k是否满足k≤N-1,若条件满足,则k递增,并继续执行步骤(5);若条件不满足,则跟踪过程结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310404989.8A CN103487800B (zh) | 2013-09-08 | 2013-09-08 | 基于残差反馈的多模型高速高机动目标跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310404989.8A CN103487800B (zh) | 2013-09-08 | 2013-09-08 | 基于残差反馈的多模型高速高机动目标跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103487800A CN103487800A (zh) | 2014-01-01 |
CN103487800B true CN103487800B (zh) | 2015-09-30 |
Family
ID=49828169
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310404989.8A Expired - Fee Related CN103487800B (zh) | 2013-09-08 | 2013-09-08 | 基于残差反馈的多模型高速高机动目标跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103487800B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104199022B (zh) * | 2014-08-30 | 2017-01-11 | 电子科技大学 | 一种基于目标模态估计的临近空间高超声速目标跟踪方法 |
CN105929391B (zh) * | 2016-04-24 | 2018-08-21 | 西安电子科技大学 | 基于后验概率的变结构交互多模型雷达机动目标跟踪方法 |
CN105954743B (zh) * | 2016-05-31 | 2018-08-21 | 西安电子科技大学 | 一种改进权值的变结构多模型机动目标跟踪方法 |
CN108037663A (zh) * | 2017-12-07 | 2018-05-15 | 电子科技大学 | 一种基于变维滤波算法的机动目标跟踪方法 |
CN111157983B (zh) * | 2020-01-06 | 2022-02-15 | 南京鹰目电子科技有限公司 | 一种雷达目标跟踪方法 |
CN111707997B (zh) * | 2020-06-03 | 2022-05-27 | 南京慧尔视智能科技有限公司 | 雷达目标跟踪方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101661104A (zh) * | 2009-09-24 | 2010-03-03 | 北京航空航天大学 | 基于雷达/红外量测数据坐标转换的目标跟踪方法 |
CN102721951A (zh) * | 2012-05-04 | 2012-10-10 | 西安电子科技大学 | 一种高机动目标跟踪方法 |
-
2013
- 2013-09-08 CN CN201310404989.8A patent/CN103487800B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101661104A (zh) * | 2009-09-24 | 2010-03-03 | 北京航空航天大学 | 基于雷达/红外量测数据坐标转换的目标跟踪方法 |
CN102721951A (zh) * | 2012-05-04 | 2012-10-10 | 西安电子科技大学 | 一种高机动目标跟踪方法 |
Non-Patent Citations (2)
Title |
---|
使用交互多模型算法的高速高机动目标跟踪;易令等;《雷达科学与技术》;20060630;第4卷(第3期);143-147 * |
基于时变马尔科夫转移概率的机动目标多模型跟踪;左东广等;《西安交通大学学报》;20030831;第37卷(第8期);824-828 * |
Also Published As
Publication number | Publication date |
---|---|
CN103487800A (zh) | 2014-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103487800B (zh) | 基于残差反馈的多模型高速高机动目标跟踪方法 | |
CN111985093B (zh) | 一种带噪声估计器的自适应无迹卡尔曼滤波状态估计方法 | |
CN104567871B (zh) | 一种基于地磁梯度张量的四元数卡尔曼滤波姿态估计方法 | |
CN103644903B (zh) | 基于分布式边缘无味粒子滤波的同步定位与地图构建方法 | |
CN105205313B (zh) | 模糊高斯和粒子滤波方法、装置及目标跟踪方法、装置 | |
CN104035095B (zh) | 基于空时最优处理器的低空风切变风速估计方法 | |
CN102622520B (zh) | 一种机动目标跟踪的分布式多模型估计融合方法 | |
CN111783307B (zh) | 一种高超声速飞行器状态估计方法 | |
CN105929391B (zh) | 基于后验概率的变结构交互多模型雷达机动目标跟踪方法 | |
CN101221238B (zh) | 基于高斯均值移动配准的动态偏差估计方法 | |
CN104793201B (zh) | 一种跟踪临近空间高超声速目标的修正变结构网格交互多模型滤波方法 | |
CN105222780B (zh) | 一种基于Stirling插值多项式逼近的椭球集员滤波方法 | |
CN104408744A (zh) | 一种用于目标跟踪的强跟踪容积卡尔曼滤波方法 | |
CN103940433B (zh) | 一种基于改进的自适应平方根ukf算法的卫星姿态确定方法 | |
CN105954743B (zh) | 一种改进权值的变结构多模型机动目标跟踪方法 | |
CN110208792A (zh) | 同时估计目标状态和轨迹参数的任意直线约束跟踪方法 | |
CN108896986A (zh) | 一种基于预测值的量测转换序贯滤波机动目标跟踪方法 | |
CN103973263B (zh) | 一种逼近滤波方法 | |
CN105445718A (zh) | 一种基于阵列重构的分布式多载舰超视距雷达的doa估计方法 | |
CN107290742A (zh) | 一种非线性目标跟踪***中平方根容积卡尔曼滤波方法 | |
CN105372653B (zh) | 面向岸基空管雷达***中一种高效转弯机动目标跟踪方法 | |
CN117272525A (zh) | 一种智能电动汽车路面附着系数估计方法 | |
CN104021285B (zh) | 一种具有最优运动模式切换参数的交互式多模型目标跟踪方法 | |
CN107621632A (zh) | 用于nshv跟踪滤波的自适应滤波方法及*** | |
CN112446010B (zh) | 自适应弱敏秩卡尔曼滤波方法及其应用 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150930 Termination date: 20200908 |
|
CF01 | Termination of patent right due to non-payment of annual fee |