太赫兹ISAR目标微动部件分离与特征参数估计方法技术

本发明专利技术提供一种太赫兹ISAR目标微动部件分离与特征参数估计方法，旨在解决在太赫兹波段由目标部件微运动引起的ISAR成像质量降低和微动特征参数估计困难等问题，提出基于Radon变换的微动部件包络分离和参数估计的技术方案，并同时实现微动目标的ISAR成像。本发明专利技术提供的技术方案有望提高现代雷达对未知非合作目标识别和解译能力，具有实际应用前景。

Separation and characteristic parameter estimation of terahertz ISAR target micromoving parts

【技术实现步骤摘要】
太赫兹ISAR目标微动部件分离与特征参数估计方法
本专利技术属于雷达信号处理
，具体涉及的是一种针对太赫兹波段的逆合成孔径雷达(ISAR)微动目标成像与特征参数估计方法。
技术介绍
作为空间、空中及海域的重要观测手段，ISAR成像技术在对非合作军事目标的探测与监视中起着不可替代的作用。针对日益提高的精细化探测需求，ISAR正朝着高分辨成像、多维特征提取和目标识别的方向发展。在对非合作目标探测与识别时，目标的形状、结构和表面材料电磁参数等非运动特征通常对先验信息的要求较高，其实用性十分受限。由于独特的表现形式，近年来，微动在目标探测与识别方面的军事价值已引起了各国的高度关注。目标整体或部件的微运动使得传统ISAR成像方法失效，同时，微动与目标或部件结构尺寸、质量分布和受力状态等物理属性密切相关，可作为目标识别的重要特征量。以喷气式飞机、旋翼飞机或螺旋桨飞机等空中目标为例，直升机主/尾旋翼的旋转、螺旋桨飞机叶片的转动、喷气式飞机发动机压缩叶片的旋转及其引擎调制(JetEngineModulation，JEM)等都是微动源。对舰船而言，则存在机械扫描雷达天线的转动等微动形式。另外，真假导弹目标的识别在一定程度上决定了反导系统的成败，而微动是真弹头与诱饵的主要差异之一。目前，美国陆军研究实验室已利用92GHz连续毫米波雷达对Mi-24/雌鹿D型直升机进行了微多普勒测量，对直升机目标的多普勒频谱构成进行了分析并测量出了不同视角的微多普勒特征。日本东芝公司共同研发中心多媒体实验室已开发出一套在线式舰船自动识别系统，该系统可快速更新ISAR获得的运动舰船目标的距离-多普勒像，可实现对包含海浪颠簸运动的航行舰船目标的识别。美国正在研制的THAADGBRX-波段雷达能够精确测量威胁目标的微动特征，弹道导弹目标的自旋、进动和翻滚等微动特征为TBMD目标识别提供了一种非常具有潜力的手段。利用精细的微多普勒幅度和频率的差异进行精确身份辨别，可以为大纵深、强机动和高度自动化的现代立体战争作战平台提供更可靠的敌我识别能力，减少战争中误击、误伤和误炸事件，更重要的是可为快速高效打击提供技术支撑。为了满足高分辨率和反隐身的需求，近几年，太赫兹雷达技术发展迅速。太赫兹频段具有亚毫米甚至更短波长以及GHz以上带宽，有利于快速目标成像和获取目标特征的细节，从而可对目标进行更精确的识别。对刚体目标能够大大提高成像分辨率；然而，由于分辨率高，微动在太赫兹频段将表现得更为明显，一方面影响成像质量，另一方面又是重要的识别依据。因此，在太赫兹波段的微动目标成像和特征提取具有重要的研究意义。
技术实现思路
为了实现在太赫兹波段对微动目标ISAR成像和特征提取，本专利技术提出了一种基于Radon变换的太赫兹ISAR目标微动部件分离与特征参数估计方法。本专利技术提供一种太赫兹ISAR目标微动部件分离与特征参数估计方法，包含以下过程：S1.太赫兹波段的雷达发射线性调频信号，通过解线频调的方式接收回波，并进行传统平动补偿；S2.距离脉压后，对其进行Radon变换；S3.在Radon变换域中寻找峰值点，并计算对应的方向角、原点到直线的距离；S4.根据方向角、原点到直线的距离，计算原距离脉压像对应的距离方程；S5.在信号域构造回波信号基，对其沿距离方向进行FFT变换，并与原距离脉压像利用最小二乘法，估计散射点对应的幅度；S6.对多个散射点分别计算回波信号基包络与幅度的乘积并进行叠加，得到刚体距离像；S7.从原距离脉压像中减去刚体距离像得到微动距离像；S8.微动部分回波信号表现为正弦曲线形式的瞬时距离，其中包括正弦曲线的振幅、频率、初始相位三个空间坐标参数；对这三个空间坐标参数分别设定搜索区间，建立离散的参数空间；S9.对参数空间中的每一个点的空间坐标参数，计算在距离和方位时刻的微动信号基；S10.将分离后的微动距离像转换至信号域，并与每个空间坐标参数下的微动信号基点乘，并进行二维FFT变换；S11.遍历参数空间，记录对应的每个熵值，得到与回波信号通过二维积分得到的ISAR图像的最小熵值对应的三个空间坐标参数，作为最优参数。本专利技术的所述太赫兹ISAR目标微动部件分离与特征参数估计方法，旨在解决在太赫兹波段由目标部件微运动引起的ISAR成像质量降低和微动特征参数估计困难等问题，提出基于Radon变换的微动部件包络分离和参数估计的技术方案，并同时实现微动目标的ISAR成像。本专利技术提供的技术方案有望提高现代雷达对未知非合作目标识别和解译能力，具有实际应用前景。本专利技术另提供一种太赫兹ISAR目标微动部件包络分离方法，利用微动在太赫兹波段距离高分辨条件下的包络起伏，基于Radon变换的直线检测估计刚体的线性包络，并通过构造信号分离刚体部件和非刚体部件。本专利技术还提供一种基于最小熵的微动参数估计方法，该方法避免了太赫兹波段微多普勒谱的混叠处理，利用Radon变换的曲线检测对微动部件回波进行曲线积分，基于最小熵的准则获得微动参数的最优估计。附图说明图1为微动目标ISAR成像几何模型；图2为技术方案总流程图；图3a～图3f为仿真结果图，图3a为X波段1GHz条件下的距离脉压像，图3b为THz波段10GHz条件下的距离脉压像，图3c为图3b的Radon变换结果，图3d为分离后的刚体距离像，图3e为分离后的微动部件距离像，图3f为三参数的图像熵序列。具体实施方式本专利技术提供一种太赫兹ISAR目标微动部件分离与特征参数估计方法，旨在解决在太赫兹波段由目标部件微运动引起的ISAR成像质量降低和微动特征参数估计困难等问题，提出基于Radon变换的微动部件包络分离和参数估计的技术方案，并同时实现微动目标的ISAR成像。(1)基于Radon直线检测的微动部件包络分离方法微动主要的形式包括振动、转动和锥旋等，其引起的瞬时距离形式可统一写为Rmic(t)＝Asin(ω·t+φ)(1)假设雷达发射线性调频信号，并通过解线频调(Dechirp)的方式进行脉冲压缩。若目标包含P个散射点，传统平动补偿后的回波模型可表示其中，n为距离向离散点序号，m为方位脉冲数，δp为第p个散射点的散射系数，c为光速，fc为载频，γ为调频率，为快时间序列，tm为慢时间序列，Rp为某个散射点的相对距离。如图1所示，设目标刚体包含P1个散射点，微动部件(非刚体)包含P2个散射点，因此(2)式可改写为刚体和非刚体的组成其中，δp是刚体部分第p个散射点的散射系数，δk是微动部分第k个散射点的散射系数，Rrig,p(tm)≈Xpθ(tm)+Yp(Xp和Yp分别是第p个散射点的方位和距离向坐标)，Rmic,k(tm)＝Aksin(ω·tm+φk)(其中，Ak和φk分别为振动幅度和初相)。对式(3)进行脉压后，目标的包络可写为其中，βp是处理后刚体部分第p个散射点的散射系数，βk是处理后微动部分第k个散射点的散射系数，fr表示距离向频率，λ为波长。刚体部分的包络Grig(n,m)在角度不大的条件下会呈现直线；由于窄带距离分辨率较低，微动部分Gmic(n,m)通常在窄带回波包络中不会体现，在常规波段和带宽，微动通常也不会有明显表现。在太赫兹波段可做到很大带宽，显然，随着带宽增大，距离分辨率的提高，微动在回波包络的起伏现象会随之凸显。在宽带本文档来自技高网...

【技术保护点】
一种太赫兹ISAR目标微动部件分离与特征参数估计方法，其特征在于，包含以下过程：S1.太赫兹波段的雷达发射线性调频信号，通过解线频调的方式接收回波，并进行传统平动补偿；S2.距离脉压后，对其进行Radon变换；S3.在Radon变换域中寻找峰值点，并计算对应的方向角、原点到直线的距离；S4.根据方向角、原点到直线的距离，计算原距离脉压像对应的距离方程；S5.在信号域构造回波信号基，对其沿距离方向进行FFT变换，并与原距离脉压像利用最小二乘法，估计散射点对应的幅度；S6.对多个散射点分别计算回波信号基包络与幅度的乘积并进行叠加，得到刚体距离像；S7.从原距离脉压像中减去刚体距离像得到微动距离像；S8.微动部分回波信号表现为正弦曲线形式的瞬时距离，其中包括正弦曲线的振幅、频率、初始相位三个空间坐标参数；对这三个空间坐标参数分别设定搜索区间，建立离散的参数空间；S9.对参数空间中的每一个点的空间坐标参数，计算在距离和方位时刻的微动信号基；S10.将分离后的微动距离像转换至信号域，并与每个空间坐标参数下的微动信号基点乘，并进行二维FFT变换；S11.遍历参数空间，记录对应的每个熵值，得到与回波信号通过二维积分得到的ISAR图像的最小熵值对应的三个空间坐标参数，作为最优参数。...

【技术特征摘要】
1.一种太赫兹ISAR目标微动部件分离与特征参数估计方法，其特征在于，包含以下过程：S1.太赫兹波段的雷达发射线性调频信号，通过解线频调的方式接收回波，并进行传统平动补偿；S2.距离脉压后，对其进行Radon变换；S3.在Radon变换域中寻找峰值点，并计算对应的方向角、原点到直线的距离；S4.根据方向角、原点到直线的距离，计算原距离脉压像对应的距离方程；S5.在信号域构造回波信号基，对其沿距离方向进行FFT变换，并与原距离脉压像利用最小二乘法，估计散射点对应的幅度；S6.对多个散射点分别计算回波信号基包络与幅度的乘积并进行叠加，得到刚体距离像；S7.从原距离脉压像中减去刚体距离像得到微动距离像；S8.微动部分回波信号表现为正弦曲线形式的瞬时距离，其中包括正弦曲线的振幅、频率、初始相位三个空间坐标参数；对这三个空间坐标参数分别设定搜索区间，建立离散的参数空间；S9.对参数空间中的每一个点的空间坐标参数，计算在距离和方位时刻的微动信号基；S10.将分离后的微动距离像转换至信号域，并与每个空间坐标参数下的微动信号基点乘，并进行二维FFT变换；S11.遍历参数空间，记录对应的每个熵值，得到与回波信号通过二维积分得到的ISAR图像的最小熵值对应的三个空间坐标参数，作为最优参数。2.如权利要求1所述太赫兹ISAR目标微动部件分离与特征参数估计方法，其特征在于，根据Radon变换域中与峰值点对应的方向角α、原点到直线的距离r，得到原距离脉压像对应的距离方程为：其中，斜率为κ＝tanα；ρr为距离分辨率；信号域构造回波信号基的包络为：其中，fr表示距离向频率，c为光速，γ为调频率，n为距离向离散点序号，m为方位脉冲数，λ为波长；刚体的包络矩阵形式为：Erig＝[Erig(n,m)]N×M利用最小二乘法估计得到刚体散射点对应幅度的表达式为：Grig＝[Grig(n,m)]N×M是目标包络矩阵；对Radon变换域中所有直线转换的亮点进行检测，提取多个刚体散射点；叠加多个刚体散射点的幅度与回波信号基包络矩阵形式的乘积，从原脉包络矩阵形式G中减去，得到的微动部件的包络矩阵形式为：3.如权利要求2所述太赫兹ISAR目标微动部件分离与特征参数估计方法，其特征在于，沿着微动部分回波信号中的瞬时距离Rmic(t)＝Asin(ω·t+φ)，对回波信号进行二维积分得到ISAR图像；基于最小熵准则进行参数估计，以使对应正弦曲线的振幅、频率、初始相位三个空间坐标参数(A,ω,φ)与目标参数匹配，ISAR图像的聚焦性能最好；其中，对参数空间的散射点计算各时刻的积分路径Rmic,k(tm)，构造匹配信号，并取ISAR图像熵值最小对应的空间坐标参数，Smic＝[Smic(n,m)]N×M是微动部分的回波矩阵，对应微动部分构造的基矩阵为Emic,k＝[Emic,k(n,m)]N×M；微动信号基为快时间序列，tm为慢时间序列；对k个散射点的微动参数估计后，从...

【专利技术属性】
技术研发人员：盛佳恋，付朝伟，倪亮，杨成山，王平，
申请(专利权)人：上海无线电设备研究所，
类型：发明
国别省市：上海,31