CN104849702B - 利用ads‑b数据的gm‑ephd滤波雷达***误差联合估计方法 - Google Patents
利用ads‑b数据的gm‑ephd滤波雷达***误差联合估计方法 Download PDFInfo
- Publication number
- CN104849702B CN104849702B CN201510218847.1A CN201510218847A CN104849702B CN 104849702 B CN104849702 B CN 104849702B CN 201510218847 A CN201510218847 A CN 201510218847A CN 104849702 B CN104849702 B CN 104849702B
- Authority
- CN
- China
- Prior art keywords
- radar
- radar system
- system error
- ads
- target
- 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
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/40—Means for monitoring or calibrating
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
一种利用ADS‑B数据的GM‑EPHD滤波雷达***误差联合估计方法。其首先将目标ADS‑B观测值转换到以雷达站位置为中心的直角坐标系下,建立目标的ADS‑B及雷达的观测方程,在线性高斯条件下建立目标及雷达***误差的状态转移方程,将雷达***误差扩展至目标状态中组成扩展状态,先使用坐标转换后的ADS‑B观测数据对目标状态进行高斯混合‑概率假设密度滤波,以获得较为准确的目标状态估计,再使用雷达观测数据对目标状态更新后的扩展状态进行高斯混合‑扩展概率假设密度滤波,进而采用两步卡尔曼滤波器联合估计出目标状态和雷达***误差,最后使用加权平均算法得到雷达***误差的融合估计结果。本方法估计精度高、估计性能好。
Description
技术领域
本发明属于传感器误差配准技术领域,特别是涉及一种利用ADS-B数据的GM-EPHD滤波雷达***误差联合估计方法。
背景技术
雷达误差一般分为两类:随机误差和***误差。在多雷达融合跟踪***中,雷达***误差估计已成为多雷达融合处理的先决条件,将会直接影响整个***的工作性能。因此需要对雷达***误差进行估计,以此对雷达测量进行相应的补偿,这个过程也称为误差配准。雷达***误差包括以下几种:雷达站点坐标标定误差、空间坐标转换误差、雷达测量误差。常用的雷达***误差估计方法一般只考虑因测距、测角引起的测量误差。现有的雷达***误差估计方法可以归纳为离线估计方法和在线估计方法。离线估计方法是通过对一段时间的雷达观测数据进行数据拟合,从而估计出雷达的***误差,如最小二乘法及最大似然法等。在线估计方法主要是利用滤波方法实现雷达***误差的递推估计,在线估计方法和离线估计方法相比具有实时误差估计的优点,因此得到国内外学者的更多关注。2006年Bar-Shalom提出利用Kalman滤波进行雷达***误差估计的方法。2007年Herrero提出一种利用Kalman滤波实现目标状态与***误差联合估计方法。
无论是基于数据拟合的离线估计方法还是基于滤波的在线估计方法,上述方法都需要满足这样一个假设:目标状态和观测之间的关联关系预先已知。现有的方法通过最近邻方法(NN)、联合概率数据关联方法(JPDA)等获得目标状态和观测之间的关联关系。而对于多目标或密集杂波场景,想要获得准确的关联关系是十分困难的,而且错误的目标状态与观测之间的关联关系将会严重影响雷达***误差估计结果。2003年Mahler在随机有限集理论框架下提出传递目标状态集合的后验概率密度的一阶统计矩的概率假设密度(Probability Hypothesis Density,PHD)滤波理论。PHD滤波器将复杂的多目标状态空间的运算转换为单目标状态空间内的运算,可有效避免多目标状态估计中复杂的数据关联过程。2006年Vo提出了高斯混合模型PHD(Gaussian Mixture PHD,GM-PHD)滤波器,给出了线性高斯条件的PHD滤波器的封闭解形式。2012年Wenling Li等将GM-PHD滤波器用于多传感器多目标跟踪中,实现了目标状态和传感器***误差的联合估计。
随着全球导航卫星***(GNSS)和空空、空地数据链通信技术的发展,一种新型的航空器运行监视技术—广播式自动相关监视技术(Automatic Dependent Surveillance-Broadcast,ADS-B)正在航空器监视中广泛应用。航空器上的ADS-B机载收发信机将本机GPS(Global Positioning System)导航设备获得的本机经度、纬度、速度、时间、高度等数据通过数据链对外广播,ADS-B地面站通过接收有效空域内航空器所发送的广播数据实现对航空器的监视,其定位精度即为机载GPS导航设备的定位精度,远远优于雷达定位精度。因此利用ADS-B进行雷达***误差估计成为该领域新的研究热点。2009年Besada提出了利用ADS-B监视数据对空中交通管理***航管雷达进行误差配准的方法。2013年He You提出利用Kalman滤波及ADS-B监视数据进行雷达***误差估计的方法。但到目前为止尚未发现利用ADS-B的概率假设密度滤波雷达***误差联合估计方面的报道。
发明内容
为了解决上述问题,本发明的目的在于提供一种利用ADS-B数据的GM-EPHD滤波雷达***误差联合估计方法。
为了达到上述目的,本发明提供的利用ADS-B数据的GM-EPHD滤波雷达***误差联合估计方法包括按顺序进行的下列步骤:
1)建立目标的ADS-B及雷达的观测方程的S1阶段;
2)建立目标及雷达***误差的状态转移方程的S2阶段;
3)利用GM-EPHD滤波器对雷达***误差进行融合估计的S3阶段。
在步骤1)中,所述的建立目标的ADS-B及雷达的观测方程的方法是首先利用坐标投影技术将目标的ADS-B观测值转换到以雷达站位置为中心的直角坐标系下,由于ADS-B监视数据的定位精度远远优于雷达的定位精度,在不考虑ADS-B定位误差的基础上,建立坐标转换后的目标的ADS-B观测方程,然后考虑雷达***误差,建立目标的雷达观测方程,为利用GM-EPHD滤波器对雷达***误差进行融合估计做准备。
在步骤2)中,所述的建立目标及雷达***误差的状态转移方程的方法是假设目标状态及雷达***误差均满足线性高斯特性,进而建立目标及雷达***误差的状态转移方程,为利用GM-EPHD滤波器对雷达***误差进行融合估计做准备。
在步骤3)中,所述的利用GM-EPHD滤波器对雷达***误差进行融合估计的方法是首先将雷达***误差扩展至目标状态中组成扩展状态,根据S2阶段中建立的目标及雷达***误差的状态转移方程分别对目标状态及雷达***误差进行预测,然后使用S1阶段中坐标转换后的ADS-B观测数据对目标状态进行高斯混合-概率假设密度滤波,以获取较为准确的目标状态估计,接着使用S1阶段中的雷达观测数据对扩展状态进行高斯混合-扩展概率假设密度滤波,进而采用两步卡尔曼滤波器联合估计出目标状态和雷达***误差,最后使用加权平均算法对雷达***误差进行融合估计。
本发明提供的利用ADS-B数据的GM-EPHD滤波雷达***误差联合估计方法首先将目标的ADS-B观测值转换到以雷达站位置为中心的直角坐标系下,由于ADS-B监视数据的定位精度远远优于雷达的定位精度,在不考虑ADS-B定位误差的基础上,建立目标的ADS-B及雷达的观测方程,然后在线性高斯条件下建立目标及雷达***误差的状态转移方程,接着将雷达***误差扩展至目标状态中组成扩展状态,先使用坐标转换后的ADS-B观测数据对目标状态进行高斯混合-概率假设密度滤波,以获得较为准确的目标状态估计,再使用雷达观测数据对目标状态更新后的扩展状态进行高斯混合-扩展概率假设密度滤波,进而采用两步卡尔曼滤波器联合估计出目标状态和雷达***误差,最后使用加权平均算法得到雷达***误差的融合估计结果。本发明方法具有估计精度高、估计性能好等优点。
附图说明
图1为本发明提供的利用ADS-B数据的GM-EPHD滤波雷达***误差联合估计方法流程图。
图2为目标真实运动轨迹及雷达和ADS-B观测数据。
图3为目标真实运动轨迹及状态估计结果。
图4为目标状态估计的OSPA距离。
图5为目标数估计结果。
图6、图7分别为采用本发明方法和文献方法获得的雷达距离、方位误差估计比较图。
图8、图9分别为采用本发明方法和文献方法获得的雷达距离、方位误差估计的均方根误差比较图。
图10为采用本发明方法对雷达***误差校正后的雷达观测图。
具体实施方式
下面结合附图和具体实施例对本发明提供的利用ADS-B数据的GM-EPHD滤波雷达***误差联合估计方法进行详细说明。
图1为本发明提供的利用ADS-B数据的GM-EPHD滤波雷达***误差联合估计方法流程图。其中的全部操作都是在计算机***中完成的,操作的主体均为计算机***。
如图1所示,本发明提供的利用ADS-B数据的GM-EPHD滤波雷达***误差联合估计方法包括按顺序进行的下列步骤:
1)建立目标的ADS-B及雷达的观测方程的S1阶段:
本阶段是首先利用坐标投影技术将目标的ADS-B观测值转换到以雷达站位置为中心的直角坐标系下,由于ADS-B监视数据的定位精度远远优于雷达的定位精度,在不考虑ADS-B定位误差的基础上,建立坐标转换后的目标的ADS-B观测方程,接着考虑雷达***误差,建立目标的雷达观测方程,然后进入下一步S2阶段。
在此阶段中,所述的建立目标的ADS-B及雷达的观测方程的具体方法如下:由于ADS-B使用WGS-84大地坐标系获得目标的经度、纬度及高度信息,而雷达使用以雷达位置为中心的极坐标系获得目标的距离及方位信息。首先将目标的ADS-B观测值转换到以雷达站位置为中心的直角坐标系下,利用ADS-B定位精度即GPS导航设备定位精度远远优于雷达定位精度的优点,在不考虑ADS-B定位误差的基础上,建立坐标转换后的目标的ADS-B观测方程,然后考虑雷达***误差,建立目标的雷达观测方程。
a、建立坐标转换后的目标的ADS-B观测方程
首先将目标的ADS-B观测值转换到地心地固(ECEF)直角坐标系下,然后再转换成以雷达站位置为中心的直角坐标系下。设ADS-B***获得的k时刻第i个目标的观测值为分别表示目标的经度、纬度及距离海平面的高度,转换后的ECEF坐标系观测值为则:
其中,e为地球的偏心率,Eq为赤道半径。
若雷达站所在位置在WGS-84坐标系下可表示为(LRS,λRS,HRS),则其ECEF坐标系下的位置(ExRS,EyRS,EzRS)同样可以由式(1)获得。
将ECEF坐标系下的目标ADS-B观测值通过式(2)转换为以雷达站位置为中心的直角坐标系下的观测值
其中,TRS为旋转矩阵。
由于ADS-B的定位精度远远优于雷达的定位精度,假设不考虑ADS-B定位误差,并且忽略坐标转换误差,则由式(2)表示的以雷达位置为中心的直角坐标系下的ADS-B二维空间测量值可以写成:
将式(3)用矩阵形式表示为:
其中,表示以雷达为中心的直角坐标系下由目标的位置、速度组成的状态矢量,表示随机噪声,且
由于监视区域内的所有目标均配备ADS-B机载设备,则k时刻可以获得Nk个目标的ADS-B观测值,即
b、建立目标的雷达观测方程
雷达使用以雷达位置为中心的极坐标系获得目标的距离和方位信息,雷达的观测模型可以用式(6)表示:
其中,分别表示k时刻第i个目标的雷达距离、方位角测量值,分别表示距离、方位对应的真实值,△r、△θ分别表示雷达距离、方位的***误差, 分别表示雷达测量距离、方位的观测噪声。
将式(6)用矩阵形式表示为:
其中,表示以雷达为中心的直角坐标系下由目标的位置、速度组成的状态矢量,bk=[△r,△θ]T表示由雷达距离、方位组成的***误差,表示雷达距离、方位的观测噪声,且
假设监视区域Nk个目标中有Mk个目标被雷达观测到(由于雷达有一定的检测概率,一般情况下Nk≥Mk),则k时刻可以获得Mk个目标的雷达观测值,即
2)建立目标及雷达***误差的状态转移方程的S2阶段:
本阶段假设目标状态及雷达***误差模型均满足线性高斯特性,进而建立目标及雷达***误差的状态转移方程,然后进入下一步S3阶段。
在此阶段中,所述的建立目标及雷达***误差的状态转移方程的具体方法如下:假设目标状态及雷达***误差模型均满足线性高斯特性,考虑目标在二维平面内做线性运动,则目标状态及雷达***误差转移方程可以分别表示为
其中,Fk-1表示目标的状态转移矩阵,G为***噪声系数矩阵,wk-1为状态过程噪声,且wk-1~N(·,0,Qk),vk-1为雷达***误差过程噪声,且vk-1~N(·,0,Qb,k)。
3)利用GM-EPHD滤波器对雷达***误差进行融合估计的S3阶段:
本阶段假设目标的状态模型及观测模型均满足线性高斯特性,目标的存活概率与目标检测概率相互独立,并且假设雷达***误差也满足线性高斯特性。为了实现雷达***误差和目标状态的联合估计,将雷达***误差扩展至目标状态中组成扩展状态,使用高斯混合-概率假设密度滤波器递归跟踪扩展状态。在此过程中,由于ADS-B的定位精度远远优于雷达的定位精度,先使用坐标转换后的ADS-B观测数据对目标状态进行高斯混合-概率假设密度滤波,以获得较为准确的目标状态估计,然后使用雷达观测数据对目标状态更新后的扩展状态进行高斯混合-扩展概率假设密度滤波。由于雷达***误差和目标状态在雷达测量更新步骤通过似然函数耦合在一起,此时采用两步卡尔曼滤波器可同时估计出雷达***误差和目标状态。
在此阶段中,所述的利用GM-EPHD滤波器对雷达***误差进行融合估计的具体方法如下:
a、假设扩展状态新生、衍生过程的强度函数可以写成高斯混合形式
其中,Jγ,k表示新生目标的个数,分别表示第j个高斯混合形式新生目标强度函数的均值、雷达***误差均值、***噪声协方差、误差协方差。Jβ,k表示衍生目标的个数,xk-1、bk-1、分别表示第i个高斯混合形式衍生目标强度函数的均值、雷达***误差均值、***噪声协方差、误差协方差。
b、预测步:
假设扩展状态后验强度函数可以写成高斯混合形式
则预测的强度函数也可以写成高斯混合形式
υk|k-1=υS,k|k-1+υβ,k|k-1+γk(xk,bk) (14)
其中,υS,k|k-1表示扩展状态存活过程的强度函数,υβ,k|k-1表示扩展状态衍生过程的强度函数,γk(xk,bk)表示扩展状态新生过程的强度函数,且由式(14)给出。
由于雷达***误差和目标状态相互独立,在预测步中,采用卡尔曼滤波器分别传递雷达***误差和目标状态的高斯混合形式强度函数的权值、均值及协方差。
c、更新步:
假设预测强度函数也可以表示成如式(15)所示的高斯混合形式
则使用坐标转换后的ADS-B测量数据对预测强度函数进行更新,得到的后验强度函数为
其中
ADS-B测量数据更新后的高斯混合形式目标强度函数的权值、均值及协方差分别为
其中,
由于ADS-B测量中不含杂波,故
当使用ADS-B测量数据对预测强度函数更新后,目标状态高斯项的权值、均值及协方差通过卡尔曼滤波器得到传递。然后使用雷达测量数据再次更新强度函数,由于雷达***误差和目标状态通过似然函数耦合在一起,采用两步卡尔曼滤波器可以依次更新目标状态和雷达***误差高斯项的权值、均值及协方差。
雷达测量数据更新后的后验强度函数为
其中
雷达测量数据更新后的目标状态高斯项的权值、均值及协方差分别为
雷达***误差高斯项的均值、协方差分别为
其中,
λ为杂波数,c(zk)为杂波分布函数,这里假设服从观测空间内的均匀分布。
d、对更新获得的扩展状态的强度函数的高斯项进行裁剪合并,获得目标状态和雷达***误差估计,方法同GM-PHD算法相同。
e、雷达***误差融合估计
雷达***误差融合估计的方法是使用加权平均算法对获得的多个雷达***误差估计值进行加权融合。
本发明提供的利用ADS-B数据的GM-EPHD滤波雷达***误差联合估计方法的效果可以通过以下仿真结果进一步说明。
仿真数据描述:实验设置在一个[-1000,1000]×[-1000,1000]的监视区域,3个目标运动,一部二坐标雷达处于监视区域的中心,观测获得目标距离、方位信息。3个目标均配备机载ADS-B收发信机,地面站接收广播信息获得目标经度、纬度及高度信息,ADS-B观测值转换到以雷达为中心的直角坐标系下。仿真雷达及ADS-B同步采样,周期均为T=1s,仿真100步,雷达检测概率为PD,k=PD=0.99。
雷达***误差真值为(50m,0.01rad)。
雷达随机测量噪声协方差矩阵
雷达观测杂波数为λ=10,杂波分布函数c(zk)服从[-1000,1000]×[-1000,1000]上的均匀分布。
***加速度扰动标准差为σv=0.1m/s2,***误差噪声协方差矩阵Bk=diag{(5m)2,(0.001rad)2}。
高斯项修剪门限T=10-5,合并门限U=4,高斯项最大个数Jmax=100。
实验运行环境为Intel Core2 Quad CPU 2.66GHz,2GB内存,仿真软件为MatlabR2010a。
为了比较本发明方法的估计结果,选取文献(Wenling Li,Jia Y,Du J,YuF.Gaussian mixture PHD filter for multi-sensor multi-target tracking withregistration errors[J].Signal Processing,2012)中的GM-EPHD联合估计方法,利用相同数据进行100次蒙特卡罗比较试验,文献方法中的“sensor1”和“sensor2”都存在***误差,利用“sensor1”和“sensor2”的融合观测对“sensor1”进行***误差估计,而本发明中使用ADS-B数据代替文献方法中的“sensor2”,使用ADS-B观测和“sensor1”观测联合估计目标状态和“sensor1”的***误差。
图2为实验所设置的目标真实运动轨迹及雷达和ADS-B观测数据,‘—*—’为雷达观测值,‘—o—’为ADS-B观测值,‘——’为目标位置的真实值。其中目标1从k=1s时刻出发,k=100s结束,其初始位置为(250m,250m);目标2从k=1s时刻出发,k=100s结束,其初始位置为(-250m,-250m);目标3在k=30s时从目标1中分开,直到k=100s结束。
图3为目标真实运动轨迹及状态估计结果,‘—o—’为目标状态估计值,‘——’为目标位置的真实值;图4为目标状态估计的OSPA距离,实验中设置p=2,c=200;图5为目标数估计结果,‘—△—’为目标数估计值,‘—o—’为目标数真实值。从图3、图4可以看出,使用本发明方法进行目标状态估计的平均误差近似为5m,估计出的目标状态和真实状态接近。从图5可以看出,采用本发明方法对目标数目估计结果和实际情况相吻合。
图6、图7分别为采用两种方法获得的雷达距离、方位误差估计比较图,‘—*—’为采用本发明方法估计结果,‘—△—’为采用文献方法的估计结果,‘——’为雷达***误差真值,对比估计结果可以看出,采用本发明方法估计稳定后雷达***误差估计结果更接近***误差真值。图8、图9分别为采用两种方法获得的雷达距离、方位误差估计的均方根误差比较图,‘—*—’为采用本发明方法估计的均方根误差,‘—△—’为采用文献方法的均方根误差,对比结果可以看出,本发明方法获得的雷达***误差估计的均方根误差要远远小于文献方法。
图10为采用本发明方法对雷达***误差校正后的雷达观测,‘—*—’为雷达***误差校正前的观测值,‘—o—’为采用本发明方法对雷达***误差校正后的观测值,‘——’为目标状态真实值,对比校正前后的雷达观测值可以看出,使用本发明方法对雷达***误差进行校正后,雷达观测与真实运动轨迹已经非常接近了,证明本发明方法对雷达***误差估计是准确有效的。
Claims (1)
1.一种利用ADS-B数据的GM-EPHD滤波雷达***误差联合估计方法,所述的利用ADS-B数据的GM-EPHD滤波雷达***误差联合估计方法包括按顺序进行的下列步骤:
1)建立目标的ADS-B及雷达的观测方程的S1阶段;
2)建立目标及雷达***误差的状态转移方程的S2阶段;
3)利用GM-EPHD滤波器对雷达***误差进行融合估计的S3阶段;
其特征在于:在步骤1)中,所述的建立目标的ADS-B及雷达的观测方程的方法是首先利用坐标投影技术将目标的ADS-B观测值转换到以雷达站位置为中心的直角坐标系下,由于ADS-B监视数据的定位精度远远优于雷达的定位精度,在不考虑ADS-B定位误差的基础上,建立坐标转换后的目标的ADS-B观测方程,然后考虑雷达***误差,建立目标的雷达观测方程,为利用GM-EPHD滤波器对雷达***误差进行融合估计做准备;
在步骤2)中,所述的建立目标及雷达***误差的状态转移方程的方法是假设目标状态及雷达***误差均满足线性高斯特性,进而建立目标及雷达***误差的状态转移方程,为利用GM-EPHD滤波器对雷达***误差进行融合估计做准备;
在步骤3)中,所述的利用GM-EPHD滤波器对雷达***误差进行融合估计的方法是首先将雷达***误差扩展至目标状态中组成扩展状态,根据S2阶段中建立的目标及雷达***误差的状态转移方程分别对目标状态及雷达***误差进行预测,然后使用S1阶段中坐标转换后的ADS-B观测数据对目标状态进行高斯混合-概率假设密度滤波,以获取较为准确的目标状态估计,接着使用S1阶段中的雷达观测数据对扩展状态进行高斯混合-扩展概率假设密度滤波,进而采用两步卡尔曼滤波器联合估计出目标状态和雷达***误差,最后使用加权平均算法对雷达***误差进行融合估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510218847.1A CN104849702B (zh) | 2015-04-30 | 2015-04-30 | 利用ads‑b数据的gm‑ephd滤波雷达***误差联合估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510218847.1A CN104849702B (zh) | 2015-04-30 | 2015-04-30 | 利用ads‑b数据的gm‑ephd滤波雷达***误差联合估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104849702A CN104849702A (zh) | 2015-08-19 |
CN104849702B true CN104849702B (zh) | 2017-10-27 |
Family
ID=53849486
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510218847.1A Expired - Fee Related CN104849702B (zh) | 2015-04-30 | 2015-04-30 | 利用ads‑b数据的gm‑ephd滤波雷达***误差联合估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104849702B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106291488B (zh) * | 2016-08-16 | 2018-08-03 | 中国人民解放军防空兵学院 | 一种雷达标定误差校正方法 |
CN106896352B (zh) * | 2017-04-17 | 2020-01-10 | 电子科技大学 | 一种基于随机集理论的多雷达异步数据分布式融合方法 |
CN107390191A (zh) * | 2017-05-23 | 2017-11-24 | 中国民航大学 | 一种ecef坐标系下概率假设密度滤波雷达空间误差配准方法 |
CN108333569B (zh) * | 2018-01-19 | 2021-01-12 | 杭州电子科技大学 | 一种基于phd滤波的异步多传感器融合多目标跟踪方法 |
CN108344981B (zh) * | 2018-01-19 | 2020-08-25 | 杭州电子科技大学 | 面向杂波的多传感器异步检测tsbf多目标跟踪方法 |
CN110703272B (zh) * | 2019-09-27 | 2022-05-17 | 江苏大学 | 一种基于车车通信和gmphd滤波的周边目标车辆状态估计方法 |
CN111580053A (zh) * | 2020-05-30 | 2020-08-25 | 中国人民解放军海军航空大学 | 大型船舶目标时变雷达散射中心测量误差联合计算及修正方法 |
CN112083409A (zh) * | 2020-09-17 | 2020-12-15 | 北京博能科技股份有限公司 | 一种航班定位方法、装置和电子设备 |
CN113257044B (zh) * | 2021-07-09 | 2022-02-11 | 中航信移动科技有限公司 | Ads-b数据的滤波方法、装置、计算机设备及存储介质 |
CN117521018B (zh) * | 2024-01-08 | 2024-03-26 | 鹏城实验室 | 基于扩展观测的融合估计方法、装置、设备及存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101894381A (zh) * | 2010-08-05 | 2010-11-24 | 上海交通大学 | 动态视频序列中多目标跟踪*** |
CN103345577A (zh) * | 2013-06-27 | 2013-10-09 | 江南大学 | 变分贝叶斯概率假设密度多目标跟踪方法 |
CN103729637A (zh) * | 2013-12-31 | 2014-04-16 | 西安工程大学 | 基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法 |
CN104156984A (zh) * | 2014-05-30 | 2014-11-19 | 西北工业大学 | 一种不均匀杂波环境下多目标跟踪的概率假设密度方法 |
CN104237862A (zh) * | 2014-09-18 | 2014-12-24 | 中国民航大学 | 基于ads-b的概率假设密度滤波雷达***误差融合估计方法 |
CN104318059A (zh) * | 2014-09-24 | 2015-01-28 | 深圳大学 | 用于非线性高斯***的目标跟踪方法和跟踪*** |
CN104501812A (zh) * | 2014-12-05 | 2015-04-08 | 江南大学 | 基于自适应新生目标强度的滤波算法 |
-
2015
- 2015-04-30 CN CN201510218847.1A patent/CN104849702B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101894381A (zh) * | 2010-08-05 | 2010-11-24 | 上海交通大学 | 动态视频序列中多目标跟踪*** |
CN103345577A (zh) * | 2013-06-27 | 2013-10-09 | 江南大学 | 变分贝叶斯概率假设密度多目标跟踪方法 |
CN103729637A (zh) * | 2013-12-31 | 2014-04-16 | 西安工程大学 | 基于容积卡尔曼滤波的扩展目标概率假设密度滤波方法 |
CN104156984A (zh) * | 2014-05-30 | 2014-11-19 | 西北工业大学 | 一种不均匀杂波环境下多目标跟踪的概率假设密度方法 |
CN104237862A (zh) * | 2014-09-18 | 2014-12-24 | 中国民航大学 | 基于ads-b的概率假设密度滤波雷达***误差融合估计方法 |
CN104318059A (zh) * | 2014-09-24 | 2015-01-28 | 深圳大学 | 用于非线性高斯***的目标跟踪方法和跟踪*** |
CN104501812A (zh) * | 2014-12-05 | 2015-04-08 | 江南大学 | 基于自适应新生目标强度的滤波算法 |
Non-Patent Citations (1)
Title |
---|
Multi-Target tracking with background discrimination using PHD filters;Roland Jonsson 等;《Information Fusion (FUSION), 2012 15th International Conference on》;20120712;正文854,856页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104849702A (zh) | 2015-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104849702B (zh) | 利用ads‑b数据的gm‑ephd滤波雷达***误差联合估计方法 | |
CN1940591B (zh) | 使用传感器融合进行目标跟踪的***和方法 | |
CN106054169B (zh) | 基于跟踪信息的多站雷达信号融合检测方法 | |
CN105761276B (zh) | 基于迭代ransac自适应新生目标强度估计的gm-phd多目标跟踪方法 | |
CN102692621A (zh) | 一种ads-b与雷达的联合***误差估计方法 | |
CN103901413A (zh) | 基于旋翼无人直升机的三坐标雷达高度动态标定设备及其方法 | |
CN108153980A (zh) | 基于ads-b与tcas数据融合的综合显示方法 | |
CN104237862B (zh) | 基于ads‑b的概率假设密度滤波雷达***误差融合估计方法 | |
CN113342059B (zh) | 基于位置和速度误差的多无人机跟踪移动辐射源方法 | |
Aernouts et al. | Combining TDoA and AoA with a particle filter in an outdoor LoRaWAN network | |
CN104297739A (zh) | 一种通航监视中对光电跟踪设备的引导方法 | |
CN110068793A (zh) | 一种定位跟踪方法 | |
Yao et al. | Integrity monitoring for Bluetooth low energy beacons RSSI based indoor positioning | |
CN113075648A (zh) | 一种无人集群目标定位信息的聚类与滤波方法 | |
Ouyang et al. | Cooperative navigation of UAVs in GNSS-denied area with colored RSSI measurements | |
Dasanayaka et al. | Analysis of vehicle location prediction errors for safety applications in cooperative-intelligent transportation systems | |
CN113933876B (zh) | 多星通讯时差定位数据融合处理方法 | |
CN104463841A (zh) | 衰减系数自适应的滤波方法及滤波*** | |
CN111816005A (zh) | 基于ads-b的远程驾驶飞机环境监视优化方法 | |
CN107390191A (zh) | 一种ecef坐标系下概率假设密度滤波雷达空间误差配准方法 | |
CN102426357B (zh) | 一种具有图像确认的多目标跟踪方法 | |
Džunda et al. | Influence of mutual position of communication network users on accuracy of positioning by telemetry method | |
Lombaerts et al. | Development and Field Test Results of Distributed Ground-Based Sensor Fusion Object Tracking | |
RU2807613C1 (ru) | Способ трассирования наземных и морских радиоизлучающих целей | |
Lombaerts et al. | Development and Field Test Results of Distributed Ground Sensor Fusion Based Object Tracking |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171027 Termination date: 20200430 |
|
CF01 | Termination of patent right due to non-payment of annual fee |