CN110456351A - 基于时变幅值lfm信号参数估计的机动目标isar成像方法 - Google Patents

基于时变幅值lfm信号参数估计的机动目标isar成像方法 Download PDF

Info

Publication number
CN110456351A
CN110456351A CN201910809413.7A CN201910809413A CN110456351A CN 110456351 A CN110456351 A CN 110456351A CN 201910809413 A CN201910809413 A CN 201910809413A CN 110456351 A CN110456351 A CN 110456351A
Authority
CN
China
Prior art keywords
signal
amplitude
variable amplitude
isar
imaging
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.)
Granted
Application number
CN201910809413.7A
Other languages
English (en)
Other versions
CN110456351B (zh
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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201910809413.7A priority Critical patent/CN110456351B/zh
Publication of CN110456351A publication Critical patent/CN110456351A/zh
Application granted granted Critical
Publication of CN110456351B publication Critical patent/CN110456351B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50—Systems of measurement based on relative movement of target
    • G01S13/505—Systems of measurement based on relative movement of target using Doppler effect for determining closest range to a target or corresponding time, e.g. miss-distance indicator
    • 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
    • G01S13/9029—SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
    • 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]

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

基于时变幅值LFM信号参数估计的机动目标ISAR成像方法，本发明涉及机动目标ISAR成像方法。本发明的目的是为了解决现有的RID方法所得的ISAR图像会有模糊、重影的问题。一、建模多分量时变幅值LFM信号；二、初始k＝1，i＝0；三、更新第k个距离门的慢时间信号并计算能量；四、若Ek＞En，i＝i+1，执行五；否则执行九；五、得到频谱；六、构造wl和wG(m)；七、得到时变幅值的估计值；八、获取剩余信号，返回三；九、k＝k+1，若k≤Nr，执行三，否则执行十；十、得到瞬时ISAR图像。本发明用于ISAR图像获取领域。

Description

基于时变幅值LFM信号参数估计的机动目标ISAR成像方法
技术领域
本发明涉及基于时变幅值参数估计的机动目标ISAR成像方法。
背景技术
逆合成孔径雷达(Inverse Synthetic Aperture Radar,ISAR)成像通常先对回波数据做运动补偿和距离压缩，将目标转换为转台目标进行成像。当目标具有机动性时，由于其多普勒频率是时变的，可将慢时间回波建模为线性调频(Linear FrequencyModulation,LFM)信号，然后采用距离-瞬时多普勒(Range Instantaneous Doppler,RID)算法进行成像。目前提出的RID算法都是基于各距离门的慢时间信号具有恒定幅值这一假设，只估计相位进行多普勒重构得到瞬时像。但实际上慢时间信号的幅值是时变的，其原因有两方面：(1)由于散射点的越距离单元徙动(Migration Through Resolution Cell,MTRC)现象，在整个成像期间散射点将从一个距离门转移到另一个距离门，导致慢时间信号具有时变的幅值；(2)在成像期间目标相对于雷达视线的姿态变化，使得雷达散射截面积发生变化，这也导致了幅值是非恒定的。在此情况下，采用现有的RID方法所得的ISAR图像将会有模糊、重影。因此，研究一种准确、有效的估计时变幅值的方法来提高ISAR图像质量是十分必要的。
发明内容
本发明的目的是为了解决现有的RID方法所得的ISAR图像会有模糊、重影的问题，而提出基于时变幅值线性调频信号参数估计的机动目标ISAR成像方法。
基于时变幅值LFM信号参数估计的机动目标ISAR成像方法具体过程为：
步骤一、获取ISAR回波，对回波信号进行运动补偿和距离压缩后，将每个距离门的慢时间信号建模为多分量时变幅值LFM信号，即单个距离门回波的形式为
式中，Ai(n)为第i个分量的时变幅值，ai1和ai2分别为第i分量的一、二阶相位参数，n为时域采样点，j为虚数单位，j2＝-1；
所述多分量为大于等于2个分量；
步骤二、初始化距离门个数k＝1，初始化散射点个数i＝0，设置能量阈值En；
步骤三、更新第k个距离门的慢时间信号sk(n)，并计算sk(n)的能量
式中，为sk(n)的共轭；
步骤四、比较Ek与En的大小，若Ek＞En，则更新i＝i+1，记录散射点i的距离向位置ri，执行步骤五；否则，执行步骤九；
步骤五、用三次相位函数法得到sk(n)的第i个分量一、二阶相位参数估计值和并构造补偿函数将sc(n)的共轭与sk(n)相乘得到含有时变幅值信息的信号sa(n)，对含有时变幅值信息的信号sa(n)做快速傅里叶变换得到sa(n)的频谱Sa(m)，m为频域采样点；
步骤六、设置门限幅度值TH，搜索Sa(m)的峰值位置记为pmax，从峰值位置pmax处向左、右搜索Sa(m)中第一个低于幅度门限值TH的索引，将索引分别记为p1和p2，得到自适应窗长wl＝2max{pmax-p1,p2-pmax}；
构造以pmax为中心、wl为窗长的自适应高斯窗函数wG(m)；
步骤七、将自适应高斯窗函数wG(m)与频谱Sa(m)相乘，提取、估计出时变幅值的频域信息对做快速傅里叶逆变换得到时变幅值的估计值
步骤八、构造窗长为wl的自适应带阻滤波器wR(m)；
用CLEAN算法获取剩余信号sk(n)，过程为：
将自适应带阻滤波器wR(m)与频谱Sa(m)相乘并做IFFT，将所得信号与补偿函数sc(n)相乘，获得剩余信号sk(n)＝IFFT[Sa(m)wR(m)]sc(n)；返回步骤三；
步骤九、令k＝k+1，比较k与快时间点数Nr的大小，若k≤Nr，则执行步骤三，否则，执行步骤十；
步骤十、选定成像时刻n＝n0，确定每个散射点的瞬时多普勒频率和瞬时幅值根据瞬时多普勒频率确定散射点i的方位向位置，根据瞬时幅值确定散射点i的强度，当得到所有散射点i的方位向位置、距离向位置ri和强度后，即得到瞬时ISAR图像。
本发明的有益效果为：
考虑到由于成像期间散射点的越距离单元徙动现象和雷达散射截面积变化将引起各距离单元慢时间回波幅度变化，本发明中提出一种自适应高斯窗函数，能够根据频谱特点自适应调整窗长，可以有效、准确地估计信号的时变幅值，弥补了经典RID算法基于回波幅值恒定这一缺陷。采用本发明方法对机动目标ISAR回波成像，可以减轻图像中的模糊、重影问题，减少图像中虚假散射点的数量，提高成像质量。
附图说明
图1为本发明基于时变幅值LFM信号参数估计的ISAR成像方法流程图；
图2为本发明相位参数和时变幅值估计流程图；
图3为实施例一中本发明方法估计时变幅值结果图；
图4a为实施例二中本发明方法第一分量估计时变幅值结果图；
图4b为实施例二中本发明方法第二分量估计时变幅值结果图；
图4c为实施例二中本发明方法第三分量估计时变幅值结果图；
图4d为实施例二中本发明方法第四分量估计时变幅值结果图；
图5为实施例三中所用的散射点模型示意图；
图6为实施例三中RD算法成像结果图；
图7为实施例三中恒定幅度RID算法成像结果图；
图8为实施例三中本发明方法成像结果图；
图9为实施例四中RD算法成像结果图；
图10a为实施例四中成像时刻0.1250s恒定幅度RID算法成像结果图；
图10b为实施例四中成像时刻0.2250s恒定幅度RID算法成像结果图；
图11a为实施例四中成像时刻0.1250s本发明方法成像结果图；
图11b为实施例四中成像时刻0.2250s本发明方法成像结果图。
具体实施方式
具体实施方式一：结合图1、图2说明本实施方式，本实施方式基于时变幅值LFM信号参数估计的机动目标ISAR成像方法具体过程为：
步骤一、获取ISAR回波，对回波信号进行运动补偿和距离压缩后，将每个距离门的慢时间信号建模为多分量时变幅值LFM信号，即单个距离门回波的形式为
式中，Ai(n)为第i个分量的时变幅值，ai1和ai2分别为第i分量的一、二阶相位参数，n为时域采样点，j为虚数单位，j2＝-1；
所述多分量为大于等于2个分量；
步骤二、初始化距离门个数k＝1，初始化散射点个数i＝0，设置能量阈值En；
步骤三、更新第k个距离门的慢时间信号sk(n)，并计算sk(n)的能量
式中，为sk(n)的共轭；
步骤四、比较Ek与En的大小，若Ek＞En，则更新i＝i+1，记录散射点i的距离向位置ri，执行步骤五；否则，执行步骤九；
步骤五、用三次相位函数(Cubic Phase Function,CPF)法得到sk(n)的第i个分量一、二阶相位参数估计值和并构造补偿函数将sc(n)的共轭与sk(n)相乘得到含有时变幅值信息的信号sa(n)，对含有时变幅值信息的信号sa(n)做快速傅里叶变换(Fast Fourier Transform,FFT)得到sa(n)的频谱Sa(m)，m为频域采样点；
步骤六、设置门限幅度值TH，搜索Sa(m)的峰值位置记为pmax，从峰值位置pmax处向左、右搜索Sa(m)中第一个低于幅度门限值TH的索引，将索引分别记为p1和p2，得到自适应窗长wl＝2max{pmax-p1,p2-pmax}；
构造以pmax为中心、wl为窗长的自适应高斯窗函数wG(m)；
步骤七、将自适应高斯窗函数wG(m)与频谱Sa(m)相乘，提取、估计出时变幅值的频域信息对做快速傅里叶逆变换(Inverse Fast FourierTransfrom,IFFT)得到时变幅值的估计值
步骤八、构造窗长为wl的自适应带阻滤波器wR(m)；
用CLEAN算法获取剩余信号sk(n)，过程为：
将自适应带阻滤波器wR(m)与频谱Sa(m)相乘并做IFFT，将所得信号与补偿函数sc(n)相乘，获得(更新)剩余信号sk(n)＝IFFT[Sa(m)wR(m)]sc(n)；返回步骤三；
步骤九、令k＝k+1，比较k与快时间点数Nr的大小，若k≤Nr，则执行步骤三，否则(得到所有散射点的相位参数和时变幅值估计值)，执行步骤十；
步骤十、选定成像时刻n＝n0，确定每个散射点的瞬时多普勒频率和瞬时幅值根据瞬时多普勒频率确定散射点i的方位向位置，根据瞬时幅值确定散射点i的强度，当得到所有散射点i的方位向位置、距离向位置ri和强度后，即得到瞬时ISAR图像。
具体实施方式二：本实施方式与具体实施方式一不同的是，所述步骤六中高斯窗函数的定义为：
式中，α为标准差，决定高斯窗函数的衰减速率，如果高斯窗函数的通带被限制在均值左右两侧三倍标准差内。
其它步骤及参数与具体实施方式一相同。
具体实施方式三：本实施方式与具体实施方式一或二不同的是，所述标准差α与信号长度N和窗长wl的关系为：
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四：本实施方式与具体实施方式一至三之一不同的是，所述步骤八中带阻滤波器的定义为：
其它步骤及参数与具体实施方式一至三之一相同。
采用以下实施例验证本发明的有益效果：
实施例一：
采用以下实施例验证本发明的有益效果：
本实施例用自适应高斯窗函数估计单分量时变幅值LFM信号的时变幅度。
构造单分量时变幅值LFM信号，信号长度N为257，LFM信号的相位参数为a1＝π/8和a2＝0.005，其时变幅值形式为：
其中，abs[·]为取绝对值或复数模函数。用具体实施方式一中所述步骤估计信号的时变幅值。
用CPF法得到相位参数a2的估计值为0.0050，用Dechirp法得到相位参数a1的估计值为0.3912，所得结果与参数真值非常接近。图3给出了利用自适应高斯窗函数估计信号时变幅值的结果，可以看到估计结果与真值十分接近。
实施例二：
本实施例用自适应高斯窗函数估计多分量时变幅值LFM信号的时变幅度。
构造一个四分量时变幅值的LFM信号，信号长度N为257，LFM信号的相位参数如表1所示，四分量LFM信号是时变幅值分别为：
A3(n)＝2abs[cos(0.002πn)]
A4(n)＝3abs[sin(0.004πn)]
表1时变幅值LFM信号的相位参数
表2为用CPF法和Dechirp法得到的参数估计值，其误差在可接受范围内。图4a、4b、4c、4d给出了利用本发明方法所得四个分量的时变幅值估计值与真值，可以看到用本发明得到的估计值与真值误差很小，估计结果精确。
表2多分量LFM信号相位参数估计结果
实施例三：
本实施例用本发明方法结合RID算法对ISAR仿真信号成像。
仿真模型为一个飞机模型如图5所示，仿真所用ISAR系统参数如表3所示。
表3 ISAR系统仿真参数
图6为RD算法成像结果，从图中可以看出所得ISAR像方位向严重散焦。图7和图8分别为恒定幅度RID算法和本发明方法的成像结果，用恒定幅度RID算法得到的ISAR瞬时像在机头、机尾处均有许多虚假散射点，这主要是由于MTRC现象使得慢时间回波的幅度时变。对比成像结果可发现本发明方法得到的ISAR像聚焦效果最佳，图像质量最好。RD、恒定幅度RID和本发明方法所得图像的熵值分别为8.2717、4.8109和3.7672。
实施例四：
本实施例用本发明方法对雅克-42飞机实测数据进行成像，雷达系统及目标参数如表4所示。
表4雷达系统及目标参数
图9为RD算法所得ISAR像，可以看到所得图像严重散焦。图10a、10b和图11a、11b分别为恒定幅度RID算法和本发明方法选定成像时刻(t＝0.1250s,0.2250s)所得ISAR瞬时像，通过对比可以看出本发明方法所得图像聚焦效果最好。三种算法所得图像熵值如表5所示，从表中可以看出本发明方法所得ISAR图像熵值最小，图像质量最佳。
表5图像熵值对比
本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (4)

1.基于时变幅值LFM信号参数估计的机动目标ISAR成像方法，其特征在于：所述方法具体过程为：
步骤一、获取ISAR回波，对回波信号进行运动补偿和距离压缩后，将每个距离门的慢时间信号建模为多分量时变幅值LFM信号，即单个距离门回波的形式为
式中，Ai(n)为第i个分量的时变幅值，ai1和ai2分别为第i分量的一、二阶相位参数，n为时域采样点，j为虚数单位，j2＝-1；
所述多分量为大于等于2个分量；
步骤二、初始化距离门个数k＝1，初始化散射点个数i＝0，设置能量阈值En；
步骤三、更新第k个距离门的慢时间信号sk(n)，并计算sk(n)的能量
式中，为sk(n)的共轭；
步骤四、比较Ek与En的大小，若Ek＞En，则更新i＝i+1，记录散射点i的距离向位置ri，执行步骤五；否则，执行步骤九；
步骤五、用三次相位函数法得到sk(n)的第i个分量一、二阶相位参数估计值和并构造补偿函数将sc(n)的共轭与sk(n)相乘得到含有时变幅值信息的信号sa(n)，对含有时变幅值信息的信号sa(n)做快速傅里叶变换得到sa(n)的频谱Sa(m)，m为频域采样点；
步骤六、设置门限幅度值TH，搜索Sa(m)的峰值位置记为pmax，从峰值位置pmax处向左、右搜索Sa(m)中第一个低于幅度门限值TH的索引，将索引分别记为p1和p2，得到自适应窗长wl＝2max{pmax-p1,p2-pmax}；
构造以pmax为中心、wl为窗长的自适应高斯窗函数wG(m)；
步骤七、将自适应高斯窗函数wG(m)与频谱Sa(m)相乘，提取、估计出时变幅值的频域信息对做快速傅里叶逆变换得到时变幅值的估计值
步骤八、构造窗长为wl的自适应带阻滤波器wR(m)；
用CLEAN算法获取剩余信号sk(n)，过程为：
将自适应带阻滤波器wR(m)与频谱Sa(m)相乘并做IFFT，将所得信号与补偿函数sc(n)相乘，获得剩余信号sk(n)＝IFFT[Sa(m)wR(m)]sc(n)；返回步骤三；
步骤九、令k＝k+1，比较k与快时间点数Nr的大小，若k≤Nr，则执行步骤三，否则，执行步骤十；
步骤十、选定成像时刻n＝n0，确定每个散射点的瞬时多普勒频率和瞬时幅值根据瞬时多普勒频率确定散射点i的方位向位置，根据瞬时幅值确定散射点i的强度，当得到所有散射点i的方位向位置、距离向位置ri和强度后，即得到瞬时ISAR图像。
2.根据权利要求1所述基于时变幅值LFM信号参数估计的机动目标ISAR成像方法，其特征在于：所述步骤六中高斯窗函数的定义为：
式中，α为标准差。
3.根据权利要求1或2所述基于时变幅值LFM信号参数估计的机动目标ISAR成像方法，其特征在于：所述标准差α与信号长度N和窗长wl的关系为：
4.根据权利要求3所述基于时变幅值LFM信号参数估计的机动目标ISAR成像方法，其特征在于：所述步骤八中带阻滤波器的定义为：
CN201910809413.7A 2019-08-29 2019-08-29 基于时变幅值lfm信号参数估计的机动目标isar成像方法 Active CN110456351B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910809413.7A CN110456351B (zh) 2019-08-29 2019-08-29 基于时变幅值lfm信号参数估计的机动目标isar成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910809413.7A CN110456351B (zh) 2019-08-29 2019-08-29 基于时变幅值lfm信号参数估计的机动目标isar成像方法

Publications (2)

Publication Number Publication Date
CN110456351A true CN110456351A (zh) 2019-11-15
CN110456351B CN110456351B (zh) 2021-06-15

Family

ID=68489843

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910809413.7A Active CN110456351B (zh) 2019-08-29 2019-08-29 基于时变幅值lfm信号参数估计的机动目标isar成像方法

Country Status (1)

Country Link
CN (1) CN110456351B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110850411A (zh) * 2019-11-20 2020-02-28 哈尔滨工业大学 基于l统计量和粒子群优化的微多普勒抑制方法
CN111142105A (zh) * 2020-01-06 2020-05-12 西南电子技术研究所（中国电子科技集团公司第十研究所） 复杂运动目标isar成像方法
CN111427039A (zh) * 2020-04-21 2020-07-17 北京航天长征飞行器研究所 基于微动特征调制的isar成像欺骗干扰方法和装置
CN112014781A (zh) * 2020-09-02 2020-12-01 无锡鸣石峻致医疗科技有限公司 一种磁共振回波信号的相位矫正方法、装置、计算机设备及计算机可读存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107085212A (zh) * 2017-04-01 2017-08-22 中国人民解放军空军工程大学 一种基于线性调频步进信号的自旋目标时变三维成像方法
CN107656274A (zh) * 2017-09-20 2018-02-02 电子科技大学 Sar动目标回波徙动校正与方位能量积累方法
CN107831473A (zh) * 2017-10-13 2018-03-23 西安电子科技大学 基于高斯过程回归的距离‑瞬时多普勒图像序列降噪方法
CN109633647A (zh) * 2019-01-21 2019-04-16 中国人民解放军陆军工程大学 一种双基地isar稀疏孔径成像方法
WO2019080521A1 (zh) * 2017-10-25 2019-05-02 清华大学 太赫兹isar成像方法及系统
CN110109108A (zh) * 2019-04-26 2019-08-09 西安电子科技大学 基于stft和frft的运动目标雷达三维成像方法
CN110133648A (zh) * 2019-05-13 2019-08-16 中国石油大学(华东) 一种选取逆合成孔径雷达船只成像时窗的方法
CN110146889A (zh) * 2019-06-17 2019-08-20 中国人民解放军国防科技大学 基于最优回波子区域选择的大转角isar成像方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107085212A (zh) * 2017-04-01 2017-08-22 中国人民解放军空军工程大学 一种基于线性调频步进信号的自旋目标时变三维成像方法
CN107656274A (zh) * 2017-09-20 2018-02-02 电子科技大学 Sar动目标回波徙动校正与方位能量积累方法
CN107831473A (zh) * 2017-10-13 2018-03-23 西安电子科技大学 基于高斯过程回归的距离‑瞬时多普勒图像序列降噪方法
WO2019080521A1 (zh) * 2017-10-25 2019-05-02 清华大学 太赫兹isar成像方法及系统
CN109633647A (zh) * 2019-01-21 2019-04-16 中国人民解放军陆军工程大学 一种双基地isar稀疏孔径成像方法
CN110109108A (zh) * 2019-04-26 2019-08-09 西安电子科技大学 基于stft和frft的运动目标雷达三维成像方法
CN110133648A (zh) * 2019-05-13 2019-08-16 中国石油大学(华东) 一种选取逆合成孔径雷达船只成像时窗的方法
CN110146889A (zh) * 2019-06-17 2019-08-20 中国人民解放军国防科技大学 基于最优回波子区域选择的大转角isar成像方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
李亚超: "ISAR目标运动参数估计及成像技术研究", 《中国博士学位论文全文数据库 信息科技辑》 *
桂小刚: "空间快速机动小目标ISAR成像算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
邢孟道 等: "基于调幅_线性调频信号参数估计的机动目标成像方法", 《现代雷达》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110850411A (zh) * 2019-11-20 2020-02-28 哈尔滨工业大学 基于l统计量和粒子群优化的微多普勒抑制方法
CN110850411B (zh) * 2019-11-20 2023-03-10 哈尔滨工业大学 基于l统计量和粒子群优化的微多普勒抑制方法
CN111142105A (zh) * 2020-01-06 2020-05-12 西南电子技术研究所（中国电子科技集团公司第十研究所） 复杂运动目标isar成像方法
CN111142105B (zh) * 2020-01-06 2023-06-13 西南电子技术研究所（中国电子科技集团公司第十研究所） 复杂运动目标isar成像方法
CN111427039A (zh) * 2020-04-21 2020-07-17 北京航天长征飞行器研究所 基于微动特征调制的isar成像欺骗干扰方法和装置
CN111427039B (zh) * 2020-04-21 2021-02-12 北京航天长征飞行器研究所 基于微动特征调制的isar成像欺骗干扰方法和装置
CN112014781A (zh) * 2020-09-02 2020-12-01 无锡鸣石峻致医疗科技有限公司 一种磁共振回波信号的相位矫正方法、装置、计算机设备及计算机可读存储介质
CN112014781B (zh) * 2020-09-02 2021-04-20 无锡鸣石峻致医疗科技有限公司 一种磁共振回波信号的相位矫正方法、装置、计算机设备及计算机可读存储介质

Also Published As

Publication number Publication date
CN110456351B (zh) 2021-06-15

Similar Documents

Publication Publication Date Title
CN110456351A (zh) 基于时变幅值lfm信号参数估计的机动目标isar成像方法
CN105137498B (zh) 一种基于特征融合的地下目标探测识别系统及方法
CN106468770B (zh) K分布杂波加噪声下的近最优雷达目标检测方法
CN108107430B (zh) 一种基于分数阶傅立叶变换的舰船目标isar成像方法
CN104155650B (zh) 一种基于熵权值法点迹质量评估的目标跟踪方法
CN109188414A (zh) 一种基于毫米波雷达的手势运动探测方法
CN111722199A (zh) 一种基于卷积神经网络的雷达信号检测方法
CN105866776B (zh) 一种地基sar的高质量动态ps点选择方法
CN109633647A (zh) 一种双基地isar稀疏孔径成像方法
CN106019256B (zh) 基于自回归模型的雷达信号自适应检测方法
CN107843875A (zh) 基于奇异值分解降噪的贝叶斯压缩感知雷达数据融合方法
CN107608935A (zh) 基于时间重排压缩变换的冲击类信号时频分析与重构方法
CN109284280A (zh) 仿真数据优化方法、装置及存储介质
CN105699952A (zh) 海杂波k分布形状参数的双分位点估计方法
CN104749564A (zh) 海杂波Weibull幅度分布参数的多分位点估计方法
CN107831473B (zh) 基于高斯过程回归的距离-瞬时多普勒图像序列降噪方法
CN108459308A (zh) 一种基于时变rcs数据的回波模拟方法及装置
CN109343060A (zh) 基于深度学习时频分析的isar成像方法及系统
CN110333489A (zh) 采用cnn与rsva组合对sar回波数据旁瓣抑制的处理方法
CN105699947B (zh) 一种sar图像旁瓣抑制方法
CN114970378A (zh) 一种基于gan网络的海杂波样本库构建方法
CN107329115A (zh) 基于gcrbf网络的lfm信号参数估计方法
CN102012510A (zh) 基于时间——相位导数分布的逆合成孔径雷达成像方法
CN109584256A (zh) 一种基于霍夫直线检测的脉冲星dm估计算法
CN111796288B (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