CN107544063B - 一种雷达跟踪状态下的目标rcs的预测方法 - Google Patents
一种雷达跟踪状态下的目标rcs的预测方法 Download PDFInfo
- Publication number
- CN107544063B CN107544063B CN201710670059.5A CN201710670059A CN107544063B CN 107544063 B CN107544063 B CN 107544063B CN 201710670059 A CN201710670059 A CN 201710670059A CN 107544063 B CN107544063 B CN 107544063B
- Authority
- CN
- China
- Prior art keywords
- target
- time
- rcs
- value
- moment
- 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
本发明属于雷达技术领域,公开了一种雷达跟踪状态下的目标RCS的预测方法,用于预测雷达跟踪过程中目标下一时刻的RCS值,包括:设定目标雷达截面积RCS的预测滤波器的阶数为M,根据目标跟踪算法得到目标自NP‑M+1时刻至NP时刻目标的距离与速度,以及自NP‑M+1时刻至NP时刻雷达对目标的观测角度,求得在自NP‑M+1时刻至NP时刻的目标距离以及运动模型,得到目标在第NP+1时刻的目标速度与目标距离的夹角预测值,由目标的观测角度求得目标RCS值的自相关矩阵和目标RCS值的互相关列向量,通过多普勒向量调节目标观测分量的协方差矩阵,求得目标RCS的预测滤波器系数,得到在第NP+1时刻的目标RCS的预测值。
Description
技术领域
本发明属于雷达技术领域,尤其涉及一种雷达跟踪状态下的目标RCS的预测方法,可用于预测雷达跟踪过程中目标下一时刻的雷达散射截面积RCS的值。
背景技术
目标雷达截面积(RCS cross section,RCS)是用来衡量目标散射特性的一个参数,一般用后向散射能量的强度来定义目标的RCS。目标RCS主要与目标的结构和表面介质、雷达频率、极化方式与目标姿态角等因素有关。
由雷达方程可知,目标RCS直接影响目标回波功率的大小,从而目标RCS直接影响目标回波信号的信噪比。而目标信噪比的大小对多站雷达中的资源调度以及相控阵雷达中的功率分配都具有重要影响。当雷达跟踪目标时,如果已知目标RCS和目标距离,则可以估计得到目标回波信号的信噪比。在目标信噪比已知的条件下,多站雷达可以调用高信噪比的雷达对目标进行观测,使资源调度更为合理;而相控阵雷达可以根据目标的信噪比分配所需的发射功率,使雷达功率得到更有效的利用。因此,雷达跟踪状态下,研究目标的RCS的预测具有重要意义。
目前,研究雷达跟踪状态下目标RCS预测的文献资料还较少,并且现有的信噪比的预测方法中考虑到RCS变化的也很少,大多数只是假设目标的RCS值保持不变。在实际应用中,随着目标与雷达相对位置的变化,目标的RCS值是不断变化的。
发明内容
针对上述现有技术的问题,本发明的目的在于提供一种雷达跟踪状态下的目标RCS的预测方法,实现了目标在下一跟踪时刻RCS值的预测。
为达到上述目的,本发明采用如下技术方案予以实现。
一种雷达跟踪状态下的目标RCS的预测方法,所述方法包括如下步骤:
步骤1,确定目标RCS的预测滤波器的阶数为M,并获取雷达接收到的从NP-M+1时刻到NP时刻的历史目标回波幅度观测向量其中,ξk表示第k时刻的目标回波幅度观测值,且NP>M,NP、k、M分别为自然数,k∈(Np-M+1,...,Np),k的初值为Np-M+1;
采用目标跟踪算法分别得到目标从NP-M+1时刻到NP时刻中,每个时刻的目标距离和目标速度,并设定初始NP-M+1时刻的目标观测角度;
步骤2,假设目标的飞行与目标姿态角度不变的假设下,根据k时刻的目标距离、目标速度以及目标观测角度,估计k+1时刻的目标观测角度;令k的值依次取Np-M+1,...,Np,直到得到当前NP+1时刻的目标观测角度;
步骤3,根据NP-M+1时刻到NP时刻的目标观测角度,得到目标RCS的自相关矩阵Φ,所述目标RCS的自相关矩阵Φ是M×M维的矩阵;
步骤4,根据NP-M+1时刻到当前NP+1时刻的目标观测角度,得到目标RCS的互相关向量η,所述目标RCS的互相关向量η是M×1维的列向量;
步骤5,根据目标RCS的自相关矩阵Φ,得到历史目标回波幅度观测向量ξ的协方差矩阵;
步骤7,根据所述预测滤波器的最优权值,所述历史目标回波幅度观测向量ξ,计算得到当前NP+1时刻的目标回波幅度预测值,从而根据所述当前NP+1时刻的目标回波幅度预测值求得当前NP+1时刻的目标RCS预测值。
本发明技术方案的特点和进一步的改进为:
(1)步骤2中,根据k时刻的目标距离、目标速度以及目标观测角度,确定k+1时刻的目标观测角度αk+1,具体采用如下公式进行计算:
其中,rk表示k时刻的目标距离,vk表示k时刻的目标速度,αk表示k时刻的目标观测角度,T表示脉冲重复间隔;
令k的值依次取Np-M+1,...,Np,直到得到当前NP+1时刻的目标观测角度。
(2)步骤3中,根据NP-M+1时刻到NP时刻的目标观测角度,得到目标RCS的自相关矩阵Φ为:
其中,ρ(·)表示相关系数函数,其表达式为:ρ(u1)=jinc(4fcd sin(u1/2)/c),u1为相关系数函数的输入,fc表示雷达载频,d表示目标尺寸的平均值,c=3×108m/s表示光速,jinc(·)表示jinc函数,jinc函数定义为:J1(·)表示第一类贝塞尔函数,u2为jinc函数的输入,αk表示k时刻的目标观测角度,k∈(Np-M+1,...,Np)。
(3)步骤4中,根据NP-M+1时刻到NP+1时刻的目标观测角度,得到目标RCS的互相关向量η为:
(4)步骤5具体包括:
根据目标RCS的自相关矩阵Φ,得到目标回波幅度观测向量的协方差矩阵C:
C=E(ξξH)=σs 2diag(ad)Φdiag(a* d);
其中,E(·)表示求期望,ξ表示NP-M+1时刻到NP时刻的历史目标回波幅度观测向量,σs 2表示目标回波幅度观测值的方差,ad表示NP-M+1时刻到NP时刻目标的多普勒导向向量,a* d表示ad的共轭,(·)H表示共轭转置,diag(·)表示向量对角化;
(5)步骤6具体包括:
其中,Rz表示干扰信号的协方差矩阵,⊙表示点乘;
(6)步骤7具体包括:
根据所述预测滤波器的最优权值wopt=[w1,…wi,…,wM]T,所述历史目标回波幅度观测向量计算得到当前NP+1时刻的目标回波幅度预测值从而根据所述当前NP+1时刻的目标回波幅度预测值求得当前NP+1时刻的目标RCS预测值其中,∑表示求和运算,wi为预测滤波器的最优权值wopt的第i个元素,i∈(1,2,...,M)。
本发明通过上述技术方案,利用目标在当前时刻与下一时刻的观测角度变化,以及预测滤波器对下一时刻目标的RCS值进行预测,得到了目标下一跟踪时刻的RCS预测值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种雷达跟踪状态下的目标RCS的预测方法的流程示意图;
图2为本发明仿真实验中飞机散射点模型示意图;
图3为本发明仿真实验中目标的运动模型示意图;
图4为本发明仿真实验中目标在不同的脉冲观测频率下回波幅度真实值、观测值与预测值的对比示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供一种雷达跟踪状态下的目标RCS的预测方法,如图1所示,所述方法包括:
步骤1,设定目标的雷达散射截面积RCS的预测滤波器的阶数为M。雷达接收到目标从NP-M+1到NP时刻的历史回波幅度观测向量ξk表示第k时刻的观测数据,(·)T表示取转置操作,NP、k为自然数。设定当前时刻的目标观测值为本实例假设所有的观测值服从均值为零,方差为σs 2的复高斯分布。
预测器阶数M的意义,利用当前时刻之前M个时刻的观测数据对当前时刻的RCS值进行预测。M的取值不超过允许的最大计算量。
步骤2,假设目标的飞行与目标姿态角度不变的假设下,利用目标跟踪算法可得到目标在任意k时刻的距离信息rk与速度信息vk,根据k时刻的目标距离rk、k时刻的目标速度vk以及k时刻的目标观测角度αk,按照目标的运动模型可以计算得到k+1时刻目标的观测角启
所述的目标跟踪算法包括,基于卡尔曼滤波的目标跟踪,基于粒子滤波的目标跟踪,基于维纳滤波的目标跟踪。本实例选用但不限于基于卡尔曼滤波的目标跟踪算法。
步骤3,(3.1)根据自NP-M+1到NP时刻的目标观测角度,可求出目标RCS的自相关矩阵Φ,自相关矩阵是一个M×M的矩阵,其形式为
其中ρ表示相关系数函数,相关系数函数的表达式为:
ρ(u1)=jinc(4fcd sin(u1/2)/c)
其中,u1表示相关系数函数的输入变量,fc表示雷达载频,d表示目标尺寸的平均值,c=3×108m/s表示光速,jinc(·)表示jinc函数,jinc函数的定义为:
J1(·)表示第一类贝塞尔函数,u2为jinc函数的输入变量;
目标的RCS值只与观测角度有关,当两个观测角度间隔很小时,这两个观测角度对应的目标RCS值具有较强的相关性。根据目标RCS值得相关性,可以使用预测器完成对目标RCS的预测;
利用相关性预测的前提条件是历史回波幅度观测向量与当前时刻的观测值是不统计独立,即它们是相关的,判断具有相关性的条件是:
C=E(ξξH)=σs 2diag(ad)Φdiag(a* d)
步骤5,(5.1)假设目标的空间响应变化是一个平稳随机过程,仅随着信号的起始角度变化。根据现有回波观测数据ξ,考虑将其以线性组合方式将预测器输出表示出来,形式为:选定度量预测器的性能准则为最小均方误差准则,目标函数可以表示为:
其中,w表示预测器的权向量,z表示不同时刻干扰分量组成的向量,将上式展开可写为
Rz=M×M,其中Rz表示干扰信号的协方差矩阵,其中干扰与信号独立。
(5.2)通过对上式w求导可解得最小均方误差对应的最佳权值形式为:
其中⊙表示点乘;
(5.3)在实际的工程应用中,不同的脉冲重复周期内干扰是统计独立的,因此Rz为一个对角矩阵,最佳权值又可以表示为:
其中,wi为预测滤波器系数w的第i个元素,i=1,2,...,M,M表示预测滤波器阶数。
如果目标RCS的相关性很强,可能会造成自相关矩阵R非奇异,导致矩阵求逆时出现错误,在这种情况下,可以使用广义矩阵求逆运算来代替矩阵求逆运算。
步骤6,根据预测器系数w以及自NP-M+1时刻至NP时刻目标回波幅度观测值,得到当前时刻的目标复幅度的预测值求取第NP+1时刻的目标RCS的预测值为其中,∑表示求和运算,i=1,2,...,M,wi为预测滤波器系数w的第i个元素;
下面结合仿真实验对本发明的效果做进一步说明:
设置仿真参数:采用一个位于坐标原点的2D雷达,设雷达载频fc=1GHz,雷达测量参数为目标的距离和方位角。目标为长度为10m的飞机目标,RCS散射特性服从散射点模型,如图2所示,各个散射点的散射系数服从均值为0,方差为1/NS的复高斯分布,且不同散射点的散射系数相互独立;其中NS表示散射点个数。通过计算可得跟踪角度差最大取值φmax=0.86(rad)。
根据仿真参数设置和目标运动模型,计算得到雷达与目标观测角度的范围是[-0.43,0.43](rad),取出[-0.43,0.43](rad)范围内的目标复幅度的真实值用于仿真,将[-0.43,0.43](rad)角度内的目标RCS的真实值画成二维曲线,并将通过本发明方法得到的预测值与加入噪声之后的真实值在不同的脉冲观测频率下对比,如图4所示(横坐标为雷达对目标的挂侧角度,单位为度,纵坐标为目标回波幅度,单位为分贝),图(4a)中,脉冲观测频率PVF=0.2Hz,图(4b)中,脉冲观测频率PVF=10Hz,图(4c)中,脉冲观测频率PVF=1Hz;由仿真结果可看出,目标RCS预测值的整体变化趋势与目标RCS真实值相同,目标RCS预测值的大小在真实值附近随机跳动,这主要是由噪声引起的,对比图4中的三个仿真结果可知,本发明方法可以实现对目标RCS的预测,能够较好的预测出目标RCS的变化趋势。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (6)
1.一种雷达跟踪状态下的目标RCS的预测方法,其特征在于,所述方法包括如下步骤:
步骤1,确定目标RCS的预测滤波器的阶数为M,并获取雷达接收到的从NP-M+1时刻到NP时刻的历史目标回波幅度观测向量其中,ξk表示第k时刻的目标回波幅度观测值,且NP>M,NP、k、M分别为自然数,k∈(Np-M+1,...,Np),k的初值为Np-M+1;
采用目标跟踪算法分别得到目标从NP-M+1时刻到NP时刻中,每个时刻的目标距离和目标速度,并设定初始NP-M+1时刻的目标观测角度;
步骤2,根据k时刻的目标距离、目标速度以及目标观测角度,确定k+1时刻的目标观测角度;令k的值依次取Np-M+1,...,Np,直到得到当前NP+1时刻的目标观测角度;
步骤3,根据NP-M+1时刻到NP时刻的目标观测角度,得到目标RCS的自相关矩阵Φ,所述目标RCS的自相关矩阵Φ是M×M维的矩阵;
步骤4,根据NP-M+1时刻到当前NP+1时刻的目标观测角度,得到目标RCS的互相关向量η,所述目标RCS的互相关向量η是M×1维的列向量;
步骤5,根据目标RCS的自相关矩阵Φ,得到历史目标回波幅度观测向量ξ的协方差矩阵;
求解所述关于预测滤波器权值的目标函数,得到预测滤波器的最优权值;
步骤7,根据所述预测滤波器的最优权值,所述历史目标回波幅度观测向量ξ,计算得到当前NP+1时刻的目标回波幅度预测值,从而根据所述当前NP+1时刻的目标回波幅度预测值求得当前NP+1时刻的目标RCS预测值;具体为:
5.根据权利要求4所述的一种雷达跟踪状态下的目标RCS的预测方法,其特征在于,步骤5具体包括:
根据目标RCS的自相关矩阵Φ,得到目标回波幅度观测向量的协方差矩阵C:
C=E(ξξH)=σs 2diag(ad)Φdiag(a* d);
其中,E(·)表示求期望,ξ表示NP-M+1时刻到NP时刻的历史目标回波幅度观测向量,σs 2表示目标回波幅度观测值的方差,ad表示NP-M+1时刻到NP时刻目标的多普勒导向向量,a* d表示ad的共轭,(·)H表示共轭转置,diag(·)表示向量对角化;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710670059.5A CN107544063B (zh) | 2017-08-08 | 2017-08-08 | 一种雷达跟踪状态下的目标rcs的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710670059.5A CN107544063B (zh) | 2017-08-08 | 2017-08-08 | 一种雷达跟踪状态下的目标rcs的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107544063A CN107544063A (zh) | 2018-01-05 |
CN107544063B true CN107544063B (zh) | 2020-05-01 |
Family
ID=60970484
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710670059.5A Active CN107544063B (zh) | 2017-08-08 | 2017-08-08 | 一种雷达跟踪状态下的目标rcs的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107544063B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108983195B (zh) * | 2018-08-17 | 2022-08-30 | 桂林电子科技大学 | 基于子阵列自适应成像的目标雷达散射截面积测量方法 |
CN111123231B (zh) * | 2018-11-01 | 2023-05-02 | 三峡大学 | 一种风电机rcs数据压缩重构方法 |
CN109802656B (zh) * | 2019-01-25 | 2022-03-04 | 西安电子科技大学 | 基于幅度信息的卡尔曼滤波方法 |
WO2021097636A1 (zh) * | 2019-11-19 | 2021-05-27 | 深圳市大疆创新科技有限公司 | 标物雷达散射截面积确定方法、装置和存储介质 |
CN111208482B (zh) * | 2020-02-28 | 2022-02-18 | 成都汇蓉国科微***技术有限公司 | 基于距离对齐的雷达精度分析方法 |
CN115356703B (zh) * | 2022-10-17 | 2023-02-03 | 中国人民解放军国防科技大学 | 基于面元分布的粗糙目标rcs缩比测量方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103760529A (zh) * | 2013-12-06 | 2014-04-30 | 河海大学 | 一种无源探测高效级联空时自适应处理方法 |
CN104076342A (zh) * | 2014-06-25 | 2014-10-01 | 西安电子科技大学 | 一种雷达跟踪状态下预测目标rcs的方法 |
CN104991240A (zh) * | 2015-07-13 | 2015-10-21 | 西安电子科技大学 | 目标散射截面积的预测方法 |
CN105842687A (zh) * | 2016-03-21 | 2016-08-10 | 西安电子科技大学 | 基于rcs预测信息的检测跟踪一体化方法 |
CN106291487A (zh) * | 2016-08-04 | 2017-01-04 | 上海无线电设备研究所 | 一种基于agc电压和回波数据的雷达接收功率和rcs估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7750842B2 (en) * | 2008-09-18 | 2010-07-06 | Raytheon Company | Parallel processing to generate radar signatures for multiple objects |
-
2017
- 2017-08-08 CN CN201710670059.5A patent/CN107544063B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103760529A (zh) * | 2013-12-06 | 2014-04-30 | 河海大学 | 一种无源探测高效级联空时自适应处理方法 |
CN104076342A (zh) * | 2014-06-25 | 2014-10-01 | 西安电子科技大学 | 一种雷达跟踪状态下预测目标rcs的方法 |
CN104991240A (zh) * | 2015-07-13 | 2015-10-21 | 西安电子科技大学 | 目标散射截面积的预测方法 |
CN105842687A (zh) * | 2016-03-21 | 2016-08-10 | 西安电子科技大学 | 基于rcs预测信息的检测跟踪一体化方法 |
CN106291487A (zh) * | 2016-08-04 | 2017-01-04 | 上海无线电设备研究所 | 一种基于agc电压和回波数据的雷达接收功率和rcs估计方法 |
Non-Patent Citations (3)
Title |
---|
一种利用幅度信息的航迹质量评估方法;刘红亮 等;《西安电子科技大学学报(自然科学版)》;20170228;第44卷(第1期);第65-70页 * |
一种用于雷达资源管理的目标雷达截面积预测算法;秦童 等;《电子与信息学报》;20150831;第37卷(第8期);第1849-1854页 * |
基于非线性机会约束规划的多基雷达***稳健功率分配算法;严俊坤 等;《电子与信息学报》;20140331;第36卷(第3期);第509-515页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107544063A (zh) | 2018-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107544063B (zh) | 一种雷达跟踪状态下的目标rcs的预测方法 | |
CN104181517B (zh) | 用于运行多输入多输出雷达的方法 | |
Connell et al. | Development of an environmental and ship motion forecasting system | |
US10571550B2 (en) | Signal correction for environmental distortion | |
Patton et al. | Efficient design of radar waveforms for optimised detection in coloured noise | |
Hamid et al. | Performance comparison of time-domain and frequency-domain beamforming techniques for sensor array processing | |
CN106443653A (zh) | 雷达装置 | |
Su et al. | Adaptive beamforming for nonstationary HF interference cancellation in skywave over-the-horizon radar | |
CN107255805B (zh) | 基于加权最小二乘的雷达目标rcs的预测方法 | |
CN104076342B (zh) | 一种雷达跟踪状态下预测目标rcs的方法 | |
Ling et al. | Real time estimation and prediction of wave excitation forces on a heaving body | |
CN105487054B (zh) | 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 | |
CN108398669A (zh) | 一种基于无需预延迟处理的空时宽带自适应单脉冲测角方法 | |
CN107576947A (zh) | 基于时间平滑的l型阵对相干信源二维波达方向估计方法 | |
CN105116408A (zh) | 一种舰船isar图像结构特征提取方法 | |
Ram et al. | Optimization of radar parameters for maximum detection probability under generalized discrete clutter conditions using stochastic geometry | |
CN103487793B (zh) | 一种基于简正波理论的宽带混响波形仿真方法 | |
US9444558B1 (en) | Synthetic robust adaptive beamforming | |
Ji et al. | The application of JDL to suppress sea clutter for shipborne HFSWR | |
CN104035078A (zh) | 一种基于阵元阶数递推的降维空时自适应权值计算方法 | |
CN110850421A (zh) | 基于混响对称谱的空时自适应处理的水下目标检测方法 | |
Dou et al. | Multi‐mode clutter suppression of multiple‐input–multiple‐output over‐the‐horizon radar based on blind source separation | |
Yan et al. | Orthogonal waveform separation in multiple‐input and multiple‐output imaging sonar with fractional Fourier filtering | |
Zhang et al. | Explicit performance limit for joint range and direction of arrival estimation in phased-array radar sensors | |
Wang et al. | Research on compound scattering modeling and imaging methods of sea surface ship target for GEO-UAV BiSAR |
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 |