CN105093198A - 一种分布式外辐射源雷达组网探测的航迹融合方法 - Google Patents
一种分布式外辐射源雷达组网探测的航迹融合方法 Download PDFInfo
- Publication number
- CN105093198A CN105093198A CN201510451060.XA CN201510451060A CN105093198A CN 105093198 A CN105093198 A CN 105093198A CN 201510451060 A CN201510451060 A CN 201510451060A CN 105093198 A CN105093198 A CN 105093198A
- Authority
- CN
- China
- Prior art keywords
- current time
- model
- constant velocity
- moment
- fast track
- 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.)
- Granted
Links
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种分布式外辐射源雷达组网探测的航迹融合方法,包括以下步骤:(1)外辐射源雷达各分布站点分别探测当前时刻的目标信息,并解析对应的当前时刻的目标状态,外辐射源雷达各分布站点获得其对应的当前时刻的局部航迹;(2)融合中心对所有来自于同一目标的外辐射源雷达分布站点的各个当前时刻的局部航迹进行加权融合,得到当前时刻的加权融合后的航迹;(3)在融合中心设置IMM滤波器,对当前时刻的加权融合后的航迹进行IMM滤波,计算得到当前时刻IMM滤波的交互输出状态向量和交互输出协方差矩阵;(4)如果外辐射源雷达各分布站点继续跟踪目标,则令当前时刻增加1,返回步骤1;否则,IMM滤波结束。
Description
技术领域
本发明属于雷达信号处理技术领域,涉及一种分布式外辐射源雷达组网探测的航迹融合方法,进一步涉及一种基于交互式多模型(IMM)算法的分布式外辐射源雷达组网探测的航迹融合方法,用于在分布式外辐射源雷达组网探测中对目标进行航迹融合。
背景技术
外辐射源雷达指的是一种本身不发射电磁波,而利用非合作照射源对目标进行探测、定位和跟踪的新体制雷达,因为其发射源的信号波形和信号功率不受控,且辐射源的空间分布和覆盖范围受限,这就导致单站外辐射源雷达的观测信息中包含较多的杂波和虚警,所探测的目标信息的置信度低。伴随着近年来通信技术的发展,各雷达站点间实时通信成为了可能,由此推进了利用多站外辐射源雷达的观测信息,通过时间和空间的信息融合,来改善多目标定位和跟踪性能的分布式外辐射源雷达组网探测的发展。目前各种体制雷达所采用的多站组网探测主要包括三种结构:第一种是雷达各分布站点分别探测目标信息(距离和多普勒等),但各站点不具备解析目标状态(目标位置和速度等)的能力,只是把探测到的目标信息传送到融合中心,由融合中心来解析目标状态,这种架构称为集中式雷达组网探测;第二种是雷达各分布站点不但能探测目标信息,同时各站点能够根据各自的探测信息解析目标状态,然后将解析的目标状态传送到融合中心,在融合中心进行目标状态的融合,这种结构称为分布式雷达组网探测;第三种是雷达各分布站点中只有部分站点同时具备目标信息探测和目标状态解析的能力,其余站点只能探测目标信息,具备目标状态解析能力的站点将目标状态传送到融合中心,只具备探测能力的站点将探测到的目标信息传送到融合中心,在融合中心站点完成目标信息和目标状态的融合,这种结构称为混合式雷达组网探测。在上述三种结构中,分布式雷达组网探测相比于集中式雷达组网探测具有传送信息量少,通信带宽窄的优势;同时,分布式雷达组网探测相比于混合式雷达组网探测具有结构简单,造价低廉的优势,因此,分布式雷达组网探测在外辐射源雷达组网探测中得到了较为广泛的应用。
分布式外辐射源雷达组网探测在融合中心进行目标状态融合时,较常见的是进行目标航迹的融合,即对各站点针对同一目标局部解析形成的航迹进行融合。目前实践中采用的航迹融合算法通常是按照一定的融合准则对各站点传送过来的目标状态进行加权,将加权的结果作为最终的目标航迹融合的结果。这样的加权融合,由于未考虑各站点对目标进行局部解析时因运动模型单一引入的模型误差,从而导致该模型误差仍然保留在目标航迹融合的结果中,造成了航迹融合精度的下降。
发明内容
针对上述已有技术的不足,本发明的目的在于提出一种分布式外辐射源雷达组网探测的航迹融合方法,该方法提出基于IMM算法对加权融合后的航迹进行滤波,即通过多模型交互转换更加精确地得到目标的运动模型;本发明方法能够克服分布式外辐射源雷达组网探测中融合中心未考虑分布站点的模型误差而造成的融合航迹精度差的缺点,因而能够提高对目标航迹融合的精度。
为达到上述技术目的,本发明采用以下技术方案予以实现。
一种分布式外辐射源雷达组网探测的航迹融合方法,其特征在于,包括以下步骤:
步骤1,外辐射源雷达各分布站点分别探测当前时刻的目标信息,并分别解析得到各分布站点对应的当前时刻的目标状态,外辐射源雷达各分布站点获得其对应的当前时刻的局部航迹,该局部航迹包含状态向量和协方差矩阵;外辐射源雷达各分布站点将获得的各个当前时刻的局部航迹传送到融合中心;
步骤2,融合中心根据外辐射源雷达每个分布站点传送过来的对应的当前时刻的局部航迹中所包含的状态向量和协方差矩阵,判断外辐射源雷达各分布站点的各个当前时刻的局部航迹是否来自于同一目标,对所有来自于同一目标的外辐射源雷达分布站点的当前时刻的局部航迹进行加权融合,得到当前时刻的加权融合后的航迹,该当前时刻的加权融合后的航迹包含当前时刻的融合状态向量和当前时刻的融合协方差矩阵;
步骤3,在融合中心设置IMM滤波器,对步骤2所述当前时刻的加权融合后的航迹进行IMM滤波,计算得到当前时刻IMM滤波的交互输出状态向量和交互输出协方差矩阵;
步骤4,如果外辐射源雷达各分布站点继续跟踪目标,则令当前时刻增加1,返回步骤1;否则,IMM滤波结束,保存步骤3所得当前时刻IMM滤波的交互输出状态向量和交互输出协方差矩阵。
本发明与现有技术相比具有以下优点:
第一,本发明在融合中心基于IMM算法对加权融合后的航迹进行滤波,在滤波过程中将目标运动模型设定为多种运动模型的叠加,从而减小了外辐射源雷达各分布站点在局部滤波时因为使用单一目标运动模型而引入的模型误差,提高了对目标航迹融合的精度;
第二,本发明在融合中心设置IMM滤波器,避免了在外辐射源雷达各分布站点同时使用IMM滤波器,减少了滤波器的个数,从而减少了运算量,并加快了运算速度,节省了存储空间;同时本发明可根据实际需要在融合中心增加或减少目标运动模型的个数,而不需要改变融合中心的IMM滤波器的结构,可移植性强;
第三,本发明方法使得分布式外辐射源雷达组网探测的整体结构更为简单,造价低廉,应用面广。
附图说明
下面结合附图说明和具体实施方式对本发明作进一步详细说明。
图1为本发明的IMM滤波的数据流向图;
图2为两种方法所得的目标航迹沿X方向的航迹精度对比图;
图3为两种方法所得的目标航迹沿X方向的航迹精度对比图;
图4为实测实验的航迹误差图。
具体实施方式
参照图1,本发明的一种分布式外辐射源雷达组网探测的航迹融合方法,包括以下具体步骤:
步骤1,外辐射源雷达各分布站点分别探测当前时刻的目标信息,并分别解析得到各分布站点对应的当前时刻的目标状态,外辐射源雷达各分布站点获得其对应的当前时刻的局部航迹,该局部航迹包含状态向量和协方差矩阵;外辐射源雷达各分布站点将获得的各个当前时刻的局部航迹传送到融合中心。
步骤2,融合中心根据外辐射源雷达每个分布站点传送过来的对应的当前时刻的局部航迹中所包含的状态向量和协方差矩阵,判断外辐射源雷达各分布站点的各个当前时刻的局部航迹是否来自于同一目标,对所有来自于同一目标的外辐射源雷达分布站点的当前时刻的局部航迹进行加权融合,得到当前时刻的加权融合后的航迹,该当前时刻的加权融合后的航迹包含当前时刻的融合状态向量和当前时刻的融合协方差矩阵;
步骤2的具体子步骤为:
2.1按照Singer-Kanyuck算法判断外辐射源雷达各分布站点的各个当前时刻的局部航迹是否来自于同一目标;
2.2对步骤2.1中所有判断为来自于同一目标的外辐射源雷达分布站点的各个当前时刻的局部航迹按照Bar-Shalom-Campo算法进行加权融合,得到当前时刻的加权融合后的航迹,该当前时刻的加权融合后的航迹包括:当前时刻的融合状态向量、当前时刻的融合协方差矩阵;将当前时刻记为k时刻,当前时刻的融合状态向量记为当前时刻的融合协方差矩阵记为
步骤3,在融合中心设置IMM滤波器,对步骤2所述当前时刻的加权融合后的航迹进行IMM滤波,计算得到当前时刻IMM滤波的交互输出状态向量和交互输出协方差矩阵;
步骤3的具体子步骤为:
3.1在融合中心设置IMM滤波器,并设定IMM滤波初始时刻,如果当前时刻是IMM滤波初始时刻,初始化目标运动模型集,该目标运动模型集包含常速度模型和常加速度模型,设定IMM滤波初始时刻目标运动模型为常速度模型的概率μ1(0)=0.5,设定IMM滤波初始时刻目标运动模型为常加速度模型的概率μ2(0)=0.5;
构建常速度模型和常加速度模型的马尔科夫转移概率矩阵p,p的维数为2×2,p中的元素pij表示模型i向模型j转移的概率,i∈{1,2},j∈{1,2},1表示常速度模型,2表示常加速度模型,p中的同模型间转移概率,即p中的元素p11和p22在区间[0.8,1]中选取;p中的不同模型间转移概率,即p中的元素p12和p21在区间[0,0.2]中选取;
根据IMM滤波初始时刻,目标运动模型分别为常速度模型和常加速度模型的概率μ1(0)和μ2(0),以及所构建的常速度模型和常加速度模型的马尔科夫转移概率矩阵p,计算IMM滤波初始时刻的归一化马尔科夫转移概率矩阵的维数为2×2,中的元素表示为i∈{1,2},j∈{1,2},1表示常速度模型,2表示常加速度模型;
将当前时刻记为k时刻;将步骤2.2中所述当前时刻的融合状态向量分别作为IMM滤波初始时刻常速度模型和常加速度模型的卡尔曼滤波输出状态向量和即 将步骤2.2中所述当前时刻的融合协方差矩阵分别作为IMM滤波初始时刻常速度模型和常加速度模型的卡尔曼滤波输出协方差矩阵和即 并令当前时刻增加1,即令k时刻增加1,返回步骤1;
如果当前时刻不是IMM滤波初始时刻,则进行下一步操作;
3.2分别计算当前时刻常速度模型和常加速度模型的交互输入条件:即根据上一时刻常速度模型和常加速度模型的卡尔曼滤波输出状态向量、上一时刻常速度模型和常加速度模型的卡尔曼滤波输出协方差矩阵、上一时刻目标运动模型分别为常速度模型和常加速度模型的概率以及上一时刻的归一化马尔科夫转移概率矩阵,分别计算当前时刻常速度模型和常加速度模型的交互输入状态向量,以及当前时刻常速度模型和常加速度模型的交互输入协方差矩阵;
具体地,将当前时刻记为k时刻;首先,计算k-1时刻模型j的归一化常数其中,j∈{1,2},1表示常速度模型,2表示常加速度模型,的表达式为:
其中,pij是常速度模型和常加速度模型的马尔科夫转移概率矩阵p的元素,pij表示模型i向模型j转移的概率,μi(k-1)表示k-1时刻目标运动模型为模型i的概率;i∈{1,2},1表示常速度模型,2表示常加速度模型;
然后,根据k-1时刻模型j的归一化常数计算k-1时刻的归一化马尔科夫转移概率矩阵中的元素的表达式为:
最后,根据k-1时刻的归一化马尔科夫转移概率矩阵计算k时刻模型j的交互输入状态向量以及k时刻模型j的的交互输入协方差矩阵其中,j∈{1,2},1表示常速度模型,2表示常加速度模型,和的表达式分别为:
其中,表示k-1时刻模型i的卡尔曼滤波输出状态向量,表示k-1时刻模型i的卡尔曼滤波输出协方差矩阵,i∈{1,2},1表示常速度模型,2表示常加速度模型,上标T表示矩阵转置;
3.3将步骤2.2中所述当前时刻的融合状态向量分别作为当前时刻常速度模型和常加速度模型的输入观测向量x1(k|k)和x2(k|k),即 将步骤2.2中所述当前时刻的融合协方差矩阵分别作为当前时刻常速度模型和常加速度模型的输入观测方差矩阵p1(k|k)和p2(k|k),即 结合步骤3.2所述当前时刻常速度模型和常加速度模型的交互输入条件,对常速度模型和常加速度模型分别进行卡尔曼滤波,即分别计算得到当前时刻常速度模型和常加速度模型的卡尔曼滤波输出状态向量和以及当前时刻常速度模型和常加速度模型的卡尔曼滤波输出协方差矩阵和
3.4根据当前时刻常速度模型和常加速度模型的输入观测向量x1(k|k)和x2(k|k)、当前时刻常速度模型和常加速度模型的输入观测方差矩阵p1(k|k)和p2(k|k)、步骤3.1所述常速度模型和常加速度模型的马尔科夫转移概率矩阵p,以及步骤3.2所述当前时刻常速度模型和常加速度模型的交互输入条件,分别计算当前时刻常速度模型和常加速度模型的似然函数Λ1(k)和Λ2(k);再根据当前时刻常速度模型和常加速度模型的似然函数Λ1(k)和Λ2(k),计算当前时刻目标运动模型分别为常速度模型和常加速度模型的概率μ1(k)和μ2(k);
具体地,将当前时刻记为k时刻;首先,计算k时刻模型j的似然函数Λj(k),其中,j∈{1,2},1表示常速度模型,2表示常加速度模型,Λj(k)的表达式为:
其中,Mj(k)表示k时刻模型j有效,xj(k|k)表示k时刻模型j的输入观测向量,表示k时刻模型j的交互输入状态向量,表示k时刻模型j的交互输入协方差矩阵;
然后,根据k时刻模型j的似然函数Λj(k),计算k时刻似然函数的归一化常数c(k),c(k)的表达式为:
其中,表示k-1时刻模型j的归一化常数;
最后,根据k时刻似然函数的归一化常数c(k),计算k时刻目标运动模型为模型j的概率μj(k),μj(k)的表达式为:
3.5根据步骤3.3所述当前时刻常速度模型和常加速度模型的卡尔曼滤波输出状态向量和以及当前时刻常速度模型和常加速度模型的卡尔曼滤波输出协方差矩阵和并结合步骤3.4所述当前时刻目标运动模型分别为常速度模型和常加速度模型的概率μ1(k)和μ2(k),加权求得当前时刻IMM滤波的交互输出结果,即加权求得当前时刻IMM滤波的交互输出状态向量和交互输出协方差矩阵
具体地,将当前时刻记为k时刻;首先,计算k时刻IMM滤波的交互输出状态向量的表达式为:
其中,表示k时刻模型j的卡尔曼滤波输出状态向量,μj(k)表示k时刻目标运动模型为模型j的概率;
然后,计算k时刻IMM滤波的交互输出状态向量的表达式为:
其中,表示k时刻模型j的卡尔曼滤波输出协方差矩阵。
步骤4,如果外辐射源雷达各分布站点继续跟踪目标,则令当前时刻增加1,返回步骤1;否则,IMM滤波结束,保存步骤3所得当前时刻IMM滤波的交互输出状态向量和交互输出协方差矩阵。
本发明的效果可由以下仿真实验和实测实验作进一步说明:
1)仿真实验参数
仿真实验的数据是在总采样帧数为200帧的采样时间内由多种运动模型叠加产生,采样周期T=1s,目标初始位置在原点,初始速度为(200m/s,200m/s),在1-50帧的采样时间内采用的运动模型是匀速直线运动叠加标准方差为0.1m/s2的高斯白噪声驱动,在51-150帧的采样时间内采用的运动模型是匀速圆周运动叠加标准方差为0.1m/s2的高斯白噪声驱动,在151-200帧的采样时间内采用的运动模型是加速度为(5m/s2,5m/s2)的匀加速直线运动叠加标准方差为0.1m/s2的高斯白噪声驱动,外辐射源雷达各分布站点的观测噪声是标准方差为0.1m/s2的高斯白噪声。
2)仿真实验内容及结果分析:
根据上述仿真实验参数,进行500次蒙特卡洛实验,在单个外辐射源雷达站点运用卡尔曼滤波计算得到目标航迹,简称卡尔曼滤波方法;同时,在多个外辐射源雷达站点运用本发明方法计算得到目标航迹,简称本发明方法;将上述两种方法所得的目标航迹分别与未叠加高斯白噪声驱动的目标的真实航迹进行比较,绘制两种方法所得的目标航迹沿X方向的航迹精度对比图,如图2所示,图2中,横坐标为采样时刻,单位为帧,纵坐标为X方向误差,单位为米(m),虚线表示在单个外辐射源雷达站点运用卡尔曼滤波计算得到的目标航迹与目标真实航迹沿X方向的误差,实线表示在多个外辐射源雷达站点运用本发明方法计算得到的目标航迹与目标真实航迹沿X方向的误差;绘制两种方法所得的目标航迹沿Y方向的航迹精度对比图,如图3所示,图3中,横坐标为采样时刻,单位为帧,纵坐标为Y方向误差,单位为米(m),虚线表示在单个外辐射源雷达站点运用卡尔曼滤波计算得到的目标航迹与目标真实航迹沿Y方向的误差,实线表示在多个外辐射源雷达站点运用本发明方法计算得到的目标航迹与目标真实航迹沿Y方向的误差;各个采样时刻的上述误差取500次蒙特卡洛实验结果的几何平均值。
从图2和图3中可以看出:在匀速直线运动的采样时间内,两种方法所得的目标航迹沿X方向和沿Y方向的误差比较一致,都能比较准确的跟踪目标;在匀速圆周运动的采样时间内,由于匀速圆周运动在X方向和Y方向的速度和加速度都不是均匀变化的,所以两种方法所得的目标航迹沿X方向和Y方向的航迹精度都开始变差,但本发明方法相比较于卡尔曼滤波方法可以得到更精确的目标航迹;在匀加速直线运动的采样时间内,本发明方法的优势更能体现出来,这是由于相比较于卡尔曼滤波方法只考虑单一目标运动模型,本发明方法将目标运动模型考虑为多种运动模型的叠加,再通过IMM滤波获得多模型交互输出的结果,减小了卡尔曼滤波方法存在的模型匹配误差,提高了对目标的跟踪精度。
3)实测实验参数:
实测实验数据在湖北省荆州市江北农场采集,两个外辐射源雷达站点分别位于东经112.4°、北纬30.1252°和东经112.414°、北纬30.1251°,两个外辐射源雷达站点相距1.3465公里,所探测的目标为正在飞行的民航飞机。两个外辐射源雷达站点的接收机具有同样的软件、硬件架构,采集速率均为1帧/秒,总采样帧数为200帧。在实测实验中,将融合中心已关联航迹设定为所选航迹,并设定所选航迹不存在错误关联和数据采集错误的情况,与所选目标航迹对应的目标的真实航迹,称为标准航迹,该标准航迹包含所探测的目标(即正在飞行的民航飞机)的真实状态信息,该标准航迹由民航ADS-B接收机获得;
4)实测实验内容及结果分析:
在实测实验中,首先,在所述两个外辐射源雷达站点中的一个外辐射源雷达站点运用卡尔曼滤波得到目标航迹;然后,在所述两个外辐射源雷达站点运用本发明方法得到目标航迹。图4为实测实验的航迹误差图,横轴是采样时刻,单位为帧,纵轴是距离误差,单位是千米(Km),实线表示在单个外辐射源雷达站点运动卡尔曼滤波所得到的目标航迹与标准航迹的距离误差,虚线表示在两个外辐射源雷达站点运用本发明方法所得到的目标航迹与标准航迹的距离误差。
从图4中可以看出,相比较于在一个外辐射源雷达站点运用卡尔曼滤波所得到的目标航迹,在两个外辐射源雷达站点运用本发明方法所得到的目标航迹具有更小的距离误差,所以本发明方法在实际应用中具有较好的优越性和实用性。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (3)
1.一种分布式外辐射源雷达组网探测的航迹融合方法,其特征在于,包括以下步骤:
步骤1,外辐射源雷达各分布站点分别探测当前时刻的目标信息,并分别解析得到各分布站点对应的当前时刻的目标状态,外辐射源雷达各分布站点获得其对应的当前时刻的局部航迹,该局部航迹包含状态向量和协方差矩阵;外辐射源雷达各分布站点将获得的各个当前时刻的局部航迹传送到融合中心;
步骤2,融合中心根据外辐射源雷达每个分布站点传送过来的对应的当前时刻的局部航迹中所包含的状态向量和协方差矩阵,判断外辐射源雷达各分布站点的各个当前时刻的局部航迹是否来自于同一目标,对所有来自于同一目标的外辐射源雷达分布站点的当前时刻的局部航迹进行加权融合,得到当前时刻的加权融合后的航迹,该当前时刻的加权融合后的航迹包含当前时刻的融合状态向量和当前时刻的融合协方差矩阵;
步骤3,在融合中心设置IMM滤波器,对步骤2所述当前时刻的加权融合后的航迹进行IMM滤波,计算得到当前时刻IMM滤波的交互输出状态向量和交互输出协方差矩阵;
步骤4,如果外辐射源雷达各分布站点继续跟踪目标,则令当前时刻增加1,返回步骤1;否则,IMM滤波结束,保存步骤3所得当前时刻IMM滤波的交互输出状态向量和交互输出协方差矩阵。
2.如权利要求1所述的一种分布式外辐射源雷达组网探测的航迹融合方法,其特征在于,所述步骤2的具体子步骤为:
2.1按照Singer-Kanyuck算法判断外辐射源雷达各分布站点的各个当前时刻的局部航迹是否来自于同一目标;
2.2对步骤2.1中所有判断为来自于同一目标的外辐射源雷达分布站点的各个当前时刻的局部航迹按照Bar-Shalom-Campo算法进行加权融合,得到当前时刻的加权融合后的航迹,该当前时刻的加权融合后的航迹包括:当前时刻的融合状态向量、当前时刻的融合协方差矩阵;将当前时刻记为k时刻,当前时刻的融合状态向量记为当前时刻的融合协方差矩阵记为
3.如权利要求1所述的一种分布式外辐射源雷达组网探测的航迹融合方法,其特征在于,所述步骤3的具体子步骤为:
3.1在融合中心设置IMM滤波器,并设定IMM滤波初始时刻,如果当前时刻是IMM滤波初始时刻,初始化目标运动模型集,该目标运动模型集包含常速度模型和常加速度模型,设定IMM滤波初始时刻目标运动模型为常速度模型的概率μ1(0)=0.5,设定IMM滤波初始时刻目标运动模型为常加速度模型的概率μ2(0)=0.5;
构建常速度模型和常加速度模型的马尔科夫转移概率矩阵p,p的维数为2×2,p中的元素pij表示模型i向模型j转移的概率,i∈{1,2},j∈{1,2},1表示常速度模型,2表示常加速度模型,p中的同模型间转移概率,即p中的元素p11和p22在区间[0.8,1]中选取;p中的不同模型间转移概率,即p中的元素p12和p21在区间[0,0.2]中选取;
根据IMM滤波初始时刻,目标运动模型分别为常速度模型和常加速度模型的概率μ1(0)和μ2(0),以及所构建的常速度模型和常加速度模型的马尔科夫转移概率矩阵p,计算IMM滤波初始时刻的归一化马尔科夫转移概率矩阵的维数为2×2,中的元素表示为i∈{1,2},j∈{1,2},1表示常速度模型,2表示常加速度模型;
将当前时刻记为k时刻;将步骤2.2中所述当前时刻的融合状态向量分别作为IMM滤波初始时刻常速度模型和常加速度模型的卡尔曼滤波输出状态向量和即 将步骤2.2中所述当前时刻的融合协方差矩阵分别作为IMM滤波初始时刻常速度模型和常加速度模型的卡尔曼滤波输出协方差矩阵和即 并令当前时刻增加1,即令k时刻增加1,返回步骤1;
如果当前时刻不是IMM滤波初始时刻,则进行下一步操作;
3.2分别计算当前时刻常速度模型和常加速度模型的交互输入条件:即根据上一时刻常速度模型和常加速度模型的卡尔曼滤波输出状态向量、上一时刻常速度模型和常加速度模型的卡尔曼滤波输出协方差矩阵、上一时刻目标运动模型分别为常速度模型和常加速度模型的概率以及上一时刻的归一化马尔科夫转移概率矩阵,分别计算当前时刻常速度模型和常加速度模型的交互输入状态向量,以及当前时刻常速度模型和常加速度模型的交互输入协方差矩阵;
具体地,将当前时刻记为k时刻;首先,计算k-1时刻模型j的归一化常数其中,j∈{1,2},1表示常速度模型,2表示常加速度模型,的表达式为:
其中,pij是常速度模型和常加速度模型的马尔科夫转移概率矩阵p的元素,pij表示模型i向模型j转移的概率,μi(k-1)表示k-1时刻目标运动模型为模型i的概率;i∈{1,2},1表示常速度模型,2表示常加速度模型;
然后,根据k-1时刻模型j的归一化常数计算k-1时刻的归一化马尔科夫转移概率矩阵中的元素的表达式为:
最后,根据k-1时刻的归一化马尔科夫转移概率矩阵计算k时刻模型j的交互输入状态向量以及k时刻模型j的的交互输入协方差矩阵其中,j∈{1,2},1表示常速度模型,2表示常加速度模型,和的表达式分别为:
其中,表示k-1时刻模型i的卡尔曼滤波输出状态向量,表示k-1时刻模型i的卡尔曼滤波输出协方差矩阵,i∈{1,2},1表示常速度模型,2表示常加速度模型,上标T表示矩阵转置;
3.3将步骤2.2中所述当前时刻的融合状态向量分别作为当前时刻常速度模型和常加速度模型的输入观测向量x1(k|k)和x2(k|k),即 将步骤2.2中所述当前时刻的融合协方差矩阵分别作为当前时刻常速度模型和常加速度模型的输入观测方差矩阵p1(k|k)和p2(k|k),即 结合步骤3.2所述当前时刻常速度模型和常加速度模型的交互输入条件,对常速度模型和常加速度模型分别进行卡尔曼滤波,即分别计算得到当前时刻常速度模型和常加速度模型的卡尔曼滤波输出状态向量和以及当前时刻常速度模型和常加速度模型的卡尔曼滤波输出协方差矩阵和
3.4根据当前时刻常速度模型和常加速度模型的输入观测向量x1(k|k)和x2(k|k)、当前时刻常速度模型和常加速度模型的输入观测方差矩阵p1(k|k)和p2(k|k)、步骤3.1所述常速度模型和常加速度模型的马尔科夫转移概率矩阵p,以及步骤3.2所述当前时刻常速度模型和常加速度模型的交互输入条件,分别计算当前时刻常速度模型和常加速度模型的似然函数Λ1(k)和Λ2(k);再根据当前时刻常速度模型和常加速度模型的似然函数Λ1(k)和Λ2(k),计算当前时刻目标运动模型分别为常速度模型和常加速度模型的概率μ1(k)和μ2(k);
具体地,将当前时刻记为k时刻;首先,计算k时刻模型j的似然函数Λj(k),其中,j∈{1,2},1表示常速度模型,2表示常加速度模型,Λj(k)的表达式为:
其中,Mj(k)表示k时刻模型j有效,xj(k|k)表示k时刻模型j的输入观测向量,表示k时刻模型j的交互输入状态向量,表示k时刻模型j的交互输入协方差矩阵;
然后,根据k时刻模型j的似然函数Λj(k),计算k时刻似然函数的归一化常数c(k),c(k)的表达式为:
其中,表示k-1时刻模型j的归一化常数;
最后,根据k时刻似然函数的归一化常数c(k),计算k时刻目标运动模型为模型j的概率μj(k),μj(k)的表达式为:
3.5根据步骤3.3所述当前时刻常速度模型和常加速度模型的卡尔曼滤波输出状态向量和以及当前时刻常速度模型和常加速度模型的卡尔曼滤波输出协方差矩阵和并结合步骤3.4所述当前时刻目标运动模型分别为常速度模型和常加速度模型的概率μ1(k)和μ2(k),加权求得当前时刻IMM滤波的交互输出结果,即加权求得当前时刻IMM滤波的交互输出状态向量和交互输出协方差矩阵
具体地,将当前时刻记为k时刻;首先,计算k时刻IMM滤波的交互输出状态向量的表达式为:
其中,表示k时刻模型j的卡尔曼滤波输出状态向量,μj(k)表示k时刻目标运动模型为模型j的概率;
然后,计算k时刻IMM滤波的交互输出状态向量的表达式为:
其中,表示k时刻模型j的卡尔曼滤波输出协方差矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510451060.XA CN105093198B (zh) | 2015-07-28 | 2015-07-28 | 一种分布式外辐射源雷达组网探测的航迹融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510451060.XA CN105093198B (zh) | 2015-07-28 | 2015-07-28 | 一种分布式外辐射源雷达组网探测的航迹融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105093198A true CN105093198A (zh) | 2015-11-25 |
CN105093198B CN105093198B (zh) | 2017-09-01 |
Family
ID=54574062
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510451060.XA Active CN105093198B (zh) | 2015-07-28 | 2015-07-28 | 一种分布式外辐射源雷达组网探测的航迹融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105093198B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105825239A (zh) * | 2016-03-31 | 2016-08-03 | 中国科学院电子学研究所 | 基于稀疏表示的多传感器航迹融合方法 |
CN105866757A (zh) * | 2016-03-30 | 2016-08-17 | 西安电子科技大学 | 基于时变定位误差的外辐射源雷达目标航迹融合方法 |
CN106990399A (zh) * | 2017-05-11 | 2017-07-28 | 西安电子科技大学 | 用于目标跟踪的组网雷达系统功率和带宽联合分配方法 |
CN109490864A (zh) * | 2018-12-18 | 2019-03-19 | 安徽四创电子股份有限公司 | 一种靶场三坐标测量雷达组网测试系统 |
CN109597104A (zh) * | 2018-10-22 | 2019-04-09 | 中国科学院光电研究院 | 一种采用多模型多系统对目标跟踪定位的方法 |
CN109633628A (zh) * | 2019-01-07 | 2019-04-16 | 西安电子科技大学 | 基于分布式组网雷达数据融合的抗rgpo干扰的方法 |
CN109752002A (zh) * | 2018-12-26 | 2019-05-14 | 北京壹氢科技有限公司 | 一种多平台无源目标跟踪方法 |
CN110542898A (zh) * | 2019-06-18 | 2019-12-06 | 同济大学 | 一种基于雷达组群的车辆行为连续跟踪探测系统及方法 |
CN111337918A (zh) * | 2020-02-17 | 2020-06-26 | 南京航空航天大学 | 一种基于神经网络的机载雷达射频隐身波形选择方法 |
CN114971436A (zh) * | 2022-08-02 | 2022-08-30 | 中国电子科技集团公司第十五研究所 | 雷达情报组网融合中的装备贡献计算方法及装置 |
CN114971283A (zh) * | 2022-05-25 | 2022-08-30 | 中国人民解放军国防科技大学 | 一种面向分布式组网雷达多目标跟踪的资源优化调度方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102622520A (zh) * | 2012-03-14 | 2012-08-01 | 北京航空航天大学 | 一种机动目标跟踪的分布式多模型估计融合方法 |
CN103743401A (zh) * | 2013-12-31 | 2014-04-23 | 电子科技大学 | 基于多模型航迹质量的异步融合方法 |
CN103853908A (zh) * | 2012-12-04 | 2014-06-11 | 中国科学院沈阳自动化研究所 | 一种自适应交互式多模型的机动目标跟踪方法 |
-
2015
- 2015-07-28 CN CN201510451060.XA patent/CN105093198B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102622520A (zh) * | 2012-03-14 | 2012-08-01 | 北京航空航天大学 | 一种机动目标跟踪的分布式多模型估计融合方法 |
CN103853908A (zh) * | 2012-12-04 | 2014-06-11 | 中国科学院沈阳自动化研究所 | 一种自适应交互式多模型的机动目标跟踪方法 |
CN103743401A (zh) * | 2013-12-31 | 2014-04-23 | 电子科技大学 | 基于多模型航迹质量的异步融合方法 |
Non-Patent Citations (2)
Title |
---|
程跃兵: "多传感器多目标航迹关联与融合算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
董亮玉: "外辐射源雷达中的定位跟踪算法和信息融合", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105866757A (zh) * | 2016-03-30 | 2016-08-17 | 西安电子科技大学 | 基于时变定位误差的外辐射源雷达目标航迹融合方法 |
CN105866757B (zh) * | 2016-03-30 | 2018-07-20 | 西安电子科技大学 | 基于时变定位误差的外辐射源雷达目标航迹融合方法 |
CN105825239A (zh) * | 2016-03-31 | 2016-08-03 | 中国科学院电子学研究所 | 基于稀疏表示的多传感器航迹融合方法 |
CN106990399A (zh) * | 2017-05-11 | 2017-07-28 | 西安电子科技大学 | 用于目标跟踪的组网雷达系统功率和带宽联合分配方法 |
CN106990399B (zh) * | 2017-05-11 | 2019-12-20 | 西安电子科技大学 | 用于目标跟踪的组网雷达系统功率和带宽联合分配方法 |
CN109597104A (zh) * | 2018-10-22 | 2019-04-09 | 中国科学院光电研究院 | 一种采用多模型多系统对目标跟踪定位的方法 |
CN109490864A (zh) * | 2018-12-18 | 2019-03-19 | 安徽四创电子股份有限公司 | 一种靶场三坐标测量雷达组网测试系统 |
CN109752002B (zh) * | 2018-12-26 | 2020-09-04 | 北京壹氢科技有限公司 | 一种多平台无源目标跟踪方法 |
CN109752002A (zh) * | 2018-12-26 | 2019-05-14 | 北京壹氢科技有限公司 | 一种多平台无源目标跟踪方法 |
CN109633628A (zh) * | 2019-01-07 | 2019-04-16 | 西安电子科技大学 | 基于分布式组网雷达数据融合的抗rgpo干扰的方法 |
CN109633628B (zh) * | 2019-01-07 | 2023-04-21 | 西安电子科技大学 | 基于分布式组网雷达数据融合的抗rgpo干扰的方法 |
CN110542898A (zh) * | 2019-06-18 | 2019-12-06 | 同济大学 | 一种基于雷达组群的车辆行为连续跟踪探测系统及方法 |
CN111337918A (zh) * | 2020-02-17 | 2020-06-26 | 南京航空航天大学 | 一种基于神经网络的机载雷达射频隐身波形选择方法 |
CN114971283A (zh) * | 2022-05-25 | 2022-08-30 | 中国人民解放军国防科技大学 | 一种面向分布式组网雷达多目标跟踪的资源优化调度方法 |
CN114971436A (zh) * | 2022-08-02 | 2022-08-30 | 中国电子科技集团公司第十五研究所 | 雷达情报组网融合中的装备贡献计算方法及装置 |
CN114971436B (zh) * | 2022-08-02 | 2022-09-30 | 中国电子科技集团公司第十五研究所 | 雷达情报组网融合中的装备贡献计算方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105093198B (zh) | 2017-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105093198A (zh) | 一种分布式外辐射源雷达组网探测的航迹融合方法 | |
CN102622520B (zh) | 一种机动目标跟踪的分布式多模型估计融合方法 | |
CN103472445B (zh) | 一种针对多目标场景的检测跟踪一体化方法 | |
CN101770024B (zh) | 多目标跟踪方法 | |
KR101628154B1 (ko) | 수신 신호 세기를 이용한 다중 표적 추적 방법 | |
CN105719312A (zh) | 基于序贯贝叶斯滤波的多目标跟踪方法及跟踪系统 | |
CN101975575A (zh) | 基于粒子滤波的被动传感器多目标跟踪方法 | |
CN103298156B (zh) | 基于无线传感器网络的无源多目标检测跟踪方法 | |
CN103152745B (zh) | 一种强自适应性移动节点定位的方法 | |
CN105549005A (zh) | 一种基于网格划分的动态目标波达方向跟踪方法 | |
CN111601253B (zh) | 无源被动智能跟踪定位方法、系统、存储介质、跟踪定位终端 | |
CN108089183A (zh) | 一种用于异步多基地雷达系统的检测跟踪一体化方法 | |
CN104535960B (zh) | 一种基于rfid的室内快速定位方法 | |
CN103152826A (zh) | 一种基于nlos状态检测补偿的移动目标追踪方法 | |
CN104199022A (zh) | 一种基于目标模态估计的临近空间高超声速目标跟踪方法 | |
CN104730510A (zh) | 一种多雷达航迹融合方法 | |
CN105954741A (zh) | 一种基于多假设拟蒙特卡罗的多目标无源协同定位方法 | |
CN105301584A (zh) | 同时解距离模糊的ipphdf机动多目标跟踪方法 | |
Liang et al. | Multi-Bernoulli filter for target tracking with multi-static Doppler only measurement | |
CN103247057A (zh) | 目标-回波-道路网三元数据关联的道路目标多假设跟踪算法 | |
CN105528796A (zh) | 一种穿墙探测雷达在多运动目标场景下对虚假目标实时剔除的方法 | |
Ronghua et al. | TDOA and track optimization of UAV swarm based on D-optimality | |
Lu et al. | Robot indoor location modeling and simulation based on Kalman filtering | |
CN103476110A (zh) | 节点自定位与目标跟踪同时进行的分布式算法 | |
CN105866757A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |