CN109917331B - 基于稀疏度阶数优化的未知杂波无源协同定位方法 - Google Patents
基于稀疏度阶数优化的未知杂波无源协同定位方法 Download PDFInfo
- Publication number
- CN109917331B CN109917331B CN201910026927.5A CN201910026927A CN109917331B CN 109917331 B CN109917331 B CN 109917331B CN 201910026927 A CN201910026927 A CN 201910026927A CN 109917331 B CN109917331 B CN 109917331B
- Authority
- CN
- China
- Prior art keywords
- clutter
- sparsity
- gaussian
- measurement
- order
- 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
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提出一种基于稀疏度阶数优化的未知杂波无源协同定位方法。该方法将稀疏度阶数优化的杂波密度估计器嵌入至高斯混合概率假设密度滤波器估计多目标状态和个数。首先通过门限技术和目标状态反馈剔除目标测量,从而获取杂波测量,其次从杂波测量中选取样本,通过GA‑SVR来拟合样本,最后通过梯度法求拟合曲线的极值点,极值点的横坐标向下取整即为最优化稀疏度阶数n。本发明能有效地提升杂波未知下的无源协同定位多目标跟踪性能,解决多目标跟踪的难题。
Description
技术领域
本发明属于目标检测跟踪领域,涉及一种基于稀疏度阶数优化的未知杂波无源协同定位方法。
背景技术
无源协同定位(Passive Coherent Location,PCL)系统是由接收站与外辐射源信号组成的系统,通常将电视或广播信号作为外辐射源信号,外辐射源的直达波与经目标反射后被接收的反射波相干处理后可对目标进行定位。由于系统本身不发射电磁能量,具有较好的隐蔽性,同时PCL具有有效对付隐身目标、低成本等优势。将传统的数据关联滤波器应用于PCL系统具有高复杂度的特点,基于随机有限集理论的概率假设密度滤波器无需复杂的数据关联,可对目标状态和个数进行估计。PCL系统环境复杂多变,尤其在海岸线附近更为明显,在复杂的环境下杂波分布将不再均匀,当假设的杂波分布模型与实际的杂波分布不一致时将会影响目标跟踪的性能。
传统的空间杂波稀疏度估计(Spatial Clutter Sparsity Estimation,SCSE)法可直接利用每帧测量数据中待计算杂波密度的测量点到另一个测量点之间的欧氏距离估计该测量位置的杂波密度,实时性好且易嵌入到不同滤波器,但存在稀疏度阶数n的选取问题,通常取n=1或n=2。本发明针对n稀疏度阶数只能取固定值,无法自适应选取n稀疏度阶数,提出基于稀疏度阶数优化的未知杂波无源协同定位方法,该方法将稀疏度阶数优化的杂波密度估计器嵌入至高斯混合概率假设密度滤波器估计多目标状态和个数。首先通过门限技术和目标状态反馈剔除目标测量,从而获取杂波测量,其次从杂波测量中选取样本,通过基于遗传算法的支持向量回归机(Genetic Algorithm-Support Vector Regression,GA-SVR)来拟合样本,然后通过梯度法求拟合曲线的极值点,极值点的横坐标向下取整即为最优化稀疏度阶数,将最优化稀疏度阶数应用杂波稀疏度中可求得杂波密度,实现在杂波分布未知下对目标进行精确地跟踪。
发明内容
本发明的目的是针对现有技术无法自适应选取稀疏度阶数,提供实现本发明的技术解决方案为:一种基于稀疏度阶数优化的高斯混合概率假设密度无源协同定位方法。
本发明方法包括以下步骤:
步骤1、利用高斯混合概率假设密度滤波器的预测方程进行多目标强度预测。根据k-1时刻多目标的后验强度的高斯混合,通过预测步骤得到k时刻预测强度的高斯混合。
步骤3、获取样本点。从杂波测量中选取样本(n,V(r(n)(y)))。
步骤4、GA-SVR拟合样本。使用GA求解θ=[C,ε,σ]T,并通过SVR来拟合样本。
步骤5、极值点的求取。利用梯度法求取决策函数的极值点,极值点的横坐标向下取整即为最优化稀疏度阶数n。
步骤6、多目标强度更新与高斯分量的修剪和合并。首先,利用估计的第k帧杂波密度ck(z),结合高斯混合概率假设密度滤波器对预测强度函数中的高斯分量进行更新,得到多目标后验强度。然后,对后验强度函数中的高斯分量进行低权重修剪和邻近分量合并的操作,得到k时刻后验强度的高斯混合
步骤7、利用多目标后验强度进行多目标状态提取。
本发明的关键技术在于针对n稀疏度阶数只能取固定值,无法自适应选取n稀疏度阶数,提出基于稀疏度阶数优化的杂波密度估计算法。首先通过门限技术和目标状态反馈剔除目标测量,从而获取杂波测量,其次从杂波测量中选取样本,通过GA-SVR来拟合样本,在杂波均匀分布时,样本点的斜率为常量,在杂波非均匀时,样本点的斜率会发生变化,通过梯度法求拟合曲线的极值点,极值点的横坐标向下取整即为最优化稀疏度阶数n。本发明能够有效提高未知杂波无源协同定位系统的多目标跟踪性能,解决未知杂波下的多目标跟踪难题。
附图说明
图1为本发明流程图。
图2为在杂波未知时的多目标跟踪效果图。
图3为传统的SCSE的杂波密度均方根误差和本发明所提算法的杂波密度均方根误差对比图。
具体实施方式
下面结合具体实施例对本发明作进一步的分析。
Nk个未知且Nk≥0的目标出现在传感器的探测范围内,接收站位置记为XRr=[xRr,yRr],考虑到外辐射源为机载外辐射源,第k时刻的发射站位置记为Xk,T=[xk,Ts,yk,Ts]。考虑单目标t(0≤t≤Nk)在k时刻的状态Xk,t为其中为k时刻目标t在直角坐标系下所处的位置,为k时刻目标t的瞬时速度,状态转换模型如下:
Xk+1.t=F·Xk,t+vk (1)
其中rk表示第k时刻的距离差,θk表示第k时刻回波路径相对于直达波路径的方位角。通常假设每时刻最多包含一个源于目标的测量,且不同时刻之间的测量相互独立。Υk为k时刻的杂波测量,测量噪声εk~N(0,R),测量噪声协方差σr为距离测量标准差,σθ为方位角标准差。||·||为欧几里得范数。
要解决的问题是在非均匀的杂波环境中估计每个时刻的目标状态,通常假设杂波的空间分布服从均匀分布,但实际环境中杂波是非均匀分布的,杂波密度估计的准确性将影响到目标状态的估计。
步骤1、利用高斯混合概率假设密度滤波器的预测方程进行多目标强度预测。
如果已知k-1时刻多目标的后验强度的高斯混合为
其中Jk-1为k-1时刻的高斯分量个数。
则k时刻预测强度的高斯混合为
步骤2、杂波测量的选取。
考虑到机载无源协同定位系统测量模型的非线性,可以利用高斯分量的均值和协方差通过无味变换得到第i个预测高斯分量的Sigma点集及其权重i=1,...,Jk|k-1。Sigma点的预测测量值个数η=0,...,L。然后,针对第i个存活高斯分量的跟踪门计算如下:
步骤3、样本选取:从步骤2中的杂波集获取样本集。
定义数集其中n为杂波稀疏度的近邻阶数,nmax为当前观测空间内的杂波测量个数。定义采样点为(n,V(r(n)(y))),其中y为维数为l的测量点,r(n)(y)为距离测量点y的第n近邻距离,V(r(n)(y))为以y为圆心,r(n)(y)为半径的超球体体积,公式如下:
V(r(n)(y))=Clr(n)(y)l (11)
其中π≈3.14,Γ(·)为Gamma函数。
步骤4、GA-SVR曲线拟合:对样本集进行GA-SVR拟合。
决策函数f(x)如下:
其中x∈[1,nmax]为连续状态量,当i=j时,K(i,j)可表示为K(i,i)。
定义推广误差Er如下:
步骤5、极值点的求取:利用梯度法公式(19)求取上述步骤4拟合曲线的极值点,极值点的横坐标向下取整即为最优化稀疏度阶数n,将自适应求取的稀疏度阶数n代入公式(13)(14)得到基于稀疏度阶数优化的杂波密度估计ck(z)。
极小值点的横坐标和极大值点的横坐标如下:
步骤6、多目标强度更新与高斯分量的修剪和合并。
利用估计的第k帧杂波密度ck(z),结合高斯混合概率假设密度滤波器对预测强度函数中的高斯分量进行更新,得到多目标后验强度:
其中对单个高斯分量使用了无味卡尔曼滤波的更新方程。具体步骤如下:
图3为传统的SCSE的杂波密度均方根误差和本发明所提算法的杂波密度均方根误差对比图,从图3可知在16s后,本发明算法相比于传统的SCSE算法更加符合真实的杂波分布,这是因为本发明算法能够自适应寻找最优稀疏度阶数,使得杂波密度更加逼近真实的杂波密度。而SCSE算法由于无法知道杂波的均匀程度,只能取固定的稀疏度阶数,在16s后所估计的杂波密度将偏离真实值,从而影响多目标跟踪的性能。图2为在杂波未知时的多目标跟踪效果图,可见使用本发明所提算法后可提升多目标跟踪的性能。
Claims (7)
1.基于稀疏度阶数优化的未知杂波无源协同定位方法,其特征在于该方法包括以下步骤:
步骤1、利用高斯混合概率假设密度滤波器的预测方程进行多目标强度预测;
步骤2、杂波测量的选取
步骤3、样本选取:从步骤2中的杂波集获取样本集(n,V(r(n)(y)));n为杂波稀疏度的近邻阶数,y为维数为l的测量点,r(n)(y)为距离测量点y的第n近邻距离,V(r(n)(y))为以y为圆心,r(n)(y)为半径的超球体体积;
步骤4、GA-SVR拟合样本:使用GA求解θ=[C,ε,σ]T,并通过SVR来拟合样本;C为惩罚因子,ε为不敏感系数,σ为高斯核带宽;
步骤5、极值点的求取:利用梯度法求取决策函数的极值点,极值点的横坐标向下取整即为最优化稀疏度阶数n,得到基于稀疏度阶数优化的杂波密度估计ck(z);
步骤6、多目标强度更新与高斯分量的修剪和合并;
首先,利用估计的第k帧杂波密度ck(z),结合高斯混合概率假设密度滤波器对预测强度函数中的高斯分量进行更新,得到多目标后验强度;然后,对后验强度函数中的高斯分量进行低权重修剪和邻近分量合并的操作,得到k时刻后验强度的高斯混合Jk为k时刻的高斯分量个数,分别为k时刻的第i个高斯分量的权重、均值和协方差;
2.如权利要求1所述的基于稀疏度阶数优化的未知杂波无源协同定位方法,其特征在于步骤1具体是:
如果已知k-1时刻多目标的后验强度的高斯混合为
则k时刻预测强度的高斯混合为
其中为新生目标强度函数; Jγ,k为k时刻新生分量的个数,分别为k时刻的第i个新生分量的权重、均值和协方差;xk|k-1表示k-1时刻在直角坐标系下所处的x轴预测位置;表示k-1时刻的第i个高斯分量的预测均值和预测协方差;
可得预测强度:
4.如权利要求3所述的基于稀疏度阶数优化的未知杂波无源协同定位方法,其特征在于步骤3具体是:
定义数集其中n为杂波稀疏度的近邻阶数,nmax为当前观测空间内的杂波测量个数;定义采样点为(n,V(r(n)(y))),其中y为维数为l的测量点,r(n)(y)为距离测量点y的第n近邻距离,V(r(n)(y))为以y为圆心,r(n)(y)为半径的超球体体积,公式如下:
V(r(n)(y))=Clr(n)(y)l (11)
其中π≈3.14,Γ(·)为Gamma函数;
决策函数f(x)如下:
其中x∈[1,nmax]为连续状态量,当i=j时,K(i,j)可表示为K(i,i);
定义推广误差Er如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910026927.5A CN109917331B (zh) | 2019-01-11 | 2019-01-11 | 基于稀疏度阶数优化的未知杂波无源协同定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910026927.5A CN109917331B (zh) | 2019-01-11 | 2019-01-11 | 基于稀疏度阶数优化的未知杂波无源协同定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109917331A CN109917331A (zh) | 2019-06-21 |
CN109917331B true CN109917331B (zh) | 2021-02-09 |
Family
ID=66960256
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910026927.5A Active CN109917331B (zh) | 2019-01-11 | 2019-01-11 | 基于稀疏度阶数优化的未知杂波无源协同定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109917331B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110596691B (zh) * | 2019-08-30 | 2021-10-22 | 杭州电子科技大学 | 一种考虑地球曲率的双基站三维无源定位方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104076348A (zh) * | 2014-07-09 | 2014-10-01 | 中国船舶重工集团公司第七二四研究所 | 一种雷达超视距基线无源协同定位方法 |
US20140372073A1 (en) * | 2011-05-04 | 2014-12-18 | Jacques Georgy | Two-stage filtering based method for multiple target tracking |
CN107064865A (zh) * | 2017-04-07 | 2017-08-18 | 杭州电子科技大学 | 基于深度聚类的极坐标动态规划无源协同定位方法 |
CN107102293A (zh) * | 2017-04-25 | 2017-08-29 | 杭州电子科技大学 | 基于滑窗累积密度估计的未知杂波无源协同定位方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106162162B (zh) * | 2016-08-01 | 2017-10-31 | 宁波大学 | 一种基于稀疏表示的重定位图像质量客观评价方法 |
CN106526559B (zh) * | 2016-09-30 | 2019-04-09 | 杭州电子科技大学 | 一种基于bcd-vsmm机动目标无源协同定位方法 |
CN109031229B (zh) * | 2017-11-27 | 2022-10-11 | 电子科技大学 | 一种杂波环境下目标跟踪的概率假设密度方法 |
-
2019
- 2019-01-11 CN CN201910026927.5A patent/CN109917331B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140372073A1 (en) * | 2011-05-04 | 2014-12-18 | Jacques Georgy | Two-stage filtering based method for multiple target tracking |
CN104076348A (zh) * | 2014-07-09 | 2014-10-01 | 中国船舶重工集团公司第七二四研究所 | 一种雷达超视距基线无源协同定位方法 |
CN107064865A (zh) * | 2017-04-07 | 2017-08-18 | 杭州电子科技大学 | 基于深度聚类的极坐标动态规划无源协同定位方法 |
CN107102293A (zh) * | 2017-04-25 | 2017-08-29 | 杭州电子科技大学 | 基于滑窗累积密度估计的未知杂波无源协同定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109917331A (zh) | 2019-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5611157B2 (ja) | 時空間適応処理を用いてレーダー信号におけるターゲットを検出する方法 | |
WO2003021289A1 (en) | Radar system and method including superresolution raid counting | |
CN110221241A (zh) | 一种基于rbf神经网络的低仰角doa估计方法 | |
CN110389327A (zh) | 接收站位置误差下多站多外辐射源雷达双基距定位方法 | |
CN108089166A (zh) | 一种基于毫米波测云雷达的自适应频域检测方法 | |
CN112597820A (zh) | 一种基于雷达信号分选的目标聚类方法 | |
Li et al. | Adaptive two-step Bayesian MIMO detectors in compound-Gaussian clutter | |
Lapierre et al. | New methods for handling the range dependence of the clutter spectrum in non-sidelooking monostatic STAP radars | |
CN107102293B (zh) | 基于滑窗累积密度估计的未知杂波无源协同定位方法 | |
CN108398678A (zh) | 一种集中式mimo雷达快速自适应目标跟踪方法 | |
CN109917331B (zh) | 基于稀疏度阶数优化的未知杂波无源协同定位方法 | |
CN112612006A (zh) | 基于深度学习的机载雷达非均匀杂波抑制方法 | |
CN111273269B (zh) | 基于ipso-bp的频率分集阵列的雷达目标定位方法 | |
Tian et al. | Performance evaluation of deception against synthetic aperture radar based on multifeature fusion | |
CN110308442B (zh) | 强杂波环境下相控阵雷达的gm-phd目标跟踪方法 | |
Zhang et al. | Joint optimisation of transmit waveform and receive filter for cognitive radar | |
CN117169818B (zh) | 一种面向海面机动目标跟踪的雷达波形设计方法 | |
CN108152796B (zh) | 一种基于灰色卡尔曼滤波的主瓣移动干扰消除方法 | |
Zhou et al. | Sea clutter distribution modeling: A kernel density estimation approach | |
Yang et al. | Multi-object Bayesian filters with amplitude information in clutter background | |
CN113093174B (zh) | 基于phd滤波雷达起伏微弱多目标的检测前跟踪方法 | |
CN114167359A (zh) | 一种弱小目标自适应关联滤波方法、系统及存储介质 | |
CN113126086A (zh) | 一种基于状态预测积累的生命探测雷达弱目标检测方法 | |
Song et al. | Box-particle cardinality balanced multi-target multi-Bernoulli filter | |
CN112333629A (zh) | 一种互耦未知条件下分布式阵列多目标定位方法 |
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 |