CN108490465B - 基于时频差与测向的地面同频多运动辐射源跟踪方法及系统 - Google Patents

基于时频差与测向的地面同频多运动辐射源跟踪方法及系统 Download PDF

Info

Publication number
CN108490465B
CN108490465B CN201810228309.4A CN201810228309A CN108490465B CN 108490465 B CN108490465 B CN 108490465B CN 201810228309 A CN201810228309 A CN 201810228309A CN 108490465 B CN108490465 B CN 108490465B
Authority
CN
China
Prior art keywords
filtering
radiation source
target
motion
frequency
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.)
Active
Application number
CN201810228309.4A
Other languages
English (en)
Other versions
CN108490465A (zh
Inventor
尤明懿
陆安南
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CETC 36 Research Institute
Original Assignee
CETC 36 Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by CETC 36 Research Institute filed Critical CETC 36 Research Institute
Priority to CN201810228309.4A priority Critical patent/CN108490465B/zh
Publication of CN108490465A publication Critical patent/CN108490465A/zh
Application granted granted Critical
Publication of CN108490465B publication Critical patent/CN108490465B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/421—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • G01S19/423—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system by combining or switching between position solutions derived from different satellite radio beacon positioning 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

本发明公开一种基于时频差与测向的地面同频多运动辐射源跟踪方法及系统，属于定位技术领域，其方法通过建立双星时频差定位系统的多运动辐射源目标的量测模型和迭代滤波模型；确定滤波模型迭代滤波的迭代停止条件；实现对未知运动状态的包含同频运动辐射源目标进行跟踪。本发明将模糊的时差、频差量测全部作为目标量测进行滤波，避免了解模糊的复杂问题；通过设置权重阈值，确定最佳运动目标滤波估计结果，实现了对同频多运动辐射源的跟踪。对于经常同时接收到多个同频信号的UHF、L/S频段目标，本发明对于提升低频段多运动辐射源的跟踪具有较大的借鉴意义。

Description

基于时频差与测向的地面同频多运动辐射源跟踪方法及系统
技术领域
本发明涉及定位技术领域，尤其是一种基于时频差与测向的地面同频多运动辐射源跟踪方法及系统。
背景技术
无源定位有着广泛的军事与商业应用，为了提供精确辐射源位置估计结果，双星(站)时频差定位体制得到了广泛应用。双星(站)时频差定位体制中的重要步骤是根据卫星接收到的辐射源信号估计其到达双星(站)的时差(TDOA)与频差(FDOA)。然而，当多个辐射源频率相近甚至相同，且信号样式一致时，在估计时频差时就无法判断双星(站)各自接收到的辐射源信号的对应关系，即产生时频差模糊，因而就无法实现辐射源的精确定位。
目前，对于UHF、L/S等较低频段信号，同时接收到多个辐射源的同频信号是很常见的，在对多运动辐射源跟踪上，由于存在频差模糊问题，跟踪效果很差，甚至无法跟踪，但，缺乏有效的解决途径。
发明内容
鉴于上述的分析，本发明旨在提供一种基于时频差与测向的地面同频多运动辐射源跟踪方法及系统，用以解决同频多运动辐射源定位时由由时频差模糊造成的无法跟踪的问题，实现对同频多运动辐射源的跟踪。
本发明的目的主要是通过以下技术方案实现的：
一种基于时频差与测向的地面同频多运动辐射源跟踪方法，包括：
基于双星时频差定位系统，建立同频多运动辐射源目标的量测模型和迭代滤波模型；
对已知运动状态的同频多运动辐射源采用上述量测模型和迭代滤波模型进行量测和迭代滤波，以确定迭代滤波模型的迭代停止条件；
采用建立的量测模型和确定迭代停止条件的滤波模型对地面同频多运动辐射源目标进行跟踪。
进一步地，所述迭代滤波模型采用的滤波算法为GM-UKF-PHD滤波算法。
进一步地，所述基于双星时频差定位系统，建立同频多运动辐射源目标的量测模型，包括：
建立多运动辐射源目标的量测模型；
在多运动辐射源目标的量测模型的基础上，结合时频差模糊，得到同频多运动辐射源目标的量测模型。
进一步地，所述多运动辐射源目标的量测模型包括地面运动辐射源状态转移方程和量测方程；
所述地面运动辐射源状态转移方程为：X(k+1)＝F·X(k)+Q；式中，
Figure BDA0001599349060000021
xe(k)为地面运动辐射源在时刻k的位置矢量，
Figure BDA0001599349060000022
为地面辐射源在时刻k的速度矢量；
Figure BDA0001599349060000023
其中ω1、ω2为位置状态转移误差，ω3、ω4为速度状态转移误差；
得到所述多运动辐射源目标的量测模型z(k)：
Figure BDA0001599349060000031
Figure BDA0001599349060000032
为辐射源信号至主星与辅星的理论时差；
Figure BDA0001599349060000033
为辐射源信号至主星与辅星的理论频差；
vt(k)为时差量测误差；
vf(k)为频差量测误差；
Figure BDA0001599349060000034
为主星对运动辐射源量测的俯仰角理论值；
Figure BDA0001599349060000035
为主星对运动辐射源量测的方位角理论值；
Figure BDA0001599349060000036
为俯仰角量测误差；
vθ(k)为方位角量测误差；
结合时频差模糊，得到同频多运动辐射源目标的量测模型zj(k)：
Figure BDA0001599349060000037
式中，
Figure BDA0001599349060000038
指辐射源ej的信号传播至主星的理论时间，
Figure BDA0001599349060000039
指辐射源ei的信号传播至辅星的理论时间，
Figure BDA00015993490600000310
指主星接收到辐射源ej的理论信号频率，
Figure BDA00015993490600000311
指辅星接收到辐射源ei的理论信号频率，
Figure BDA00015993490600000312
分别为主星对辐射源ei和ej测向得到的理论俯仰角与方位角，N为辐射源数。
进一步地，所述迭代滤波运算包括：
1)对滤波器进行初始化；
2)采用GM-UKF-PHD滤波算法进行滤波，并迭代计算每一个目标滤波估计结果的权重值；
3)对获取的目标权重，进行归一化处理。
进一步地，所述迭代停止是通过设置权重阈值实现的，所述权重阈值的确定方法为：
1)对已知运动状态的地面同频多静止辐射源依据建立多运动辐射源目标的量测模型和滤波模型，进行量测和滤波，得到不同权重值的目标滤波估计结果；
2)将已知的运动状态与不同权重值的目标滤波估计结果进行运动位置比对；
3)找出位置误差最小的目标滤波估计结果，其对应的权重值即为权重阈值。
进一步地，采用建立的量测模型和确定迭代停止条件的迭代滤波模型对未知运动状态包含同频运动辐射源的目标进行跟踪，当滤波模型输出的目标滤波估计结果与权重阈值最接近时，即认为达到迭代停止条件，输出跟踪结果。
一种基于时频差与测向的地面同频多运动辐射源跟踪系统，包括多目标量测模块、多目标滤波模块和多目标跟踪状态提取模块；
所述多目标量测模块根据建立双星时频差定位系统的量测模型；对多运动辐射源目标进行量测，估算目标的时频差和对目标测向，得到量测矢量，输出到所述多目标滤波模块；
所述多目标滤波模块对所述多目标量测模块输出的量测矢量进行GM-UKF-PHD滤波，得到不同权重的运动辐射源目标状态；
所述多目标跟踪状态提取模块与所述多目标滤波模块相连，根据设置的权重阈值，停止多目标滤波模块的迭代运算，提取多运动辐射源的跟踪状态。
进一步地，所述多目标滤波模块包括初始值估计模块、滤波模块和更新模块；
所述初始值估计模块连接滤波模块，为滤波模块的滤波提供初始输入信息；
所述滤波模块连接所述初始值估计模和更新模块，接收所述初始值估计模块的初始输入信息，开始进行GM-UKF-PHD滤波；将每次滤波的结果存储到更新模块；在初始滤波处理后，接收更新模块输出的上一次滤波结果，进行迭代GM-UKF-PHD滤波；
所述更新模块的输入与输出与所述预测模块连接，所述更新模块存储上一次滤波模块的滤波结果，并将存储的滤波结果输出到滤波模块进行当前的迭代滤波。
进一步地，当多目标滤波模块输出的目标滤波估计结果对应的权重与权重阈值最接近时，输出跟踪结果。
本发明有益效果如下：
将模糊的时差、频差量测全部作为目标量测进行滤波，避免了解模糊的复杂问题。通过多步滤波，解除了时差频差模糊，同时取得更高精度的辐射源定位结果，通过设置权重阈值，确定最佳运动目标滤波估计结果，作为跟踪结果输出，实现了对同频多运动辐射源的跟踪。
附图说明
附图仅用于示出具体实施例的目的，而并不认为是对本发明的限制，在整个附图中，相同的参考符号表示相同的部件。
图1为基于时频差与测向的地面同频多运动辐射源跟踪方法流程图；
图2为双星定位系统坐标体系图；
图3为基于时频差与测向的地面同频多运动辐射源跟踪系统组成示意图；
图4权重阈值为0.5时的目标跟踪情况图；
图5为权重阈值为0.25时的目标跟踪情况图；
图6为权重阈值为0.1时的目标跟踪情况图；
图7为权重阈值为0.01时的目标跟踪情况图。
具体实施方式
下面结合附图来具体描述本发明的优选实施例，其中，附图构成本申请一部分，并与本发明的实施例一起用于阐释本发明的原理。
本发明的一个具体实施例，公开了基于时频差与测向的地面同频多运动辐射源跟踪方法，如图1所示，包括以下步骤：
步骤S1、基于双星时频差定位系统，建立多运动辐射源目标的量测模型和迭代滤波模型；
所述量测模型包括地面运动辐射源状态转移方程和量测方程；
所述双星时频差定位系统的坐标体系如图2所示，
其中，
地面辐射源E在时刻k的位置矢量为：xe＝(xe(k),ye(k),0)T，
地面辐射源E在时刻k的速度矢量为：
Figure BDA0001599349060000061
双星(站)系统中主星(站)在时刻k的位置矢量为：xs1＝(xs1(k),ys1(k),zs1(k))T；
双星(站)系统中主星(站)在时刻k的速度矢量为：
Figure BDA0001599349060000071
双星(站)系统中辅星(站)在时刻k已知的位置矢量为：xs2＝(xs2(k),ys2(k),zs2(k))T；
双星(站)系统中辅星(站)在时刻k已知的速度矢量为：
Figure BDA0001599349060000072
对于地面运动辐射源E，其状态转移方程为：
Figure BDA0001599349060000073
式(3)中，
Figure BDA0001599349060000074
其中ωi为符合均值为0方差为
Figure BDA0001599349060000075
的状态转移误差，根据辐射源运动特性确定。通常根根类似辐射源的运动轨迹，采用轨迹曲线拟合结果残差的方差作为状态转移误差参数ω1、ω2的方差估计，而采用轨迹变化率(即速度)曲线拟合结果残差的方差作为状态转移误差参数ω3、ω4的方差估计。
所述量测方程中运动辐射源的量测由时频差量测方程和测向方程构成，具体为：
Figure BDA0001599349060000076
z(k)为在时刻k对地面辐射源的量测矢量；
Figure BDA0001599349060000077
为辐射源信号至主星与辅星的理论时差；
Figure BDA0001599349060000081
为辐射源信号至主星与辅星的理论频差；
vt(k)为时差量测误差；
vf(k)为频差量测误差；
Figure BDA0001599349060000082
为主星对运动辐射源量测的俯仰角理论值；
Figure BDA0001599349060000083
为主星对运动辐射源量测的方位角理论值；
为俯仰角量测误差；
vθ(k)为方位角量测误差。
其中：
所述时频差量测方程为：
Figure BDA0001599349060000085
式中，Δt(k)、Δf(k)分别为时刻k辐射源E辐射的信号传播至主星(站)与辅星(站)的时间差与频率差，c＝300000km/s为光速，||·||为矢量的模，vt(k)为时差量测误差，一般有
Figure BDA0001599349060000086
即时差量测误差符合均值为零，方差为
Figure BDA0001599349060000087
的正态分布，fe为辐射源信号频率，vf(k)为时差量测误差，一般有
Figure BDA0001599349060000088
即时差量测误差符合均值为零，方差为
Figure BDA0001599349060000089
的正态分布。te-s1、te-s2分别为辐射源信号传播至主星(站)与辅星(站)的时间。
所述测向方程为由主星完成对运动辐射源测向，即量测辐射源的俯仰角
Figure BDA00015993490600000810
与方位角
Figure BDA00015993490600000811
有：
Figure BDA0001599349060000091
式中，
Figure BDA0001599349060000092
为俯仰角量测误差，一般有
Figure BDA0001599349060000093
即时差量测误差符合均值为零，方差为
Figure BDA0001599349060000094
的正态分布，vθ(k)为方位角量测误差，一般有
Figure BDA0001599349060000095
即时差量测误差符合均值为零，方差为
Figure BDA0001599349060000096
的正态分布。
由于同频多运动辐射源的量测中存在时频差模糊，结合时频差模糊，得到同频多运动辐射源目标的量测模型zj(k)：
Figure BDA0001599349060000097
式中，
Figure BDA0001599349060000098
指辐射源ej的信号传播至主星的理论时间，
Figure BDA0001599349060000099
指辐射源ei的信号传播至辅星的理论时间，
Figure BDA00015993490600000910
指主星接收到辐射源ej的理论信号频率，
Figure BDA00015993490600000911
指辅星接收到辐射源ei的理论信号频率，
Figure BDA00015993490600000912
分别为主星对辐射源ei和ej测向得到的理论俯仰角与方位角，N为辐射源数；式中，只有当i＝j时，时频差量测为正确量测，而当i≠j时的时频差量测为虚假量测。需要注意的是，式中，仅时频差量测存在模糊，而测向结果不存在模糊。
对于N个运动辐射源的跟踪，将N个辐射源在时刻k的状态组合可视为目标状态集X(k)＝[xe1(k),...,xei(k),...,xeN(k)]T，i＝1,...,N，xei(k)为第i个辐射源的目标状态；对N个辐射源的所有量测结果进行组合，得到一个量测集：Z(k)＝[z1(k),...,zi×j(k),...,zN×N(k)]T，j＝1,...,N,i＝1,...,N；
由于，存在时频差模糊，难以识别主、副星(站)接收到的辐射源信号的配对关系，影响对辐射源的定位，无法对运动辐射源跟踪，如果在滤波前采用去模糊方法去除存在时频差模糊，其去模糊方法本身存在复杂性；本实施例，为了减少运算的复杂性，建立迭代滤波模型，直接将包含时差频差模糊的全部量测结果作为目标量测进行迭代滤波，避免了解模糊的复杂问题。
特殊的，所述迭代滤波模型采用高斯混合无迹卡尔曼滤波概率假设密度函数(GM-UKF-PHD)滤波算法进行；
在滤波算法中，定义
Figure BDA0001599349060000101
为(k-1)时刻滤波算法的混合高斯分布集，其中
Figure BDA0001599349060000102
为分布i的权重，
Figure BDA0001599349060000103
为分布i的均值矢量，
Figure BDA0001599349060000104
为分布i的协方差矩阵；Jk-1为(k-1)时刻的进入滤波目标数，i＝1,…,Jk-1。
所述滤波算法具体包括以下步骤：
1)、对滤波器进行初始化
所述初始化包括：
分布i的初始权重
Figure BDA0001599349060000105
分布i的各个辐射源目标状态的初始估计
Figure BDA0001599349060000106
由其他定位手段在先引导获取或者根据在先掌握的情报信息获得，所述其他定位手段包括测向定位，光学定位等；
分布i的协方差矩阵
Figure BDA0001599349060000111
的初始值根据目标位置初始估计精度设置，为了避免目标位置初始估计精度可能较差，设置较大的
Figure BDA0001599349060000112
遵循的原理是协方差矩阵
Figure BDA0001599349060000113
参数的设置尽可能覆盖目标位置的全部可能区域。其效果即是避免因
Figure BDA0001599349060000114
参数设置过小导致滤波过程发散。如：可设置为
Figure BDA0001599349060000115
辐射源目标数初始估计J0，J0是根据先验信息或起始时刻量测数量进行估计的；
辐射源目标状态转移的协方差矩阵Qk-1为目标状态转移过程噪声，且Qk-1＝Q；
量测矢量的协方差矩阵
Figure BDA0001599349060000116
2)采用GM-UKF-PHD滤波方法对每一个目标，基于每一组量测进行滤波，并迭代计算每一个目标滤波估计结果的权重值。
在第k次迭代滤波过程中，根据上一次滤波的归一化权重，对本次权重值进行赋值，即
Figure BDA0001599349060000117
式中，
l为第k次滤波时的量测序号变量，l＝1,...,Lk，Lk为第k次滤波时的量测数；
j为第k次滤波的目标序号变量，j＝1,...,Jk，Jk为进入第k次滤波的目标数；
Figure BDA0001599349060000125
为目标j在第k-1次滤波的归一化权重；
N(A；B,C)的含义为对于均值为B，方差为C的多元正态分布，矢量A的概率密度；
zkl为在第k次滤波时得到的第l组量测；
Figure BDA0001599349060000122
为基于GM-UKF-PHD滤波方法，根据目标j在第k-1次的滤波结果，其在第k次滤波时的量测的预测值；
Figure BDA0001599349060000123
为基于GM-UKF-PHD滤波方法，根据目标j在第k-1次的滤波结果，其在第k次滤波时的量测的协方差矩阵预测值。
3)对获取的目标权重，进行归一化处理；
为统一选取权重预置，根据公式
Figure BDA0001599349060000124
对获取的每一个进入第k次滤波的目标权重，进行归一化处理。
步骤S2、对已知运动状态的同频多运动辐射源采用上述量测模型和迭代滤波模型进行量测和迭代滤波，以确定滤波模型的迭代滤波的迭代停止条件；
所述迭代停止是通过设置权重阈值实现的，所述权重阈值的确定方法为：
1)对已知运动状态的地面同频多静止辐射源依据建立多运动辐射源目标的量测模型和滤波模型，进行量测和滤波，得到不同权重值的目标滤波估计结果；
2)将已知的运动状态与不同权重值的目标滤波估计结果进行运动位置比对；
3)找出位置误差最小的目标滤波估计结果，其对应的权重值即为权重阈值wT。
步骤S3、对包含同频的多运动辐射源目标进行跟踪。
采用建立的量测模型和确定迭代停止条件的迭代滤波模型对包含同频的多运动辐射源目标进行跟踪，通过量测模型的量测得到包含时差频差模糊的全部量测结果，将全部量测结果送入迭代滤波模型进行迭代滤波，当滤波模型输出的目标滤波估计结果的权重值与权重阈值最接近时，即认为达到迭代停止条件，输出跟踪结果。
一种基于时频差与测向的地面同频多运动辐射源跟踪系统，如图3所示，包括多目标量测模块、多目标滤波模块和多目标跟踪状态提取模块；
所述多目标量测模块根据建立双星时频差定位系统的量测模型；对多个辐射源目标进行量测，估算目标的时频差，得到量测矢量，对多个辐射源目标的运动俯仰角和方位角进行量测；
特殊的，由于存在多个同频辐射源，在进行量测估算时频差时，会得到含模糊时频差信息的量测矢量。
所述多目标滤波模块对所述多目标量测模块输出的量测矢量进行迭代滤波运算，得到不同权重的辐射源目标状态；
所述多目标滤波模块由初始值估计模块、滤波模块与更新模块组成。
所述初始值估计模块连接滤波模块，为滤波模块的滤波提供初始输入信息；所述初始值估计模块的输入项为其他定位手段在先引导获取或者根据在先掌握的情报信息，输出为目标数初始估计J0、各个目标状态的初始估计
Figure BDA0001599349060000131
各个目标状态估计的初始协方差矩阵
Figure BDA0001599349060000132
目标状态转移的协方差矩阵Q，量测矢量的协方差矩阵R。
所述滤波模块连接所述初始值估计模和更新模块，接收所述初始值估计模块的初始输入信息，开始进行GM-UKF-PHD滤波；将每次的滤波的结果存储到更新模块；并在除初始滤波外，接收更新模块输出的上一次滤波结果，进行当前的GM-UKF-PHD滤波。
所述更新模块的输入与输出与所述预测模块连接，所述更新模块存储上一次滤波模块的滤波结果，并将存储的滤波结果输出到滤波模块进行当前的迭代滤波；
所述更新模块的输出为滤波模块上一次滤波所得的进入滤波目标数Jk-1，上一次滤波所得的目标权重
Figure BDA0001599349060000141
上一次滤波所得的目标状态估计结果
Figure BDA0001599349060000142
上一次滤波所得的目标状态估计协方差矩阵
Figure BDA0001599349060000143
所述更新模块的输入项目包括：滤波模块滤波更新后输出的进入下一次滤波目标数Jk，滤波更新后的目标权重
Figure BDA0001599349060000144
滤波更新后的目标状态估计结果
Figure BDA0001599349060000145
滤波更新后的目标状态估计协方差矩阵
Figure BDA0001599349060000146
所述多目标跟踪状态提取模块与所述多目标滤波模块的更新模块相连，根据设置权重阈值wT，停止多目标滤波模块的迭代运算，提取多运动辐射源的跟踪状态；具体：根据设置的权重阈值wT，当滤波结果
Figure BDA0001599349060000147
最接近时wT，停止多目标滤波模块的迭代运算，提取更新模块存储的对应
Figure BDA0001599349060000148
作为输出的多目标状态；所述权重阈值wT根据上述权重阈值的确定方法确定的。
考虑在图2所示的参考坐标系内，以两个同频运动辐射源为例：
两个辐射源的实际起始位置分别为(250，250，0)与(-250，-250，0)，匀速直线运动速度分别为(2.5，11.5，0)与(-11.5，-2.5，0)。主星(站)在各个时刻的位置为：(4k,4k,500)，辅星(站)在各个时刻的位置为：(500-10k,500+5k,1000)。假设两辐射源发射同频信号，且频率fe＝2×107Hz。基于工程可实现性的考虑，设双星定位系统测时误差σt＝50ns，测频误差σf＝10Hz，俯仰、方位角量测误差均为1°。假设状态转移协方差矩阵为：Q＝0.52×[11 0 1 1 0]T。图4～图7给出了不同状态提取阈值情况下的目标跟踪结果，其中·代表目标实际位置，○代表目标位置估计结果。由图3～图6的结果可见，状态提取阈值的选取对目标跟踪结果影响较大，当合理选取状态提取阈值(如wT＝0.25)时，本专利公开的跟踪系统能够较好地跟踪两同频运动辐射源的状态。
综上所述，本发明实施例提供的基于时频差与测向的地面同频多运动辐射源跟踪方法及系统将模糊的时差、频差量测全部作为目标量测进行滤波，避免了解模糊的复杂问题。通过多步滤波，解除了时差频差模糊，同时取得更高精度的辐射源定位结果，通过设置权重阈值，确定最佳运动目标滤波估计结果，作为跟踪结果输出，实现了对同频多运动辐射源的跟踪。
本领域技术人员可以理解，实现上述实施例方法的全部或部分流程，可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于计算机可读存储介质中。其中，所述计算机可读存储介质为磁盘、光盘、只读存储记忆体或随机存储记忆体等。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。

Claims (9)

1.一种基于时频差与测向的地面同频多运动辐射源跟踪方法，其特征在于，包括：
基于双星时频差定位系统，建立同频多运动辐射源目标的量测模型和迭代滤波模型；
多运动辐射源目标的量测模型包括地面运动辐射源状态转移方程和量测方程；
所述地面运动辐射源状态转移方程为：X(k+1)＝F·X(k)+Q；式中，
Figure FDA0002459157220000011
xe(k)为地面运动辐射源在时刻k的位置矢量，
Figure FDA0002459157220000012
为地面辐射源在时刻k的速度矢量；
Figure FDA0002459157220000013
其中ω1、ω2为位置状态转移误差，ω3、ω4为速度状态转移误差；
所述多运动辐射源目标的量测模型z(k)：
Figure FDA0002459157220000014
Figure FDA0002459157220000015
为辐射源信号至主星与辅星的理论时差；
Figure FDA0002459157220000016
为辐射源信号至主星与辅星的理论频差；
vt(k)为时差量测误差；
vf(k)为频差量测误差；
Figure FDA0002459157220000021
为主星对运动辐射源量测的俯仰角理论值；
Figure FDA0002459157220000022
为主星对运动辐射源量测的方位角理论值；
Figure FDA0002459157220000023
为俯仰角量测误差；
vθ(k)为方位角量测误差；
结合时频差模糊，得到同频多运动辐射源目标的量测模型zj(k)：
Figure FDA0002459157220000024
式中，
Figure FDA0002459157220000025
指辐射源ej的信号传播至主星的理论时间，
Figure FDA0002459157220000026
指辐射源ei的信号传播至辅星的理论时间，
Figure FDA0002459157220000027
指主星接收到辐射源ej的理论信号频率，
Figure FDA0002459157220000028
指辅星接收到辐射源ei的理论信号频率，
Figure FDA0002459157220000029
分别为主星对辐射源ei和ej测向得到的理论俯仰角与方位角，N为辐射源数；
对已知运动状态的同频多运动辐射源采用上述量测模型和迭代滤波模型进行量测和迭代滤波，以确定迭代滤波模型的迭代停止条件；
采用建立的量测模型和确定迭代停止条件的滤波模型对地面同频多运动辐射源目标进行跟踪。
2.根据权利要求1所述的地面同频多运动辐射源跟踪方法，其特征在于，所述迭代滤波模型采用的滤波算法为GM-UKF-PHD滤波算法。
3.根据权利要求1或2所述的地面同频多运动辐射源跟踪方法，其特征在于，所述基于双星时频差定位系统，建立同频多运动辐射源目标的量测模型，包括：
建立多运动辐射源目标的量测模型；
在多运动辐射源目标的量测模型的基础上，结合时频差模糊，得到同频多运动辐射源目标的量测模型。
4.根据权利要求1或2所述的地面同频多运动辐射源跟踪方法，其特征在于，
所述迭代滤波运算包括：
1)对滤波器进行初始化；
2)采用GM-UKF-PHD滤波算法进行滤波，并迭代计算每一个目标滤波估计结果的权重值；
3)对获取的目标权重，进行归一化处理。
5.根据权利要求1所述的地面同频多运动辐射源跟踪方法，其特征在于，所述迭代停止是通过设置权重阈值实现的，所述权重阈值的确定方法为：
1)对已知运动状态的地面同频多静止辐射源依据建立多运动辐射源目标的量测模型和滤波模型，进行量测和滤波，得到不同权重值的目标滤波估计结果；
2)将已知的运动状态与不同权重值的目标滤波估计结果进行运动位置比对；
3)找出位置误差最小的目标滤波估计结果，其对应的权重值即为权重阈值。
6.根据权利要求5所述的地面同频多运动辐射源跟踪方法，其特征在于，采用建立的量测模型和确定迭代停止条件的迭代滤波模型对未知运动状态多同频运动辐射源目标进行跟踪，当滤波模型输出目标滤波估计结果的权重值与权重阈值最接近时，即认为达到迭代停止条件，输出跟踪结果。
7.一种采用权利要求1-6任一项所述跟踪方法的地面同频多运动辐射源跟踪系统，其特征在于，包括多目标量测模块、多目标滤波模块和多目标跟踪状态提取模块；
所述多目标量测模块，用于根据建立双星时频差定位系统的量测模型，对多运动辐射源目标进行量测，估算目标的时频差和对目标测向，得到量测矢量；
所述多目标滤波模块，用于对所述多目标量测模块输出的量测矢量进行GM-UKF-PHD滤波，得到不同权重的运动辐射源目标状态；
所述多目标跟踪状态提取模块，用于根据设置的权重阈值，停止多目标滤波模块的迭代运算，提取多运动辐射源的跟踪状态。
8.根据权利要求7所述的地面同频多运动辐射源跟踪系统，其特征在于，所述多目标滤波模块包括初始值估计模块、滤波模块和更新模块；
所述初始值估计模块连接滤波模块，为滤波模块的滤波提供初始输入信息；
所述滤波模块连接所述初始值估计模和更新模块，接收所述初始值估计模块的初始输入信息，开始进行GM-UKF-PHD滤波；将每次滤波的结果存储到更新模块；在初始滤波处理后，接收更新模块输出的上一次滤波结果，进行迭代GM-UKF-PHD滤波；
所述更新模块的输入与输出与预测模块连接，所述更新模块存储上一次滤波模块的滤波结果，并将存储的滤波结果输出到滤波模块进行当前的迭代滤波。
9.根据权利要求7或8所述的地面同频多运动辐射源跟踪系统，其特征在于，当多目标滤波模块输出的目标滤波估计结果对应的权重与权重阈值最接近时，输出跟踪结果。
CN201810228309.4A 2018-03-16 2018-03-16 基于时频差与测向的地面同频多运动辐射源跟踪方法及系统 Active CN108490465B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810228309.4A CN108490465B (zh) 2018-03-16 2018-03-16 基于时频差与测向的地面同频多运动辐射源跟踪方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810228309.4A CN108490465B (zh) 2018-03-16 2018-03-16 基于时频差与测向的地面同频多运动辐射源跟踪方法及系统

Publications (2)

Publication Number Publication Date
CN108490465A CN108490465A (zh) 2018-09-04
CN108490465B true CN108490465B (zh) 2020-06-16

Family

ID=63318553

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810228309.4A Active CN108490465B (zh) 2018-03-16 2018-03-16 基于时频差与测向的地面同频多运动辐射源跟踪方法及系统

Country Status (1)

Country Link
CN (1) CN108490465B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109255353B (zh) * 2018-09-12 2022-06-28 首都师范大学 一种动目标检测方法、装置、电子设备及存储介质
CN110133609B (zh) * 2019-04-04 2020-10-16 杭州电子科技大学 一种机载外辐射源状态未知下移动接收源路径优化方法
CN110471029B (zh) * 2019-08-16 2023-11-24 天津大学 一种基于扩展卡尔曼滤波的单站无源定位方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103744052A (zh) * 2013-12-23 2014-04-23 中国电子科技集团公司第三十六研究所 一种应用于空中目标定位的双星测时差测向方法及装置
CN105044667A (zh) * 2015-07-29 2015-11-11 中国电子科技集团公司第三十六研究所 一种运动目标的双星跟踪方法、装置和系统
CN105699964A (zh) * 2016-02-29 2016-06-22 无锡南理工科技发展有限公司 一种基于汽车防撞雷达的道路多目标跟踪方法
CN107526073A (zh) * 2017-08-22 2017-12-29 哈尔滨工程大学 一种运动多站无源时差频差联合定位方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8866672B2 (en) * 2012-04-05 2014-10-21 L-3 Communications Integrated Systems Lp Cooperative systems and methods for TDOA-based emitter location

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103744052A (zh) * 2013-12-23 2014-04-23 中国电子科技集团公司第三十六研究所 一种应用于空中目标定位的双星测时差测向方法及装置
CN105044667A (zh) * 2015-07-29 2015-11-11 中国电子科技集团公司第三十六研究所 一种运动目标的双星跟踪方法、装置和系统
CN105699964A (zh) * 2016-02-29 2016-06-22 无锡南理工科技发展有限公司 一种基于汽车防撞雷达的道路多目标跟踪方法
CN107526073A (zh) * 2017-08-22 2017-12-29 哈尔滨工程大学 一种运动多站无源时差频差联合定位方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种基于TDOA、FDOA、DOA 量测与粒子滤波的动目标双星跟踪方法;尤明懿;《通信对抗》;20170315;第36卷(第1期);第29-32页 *
基于IMM的UK-GMPHDF算法在多机动目标跟踪中的应用;郝燕玲 等;《系统工程理论与实践》;20111115;第31卷(第11期);第2225-2233页 *

Also Published As

Publication number Publication date
CN108490465A (zh) 2018-09-04

Similar Documents

Publication Publication Date Title
CN110503071B (zh) 基于变分贝叶斯标签多伯努利叠加模型的多目标跟踪方法
CN108490465B (zh) 基于时频差与测向的地面同频多运动辐射源跟踪方法及系统
Niedfeldt et al. Recursive RANSAC: Multiple signal estimation with outliers
CN109581281B (zh) 基于到达时间差和到达频率差的移动目标定位方法
CN107346020B (zh) 一种用于异步多基地雷达系统的分布式批估计融合方法
CN105182311A (zh) 全向雷达数据处理方法及系统
WO2016112758A1 (zh) 终端的定位方法及装置
CN115204212A (zh) 一种基于stm-pmbm滤波算法的多目标跟踪方法
Anderson et al. Track association for over-the-horizon radar with a statistical ionospheric model
CN108519595A (zh) 联合多传感器配准与多目标跟踪方法
CN108469627B (zh) 基于时频差的地面同频多静止辐射源定位方法及系统
CN105866757B (zh) 基于时变定位误差的外辐射源雷达目标航迹融合方法
CN106500588A (zh) 一种相位干涉仪通道间相位差噪声协方差估计方法
Kim et al. Joint CKF-PHD filter and map fusion for 5G multi-cell SLAM
EP3792651B1 (en) Path estimating device and portable information terminal
CN206292389U (zh) 一种移动目标的定位系统
Son et al. SIMM method based on acceleration extraction for nonlinear maneuvering target tracking
US11408964B2 (en) Location estimating apparatus, location estimating method and program storing recording medium, and location estimating system
Li et al. Simultaneous target and multipath positioning via multi-hypothesis single-cluster PHD filtering
CN106066471A (zh) 一种移动目标的定位方法及系统
Baek et al. Grid-based target estimation scheme for passive coherent location system
Crouse et al. Sensor bias estimation in the presence of data association uncertainty
CN108828583A (zh) 一种基于模糊c均值点迹分簇方法
Wang et al. WiFi Realtime Localization for Smartphones by MIMO-OFDM Tracking Algorithm
Huang et al. Helicopter tracking and classification with multiple interacting multiple model estimator with out-of-sequence acoustic and EO measurements

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant