CN103941233A - 一种多平台主被动传感器协同跟踪的雷达间歇交替辐射控制方法 - Google Patents
一种多平台主被动传感器协同跟踪的雷达间歇交替辐射控制方法 Download PDFInfo
- Publication number
- CN103941233A CN103941233A CN201410082846.4A CN201410082846A CN103941233A CN 103941233 A CN103941233 A CN 103941233A CN 201410082846 A CN201410082846 A CN 201410082846A CN 103941233 A CN103941233 A CN 103941233A
- Authority
- CN
- China
- Prior art keywords
- radar
- measurement
- ned
- passive
- pas
- 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/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
-
- 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
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
本发明公开了一种多平台主被动传感器协同跟踪的雷达辐射控制方法,对多个搭载主被动传感器的平台进行雷达辐射协同管控,在满足跟踪任务需求的条件下,控制雷达在时间上间歇开机、在空间上交替开机,并利用基于坐标交替转换的多平台主被动序贯滤波实现多平台辐射管控下的融合跟踪问题,该跟踪方法不需要多平台被动传感器量测时间同步,也不需要多平台被动传感器的交叉定位,能够适应被动量测不连续和无规律出现的情况,具有良好的适应性,适合工程应用,解决雷达长时间开机辐射易被侦察的问题,进一步提高雷达的隐蔽能力,进而提高多平台***整体生存能力。
Description
一、技术领域
本发明隶属于传感器管理领域,适用于解决多平台主被动传感器协同跟踪的雷达辐射控制问题。
二、背景技术
现代复杂电磁环境下,各种侦察和干扰设备日趋完善,为了尽可能提高作战平台的隐蔽性,防止已方雷达被对方侦察到进而受到干扰和攻击,就必须在满足跟踪精度的前提下尽可能减少雷达的电磁辐射,提高作战平台的隐蔽性和反侦察能力。
现有的雷达辐射控制以及辐射控制下的主被动传感器融合跟踪技术主要集中在单平台,如单机载平台雷达、红外、电子支援措施的协同跟踪与雷达开关机管理等,这些方法只对单平台上的雷达在时间上进行间歇开机管理,而多平台环境下,不同平台被动传感器只能提供目标方位角、俯仰角,不能提供目标距离,无法将被动量测点坐标转换到同一平台坐标系下进行滤波跟踪,而双平台交叉定位方法虽然在一定程度上解决目标距离未知的问题,但一般需要被动量测时间同步,难以解决被动量测不同步或间歇出现的情况。可见,现有雷达辐射控制方法不能适用于多平台辐射管理下多平台主被动协同跟踪。
三、发明内容
本发明所要解决的问题就是,针对现有方法不适用多平台协同雷达辐射管控的问题,提供一种多平台主被动传感器协同跟踪的雷达辐射控制方法,对多个搭载主被动传感器的平台进行雷达辐射协同管控,在满足跟踪任务需求的条件下,控制雷达时间上间歇开机,空间上交替开机,并实现了多平台雷达间歇交替辐射时的主被动融合跟踪,提高多平台探测***的隐蔽性。
本发明解决所述技术问题,采用技术方案步骤如下:
(1)被动传感器引导雷达开机完成目标航迹起始;
(2)利用协方差控制方法在时间上确定下一次需要雷达辐射的时刻;
(3a)当预定开机时间到,没有得到新的被动量测,则跳到第(4a)步;
(3b)当预定开机时间未到,有被动量测时,首先用所得到的被动量测点迹与融合航迹进行关联,然后将关联上的被动量测利用基于坐标交替转换的多平台主被动序贯滤波进行跟踪,回到第(2)步;
(4a)在空间上选择不同平台交替开机辐射;
(4b)雷达开机得到新的雷达量测后,先利用所得雷达量测点迹与融合航迹进行关联,然后将关联上的雷达量测利用基于坐标交替转换的多平台主被动序贯滤波算法进行跟踪,跳回到第(2)步。
具体的,所述步骤(2)具体为:
a1)假设雷达最小和最大辐射间隔分别用Tmin、Tmax表示,并令Tleft=Tmin,Tright=Tmax,进入下一步;
a2)如果Tright-Tleft<ε(ε为一小的常数),则令下次雷达开机辐射间隔的时间T1=Tleft或T1=Tright,跳到第a5)步,否则令Ttest=Tleft+[(Tright-Tleft)/2];
a3)令t1=Ttest,计算间隔时间为t1时扩展卡尔曼滤波的预测协方差P(t1|0)
P(t1|0)=F(t1)·P(0|0)·F′(t1)+Γ(t1)·σv 2·Γ′(t1)
其中,P(0|0)表示最近一次滤波协方差估计,F(t1)表示状态转移矩阵,Γ(t1)表示过程噪声分布矩阵,σv 2为零均值高斯白噪声的协方差;
a4)判断预测协方差是否超过门限。假设跟踪精度门限为Rth,如果 (其中分别是预测协方差矩阵P(t1|0)的对角线上X、Y、Z位置误差方差项),则Tleft=Ttest,否则Tright=Tlest,跳回到第a2)步;
a5)得到下次预定雷达开机辐射间隔时间为T1。
具体的,所述步骤(3b)具体为:
b1)得到与融合航迹关联的被动量测后,首先将该被动量测到来前最近的融合航迹点 以及对应的协方差P(k|k)转换到获得该量测的传感器所在平台的载体坐标系下,得到对应的航迹点P1 NED(k|k);
b2)利用P1 NED(k|k)进行外推,得到状态预测
预测协方差
P1 NED(k+1|k)=F·P1 NED(k|k)·FT
b3)然后将载体坐标系下的预测航迹和对应的协方差P1 NED(k+1|k)从X-Y-Z直角坐标下转换到ρ-θ-ε(距离-方位-俯仰)的极坐标系
Pi SP(k+1|k)=Uk TP1 NED(k+1|k)|(x,y,z)Uk
式中,Uk为直角坐标到球极坐标的转换矩阵。
构造关联统计量
g1(k+1)=d1(k)T[P1+R1]-1d1(k)
其中,P1为Pi SP(k|k+1)在方位和俯仰上的误差协方差矩阵,R1为被动量测在方位俯仰上的量测误差矩阵,d1(k)为在方位和俯仰上的新息。
利用如下公式进行量测与航迹的关联判决
g1(k+1)<th 关联
g1(k+1)>th 不关联
其中th是二自由度χ2分布的统计量g1(k)的假设检验门限。
b4)得到与融合航迹关联的被动量测Z(k+1)=[θPAS(k+1),εPAS(k+1)]′后,先利用 计算被动量测方位和俯仰的预测值
其中坐标系以载体为坐标原点的北-天-东为X-Y-Z方向(以下均以此为坐标系),arctan2表示matlab软件中对应的反正切函数;
利用计算被动量测时扩展卡尔曼滤波的量测矩阵
其中,HPAS(k+1)对应的状态向量表示形式为
利用雷达测量精度性能指标,得到量测误差矩阵
其中,分别表示被动量测的方位和俯仰的误差方差。
将被动量测Z(k+1)=[θPAS(k+1),εPAS(k+1)]′、被动量测时扩展卡尔曼滤波的量测矩阵HPAS、量测误差矩阵R(k+1)、被动量测方位和俯仰的预测值h(k+1|k)分别带入扩展卡尔曼滤波公式:
P1 NED(k+1|k+1)=[I-K(k+1)HPAS(k+1)]P1 NED(k+1|k)
其中,I表示与状态向量维数相同的单位矩阵。
b5)然后将和P1 NED(k+1|k+1)分别转换回原融合航迹的坐标系下,成为 原航迹的下一个滤波状态估计点及对应的协方差。
本发明的有益效果是,可以解决现代战场复杂电场环境下雷达长时间开机辐射易被侦察的问题,控制雷达在时间上间歇开机,在空间上交替开机,提高雷达的隐蔽能力,进一步提高多平台***整体生存能力。本发明利用基于坐标交替转换的多平台主被动序贯滤波实现多平台辐射管控下的融合跟踪问题,该方法不需要多平台被动传感器量测时间同步,也不需要多平台被动传感器的交叉定位,能够适应被动量测不连续和无规律出现的情况,适应性强,适合工程应用。
四、附图说明
附图1本发明的方法步骤流程图;
附图2是用本发明进行仿真实验时己方飞机和目标飞机的飞行轨迹图
附图3是用本发明进行仿真实验时被动传感器量测数据率为1s一个时2架己方飞机雷达辐射情况
附图4是用本发明进行仿真实验时被动传感器量测数据率为1s一个时目标跟踪误差
附图5是用本发明进行仿真实验时被动传感器量测数据率为10s一个时雷达辐射情况
附图6是用本发明进行仿真实验时被动传感器量测数据率为10s一个时跟踪误差
附图7是用现有单平台方法进行仿真实验时被动传感器量测数据率为1s一个时己方飞机1雷达辐射情况
附图8是用现有单平台方法进行仿真实验时被动传感器量测数据率为1s一个时目标跟踪误差
附图9是用现有单平台方法进行仿真实验时被动传感器量测数据率为10s一个时己方飞机1雷达辐射情况
附图10是用现有单平台方法进行仿真实验时被动传感器量测数据率为10s一个时目标跟踪误差
五、具体实施方法
下面参照附图对本发明创造做进一步详细说明
参照图1,本发明的具体步骤如下
步骤一:首先使多个平台的被动传感器一直处于探测状态,并利用被动传感器引导两个平台中的某一雷达进行开机探测,得到雷达测量后,将雷达量测利用逻辑法进行航迹起始,利用三点法完成滤波初始化,得到初始化的状态估计和协方差P(0|0),并通过坐标换转,将状态估计和协方差转换到融合中心坐标系下,得到融合中心的初始状态和初 始协方差P1 NED(0|0);
步骤二:根据目标跟踪精度的任务需求,预测下一次雷达开机的时间,具体步骤为:
a1)假设雷达最小和最大辐射间隔分别用Tmin、Tmax表示,并令Tleft=Tmin,Tright=Tmax,进入下一步;
a2)如果Tright-Tleft<ε(ε为一小的常数),则令下次雷达开机辐射间隔的时间T1=Tleft或T1=Tright,跳到第a5)步,否则令Ttest=Tleft+[(Tright-Tleft)/2];
a3)令t1=Ttest,计算间隔时间为t1时扩展卡尔曼滤波的预测协方差P(t1|0)
P(t1|0)=F(t1)·P(0|0)·F′(t1)+Γ(t1)·σv 2·Γ′(t1)
其中,P(0|0)表示最近一次滤波协方差估计,F(t1)表示状态转移矩阵,Γ(t1)表示过程噪声分布矩阵,σv 2为零均值高斯白噪声的协方差;
a4)判断预测协方差是否超过门限。假设跟踪精度门限为Rth,如果 (其中分别是预测协方差矩阵P(t1|0)的对角线上X、Y、Z位置误差方差项),则Tleft=Ttest,否则Tright=Ttest,跳回到第a2)步;
a5)得到下次预定雷达开机辐射间隔时间为T1。
步骤三:1)当雷达预定开机时间到,没有得到新的被动量测,则跳到步骤四;否则,当预定开机时间未到,有被动量测时,先将被动量测点迹与融合航迹进行关联判决,关联判决的具体步骤为:
b1)如果得到新的被动量测,首先将融合航迹和对应的协方差转换到该被动量测所在平台的载体坐标系下;
b2)然后将载体坐标系下的航迹和对应的协方差从以该载体为原点的X-Y-Z直角坐标下转换到对应的ρ-θ-ε(距离-方位-俯仰)的极坐标系
Xi SP(k+1|k)=Uk TXi(k+1|k)|(x,y,z)Uk
Pi SP(k+1|k)=Uk TPi(k+1|k)|(x,y,z)Uk
式中,Uk为直角坐标到极坐标的转换矩阵。
b3)构造关联统计量
g1(k+1)=d1(k)T[P1+R1]-1d1(k)
其中,P1为Pi SP(k|k+1)在方位和俯仰上的误差协方差矩阵,R1为被动量测在方位俯仰上的量测误差矩阵,d1(k)为被动量测在方位和俯仰上的新息。
b4)利用如下公式进行量测与航迹的关联判决
g1(k+1)<th1 关联
g1(k+1)>th1 不关联
其中th1是二自由度χ2分布的统计量g1(k)的假设检验门限。
2)将关联上的被动量测利用基于坐标交替转换的多平台主被动序贯滤波进行跟踪,然后跳回到步骤二。其中基于坐标交替转换的多平台主被动序贯滤波进行跟踪的具体步骤为:
c1)得到新的被动量测后,首先将该被动量测到来前最近的融合航迹点以及对应的协方差P(k|k)转换到获得该量测的传感器所在平台的载体坐标系下,得到对应的航迹点 P1 NED(k|k);
c2)利用P1 NED(k|k)进行外推,得到状态预测
预测协方差
P1 NED(k+1|k)=F·P1 NED(k|k)·FT
利用计算被动量测方位和俯仰的预测值
其中坐标系以载体为坐标原点的北-天-东为X-Y-Z方向(以下均以此为坐标系),arctan2表示matlab软件中对应的反正切函数;
利用计算被动量测时扩展卡尔曼滤波的量测矩阵
其中,状态向量表示形式为
利用雷达测量精度性能指标,得到量测误差矩阵
其中,分别表示被动量测的方位和俯仰的误差方差。
将被动量测Z(k+1)=[θPAS(k+1),εPAS(k+1)]′、HPAS、R(k+1)、HPAS、h(k+1|k)带入扩展卡尔曼滤波公式:
P1 NED(k+1|k+1)=[I-K(k+1)HPAS(k+1)]P1 NED(k+1|k)
其中,I表示与状态向量维数相同的单位矩阵。
c3)然后将和P1 NED(k+1|k+1)分别转换回原融合航迹的坐标系下,成为原航迹的下一个滤波状态估计点及对应的协方差。
步骤四:1)在空间上选择不同的平台上的雷达进行交替辐射,选择方法可采用顺序交替或按平台的优先等级交替;
2)当雷达开机辐射,得到新的雷达量测,先将雷达量测点迹与融合航迹进行关联判决,具体步骤为:
d1)如果得到新的雷达量测,首先将融合航迹和对应的协方差转换到该雷达量测所在平台的载体坐标系下;
d2)然后将载体坐标系下的航迹和对应的协方差从以载体为原点的X-Y-Z直角坐标下转换到ρ-θ-ε(距离-方位-俯仰)的极坐标系
Xi SP(k+1|k)=Uk TXi(k+1|k)|(x,y,z)Uk
Pi SP(k+1|k)=Uk TPi(k+1|k)|(x,y,z)Uk
式中,Uk为直角坐标到极坐标的转换矩阵。
d3)构造关联统计量
g2(k+1)=d2(k)T[P2+R2]-1d2(k)
其中,P2为Pi SP(k|k+1)在距离、方位和俯仰上的误差协方差矩阵,R2为雷达量测在距离、方位、俯仰上的量测误差矩阵,d2(k)为雷达量测在距离、方位和俯仰上的新息。
d4)利用如下公式进行量测与航迹的关联判决
g2(k+1)<th2 关联
g2(k+1)>th2 不关联
其中th2是三自由度的χ2分布的统计量g2(k)的假设检验门限。
3)将关联上的雷达量测利用基于坐标交替转换的多平台主被动序贯滤波算法进行跟踪,跳回到步骤二,其中雷达量测利用基于坐标交替转换的多平台主被动序贯滤波具体步骤如下:
e1)得到新的雷达量测后,首先将该被动量测到来前最近的融合航迹点以及对应的滤波协方差P(k|k)转换到获得该量测的传感器所在平台的载体坐标系下,得到对应的航迹点P1 NED(k|k);
e2)利用P1 NED(k|k)进行外推,得到一步预测
预测协方差
P1 NED(k+1|k)=F·P1 NED(k|k)·FT
利用计算雷达量测距离、方位和俯仰的预测值
其中坐标系以载体为坐标原点的北-天-东为X-Y-Z方向(以下均以此为坐标系),arctan2表示matlab软件中对应的反正切函数;
利用计算雷达量测时扩展卡尔曼滤波的量测矩阵
其中,状态向量表示形式为
利用雷达测量精度性能指标,得到雷达量测误差矩阵
其中,分别表示雷达量测的距离、方位和俯仰的误差方差。
将雷达量测Z(k+1)=[ρradar(k+1),θradar(k+1),εradar(k+1)]′、Hradar、R(k+1)、HPAS、h(k+1|k)带入扩展卡尔曼滤波公式:
P1 NED(k+1|k+1)=[I-K(k+1)HPAS(k+1)]P1 NED(k+1|k)
其中,I表示与状态向量维数相同的单位矩阵。
e3)然后将和P1 NED(k+1|k+1)分别坐标转换回到原融合航迹的坐标系下,成为原航迹的下一个滤波状态估计点及对应的协方差。
本发明的效果可以通过以下matlab仿真结果进一步说明:
仿真实验条件:
设有2架己方作战飞机,初始地理坐标纬度、经度、高度分别为(36.14°,120°,10000m)、(36.23°,120°,10000m),其纬度、经度、高度方向上初始速度分别为(0,-150m/s,0)、(0,-150m/s,0),且均作匀速直线运动,仿真总时间110s,飞行的轨迹图如图2所示。每架飞机上分别搭载一个雷达、一个红外传感器和一个电子支援措施,雷达能测量目标的距离、方位、俯仰,红外传感器和电子支援措施能测量目标的方位、俯仰,目标飞机1架,其地理坐标纬度、经度、高度为(36.18°,119.5°,10000m),其纬度、经度、高度方向上初始速度为(0,200m/s,0);机载雷达最小辐射间隔为1秒,被动量测的数据率分两种情况:1)机载被动传感器量测每1秒有一个量测,2)机载被动传感器量测每10秒有一个量测;机载雷达测距误差为100m,方位角均方根误差为0.2°,俯仰角均方根误差为0.2°;机载红外传感器方位角均方根误差为0.1°,俯仰角均方根误差为0.1°;机载电子资源措施的方位角均方根误差为0.3°,俯仰角均方根误差为0.3°。跟踪任务需求精度为100米以内。
实验一:
用本发明方法对2架己方飞机探测目标过程中的雷达辐射进行管控,在被动传感器量测数据率为每1秒一个时,2架己方飞机雷达辐射情况如图3所示,融合跟踪误差如图4所示。在被动传感器量测数据率为每10秒一个时,2架己方飞机雷达辐射情况如图5所示,融合跟踪误差如图6所示。
从图3可以看出,在需求精度限定条件下,己方飞机1和飞机2的雷达在时间上间歇开机,随着跟踪精度的提高开机时间间隔越来越大,并在空间上交替开机。同时,由于2架飞机的红外和电子支援措施等被动量测的辅助跟踪,雷达在一段时间后完全不需要开机就能满足跟踪需求。
比较图3和图5可以看出,雷达辐射次数除了与跟踪精度需求有关,还与被动量测的数据率有关,即被动量测数据高时雷达辐射次数少,被动量测数据率低时雷达辐射次数多,雷达能自适应根据被动量测的数据率来调整辐射的时间。
实验二:
在同一仿真场景下,用现有的单平台方法来让己方飞机1完成对目标的跟踪和雷达辐射管控,对应的在被动传感器量测数据率为每1秒一个时,己方飞机1雷达辐射情况如图7所示,融合跟踪误差如图8所示。在被动传感器量测数据率为每10秒一个时,1架飞机搭载的 雷达辐射情况如图9所示,融合跟踪误差如图10所示。本发明的方法与传统的方法在相同的仿真场景下对应的雷达辐射次数比较如表1所示
表1 现有方法与本发明方法对应的雷达辐射次数比较
从表1可以看出现有单平台的雷达辐射方法中在被动传感器量测间隔为1s,本发明方法比现有方法飞机1雷达辐射次数减少约70%,在被动传感器量测间隔为10s时,飞机1雷达辐射次数减少约53%,极大的减少了飞机的雷达辐射时间,提高了平台的隐蔽性。
Claims (3)
1.一种多平台主被动传感器协同跟踪的雷达间歇交替辐射控制方法,其特征在于包括以下步骤:
(1)被动传感器引导雷达开机完成目标航迹起始;
(2)利用协方差控制方法在时间上确定下一次需要雷达辐射的时刻;
(3a)当预定开机时间到,没有得到新的被动量测,则跳到第(4a)步;
(3b)当预定开机时间未到,有被动量测时,首先用所得到的被动量测点迹与融合航迹进行关联,然后将关联上的被动量测利用基于坐标交替转换的多平台主被动序贯滤波进行跟踪,回到第(2)步:
(4a)在空间上选择不同平台交替开机辐射;
(4b)雷达开机得到新的雷达量测后,先利用所得雷达量测点迹与融合航迹进行关联,然后将关联上的雷达量测利用基于坐标交替转换的多平台主被动序贯滤波算法进行跟踪,跳回到第(2)步。
2.根据权利要求1所述的一种多平台主被动传感器协同跟踪的雷达间歇交替辐射控制方法,其特征在于,所述步骤(2)具体为:
a1)假设雷达最小和最大辐射间隔分别用Tmin、Tmax表示,并令Tleft=Tmin,Tright=Tmax,进入下一步;
a2)如果Tright-Tleft<ε(ε为一小的常数),则令下次雷达开机辐射间隔的时间T1=Tleft或T1=Tright,跳到第a5)步,否则令Ttest=Tlest+[(Tright-Tleft)/2];
a3)令t1=Tlest,计算间隔时间为t1时扩展卡尔曼滤波的预测协方差P(t1|0)
P(t1|0)=F(t1)·P(0|0)·F′(t1)+Γ(t1)·σv 2·Γ′(t1)
其中,P(0|0)表示最近一次滤波协方差估计,F(t1)表示状态转移矩阵,Γ(t1)表示过程噪声分布矩阵,σv 2为零均值高斯白噪声的协方差;
a4)判断预测协方差是否超过门限。假设跟踪精度门限为Rth,如果 (其中分别是预测协方差矩阵P(t1|0)的对角线上X、Y、Z位置误差方差项),则Tleft=Ttest,否则Tright=Tlest,跳回到第a2)步;
a5)得到下次预定雷达开机辐射间隔时间为T1。
3.根据权利要求1所述的一种多平台主被动传感器协同跟踪的雷达间歇交替辐射控制方法,其特征在于,所述步骤(3b)具体为:
b1)得到与融合航迹关联的被动量测后,首先将该被动量测到来前最近的融合航迹点以及对应的协方差P(k|k)转换到获得该量测的传感器所在平台的载体坐标系下,得到对应的航迹点P1 NED(k|k);
b2)利用P1 NED(k|k)进行外推,得到状态预测
预测协方差
P1 NED(k+1|k)=F·P1 NED(k|k)·FT
b3)然后将载体坐标系下的预测航迹和对应的协方差P1 NED(k+1|k)从X-Y-Z直角坐标下转换到ρ-θ-ε(距离-方位-俯仰)的极坐标系
Pi SP(k+1|k)=Uk TP1 NED(k+1|k)|(x,y,z)Uk
式中,Uk为直角坐标到球极坐标的转换矩阵。
构造关联统计量
g1(k+1)=d1(k)T[P1+R1]-1d1(k)
其中,P1为Pi SP(k|k+1)在方位和俯仰上的误差协方差矩阵,R1为被动量测在方位俯仰上的量测误差矩阵,d1(k)为在方位和俯仰上的新息。
利用如下公式进行量测与航迹的关联判决
g1(k+1)<th 关联
g1(k+1)>th 不关联
其中th是二自由度χ2分布的统计量g1(k)的假设检验门限。
b4)得到与融合航迹关联的被动量测Z(k+1)=[θPAS(k+1),εPAS(k+1)]′后,先利用计算被动量测方位和俯仰的预测值
其中坐标系以载体为坐标原点的北-天-东为X-Y-Z方向(以下均以此为坐标系),arctan2表示matlab软件中对应的反正切函数。
利用计算被动量测时扩展卡尔曼滤波的量测矩阵
其中,HPAS(k+1)对应的状态向量表示形式为
利用雷达测量精度性能指标,得到量测误差矩阵
其中,分别表示被动量测的方位和俯仰的误差方差。
将被动量测Z(k+1)=[θPAS(k+1),εPAS(k+1)]′、被动量测时扩展卡尔曼滤波的量测矩阵HPAS、量测误差矩阵R(k+1)、被动量测方位和俯仰的预测值h(k+1|k)分别带入扩展卡尔曼滤波公式:
P1 NED(k+1|k+1)=[I-K(k+1)HPAS(k+1)]P1 NED(k+1|k)
其中,I表示与状态向量维数相同的单位矩阵。
b5)然后将和P1 NED(k+1|k+1)分别转换回原融合航迹的坐标系下,成为原航迹的下一个滤波状态估计点及对应的协方差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410082846.4A CN103941233B (zh) | 2014-03-04 | 2014-03-04 | 一种多平台主被动传感器协同跟踪的雷达间歇交替辐射控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410082846.4A CN103941233B (zh) | 2014-03-04 | 2014-03-04 | 一种多平台主被动传感器协同跟踪的雷达间歇交替辐射控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103941233A true CN103941233A (zh) | 2014-07-23 |
CN103941233B CN103941233B (zh) | 2016-06-29 |
Family
ID=51188979
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410082846.4A Expired - Fee Related CN103941233B (zh) | 2014-03-04 | 2014-03-04 | 一种多平台主被动传感器协同跟踪的雷达间歇交替辐射控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103941233B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105116916A (zh) * | 2015-09-25 | 2015-12-02 | 北京机械设备研究所 | 一种分布式光电跟踪***协同跟踪方法 |
CN105676180A (zh) * | 2016-01-05 | 2016-06-15 | 中国航空无线电电子研究所 | 一种基于多级拍卖算法的多平台点迹融合方法 |
CN106443665A (zh) * | 2016-09-13 | 2017-02-22 | 中国人民解放军海军航空工程学院 | 基于辅助距离的雷达与esm航迹关联方法 |
CN106789740A (zh) * | 2016-11-11 | 2017-05-31 | 中国人民解放军海军航空工程学院 | 按任务优先级序贯拍卖的多平台传感器协同管理方法 |
CN106932759A (zh) * | 2017-01-17 | 2017-07-07 | 电子科技大学 | 一种用于主被动雷达的协同定位方法 |
CN106980114A (zh) * | 2017-03-31 | 2017-07-25 | 电子科技大学 | 无源雷达目标跟踪方法 |
CN109116349A (zh) * | 2018-07-26 | 2019-01-01 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 多传感器协同跟踪联合优化决策方法 |
CN109164419A (zh) * | 2018-08-15 | 2019-01-08 | 中国电子科技集团公司第二十研究所 | 基于交互式多模型的多平台无序量测处理方法 |
CN109239704A (zh) * | 2018-08-21 | 2019-01-18 | 电子科技大学 | 一种基于序贯滤波交互式多模型的自适应采样方法 |
CN111624594A (zh) * | 2020-05-12 | 2020-09-04 | 中国电子科技集团公司第三十八研究所 | 一种基于转换量测重构的组网雷达跟踪方法 |
CN111829506A (zh) * | 2020-06-12 | 2020-10-27 | 中国船舶重工集团公司第七二四研究所 | 一种基于序贯处理的主被动航迹运动特征关联方法 |
CN111985379A (zh) * | 2020-08-13 | 2020-11-24 | 中国第一汽车股份有限公司 | 基于车载雷达的目标追踪方法、装置、设备及车辆 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7719461B1 (en) * | 2008-08-05 | 2010-05-18 | Lockheed Martin Corporation | Track fusion by optimal reduced state estimation in multi-sensor environment with limited-bandwidth communication path |
CN101770024A (zh) * | 2010-01-25 | 2010-07-07 | 上海交通大学 | 多目标跟踪方法 |
US8416122B1 (en) * | 2010-05-13 | 2013-04-09 | Lockheed Martin Corporation | Point-in-polygon target location |
CN103424742A (zh) * | 2013-07-24 | 2013-12-04 | 深圳大学 | 一种序贯处理测量数据的目标跟踪方法与目标跟踪*** |
-
2014
- 2014-03-04 CN CN201410082846.4A patent/CN103941233B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7719461B1 (en) * | 2008-08-05 | 2010-05-18 | Lockheed Martin Corporation | Track fusion by optimal reduced state estimation in multi-sensor environment with limited-bandwidth communication path |
CN101770024A (zh) * | 2010-01-25 | 2010-07-07 | 上海交通大学 | 多目标跟踪方法 |
US8416122B1 (en) * | 2010-05-13 | 2013-04-09 | Lockheed Martin Corporation | Point-in-polygon target location |
CN103424742A (zh) * | 2013-07-24 | 2013-12-04 | 深圳大学 | 一种序贯处理测量数据的目标跟踪方法与目标跟踪*** |
Non-Patent Citations (2)
Title |
---|
吴巍 等: ""多机载平台多目标跟踪与辐射控制"", 《***工程与电子技术》, vol. 34, no. 3, 31 March 2012 (2012-03-31) * |
吴巍 等: ""给定辐射时间限制下的雷达辐射时机最优规划研究"", 《兵工学报》, vol. 33, no. 1, 31 January 2012 (2012-01-31) * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105116916A (zh) * | 2015-09-25 | 2015-12-02 | 北京机械设备研究所 | 一种分布式光电跟踪***协同跟踪方法 |
CN105676180A (zh) * | 2016-01-05 | 2016-06-15 | 中国航空无线电电子研究所 | 一种基于多级拍卖算法的多平台点迹融合方法 |
CN105676180B (zh) * | 2016-01-05 | 2018-06-01 | 中国航空无线电电子研究所 | 一种基于多级拍卖算法的多平台点迹融合方法 |
CN106443665A (zh) * | 2016-09-13 | 2017-02-22 | 中国人民解放军海军航空工程学院 | 基于辅助距离的雷达与esm航迹关联方法 |
CN106443665B (zh) * | 2016-09-13 | 2018-09-14 | 中国人民解放军海军航空大学 | 基于辅助距离的雷达与esm航迹关联方法 |
CN106789740B (zh) * | 2016-11-11 | 2019-12-10 | 中国人民解放军海军航空大学 | 按任务优先级序贯拍卖的多平台传感器协同管理方法 |
CN106789740A (zh) * | 2016-11-11 | 2017-05-31 | 中国人民解放军海军航空工程学院 | 按任务优先级序贯拍卖的多平台传感器协同管理方法 |
CN106932759A (zh) * | 2017-01-17 | 2017-07-07 | 电子科技大学 | 一种用于主被动雷达的协同定位方法 |
CN106980114A (zh) * | 2017-03-31 | 2017-07-25 | 电子科技大学 | 无源雷达目标跟踪方法 |
CN109116349A (zh) * | 2018-07-26 | 2019-01-01 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 多传感器协同跟踪联合优化决策方法 |
CN109116349B (zh) * | 2018-07-26 | 2022-12-13 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 多传感器协同跟踪联合优化决策方法 |
CN109164419A (zh) * | 2018-08-15 | 2019-01-08 | 中国电子科技集团公司第二十研究所 | 基于交互式多模型的多平台无序量测处理方法 |
CN109164419B (zh) * | 2018-08-15 | 2022-07-05 | 中国电子科技集团公司第二十研究所 | 基于交互式多模型的多平台无序量测处理方法 |
CN109239704A (zh) * | 2018-08-21 | 2019-01-18 | 电子科技大学 | 一种基于序贯滤波交互式多模型的自适应采样方法 |
CN109239704B (zh) * | 2018-08-21 | 2023-03-10 | 电子科技大学 | 一种基于序贯滤波交互式多模型的自适应采样方法 |
CN111624594A (zh) * | 2020-05-12 | 2020-09-04 | 中国电子科技集团公司第三十八研究所 | 一种基于转换量测重构的组网雷达跟踪方法 |
CN111624594B (zh) * | 2020-05-12 | 2022-09-23 | 中国电子科技集团公司第三十八研究所 | 一种基于转换量测重构的组网雷达跟踪方法及*** |
CN111829506A (zh) * | 2020-06-12 | 2020-10-27 | 中国船舶重工集团公司第七二四研究所 | 一种基于序贯处理的主被动航迹运动特征关联方法 |
CN111985379A (zh) * | 2020-08-13 | 2020-11-24 | 中国第一汽车股份有限公司 | 基于车载雷达的目标追踪方法、装置、设备及车辆 |
Also Published As
Publication number | Publication date |
---|---|
CN103941233B (zh) | 2016-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103941233A (zh) | 一种多平台主被动传感器协同跟踪的雷达间歇交替辐射控制方法 | |
CN106646450B (zh) | 基于距离分步聚类的雷达航迹抗差关联方法 | |
CN109116349A (zh) | 多传感器协同跟踪联合优化决策方法 | |
CN104484717B (zh) | 基于当前航路的一种机载航迹预测方法 | |
CN103759732B (zh) | 一种角度信息辅助的集中式多传感器多假设跟踪方法 | |
CN106597369B (zh) | 一种无人机的控制方法、控制平台、控制*** | |
CN102707276A (zh) | 基于目标拓扑信息的ais与雷达航迹抗差关联算法 | |
CN107765226A (zh) | 一种sar卫星雷达回波模拟方法、***和介质 | |
CN108089183A (zh) | 一种用于异步多基地雷达***的检测跟踪一体化方法 | |
CN103247057A (zh) | 目标-回波-道路网三元数据关联的道路目标多假设跟踪算法 | |
CN106323272B (zh) | 一种获取航迹起始轨迹的方法及电子设备 | |
CN107153186A (zh) | 激光雷达标定方法及激光雷达 | |
CN111601253A (zh) | 无源被动智能跟踪定位方法、***、存储介质、跟踪定位终端 | |
CN103926569B (zh) | 三坐标雷达网基于交叉定位点关联的集中式与分布式压制干扰鉴别方法 | |
CN101308206B (zh) | 一种白噪声背景下的圆周轨迹机动目标跟踪方法 | |
CN105676170B (zh) | 一种用于多目标无源定位与跟踪的虚假点去除方法 | |
CN103499809B (zh) | 一种纯方位双机协同目标跟踪定位路径规划方法 | |
CN110196601A (zh) | 无人机控制方法、装置、***和计算机可读存储介质 | |
CN104237862A (zh) | 基于ads-b的概率假设密度滤波雷达***误差融合估计方法 | |
CN103698767A (zh) | 获取空间目标与观测站距离的方法及装置 | |
CN105607054A (zh) | 收发一体式超宽带飞行雷达定时采样探测方法 | |
CN110426689A (zh) | 一种基于em-cks的机载多平台多传感器***误差配准算法 | |
CN104457756B (zh) | 一种基于双机测距的海面物体定位方法 | |
CN105242248A (zh) | 一种基于测控设备的雷达挂飞试验位置参数自动装订方法 | |
CN108957438A (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: 20160629 Termination date: 20200304 |
|
CF01 | Termination of patent right due to non-payment of annual fee |