CN105738897B - 基于联合稀疏特性的信号距离像重构方法 - Google Patents

基于联合稀疏特性的信号距离像重构方法 Download PDF

Info

Publication number
CN105738897B
CN105738897B CN201610108222.4A CN201610108222A CN105738897B CN 105738897 B CN105738897 B CN 105738897B CN 201610108222 A CN201610108222 A CN 201610108222A CN 105738897 B CN105738897 B CN 105738897B
Authority
CN
China
Prior art keywords
signal
random
echo
subpulse
distance
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.)
Expired - Fee Related
Application number
CN201610108222.4A
Other languages
English (en)
Other versions
CN105738897A (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201610108222.4A priority Critical patent/CN105738897B/zh
Publication of CN105738897A publication Critical patent/CN105738897A/zh
Application granted granted Critical
Publication of CN105738897B publication Critical patent/CN105738897B/zh
Expired - Fee Related 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/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
    • 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]

Abstract

本发明公开了一种基于联合稀疏特征的信号距离像重构方法，该方法包括以下主要步骤：(1)首先生成随机稀疏频率步进信号回波；(2)基于回波信号联合稀疏信息建立随机稀疏频率步进信号距离像联合重构模型；(3)按照子脉冲随机方式不同对每组子脉冲信号构建不同的量测矩阵，并利用傅里叶基作为回波信号稀疏基矩阵；(4)在寻找回波信号支撑集时利用所有信号支撑集相同的信息，并利用DC－OMP算法对距离像进行联合重构。

Description

基于联合稀疏特性的信号距离像重构方法
技术领域
本发明是一种距离像重构方法，在随机频率步进信号的距离像重构过程中充分利用回波信号联合稀疏的特性提高距离像的重构精度。
背景技术
逆合成孔径雷达作为一种可以对非合作目标进行观测、成像、跟踪的雷达系统，主要通过大带宽信号得到距离向的高分辨距离像。高分辨距离像也是ISAR成像的关键。步进频信号作为一种高效的大带宽信号常被用于成像系统中。传统的步进频率信号通过发射载频均匀步进的脉冲信号实现距离向上的高分辨，但是必须连续发射多个载频连续变化的子脉冲，所需成像时间较长，当要得到二维图像时，所需的时间又会延长。稀疏步进频(sparsestepped frequency，SSF)信号通过随机发射传统步进频信号的部分子脉冲来实现传统信号的高分辨能力，极大的缩短了成像时间。为增加信号的抗干扰性能，通常每组信号之间载频变化规律都不相同。然而，传统的步进频率信号通过IFFT就能得到目标的高分辨率距离像，但是对这种随机稀疏步进频(Random sparse stepped frequency，RSSF)信号来说并不适用，怎样处理此类信号得到了学者的广泛关注。
压缩感知理论的引入为处理此类信号提供了一条有效的解决方法。然而，目前的大多数压缩感知距离像成像方法都是对各组脉冲进行分别重构，没有考虑各组脉冲信号之间的联系。实际上，在雷达观测目标时，由于目标在观测区域具有稀疏性，因此雷达回波都具有相同的稀疏结构。本文基于回波信号的联合稀疏模型(Joint Sparsity Model，JSM)，在重构目标高分辨距离像时，联合利用所有脉冲支撑集相同的信息，提高了一维距离像的重构精度，同时降低了量测值的数量，增强了抗噪性能。理论分析与仿真验证了本文算法的有效性。
发明内容
(1)首先生成随机稀疏频率步进信号回波；
(2)基于回波信号联合稀疏信息建立随机稀疏频率步进信号距离像联合重构模型；
(3)按照子脉冲随机方式不同对每组子脉冲信号构建不同的量测矩阵，并利用傅里叶基作为回波信号稀疏基矩阵；
(4)在寻找回波信号支撑集时利用所有信号支撑集相同的信息，并利用DC－OMP算法对距离像进行联合重构。
本发明相比背景技术具有如下优点：
(1)本发明可以在发射较少子脉冲的情况下实现对距离像的高分辨重构，减少了成像时间，有利于节省雷达系统的时间资源。
(2)本发明可以在较低信噪比条件下实现对距离像的高精度重构，具有较强的鲁棒性。
附图说明
附图1是本发明的流程图。参照附图1，本发明的流程由建立随机稀疏频率步进联合重构回波模型、依据子脉冲随机方式构建量测矩阵、利用子脉冲支撑集相同的特征对支撑集进行估计、利用最小二乘法对对应位置的系数进行重构、将重构得到的系数赋予距离像对应的位置以及获得最终的距离象重构结果6个部分组成。其中1用于建立随机稀疏频率步进联合重构回波模型；2用于依据子脉冲随机方式构建量测矩阵；3用于利用子脉冲支撑集相同的特征对支撑集进行估计；4用于利用最小二乘法对对应位置的系数进行重构；5用于将重构得到的系数赋予距离像对应的位置；6用于获得最终的距离象重构结果。
附图2是信噪比为－5dB条件下的距离像重构图，图中横坐标为距离，纵坐标为散射点强度。
附图3为不同信噪比条件下重构误差对比曲线，图中横坐标为信噪比，纵坐标为误差。
附图4为随机发射32个子脉冲条件下的距离像重构结果，图中横坐标为距离，纵坐标为散射点强度。
附图5为随机发射不同个数子脉冲条件下的重构误差曲线，图中横坐标为距离像量测数，纵坐标为误差。
具体实施方式
实施本发明的原理如下：基于压缩感知理论构建随机稀疏频率步进信号的距离像联合重构模型，然后根据子脉冲信号随机稀疏方式的不同设计出不同的量测矩阵，最后基于随机稀疏频率步进信号的联合稀疏重构模型，利用DC－OMP算法实现对距离像的联合稀疏重构。本发明可以在较低的信噪比，较少的距离向量测数的条件下实现对距离像的高精度重构。
(1)构建随机稀疏频率步进信号联合稀疏回波模型。假设随机稀疏频率步进信号其回波基频信号为：
式中，σi,i＝1,......I为目标散射点强度；M为随机发射的子脉冲个数；Tr为脉冲重复间隔；T为脉冲脉宽；；Ri(t,na)表示雷达与目标散射点的距离；Na为脉组数；为其载频变化规律，f0为载频，Δf为频率步进量。
对上式的子脉冲进行采样，得到的采样信号为：
Δθ为转动步长；(xi,yi)为目标坐标。
对于由于观测的是同一目标，因此所有回波信号其稀疏特性相同，即具有相同的稀疏结构。上述性质可以表示为：
式中，为第na组回波的支撑集。
因此，基于压缩感知理论的目标高精度距离像联合重构模型为：
其中U(m；na)为传统步进频回波数据；X(na)表示高分辨率距离像； 为量测矩阵；Ψ为距离向稀疏基。
在求解上述模型时，可以充分利用回波信号联合稀疏这一先验信息，提高重构精度。
(2)设计距离向量测矩阵距离向稀疏基Ψ。按照子脉冲随机规律设计距离向量测矩阵
由于回波在频域时是稀疏的，因此，设计距离向稀疏基Ψ为傅里叶基。
(3)利用DC－OMP算法进行联合稀疏重构。具体为：首先，利用回波信号之间联合稀疏特性，寻找并更新支撑集位置
式中，为残差信号。
为估计的支撑集集合。
其次，利用最小二乘法重构出目标支撑集位置的系数
(4)将得到的值赋予X(na)中相应的位置，得到Na组重构出的高精度距离像。
下面结合实例说明一下整个发明的优势。
假设发射信号参数如下：f0＝10GHz，Δf＝3MHz，Na＝300，M＝128，N＝256，散射点相对位置[022.5]m，强度为[10.50.1]，图2为信噪比为－5dB时的重构图，图3为不同信噪比条件下的重构误差曲线。可以看出：在信噪比为－5dB时，传统压缩感知方法(OMP算法)已经出现了重构误差(弱散射点位置重构不准)，而本发明方法可以准确的重构出弱散射点，显示出较强的抗噪性能。再考虑量测值的影响。假设其他条件不变，信噪比为10dB，图4为距离向随机发射M＝32个子脉冲时的重构结果，图5为不同量测值条件下的重构误差曲线图。当量测数M降为32时，传统压缩感知方法(OMP算法)估计结果出错，而本发明仍能准确重构出目标的高分辨距离像。从误差曲线图也可以看出，在不同的量测值条件下，本发明重构误差较小，证明了该方法可以在更低的量测值条件下得到精确的重构结果。

Claims (1)

1.基于联合稀疏特性的信号距离像重构方法，其特征在于包括以下主要技术措施：
(1)雷达发射Na组随机稀疏频率步进脉冲信号，随机稀疏即每组脉冲只随机发射传统频率步进N个子脉冲中的M个，得到随机稀疏频率步进回波信号的压缩感知重构模型为
其中，X(na)表示高分辨率距离像； 为量测矩阵；Ψ为距离向稀疏基；U(m；na)为传统步进频回波数据
(m,na)为随机正整数，表示第na组脉冲子脉冲步进倍数；σi,i＝1,......I为目标散射点强度；Δθ为转动步长；(xi,yi)为目标坐标；
(2)按照子脉冲随机规律设计
(3)设计距离向稀疏基Ψ为傅里叶基；
(4)基于回波信号之间联合稀疏特性并利用DC－OMP算法进行距离像联合重构，首先寻找并更新支撑集位置
式中，为残差信号，
为估计的支撑集集合；
(5)利用最小二乘法重构出目标支撑集位置的系数
(6)将得到的值赋予X(na)中相应的位置，得到Na组重构出的高精度距离像。
CN201610108222.4A 2016-02-26 2016-02-26 基于联合稀疏特性的信号距离像重构方法 Expired - Fee Related CN105738897B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610108222.4A CN105738897B (zh) 2016-02-26 2016-02-26 基于联合稀疏特性的信号距离像重构方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610108222.4A CN105738897B (zh) 2016-02-26 2016-02-26 基于联合稀疏特性的信号距离像重构方法

Publications (2)

Publication Number Publication Date
CN105738897A CN105738897A (zh) 2016-07-06
CN105738897B true CN105738897B (zh) 2018-09-25

Family

ID=56249511

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610108222.4A Expired - Fee Related CN105738897B (zh) 2016-02-26 2016-02-26 基于联合稀疏特性的信号距离像重构方法

Country Status (1)

Country Link
CN (1) CN105738897B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106990392B (zh) * 2017-02-15 2019-06-25 中国人民解放军国防科学技术大学 一种基于随机步进频信号的空间目标微动信息获取方法
CN107037409B (zh) * 2017-04-28 2021-09-14 中国人民解放军空军工程大学 基于压缩感知的mimo雷达波形分离方法
CN107192991B (zh) * 2017-07-18 2019-10-11 中国人民解放军国防科学技术大学 一种射频辐射式仿真宽带脉冲雷达成像等效模拟方法
CN109085589B (zh) * 2018-10-16 2019-04-30 中国人民解放军国防科技大学 基于图像质量引导的稀疏孔径isar成像相位自聚焦方法
CN112130149B (zh) * 2020-09-24 2022-07-12 哈尔滨工程大学 结合正交匹配追踪算法的稀疏发射的isar成像方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5424667B2 (ja) * 2008-06-02 2014-02-26 三菱電機株式会社 画像レーダ装置
KR101103493B1 (ko) * 2010-07-13 2012-01-10 영남대학교 산학협력단 역합성 개구면 레이더 영상을 이용한 표적 식별 방법 및 그 장치
CN102135618B (zh) * 2010-12-23 2012-09-05 北京航空航天大学 一种针对空域稀疏sar图像的频域压缩感知方法
CN102928838A (zh) * 2012-10-08 2013-02-13 西北工业大学 基于稀疏表示与时频变换的距离-瞬时多普勒成像方法
CN104749573B (zh) * 2013-12-31 2017-05-03 中国科学院电子学研究所 空频域二维稀疏的步进频率sar成像方法

Also Published As

Publication number Publication date
CN105738897A (zh) 2016-07-06

Similar Documents

Publication Publication Date Title
CN105738897B (zh) 基于联合稀疏特性的信号距离像重构方法
Johansson et al. Three-dimensional ground-penetrating radar imaging using synthetic aperture time-domain focusing
CN107167781B (zh) 海杂波幅度对数正态分布参数的分位点估计方法
Sun et al. Advanced signal processing methods for ground-penetrating radar: Applications to civil engineering
CN106872974B (zh) 基于高超声速平台双通道雷达的高精度运动目标成像方法
CN106405552B (zh) 基于wvd—pga算法的sar雷达目标聚焦方法
CN110988834B (zh) 一种基于自适应阈值双门限的脉冲到达时间测量方法
Zhang et al. A novel moving target imaging algorithm for HRWS SAR based on local maximum-likelihood minimum entropy
CN107255814A (zh) 一种基于lfmsk波形的雷达目标检测方法
CN108318879A (zh) 基于iaa谱估计技术的isar图像横向定标方法
CN108535726A (zh) 基于幂权Fourier变换的ISAR成像方法
RU2416105C1 (ru) Способ определения параметров движения воздушных объектов в обзорных радиолокаторах за счет использования когерентных свойств отраженных сигналов
CN110596706B (zh) 一种基于三维图像域投射变换的雷达散射截面积外推方法
CN110031838A (zh) 一种不需要参考发射波形的穿墙雷达墙体参数估计方法
CN107390213A (zh) 一种基于滑动时窗的探地雷达记录剖面的时延曲线提取方法
CN107741581B (zh) 基于截断矩的广义帕累托分布参数估计方法
EP2817655B1 (en) Systems and methods for image sharpening
CN112505719A (zh) 多普勒频率校正二次补偿的激光测风雷达测风方法及系统
CN109856636B (zh) 曲线合成孔径雷达自适应三维成像方法
CN106707284A (zh) 一种用于双基地逆合成孔径雷达的成像方法
CN108415017A (zh) 复杂目标雷达散射特性稀疏表征的一维增广状态空间方法
CN106093927B (zh) 基于雷达大脉宽信号的目标测速测距方法
CN116430384A (zh) 一种基于多主影像处理的层析sar三维成像方法
CN114252878A (zh) 一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法
CN114428235A (zh) 基于决策级融合的空间微动目标识别方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180925

Termination date: 20190226