CN113030967A - 一种用于逆合成孔径雷达成像的距离对准方法及系统 - Google Patents

一种用于逆合成孔径雷达成像的距离对准方法及系统 Download PDF

Info

Publication number
CN113030967A
CN113030967A CN202110263357.9A CN202110263357A CN113030967A CN 113030967 A CN113030967 A CN 113030967A CN 202110263357 A CN202110263357 A CN 202110263357A CN 113030967 A CN113030967 A CN 113030967A
Authority
CN
China
Prior art keywords
recursion
pulse
distance
signal
alignment
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.)
Pending
Application number
CN202110263357.9A
Other languages
English (en)
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.)
Chinese People's Liberation Army 91604
Original Assignee
Chinese People's Liberation Army 91604
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 Chinese People's Liberation Army 91604 filed Critical Chinese People's Liberation Army 91604
Priority to CN202110263357.9A priority Critical patent/CN113030967A/zh
Publication of CN113030967A publication Critical patent/CN113030967A/zh
Pending legal-status Critical Current

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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904—SAR modes
    • G01S13/9064—Inverse SAR [ISAR]
    • 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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021—SAR image post-processing techniques
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40—Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种用于逆合成孔径雷达成像的距离对准方法及系统。该方法包括：保存脉压后的ISAR原始数据，然后将每个脉冲回波脉压进行傅里叶变换获得每个脉冲的频域数据；对每个脉冲的频域数据通过递推算法获得递推参数并递推出雷达发射信号带宽之外的信息；对递推后的脉冲频域数据通过傅立叶逆变换成像，得到超分辨图像；通过全局最小熵算法对超分辨图像进行距离对准处理，得到每个脉冲压缩图像距离对准误差；原始脉压后的回波信号按照Δr(n)进行移动，该方法具有计算效率高，对准效果好的基础上，实现亚距离单元对准，提高了ISAR距离对准的精度，提高ISAR成像质量。

Description

一种用于逆合成孔径雷达成像的距离对准方法及系统
技术领域
本发明涉及雷达成像技术领域，特别是涉及一种用于逆合成孔径雷达成像的距离对准方法及系统。
背景技术
逆合成孔径雷达(Inverse Synthetic Aperture Radar，简称ISAR)的成像目标一般为非合作目标，目标相对于雷达的运动可以分解为两部分，即平动和转动。其中，对成像有贡献的只有转动部分，平动必须要被补偿掉。经过平动补偿以后，ISAR成像就转变成转台成像。平动补偿是ISAR成像的关键技术，它可以分为两步。第一步是补偿复包络的时间延迟，也就是所谓的距离对准；第二步是找到一个参考点，然后根据这个参考点来调整初相，通常把这一步叫相位补偿。其中，距离对准作为相位补偿的基础直接影响到ISAR成像的质量，是ISAR成像的关键技之一。
随着对ISAR成像分辨要求的不断提高，非合作目标高分辨率成像逐渐成为一个重要的研究方向。对于ISAR成像来说，分辨率的提高意味着需要更大的积累角和更长的孔径时间，同时要求距离对准算法具有更好的精度和鲁棒性。现有的距离对准算法主要有：相邻包络互相关法、包络最大修正峰度法、全局对准算法、全局最小熵算法等。以上算法都具有比较好的鲁棒性，对提高空间目标的成像质量具有一定作用，但是以上算法的距离对准精度是一个距离单元，无法满足非合作目标日益增加的高分辨率成像要求。
发明内容
本发明旨在至少解决现有技术中存在的技术问题，特别创新地提出了一种用于逆合成孔径雷达成像的距离对准方法及系统。
为了实现本发明的上述目的，根据本发明的第一个方面，本发明提供了一种用于逆合成孔径雷达成像的距离对准方法，包括：步骤S1，保存脉压后的ISAR原始数据，然后将每个脉冲回波脉压进行傅里叶变换获得每个脉冲的频域数据；步骤S2，对每个脉冲的频域数据通过递推算法获得递推参数，利用所述递推参数基于公式
Figure BDA0002971008460000021
递推出雷达发射信号带宽之外的信息，其中
Figure BDA0002971008460000022
为当前位置的递推信息；n为每个脉冲频域数据之外的递推数据点序号，且N≤n≤4N，N表示脉冲的频域数据的点数；x[n-k]为当前位置之前第k个位置的有效数据；
Figure BDA0002971008460000023
为前位置之前第k个位置的p阶递推参数，p为所选取的递推阶数，且
Figure BDA0002971008460000024
步骤S3，对递推后的脉冲频域数据通过傅立叶逆变换成像，得到超分辨图像；步骤S4，通过全局最小熵算法对超分辨图像进行距离对准处理，得到每个脉冲压缩图像距离对准的Δr(m),0≤m≤M-1，其中，Δr(m)为通过全局最小熵算法得到的第m个脉冲的距离对准误差，M为脉冲数；步骤S5，原始脉压后的回波信号按照Δr(m)进行移动完成亚距离单元精度的距离对准。
在本发明的一种优选实施方式中，在所述步骤S2中，通过递推算法获得递推参数的过程为：步骤S21，初始化预测参数：
Figure BDA0002971008460000031
Figure BDA00029710084600000325
表示前向预测误差的初始值，
Figure BDA0002971008460000032
表示后向预测误差的初始值，
Figure BDA0002971008460000033
表示信号的自相关函数；步骤S22，通过Burg算法估计反射系数，使预测值和真实值之间的均方误差之和最小来求取前位置之前第k个位置第k阶递推的反射系数
Figure BDA0002971008460000034
Figure BDA0002971008460000035
Figure BDA00029710084600000326
中下标k表示递推阶数，k＝1,2,...,p；
Figure BDA0002971008460000036
表示对
Figure BDA0002971008460000037
信号取共轭；
Figure BDA0002971008460000038
表示第n－1个递推数据点的k－1阶递推数据点的后向预测误差；
Figure BDA0002971008460000039
表示第n个递推数据点的k－1阶递推数据点的前向预测误差；利用Levinson－Durbin递推得到AR参数估计值
Figure BDA00029710084600000310
Figure BDA00029710084600000327
表示前位置之前第1个位置的p阶递推参数，
Figure BDA00029710084600000311
表示前位置之前第2个位置的p阶递推参数，
Figure BDA00029710084600000312
表示前位置之前第p个位置的p阶递推参数；其中，
Figure BDA00029710084600000313
的递推过程为：按照k＝1,2,...,p的取值顺序，递推阶数k从1到p逐渐增加，第k阶递推公式为：
Figure BDA00029710084600000314
上述递推公式中，
Figure BDA00029710084600000315
中存在如下关系：
Figure BDA00029710084600000316
其中，
Figure BDA00029710084600000317
表示对
Figure BDA00029710084600000318
信号取共轭；
Figure BDA00029710084600000319
表示对
Figure BDA00029710084600000320
信号取共轭，如果k＝1，则
Figure BDA00029710084600000321
利用公式
Figure BDA00029710084600000322
递推得到末尾参数
Figure BDA00029710084600000323
将数组
Figure BDA00029710084600000324
作为递推参数。
在本发明的一种优选实施方式中，在所述步骤S1中，脉冲回波脉压为：
Figure BDA0002971008460000041
其中，B＝K·Tr表示距离向带宽，K为信号调频斜率，Tr为信号脉冲宽度，fτ为每个脉冲傅里叶变换后的距离向频率，fc为雷达发射信号频率，Rt为目标与雷达之间的斜距，c为光速，rect(·)表示矩形窗函数；对脉冲回波脉压做傅里叶变换，得到：
Figure BDA0002971008460000042
其中，令
Figure BDA0002971008460000043
τ为距离向时间。
上述技术方案的有益效果为：上述方法通过递推雷达发射信号带宽之外的信息来获得超分辨图像，且利用全局最小熵算法对超分辨图像进行距离对准处理，进而实现了亚距离单元精度的距离对准，该方法除了具有全局最小熵算法计算效率高，对准效果好的有益效果外，还实现亚距离单元对准，提高了ISAR距离对准的精度，提高ISAR成像质量。
为了实现本发明的上述目的，根据本发明的第二个方面，本发明提供了一种雷达成像系统，包括雷达和成像处理器，所述成像处理器接收雷达发射信号在非合作目标表面反射回来的回波信号，并按照本发明所述的距离对准方法进行回波信号距离对准。
上述技术方案的有益效果为：该系统通过递推雷达发射信号带宽之外的信息来获得超分辨图像，且利用全局最小熵算法对超分辨图像进行距离对准处理，进而实现了亚距离单元精度的距离对准，除了具有全局最小熵算法计算效率高，对准效果好的有益效果外，还实现亚距离单元对准，提高了ISAR距离对准的精度，提高ISAR成像质量。
附图说明
图1是本发明一具体实施方式中距离对准方法的流程示意图；
图2是本发明在一种应用场景中的距离对准效果对比图，其中，图2(a)为雷达回波信号未进行距离对准前的数据图形，图2(b)为雷达回波信号按照现有的全局最小熵距离对准算法对准后的数据图像，图2(c)为雷达回波信号按照本发明提供的距离对准方法对准后的数据图形；
图3是本发明在一种应用场景中各距离单元所有脉冲幅度累加对比图，其中，点A为使用本发明提供的距离对准方法后的脉冲幅度累加值最高点，点B为使用现有的全局最小熵距离对准算法后的脉冲幅度累加值最高点；
图4是本发明在一种应用场景中成像结果对比图，其中，图4(a)为使用现有的全局最小熵距离对准算法对准后的ISAR成像结果，图4(b)为按照本发明提供的距离对准方法对准后的ISAR成像结果；
图5是本发明在一种应用场景中方位向点目标剖面图，其中，点C为使用本发明提供的距离对准方法后的方位向点目标剖面图中的幅度最高点，点D为使用现有的全局最小熵距离对准算法后的方位向点目标剖面图中的幅度最高点；
图6是本发明在一种应用场景中距离向点目标剖面图，其中，点E为使用本发明提供的距离对准方法后的距离向点目标剖面图中的幅度最高点，点F为使用现有的全局最小熵距离对准算法后的距离向点目标剖面图中的幅度最高点。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
在本发明的描述中，需要理解的是，术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
在本发明的描述中，除非另有规定和限定，需要说明的是，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是机械连接或电连接，也可以是两个元件内部的连通，可以是直接相连，也可以通过中间媒介间接相连，对于本领域的普通技术人员而言，可以根据具体情况理解上述术语的具体含义。
本发明提供了一种用于逆合成孔径雷达成像的的距离对准方法，在一种优选实施方式中，该方法的流程示意图如图1所示，包括：
步骤S1，保存脉压后的ISAR原始数据，然后将每个脉冲回波脉压进行傅里叶变换获得每个脉冲的频域数据；
步骤S2，对每个脉冲的频域数据通过递推算法获得递推参数，利用所述递推参数基于公式
Figure BDA0002971008460000061
递推出雷达发射信号带宽之外的信息，其中
Figure BDA0002971008460000062
为当前位置的递推信息；n为每个脉冲频域数据之外的递推数据点序号，且N≤n≤4N，N表示脉冲的频域数据的点数；x[n-k]为当前位置之前第k个位置的有效数据；
Figure BDA0002971008460000071
为前位置之前第k个位置的p阶递推参数，p为所选取的递推阶数，且
Figure BDA0002971008460000072
步骤S3，对递推后的脉冲频域数据通过傅立叶逆变换成像，得到超分辨图像；
步骤S4，通过全局最小熵算法对超分辨图像进行距离对准处理，得到每个脉冲压缩图像距离对准的Δr(m),0≤m≤M-1，其中，Δr(m)为通过全局最小熵算法得到的第m个脉冲的距离对准误差，M为脉冲数；
步骤S5，原始脉压后的回波信号按照Δr(m)进行移动完成亚距离单元精度的距离对准。具体的，对脉压后的每个距离向信号由步骤S4得到的Δr(m)进行移动，实现各个脉冲回波信号有用回波对准成一条线线，完成距离对准。
在本实施方式中，全局最小熵算法为现有的距离对准算法中的一种，全局最小熵法算法通过迭代，得到使所有ISAR目标回波和的熵达到最小时各目标回波需要补偿的相移量，完成最后距离对准。该算法计算效率高，对准效果好，是目前最常用的一种距离对准算法。该算法的对准精度通常为一个距离单元。
在本实施方式中，优选的，脉冲回波脉压为：
Figure BDA0002971008460000073
其中，B＝K·Tr表示距离向带宽，K为信号调频斜率，Tr为信号脉冲宽度，fτ为每个脉冲傅里叶变换后的距离向频率，fc为雷达发射信号频率，Rt为目标与雷达之间的斜距，c为光速，rect(·)表示矩形窗函数；
对脉冲回波脉压做傅里叶变换，得到：
Figure BDA0002971008460000074
其中，令
Figure BDA0002971008460000075
Sinc(x)＝sin(πx)/πx，τ为距离向时间。
由上式可得，距离向脉冲回波脉压信号的分辨率取决于距离向信号的带宽B。要想得到距离向超分辨的成像结果，用Burg算法线性外推发射信号带宽之外的信息，使其能在不扩大发射信号带宽的情况下，有效增加目标的信号带宽，得到距离向亚距离单元信息。
在一种优选实施方式中，在步骤S2中，通过递推算法获得递推参数的过程为：
步骤S21，初始化预测参数：
Figure BDA0002971008460000081
Figure BDA0002971008460000082
Figure BDA0002971008460000083
Figure BDA0002971008460000084
Figure BDA0002971008460000085
表示前向预测误差的初始值，
Figure BDA0002971008460000086
表示后向预测误差的初始值，
Figure BDA0002971008460000087
表示信号的自相关函数；
步骤S22，通过Burg算法估计反射系数，使预测值和真实值之间的均方误差之和最小来求取前位置之前第k个位置第k阶递推的反射系数
Figure BDA0002971008460000088
Figure BDA0002971008460000089
Figure BDA00029710084600000810
中下标k表示递推阶数，k＝1,2,...,p；
Figure BDA00029710084600000811
表示对
Figure BDA00029710084600000812
信号取共轭；
Figure BDA00029710084600000813
表示第n－1个递推数据点的k－1阶递推数据点的后向预测误差；
Figure BDA00029710084600000814
表示第n个递推数据点的k－1阶递推数据点的前向预测误差；
利用Levinson－Durbin递推得到AR参数估计值
Figure BDA00029710084600000815
Figure BDA00029710084600000817
表示前位置之前第1个位置的p阶递推参数，
Figure BDA00029710084600000816
表示前位置之前第2个位置的p阶递推参数，
Figure BDA0002971008460000091
表示前位置之前第p个位置的p阶递推参数；其中，
Figure BDA0002971008460000092
的递推过程为：按照k＝1,2,...,p的取值顺序，递推阶数k从1到p逐渐增加，第k阶递推公式为：
Figure BDA0002971008460000093
上述递推公式中，
Figure BDA0002971008460000094
中存在如下关系：
Figure BDA0002971008460000095
其中，
Figure BDA0002971008460000096
表示对
Figure BDA0002971008460000097
信号取共轭；
Figure BDA0002971008460000098
表示对
Figure BDA0002971008460000099
信号取共轭，如果k＝1，则
Figure BDA00029710084600000910
利用公式
Figure BDA00029710084600000911
递推得到末尾参数
Figure BDA00029710084600000912
将数组
Figure BDA00029710084600000913
作为递推参数。
在本实施方式中，Burg算法先估计反射系数，然后利用Levinson－Durbin递推得到AR参数估计值。基于使预测误差功率的估计值达到极小利用递推方法对不同阶反射系数预测，得到反射系数的估计值。
在本发明提供的距离对准方法的一种应用场景中，获取了一组实测的雷达回波脉压数据，该组数据方位向为512个脉冲和距离向为256个距离波门。
在该应用场景中，如图2所示为距离对准效果对比图，其中，图2(a)为雷达回波信号未进行距离对准前的数据图形，图2(b)为雷达回波信号按照现有的全局最小熵距离对准算法对准后的数据图像，图2(c)为雷达回波信号按照本发明提供的距离对准方法对准后的数据图形；通过图2可以看出本发明提出的距离对准算法和现有的全局最小熵距离对准算法都能够实现距离对准。
在该应用场景中，为了对本发明提出的距离对准算法和现有的全局最小熵距离对准算法进行比较，下面对每个距离单元上各个脉冲回波幅度进行累加，其结果如图3所示，图3中，横坐标表示距离向采样像素点，纵坐标表示幅度，点A为使用本发明提供的距离对准方法后的脉冲幅度累加值最高点，点B为使用现有的全局最小熵距离对准算法后的脉冲幅度累加值最高点。通过图3可以看出，本发明提供的距离对准方法的幅度明显高于现有的全局最小熵距离对准算法幅度，也就是证明了本发明提供的距离对准方法的对准结果优于现有的全局最小熵距离对准算法的结果。
在该应用场景中，为了评估本发明提供的距离对准方法的性能，该组数据通过本发明提供的距离对准算法和现有的全局最小熵距离对准算法的成像性能进行比较，其结果如图4、图5和图6所示，图4是本发明在一种应用场景中成像结果对比图，其中，图4中横坐标表示方位向像素点，纵坐标表示距离向像素点，图4(a)为使用现有的全局最小熵距离对准算法对准后的ISAR成像结果，图4(b)为按照本发明提供的距离对准方法对准后的ISAR成像结果。图5是本发明在一种应用场景中方位向点目标剖面图，图5中横坐标表示方位向采样像素点，纵坐标表示幅度，其中，点C为使用本发明提供的距离对准方法后的方位向点目标剖面图中的幅度最高点，点D为使用现有的全局最小熵距离对准算法后的方位向点目标剖面图中的幅度最高点；图6是本发明在一种应用场景中距离向点目标剖面图，图6中横坐标表示距离向采样像素点，纵坐标表示幅度，其中，点E为使用本发明提供的距离对准方法后的距离向点目标剖面图中的幅度最高点，点F为使用现有的全局最小熵距离对准算法后的距离向点目标剖面图中的幅度最高点。
通过观察图4、图5和图6可见，通过对比数据的成像结果和点目标剖面图，可以得到本发明提供的距离对准方法的成像结果优于现有的全局最小熵的成像结果，通过数据的处理证明本发明提出距离对准算法的有效性，具有更优的距离对准精度和成像质量。
本发明还提供了一种雷达成像系统，在一种优选实施方式中，该系统包括雷达和成像处理器，成像处理器接收雷达发射信号在非合作目标表面反射回来的回波信号，并按照上述距离对准方法进行回波信号距离对准。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。

Claims (4)

1.一种用于逆合成孔径雷达成像的距离对准方法，其特征在于，包括：
步骤S1，保存脉压后的ISAR原始数据，然后将每个脉冲回波脉压进行傅里叶变换获得每个脉冲的频域数据；
步骤S2，对每个脉冲的频域数据通过递推算法获得递推参数，利用所述递推参数基于公式
Figure FDA0002971008450000011
递推出雷达发射信号带宽之外的信息，其中
Figure FDA0002971008450000012
为当前位置的递推信息；n为每个脉冲频域数据之外的递推数据点序号，且N≤n≤4N，N表示脉冲的频域数据的点数；x[n-k]为当前位置之前第k个位置的有效数据；
Figure FDA0002971008450000013
为前位置之前第k个位置的p阶递推参数，p为所选取的递推阶数，且
Figure FDA0002971008450000014
步骤S3，对递推后的脉冲频域数据通过傅立叶逆变换成像，得到超分辨图像；
步骤S4，通过全局最小熵算法对超分辨图像进行距离对准处理，得到每个脉冲压缩图像距离对准的Δr(m),0≤m≤M-1，其中，Δr(m)为通过全局最小熵算法得到的第m个脉冲的距离对准误差，M为脉冲数；
步骤S5，原始脉压后的回波信号按照Δr(m)进行移动完成亚距离单元精度的距离对准。
2.如权利要求1所述的用于逆合成孔径雷达成像的距离对准方法，其特征在于，在所述步骤S2中，通过递推算法获得递推参数的过程为：
步骤S21，初始化预测参数：
Figure FDA0002971008450000021
Figure FDA0002971008450000022
Figure FDA0002971008450000023
Figure FDA0002971008450000024
Figure FDA0002971008450000025
表示前向预测误差的初始值，
Figure FDA0002971008450000026
表示后向预测误差的初始值，
Figure FDA0002971008450000027
表示信号的自相关函数；
步骤S22，通过Burg算法估计反射系数，使预测值和真实值之间的均方误差之和最小来求取前位置之前第k个位置第k阶递推的反射系数
Figure FDA0002971008450000028
Figure FDA0002971008450000029
Figure FDA00029710084500000210
中下标k表示递推阶数，k＝1,2,...,p；
Figure FDA00029710084500000211
表示对
Figure FDA00029710084500000212
信号取共轭；
Figure FDA00029710084500000213
表示第n－1个递推数据点的k－1阶递推数据点的后向预测误差；
Figure FDA00029710084500000214
表示第n个递推数据点的k－1阶递推数据点的前向预测误差；
利用Levinson－Durbin递推得到AR参数估计值
Figure FDA00029710084500000215
Figure FDA00029710084500000216
表示前位置之前第1个位置的p阶递推参数，
Figure FDA00029710084500000217
表示前位置之前第2个位置的p阶递推参数，
Figure FDA00029710084500000218
表示前位置之前第p个位置的p阶递推参数；其中，
Figure FDA00029710084500000219
的递推过程为：按照k＝1,2,...,p的取值顺序，递推阶数k从1到p逐渐增加，第k阶递推公式为：
Figure FDA00029710084500000220
上述递推公式中，
Figure FDA00029710084500000221
中存在如下关系：
Figure FDA0002971008450000031
Figure FDA0002971008450000032
其中，
Figure FDA0002971008450000033
表示对
Figure FDA0002971008450000034
信号取共轭；
Figure FDA0002971008450000035
表示对
Figure FDA0002971008450000036
信号取共轭，如果k＝1，则
Figure FDA0002971008450000037
利用公式
Figure FDA0002971008450000038
递推得到末尾参数
Figure FDA0002971008450000039
将数组
Figure FDA00029710084500000310
作为递推参数。
3.如权利要求1所述的用于逆合成孔径雷达成像的距离对准方法，其特征在于，在所述步骤S1中，脉冲回波脉压为：
Figure FDA00029710084500000311
其中，B＝K·Tr表示距离向带宽，K为信号调频斜率，Tr为信号脉冲宽度，fτ为每个脉冲傅里叶变换后的距离向频率，fc为雷达发射信号频率，Rt为目标与雷达之间的斜距，c为光速，rect(·)表示矩形窗函数；
对脉冲回波脉压做傅里叶变换，得到：
Figure FDA00029710084500000312
其中，令
Figure FDA00029710084500000313
Sinc(x)＝sin(πx)/πx，τ为距离向时间。
4.一种雷达成像系统，其特征在于，包括雷达和成像处理器，所述成像处理器接收雷达发射信号在非合作目标表面反射回来的回波信号，并按照权利要求1－3之一所述的距离对准方法进行回波信号距离对准。
CN202110263357.9A 2021-03-11 2021-03-11 一种用于逆合成孔径雷达成像的距离对准方法及系统 Pending CN113030967A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110263357.9A CN113030967A (zh) 2021-03-11 2021-03-11 一种用于逆合成孔径雷达成像的距离对准方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110263357.9A CN113030967A (zh) 2021-03-11 2021-03-11 一种用于逆合成孔径雷达成像的距离对准方法及系统

Publications (1)

Publication Number Publication Date
CN113030967A true CN113030967A (zh) 2021-06-25

Family

ID=76469571

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110263357.9A Pending CN113030967A (zh) 2021-03-11 2021-03-11 一种用于逆合成孔径雷达成像的距离对准方法及系统

Country Status (1)

Country Link
CN (1) CN113030967A (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6255981B1 (en) * 1999-08-04 2001-07-03 Raytheon Company Method for range alignment and rotation correction of a high resolution image in an inverse synthetic aperture radar system
US20090121926A1 (en) * 2007-11-09 2009-05-14 The Boeing Company Multi-spot inverse synthetic aperture radar imaging
CN103901429A (zh) * 2014-04-09 2014-07-02 西安电子科技大学 基于稀疏孔径的机动目标逆合成孔径雷达成像方法
CN106291552A (zh) * 2016-07-28 2017-01-04 南京航空航天大学 一种能够去除突跳和漂移误差的改进isar距离对准方法
KR101738811B1 (ko) * 2016-02-26 2017-05-22 포항공과대학교 산학협력단 살리스 엔트로피를 이용한 역합성 개구면 레이더 영상 생성 장치 및 그것을 이용한 자동 초점 방법
KR20180046966A (ko) * 2016-10-28 2018-05-10 국방과학연구소 효율적인 회전 운동 보상을 위한 isar 병진 운동 보상 방법 및 장치
CN108020833A (zh) * 2017-10-25 2018-05-11 清华大学 太赫兹isar成像方法及系统
CN209673989U (zh) * 2019-02-13 2019-11-22 中国人民解放军91604部队 一种测控雷达角度零值激光标校系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6255981B1 (en) * 1999-08-04 2001-07-03 Raytheon Company Method for range alignment and rotation correction of a high resolution image in an inverse synthetic aperture radar system
US20090121926A1 (en) * 2007-11-09 2009-05-14 The Boeing Company Multi-spot inverse synthetic aperture radar imaging
CN103901429A (zh) * 2014-04-09 2014-07-02 西安电子科技大学 基于稀疏孔径的机动目标逆合成孔径雷达成像方法
KR101738811B1 (ko) * 2016-02-26 2017-05-22 포항공과대학교 산학협력단 살리스 엔트로피를 이용한 역합성 개구면 레이더 영상 생성 장치 및 그것을 이용한 자동 초점 방법
CN106291552A (zh) * 2016-07-28 2017-01-04 南京航空航天大学 一种能够去除突跳和漂移误差的改进isar距离对准方法
KR20180046966A (ko) * 2016-10-28 2018-05-10 국방과학연구소 효율적인 회전 운동 보상을 위한 isar 병진 운동 보상 방법 및 장치
CN108020833A (zh) * 2017-10-25 2018-05-11 清华大学 太赫兹isar成像方法及系统
CN209673989U (zh) * 2019-02-13 2019-11-22 中国人民解放军91604部队 一种测控雷达角度零值激光标校系统

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
SUSU0203: "Burg法求解AR(p)模型参数及MATLAB实现", 《CSDN》 *
WANG RUI: "Modified sub-integer range alignment based on minimum entropy for ISAR", 《 IET INTERNATIONAL RADAR CONFERENCE 2013》 *
YUAN LIU: "Novel ISAR Range Alignment via Minimizing the Entropy of the Sum Range Profile", 《2020 21ST INTERNATIONAL RADAR SYMPOSIUM (IRS)》 *
俞翔等: "一种改进型全局最小熵ISAR距离对准算法", 《数据采集与处理》 *
姚红梅等: "最小熵谱外推技术及其在雷达超分辨中的应用", 《现代雷达》 *
李明等著: "《战略预警雷达信号处理新技术》", 31 December 2017, 国防工业出版社 *
李源著: "《逆合成孔径雷达理论与对抗》", 31 December 2013, 国防工业出版社 *
田彪等: "宽带逆合成孔径雷达高分辨成像技术综述", 《雷达学报》 *
郭少南等: "利用超分辨技术实现ISAR包络对齐", 《第十届全国雷达学术年会论文集》 *
金胜等: "ISAR高分辨率成像方法综述", 《雷达科学与技术》 *

Similar Documents

Publication Publication Date Title
CN107229048B (zh) 一种高分宽幅sar动目标速度估计与成像方法
CN109799499B (zh) 一种穿墙雷达墙体参数估计方法
CN110133654B (zh) 一种高轨卫星sar动目标检测方法
CN110161472B (zh) 一种基于信号复用的宽带车载毫米波雷达解速度模糊方法
CN111239721B (zh) 车载mimo雷达求熵解速度模糊的方法
CN105301589B (zh) 高分辨宽测绘带sar地面运动目标成像方法
JP2012122743A (ja) レーダ装置
CN115079171A (zh) 基于混合范数的自适应稀疏约束逆合成孔径雷达成像方法
CN111060909B (zh) 一种机载雷达斜前视超分辨成像方法
CN116299551A (zh) 一种太赫兹sar二维自聚焦成像算法
CN118259252A (zh) 一种基于图信号处理的频控阵雷达多普勒估计方法
CN109143235B (zh) 一种双基前视合成孔径雷达地面运动目标检测方法
Callow et al. Motion-compensation improvement for widebeam, multiple-receiver SAS systems
CN113030967A (zh) 一种用于逆合成孔径雷达成像的距离对准方法及系统
JP3740657B2 (ja) 合成開口レーダ装置及びその像再生方法
Li et al. Stepped-frequency inverse synthetic aperture radar imaging based on adjacent pulse correlation integration and coherent processing
CN114488067B (zh) 基于Staggered SAR体制低过采样率下数据重建方法
Liu et al. Novel isar range alignment via minimizing the entropy of the sum range profile
CN114488149A (zh) 基于Tsallis熵最小化的ISAR全局包络对齐方法
CN109765540B (zh) 一种频率步进体制米波雷达目标提取方法
CN116736297B (zh) 一种异构多帧联合相参积累方法
CN115144833B (zh) 一种宽带机载雷达二维ckt-stap方法
CN113608188B (zh) 一种基于天基tdm-mimo雷达空间碎片的测角方法及系统
CN114488054B (zh) 一种计算高效的合成孔径雷达地面运动目标聚焦方法
CN116482685B (zh) 一种基于波束域相位中心互相关法的自适应dbf方法

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210625