CN109655819A - 一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法 - Google Patents

Abstract

本发明公开一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法，包含：导引头依次发射多个线性调频脉冲，并接收目标回波；每个脉冲进行匹配滤波获得目标的一维距离像；利用多个一维距离像相参处理获得目标的距离多普勒图像；利用恒虚警检测将目标不同散射点置于不同的距离多普勒单元中；利用相控阵比幅单脉冲测角获得目标不同散射点的方位角和俯仰角信息；利用“多径杂波抑制算法”抑制多径杂波；利用坐标系欧拉转换获得目标坐标系中的三维图像。本发明相较于二维图像能够获得目标的高度信息，利用“多径杂波抑制算法”抑制多径杂波，获得打击目标的精确三维图像；本发明在准确获得打击目标结构特征的情况下，为后续精确打击提供图像基础。

Description

一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法

技术领域

本发明涉及弹载脉冲多普勒雷达领域，特别涉及一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法。

背景技术

随着现代作战环境的日趋复杂，多任务导弹已经成为未来空射导弹的主要发展方向。空射导弹在有效地打击空中机动目标的同时，具有打击地面军事基地的多重打击能力。但是空射导弹受到战斗部威力限制，打击地面军事基地时无法获得最大打击效能。在制导末端对地面军事基地进行三维成像能够获得基地的结构特征，结合结构特征进行易损部位的打击可以有效地提高打击效能。同时制导末端打击地面军事基地时需要考虑多径杂波问题，成像时多径杂波会造成雷达目标检测困难和目标图像失真。

合成孔径成像和实孔径成像是目前较为成熟的成像体制，合成孔径雷达无法完成前视成像，因此本文将基于实孔径多普勒波束锐化成像，结合相控阵比幅单脉冲测角获得目标的三维图像。同时利用“多径杂波抑制算法”在高度域进行多径杂波抑制，解决制导末端多径杂波对目标图像的干扰。获得打击目标的精确三维图像，为易损部位的选择提供基础。

现有技术中公开有一种弹载/机载雷达单脉冲前视成像处理方法，利用匹配滤波完成距离向聚焦，距离徙动校正后利用单脉冲测角完成方位能量积累，最终获得目标图像。该方法可以减小雷达与目标之间相对运动对成像质量的影响，也可以提高测角精度，但是不具备频率分辨能力，同时无法完成目标的三维成像和多径杂波的抑制，不能准确的获得目标的结构特征。

现有技术中还记载有一种用于机载雷达前视成像的自聚焦方法，通过提取距离门样本数据，对鉴角曲线进行估计，经过多次迭代获得准确的鉴角曲线，然后进行单脉冲成像。该方法可解决测角误差导致的单脉冲图像分辨率下降问题，但是不具备频率分辨能力，同时无法完成目标的三维成像和多径杂波的抑制，不能准确的获得目标的结构特征。

目前，学者还研究了一种基于天线方向图的机载雷达单脉冲前视成像方法，该方法利用测试的和通道数据和差通道数据计算方向图，修正理想的天线方向图，将同距离单元内相同角度的散射点进行非相参累加，实现前视成像。该方法能解决天线方向图失配的问题，但是不具备频率分辨能力，同时无法完成目标的三维成像和多径杂波的抑制，不能准确地获得目标的结构特征。

现有技术中的以上方法均不能实现目标的三维成像，获得目标准确的特征结构，为后续精确打击提供基础，同时都不能完成多径杂波的抑制，影响了雷达的检测性能，导致目标图像的失真。

发明内容

本发明的目的是提供一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法，能够利用实孔径多普勒波束锐化和相控阵比幅单脉冲测角获得目标的三维图像，利用“多径杂波抑制算法”抑制多径杂波。

为了达到上述发明目的，本发明公开了一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法，该方法包含以下步骤：

S1、导引头依次发射多个线性调频脉冲；

S2、接收目标回波；

S3、每个脉冲进行匹配滤波获得目标的一维距离像；

S4、利用多个一维距离像相参处理获得目标的距离多普勒图像；

S5、利用恒虚警检测将目标不同散射点置于不同的距离多普勒单元中；

S6、利用相控阵比幅单脉冲测角获得目标不同散射点的方位角和俯仰角信息；

S7、利用多径杂波抑制算法抑制多径杂波；

S8、利用坐标系欧拉转换获得目标坐标系中的三维图像。、

优选地，所述步骤S1中，进一步包含：导引头依次发射N个线性调频脉冲，发射的线性调频脉冲表达式为：A·exp[j(2πf₀t+πkt²)]，其中，A为发射信号的幅度，t为时间，f₀为载频，k为调频率，j为虚数单位。

优选地，所述步骤S2中，进一步包含：所述发射的线性调频脉冲经过距离为R₀处的目标反射后，接收到的目标回波为：A·exp{j[2πf₀(t-2R₀/c)+πk(t-2R₀/c)²]}。

优选地，所述步骤S3中，进一步包含：经过下变频后获得的混频信号为A’·exp{j[2πf₀(-2R₀/c)+πk(t-2R₀/c)²]}，利用线性调频信号的匹配滤波器对每个线性调频信号的脉冲做脉冲压缩后得到的一维距离像可表示为A”·sinc(R-R₀)；其中，N个线性调频信号脉冲压缩后获得N个一维距离像。

优选地，所述步骤S4中，进一步包含：将N个一维距离像按照距离门进行重排，每个距离门内有N个脉冲压缩数据，对每个距离门内的N个数据进行离散傅里叶变换处理，获得每个距离门内的多普勒信息，所有的距离门内N个数据处理完毕后获得目标散射点的距离多普勒二维图像。

优选地，所述步骤S5中，进一步包含：对距离多普勒图像进行二维恒虚警检测，获得每个目标散射点的距离多普勒单元(R_i,f_i)。

优选地，所述步骤S6中，进一步包含：对存在目标散射点的每个距离多普勒单元(R_i,f_i)进行相控阵比幅单脉冲测角，利用相控阵比幅单脉冲测角获得每个目标散射点的方位角信息α_i和俯仰角信息β_i；测角公式为：

其中，F₁(θ)表示子阵1的方向图，F₂(θ)表示子阵2的方向图，M表示阵元数，D表示阵元间隔，ξ表示子阵与轴线的夹角，ε表示待测角度。

优选地，所述步骤S7中，进一步包含：利用每个散射点的距离R_i、俯仰角β_i计算出其高度h_i＝R_i sinβ_i，通过每个散射点的高度统计得到地平面的高度估计H，然后剔除掉h_i<H的散射点，完成多径杂波的抑制；其中，多径杂波是指雷达回波在目标不同散射点之间多次反射或者在目标散射点与地面之间多次反射，该多径杂波的时延比真实散射点的时延长，多径杂波成像时位于地平面以下。

优选地，所述步骤S8中，进一步包含：将每个散射点的坐标(R_i,α_i,β_i)进行转换获得坐标(x_i,y_i,z_i)，转换公式为：x_i＝R_i cosβ_i cosα_i，y_i＝R_i cosβ_i sinα_i，z_i＝R_i sinβ_i，利用转换公式完成散射点三维成像。

与现有技术相比，本发明的有益效果为：

本发明应用于脉冲多普勒雷达领域，能够在制导末端利用实孔径多普勒波束锐化和相控阵比幅单脉冲测角对打击目标进行三维成像处理，相较于二维图像能够获得目标的高度信息，利用“多径杂波抑制算法”抑制多径杂波，获得打击目标的精确三维图像；本发明在准确获得打击目标结构特征的情况下，为后续精确打击提供图像基础。

附图说明

图1本发明的基于实孔径多普勒波束锐化的杂波抑制三维成像所算法流程示意图；

图2本发明的理论目标图像示意图；

图3本发明的采样后时域图像示意图；

图4本发明的脉冲压缩后的一维距离像示意图；

图5本发明的时频二维图像示意图；

图6本发明的未剔除成像结果示意图；

图7本发明的最终成像结果示意图。

具体实施方式

本发明公开了一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法，为了使本发明更加明显易懂，以下结合附图和具体实施方式对本发明做进一步说明。

如图1所示，本发明的基于实孔径多普勒波束锐化的杂波抑制三维成像方法包含以下步骤：

S1、导引头依次发射多个线性调频脉冲；

S2、接收目标回波；

S3、每个脉冲进行匹配滤波获得目标的一维距离像；

S4、利用多个一维距离像相参处理获得目标的距离多普勒图像；

S5、利用恒虚警检测将目标不同散射点置于不同的距离多普勒单元中；

S6、利用相控阵比幅单脉冲测角获得目标不同散射点的方位角和俯仰角信息；

S7、利用“多径杂波抑制算法”抑制多径杂波；

S8、利用坐标系欧拉转换获得目标坐标系中的三维图像。

所述步骤S1中，进一步包含：

导引头依次发射N个线性调频脉冲，发射的线性调频脉冲表示式为A·exp[j(2πf₀t+πkt²)]，其中，A为发射信号的幅度，t为时间，f₀为载频，k为调频率，j为虚数单位；

所述步骤S2中，进一步包含：

发射的线性调频脉冲经过距离为R₀处的目标反射后，接收到的目标回波为A·exp{j[2πf₀(t-2R₀/c)+πk(t-2R₀/c)²]}，其中，c表示光速。

所述步骤S3中，进一步包含：

经过下变频后获得的混频信号为A’·exp{j[2πf₀(-2R₀/c)+πk(t-2R₀/c)²]}，利用线性调频信号的匹配滤波器对每个线性调频信号的脉冲做脉冲压缩后得到的一维距离像可表示为A”·sinc(R-R₀)，其中，A’表示回波信号的幅度，受到目标特性和传播距离等因素的影响，与A不相同且无确定关系；另，A”表示经过脉冲压缩之后的信号，受滤波器等因素的影响，与之前A和A’不同，三者都不相同且无确定关系；R表示距离变量。N个线性调频信号脉冲压缩后可获得N个一维距离像。

所述步骤S4中，进一步包含：

将N个一维距离像按照距离门进行重排，每个距离门内有N个脉冲压缩数据，对每个距离门内的N个数据进行DFT(Discrete Fourier Transform，离散傅里叶变换)处理，获得每个距离门内的多普勒信息，所有的距离门内N个数据处理完毕后获得目标散射点的距离多普勒二维图像。

所述步骤S5中，进一步包含：

对距离多普勒图像进行二维恒虚警检测，获得每个目标散射点的距离多普勒单元(R_i,f_i)。

所述步骤S6中，进一步包含：

对存在目标散射点的每个距离多普勒单元(R_i,f_i)进行相控阵比幅单脉冲测角，F₁(θ)表示子阵1的方向图，F₂(θ)表示子阵2的方向图，M表示阵元数，D表示阵元间隔，ξ表示子阵与轴线的夹角，ε表示待测角度，测角公式为：，利用相控阵比幅单脉冲测角获得每个目标散射点的方位角信息α_i和俯仰角信息β_i。

所述步骤S7中，进一步包含：

多径杂波指的是雷达回波在目标不同散射点之间多次反射或者在目标散射点与地面之间多次反射，这种杂波的时延一般比真实散射点的时延要长，因此其成像时往往位于地平面以下。利用每个散射点的距离R_i、俯仰角β_i计算出其高度h_i＝R_i sinβ_i，通过每个散射点的高度统计得到地平面的高度估计H，然后剔除掉h_i<H的散射点，即可完成多径杂波的抑制。

所述步骤S8中，进一步包含：

将每个散射点的坐标(R_i,α_i,β_i)进行转换获得坐标(x_i,y_i,z_i)，具体的转换公式为x_i＝R_i cosβ_i cosα_i，y_i＝R_i cosβ_i sinα_i，z_i＝R_i sinβ_i，利用转换公式完成散射点三维成像。

如图2所示，本发明的模拟打击目标为具有五个强散射点的地面军事目标，其在目标坐标系中所成的三维目标像如图所示。

其中，五个强散射点的x轴坐标分别为[-100，-100，100，100，0](单位：m)，y轴坐标分别为[-50，50，-50，50，0](单位：m)，z轴坐标分别为[0，0，0，0，100](单位：m)，弹目相对速度为[1103，1110，1106，1097，1094](单位：m/s)，导引头与打击目标中心的相对距离为1385m，波束方位角指向为45.00°，波束俯仰角指向为35.00°。导引头发射32个线性调频脉冲，每个线性调频脉冲的脉冲宽度为0.5μs，脉冲重复周期为80μs，信号带宽为40MHz。接收到的信号经过下变频、低通滤波后得到的图像如图3所示。

对32个线性调频回波脉冲分别做脉冲压缩，脉冲压缩所使用的滤波器为exp(jπkt²)，第一个线性调频脉冲做脉冲压缩获得的一维距离图像如图4所示，强散射点集中在第347，351，363，367，378，386，393个距离单元内，对每个脉冲回波做脉冲压缩可得到32个一维距离像。

将32个一维距离像按照距离门进行重排，然后对每个距离门内的数据做32点DFT，对所有距离门DFT后获得时频二维图像，由于PRF(pulse repetition frequency，脉冲重复频率)为12.5kHz，因此多普勒维的分辨率为390.6Hz，对于目标模式所给出的五个强散射点，其在多普勒维应该位于不同的多普勒分辨单元。所获得的距离多普勒二维图像如图5所示。强散射点所在的距离单元对应的多普勒单元分别为18，17，20，17，19，15，16。

对所获得距离多普勒图像进行二维恒虚警检测，恒虚警门限设置为20dB，则获得的强散射点的距离多普勒单元分别为(347，18)、(351，17)、(363，20)、(367，17)、(378，19)、(386，15)、(393，16)。

相控阵雷达的阵元数M＝32，表示阵元间隔D＝λ/2，子阵与轴线的夹角ξ＝2.5°，对每个距离多普勒单元的散射点进行相控阵比幅单脉冲测角获得每个散射点的方位角信息分别为46.97°、43.12°、50.53°、46.66°、39.81°、45.00°、43.36°，每个散射点的俯仰角信息分别为37.95°、33.61°、36.00°、32.08°、34.33°、38.49°、32.85°。

由散射点的距离单元和多普勒单元计算得到的散射点的距离信息分别为1301m、1316m、1361m、1376m、1417m、1447m、1474m，此时目标的三维图像如图6所示，F点和G点为多径杂波。结合散射点的俯仰角信息进行“多径杂波抑制算法”，由公式h_i＝R_isinβ_i计算得到每个散射点的高度分别为800.18m、728.60m、800.11m、730.99m、799.38m、900.85m、799.36m，经过散射点个数与散射点高度的统计得到的地平面高度估计为800.00m，由于成像场景设置为空射导弹打击地面军事目标，因此直接用公式计算得到的高度正方向垂直向下，即高度小于800.00m的散射点实际上是位于地平面以下的散射点。因此高度为728.60m和730.99m的散射点被剔除。

坐标系欧拉转换：首先计算得到目标点在与弹体固联的大地坐标系中的x轴坐标，计算公式为x_i＝R_icosβ_icosα_i，得到的各散射点的x轴坐标分别为700.17m、700.06m、899.28m、801.17m、900.13m，然后计算各散射点在与弹体固联的大地坐标系中的y轴坐标，计算公式为y_i＝R_icosβ_isinα_i，得到的各散射点的y轴坐标分别为750.15m、850.14m、749.40m、801.14m、850.14m，最后计算各散射点在与弹体固联的大地坐标系中的z轴坐标，计算公式为z_i＝R_isinβ_i，得到的各散射点的z轴坐标分别为800.18m、800.11m、799.38m、900.85m、799.36m。

将与弹体固联的大地坐标系转换至与目标固联的大地坐标系中，即将与弹体固联的大地坐标系中各散射点的三维坐标减去弹目距离在各坐标轴上的投影。导引头与打击目标中心的相对距离R₀为1385m，波束方位角指向α₀为45.00°，波束俯仰角指向β₀为45.00°。则在x轴上的投影距离x₀的计算公式为x₀＝R₀cosβ₀cosα₀，计算结果为800m，在y轴上的投影距离y₀的计算公式为y₀＝R₀cosβ₀sinα₀，计算结果为800m，在z轴上的投影距离z₀的计算公式为z₀＝R₀sinβ₀，计算结果为800m。由此计算得到的目标坐标系中各散射点的x轴坐标分别为-99.83m、-99.94m、99.28m、1.17m、100.13m，y轴坐标分别为49.85m、50.14m、50.60m、1.14m、50.14m，z轴坐标分别为0.18m、0.11m、-0.62m、100.85m、-0.67m，得到目标的三维图像如图7所示。

尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。

Claims (9)

1.一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法，其特征在于，该方法包含以下步骤：

S1、导引头依次发射多个线性调频脉冲；

S2、接收目标回波；

S3、每个脉冲进行匹配滤波获得目标的一维距离像；

S4、利用多个一维距离像相参处理获得目标的距离多普勒图像；

S5、利用恒虚警检测将目标不同散射点置于不同的距离多普勒单元中；

S6、利用相控阵比幅单脉冲测角获得目标不同散射点的方位角和俯仰角信息；

S7、利用多径杂波抑制算法抑制多径杂波；

S8、利用坐标系欧拉转换获得目标坐标系中的三维图像。

2.如权利要求1所述的基于实孔径多普勒波束锐化的杂波抑制三维成像方法，其特征在于，

所述步骤S1中，进一步包含：

导引头依次发射N个线性调频脉冲，发射的线性调频脉冲表达式为：A·exp[j(2πf₀t+πkt²)]，其中，A为发射信号的幅度，t为时间，f₀为载频，k为调频率，j为虚数单位。

3.如权利要求2所述的基于实孔径多普勒波束锐化的杂波抑制三维成像方法，其特征在于，

所述步骤S2中，进一步包含：

所述发射的线性调频脉冲经过距离为R₀处的目标反射后，接收到的目标回波为：A·exp{j[2πf₀(t-2R₀/c)+πk(t-2R₀/c)²]}，其中，c表示光速。

4.如权利要求3所述的基于实孔径多普勒波束锐化的杂波抑制三维成像方法，其特征在于，

所述步骤S3中，进一步包含：

经过下变频后获得的混频信号为A’·exp{j[2πf₀(-2R₀/c)+πk(t-2R₀/c)²]}，利用线性调频信号的匹配滤波器对每个线性调频信号的脉冲做脉冲压缩后得到的一维距离像可表示为A”·sinc(R-R₀)；其中，N个线性调频信号脉冲压缩后获得N个一维距离像。

5.如权利要求4所述的基于实孔径多普勒波束锐化的杂波抑制三维成像方法，其特征在于，

所述步骤S4中，进一步包含：

将N个一维距离像按照距离门进行重排，每个距离门内有N个脉冲压缩数据，对每个距离门内的N个数据进行离散傅里叶变换处理，获得每个距离门内的多普勒信息，所有的距离门内N个数据处理完毕后获得目标散射点的距离多普勒二维图像。

6.如权利要求5所述的基于实孔径多普勒波束锐化的杂波抑制三维成像方法，其特征在于，

所述步骤S5中，进一步包含：

对距离多普勒图像进行二维恒虚警检测，获得每个目标散射点的距离多普勒单元(R_i,f_i)。

7.如权利要求6所述的基于实孔径多普勒波束锐化的杂波抑制三维成像方法，其特征在于，

所述步骤S6中，进一步包含：

对存在目标散射点的每个距离多普勒单元(R_i,f_i)进行相控阵比幅单脉冲测角，利用相控阵比幅单脉冲测角获得每个目标散射点的方位角信息α_i和俯仰角信息β_i；测角公式为：

其中，F₁(θ)表示子阵1的方向图，F₂(θ)表示子阵2的方向图，M表示阵元数，D表示阵元间隔，ξ表示子阵与轴线的夹角，ε表示待测角度。

8.如权利要求7所述的基于实孔径多普勒波束锐化的杂波抑制三维成像方法，其特征在于，

所述步骤S7中，进一步包含：

利用每个散射点的距离R_i、俯仰角β_i计算出其高度h_i＝R_i sinβ_i，通过每个散射点的高度统计得到地平面的高度估计H，然后剔除掉h_i<H的散射点，完成多径杂波的抑制；其中，多径杂波是指雷达回波在目标不同散射点之间多次反射或者在目标散射点与地面之间多次反射，该多径杂波的时延比真实散射点的时延长，多径杂波成像时位于地平面以下。

9.如权利要求8所述的基于实孔径多普勒波束锐化的杂波抑制三维成像方法，其特征在于，

所述步骤S8中，进一步包含：

将每个散射点的坐标(R_i,α_i,β_i)进行转换获得坐标(x_i,y_i,z_i)，转换公式为：x_i＝R_i cosβ_i cosα_i，y_i＝R_i cosβ_i sinα_i，z_i＝R_i sinβ_i，利用转换公式完成散射点三维成像。