CN105738897A - 基于联合稀疏特性的信号距离像重构方法 - Google Patents
基于联合稀疏特性的信号距离像重构方法 Download PDFInfo
- Publication number
- CN105738897A CN105738897A CN201610108222.4A CN201610108222A CN105738897A CN 105738897 A CN105738897 A CN 105738897A CN 201610108222 A CN201610108222 A CN 201610108222A CN 105738897 A CN105738897 A CN 105738897A
- Authority
- CN
- China
- Prior art keywords
- signal
- random
- theta
- sparse
- echo
- 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
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]
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于联合稀疏特征的信号距离像重构方法,该方法包括以下主要步骤:(1)首先生成随机稀疏频率步进信号回波;(2)基于回波信号联合稀疏信息建立随机稀疏频率步进信号距离像联合重构模型;(3)按照子脉冲随机方式不同对每组子脉冲信号构建不同的量测矩阵,并利用傅里叶基作为回波信号稀疏基矩阵;(4)在寻找回波信号支撑集时利用所有信号支撑集相同的信息,并利用DC-OMP算法对距离像进行联合重构。
Description
技术领域
本发明是一种距离像重构方法,在随机频率步进信号的距离像重构过程中充分利用回波信号联合稀疏的特性提高距离像的重构精度。
背景技术
逆合成孔径雷达作为一种可以对非合作目标进行观测、成像、跟踪的雷达系统,主要通过大带宽信号得到距离向的高分辨距离像。高分辨距离像也是ISAR成像的关键。步进频信号作为一种高效的大带宽信号常被用于成像系统中。传统的步进频率信号通过发射载频均匀步进的脉冲信号实现距离向上的高分辨,但是必须连续发射多个载频连续变化的子脉冲,所需成像时间较长,当要得到二维图像时,所需的时间又会延长。稀疏步进频(sparsesteppedfrequency,SSF)信号通过随机发射传统步进频信号的部分子脉冲来实现传统信号的高分辨能力,极大的缩短了成像时间。为增加信号的抗干扰性能,通常每组信号之间载频变化规律都不相同。然而,传统的步进频率信号通过IFFT就能得到目标的高分辨率距离像,但是对这种随机稀疏步进频(Randomsparsesteppedfrequency,RSSF)信号来说并不适用,怎样处理此类信号得到了学者的广泛关注。
压缩感知理论的引入为处理此类信号提供了一条有效的解决方法。然而,目前的大多数压缩感知距离像成像方法都是对各组脉冲进行分别重构,没有考虑各组脉冲信号之间的联系。实际上,在雷达观测目标时,由于目标在观测区域具有稀疏性,因此雷达回波都具有相同的稀疏结构。本文基于回波信号的联合稀疏模型(JointSparsityModel,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组重构出的高精度距离像。
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 true CN105738897A (zh) | 2016-07-06 |
CN105738897B 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) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106990392A (zh) * | 2017-02-15 | 2017-07-28 | 中国人民解放军国防科学技术大学 | 一种基于随机步进频信号的空间目标微动信息获取方法 |
CN107037409A (zh) * | 2017-04-28 | 2017-08-11 | 中国人民解放军空军工程大学 | 基于压缩感知的mimo雷达波形分离方法 |
CN107192991A (zh) * | 2017-07-18 | 2017-09-22 | 中国人民解放军国防科学技术大学 | 一种射频辐射式仿真宽带脉冲雷达成像等效模拟方法 |
CN109085589A (zh) * | 2018-10-16 | 2018-12-25 | 中国人民解放军国防科技大学 | 基于图像质量引导的稀疏孔径isar成像相位自聚焦方法 |
CN112130149A (zh) * | 2020-09-24 | 2020-12-25 | 哈尔滨工程大学 | 结合正交匹配追踪算法的稀疏发射的isar成像方法 |
CN113985412A (zh) * | 2021-11-04 | 2022-01-28 | 西安电子科技大学 | 矢量建模优化反演的运动目标三维成像方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010014700A (ja) * | 2008-06-02 | 2010-01-21 | Mitsubishi Electric Corp | 画像レーダ装置 |
CN102135618A (zh) * | 2010-12-23 | 2011-07-27 | 北京航空航天大学 | 一种针对空域稀疏sar图像的频域压缩感知方法 |
KR101103493B1 (ko) * | 2010-07-13 | 2012-01-10 | 영남대학교 산학협력단 | 역합성 개구면 레이더 영상을 이용한 표적 식별 방법 및 그 장치 |
CN102928838A (zh) * | 2012-10-08 | 2013-02-13 | 西北工业大学 | 基于稀疏表示与时频变换的距离-瞬时多普勒成像方法 |
CN104749573A (zh) * | 2013-12-31 | 2015-07-01 | 中国科学院电子学研究所 | 空频域二维稀疏的步进频率sar成像方法 |
-
2016
- 2016-02-26 CN CN201610108222.4A patent/CN105738897B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010014700A (ja) * | 2008-06-02 | 2010-01-21 | Mitsubishi Electric Corp | 画像レーダ装置 |
KR101103493B1 (ko) * | 2010-07-13 | 2012-01-10 | 영남대학교 산학협력단 | 역합성 개구면 레이더 영상을 이용한 표적 식별 방법 및 그 장치 |
CN102135618A (zh) * | 2010-12-23 | 2011-07-27 | 北京航空航天大学 | 一种针对空域稀疏sar图像的频域压缩感知方法 |
CN102928838A (zh) * | 2012-10-08 | 2013-02-13 | 西北工业大学 | 基于稀疏表示与时频变换的距离-瞬时多普勒成像方法 |
CN104749573A (zh) * | 2013-12-31 | 2015-07-01 | 中国科学院电子学研究所 | 空频域二维稀疏的步进频率sar成像方法 |
Non-Patent Citations (2)
Title |
---|
WANG HONGXIAN ET AL.: "ISAR imaging via sparse frequency-stepped chirp signal", 《SCIENCE CHINA(INFORMATION SCIENCES)》 * |
顾福飞等: "一种孔径和频率二维稀疏的步进频SAR成像方法", 《航空学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106990392A (zh) * | 2017-02-15 | 2017-07-28 | 中国人民解放军国防科学技术大学 | 一种基于随机步进频信号的空间目标微动信息获取方法 |
CN106990392B (zh) * | 2017-02-15 | 2019-06-25 | 中国人民解放军国防科学技术大学 | 一种基于随机步进频信号的空间目标微动信息获取方法 |
CN107037409A (zh) * | 2017-04-28 | 2017-08-11 | 中国人民解放军空军工程大学 | 基于压缩感知的mimo雷达波形分离方法 |
CN107037409B (zh) * | 2017-04-28 | 2021-09-14 | 中国人民解放军空军工程大学 | 基于压缩感知的mimo雷达波形分离方法 |
CN107192991A (zh) * | 2017-07-18 | 2017-09-22 | 中国人民解放军国防科学技术大学 | 一种射频辐射式仿真宽带脉冲雷达成像等效模拟方法 |
CN107192991B (zh) * | 2017-07-18 | 2019-10-11 | 中国人民解放军国防科学技术大学 | 一种射频辐射式仿真宽带脉冲雷达成像等效模拟方法 |
CN109085589A (zh) * | 2018-10-16 | 2018-12-25 | 中国人民解放军国防科技大学 | 基于图像质量引导的稀疏孔径isar成像相位自聚焦方法 |
CN109085589B (zh) * | 2018-10-16 | 2019-04-30 | 中国人民解放军国防科技大学 | 基于图像质量引导的稀疏孔径isar成像相位自聚焦方法 |
CN112130149A (zh) * | 2020-09-24 | 2020-12-25 | 哈尔滨工程大学 | 结合正交匹配追踪算法的稀疏发射的isar成像方法 |
CN113985412A (zh) * | 2021-11-04 | 2022-01-28 | 西安电子科技大学 | 矢量建模优化反演的运动目标三维成像方法 |
CN113985412B (zh) * | 2021-11-04 | 2024-05-14 | 西安电子科技大学 | 矢量建模优化反演的运动目标三维成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105738897B (zh) | 2018-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105738897A (zh) | 基于联合稀疏特性的信号距离像重构方法 | |
CN103698763B (zh) | 基于硬阈值正交匹配追踪的线阵sar稀疏成像方法 | |
CN106772365B (zh) | 一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法 | |
Rao et al. | Adaptive sparse recovery by parametric weighted L $ _ {1} $ minimization for ISAR imaging of uniformly rotating targets | |
Gurbuz et al. | A compressive sensing data acquisition and imaging method for stepped frequency GPRs | |
CN107462887B (zh) | 基于压缩感知的宽幅星载合成孔径雷达成像方法 | |
CN103713288B (zh) | 基于迭代最小化稀疏贝叶斯重构线阵sar成像方法 | |
CN104101864B (zh) | 基于eof分解的导航x波段雷达海浪参数反演算法 | |
CN109061589A (zh) | 随机跳频雷达的目标运动参数估计方法 | |
CN106291524A (zh) | 基于人体模型的lfmcw雷达探测运动人体目标的方法 | |
CN109507666B (zh) | 基于离网变分贝叶斯算法的isar稀疏频带成像方法 | |
CN102854505B (zh) | 一种加权稀疏驱动自聚焦sar成像方法 | |
CN106569191A (zh) | 一种利用高分辨率成像获取目标rcs的方法 | |
CN108226928A (zh) | 基于期望传播算法的逆合成孔径雷达成像方法 | |
CN103235295A (zh) | 基于压缩卡尔曼滤波的小场景雷达目标距离像估计方法 | |
CN109884627B (zh) | 任意线阵构型的近程毫米波快速三维成像方法 | |
CN109738878A (zh) | 基于压缩感知和频率步进波形的雷达一维距离像识别方法 | |
CN110879391B (zh) | 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法 | |
CN105116408A (zh) | 一种舰船isar图像结构特征提取方法 | |
CN106054154A (zh) | 一种基于icpf的机动目标的步进频率信号高分辨率成像方法 | |
CN106707284A (zh) | 一种用于双基地逆合成孔径雷达的成像方法 | |
CN113484859B (zh) | 一种基于融合技术的二维超分辨雷达成像方法 | |
CN104977570B (zh) | 改进基于零空间调整的双通道稀疏sar动目标检测方法 | |
CN107024684B (zh) | 一种空间高速运动目标干涉式三维成像方法 | |
CN110348169B (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 |