基于压缩感知重构的太赫兹雷达成像算法的制作方法

本发明涉及太赫兹调频连续波雷达成像技术，特别是一种基于压缩感知重构的太赫兹雷达成像算法。

背景技术：

太赫兹波(thz)是指频率在0.1thz—10thz之间的电磁波，是最后一个尚未被人类完全认知和利用的电磁波段。thz雷达具有大的带宽，能够实现高分辨率成像，同时可以穿透其他频段的电磁波不能穿透的物质，适用于无损检测、安检等一系列实际应用当中。

压缩感知理论利用信号的稀疏性，通过较少的测量数据即可准确恢复出原始信号。在雷达成像领域，尽管thz成像技术可极大提高成像分辨率，但同时对采样率、方位采样间隔提出了更高要求，由此带来数据获取时间长、数据传输困难等问题。依据压缩感知理论，如果雷达场景稀疏或者可通过某种变换稀疏表示，则可通过采集远低于传统奈奎斯特采样方法要求的数据量进行场景重构，即利用较少的采样数据来实现对原始完整信号的近似逼近，获得较理想的恢复重建。由于压缩采样方法的采用，压缩感知雷达成像技术可以降低系统采样率、增大测绘带宽度、降低采样数据量。所以本发明中将压缩感知理论应用到thz雷达成像中以减少采样的数据量。

目前已经有了一些压缩感知应用于thz雷达成像的算法，“基于压缩感知的太赫兹调频连续波雷达成像算法研究”一文中提出了基于sl0压缩感知重构算法的thz-fmcw雷达成像算法。该成像算法是将压缩感知中的多维优化问题转化为一维优化问题进行研究，但是这样处理存在一个问题：对于大场景目标，测量矩阵尺寸会过于庞大，需要占用大量存储空间，增加了成像处理时间。

技术实现要素：

本发明的目的在于提供一种基于压缩感知重构的太赫兹雷达成像算法。

实现本发明目的的技术方案为：一种基于压缩感知重构的太赫兹雷达成像算法，包括以下步骤：

步骤1，根据回波信号的推导公式分离出两个隐藏式的测量矩阵；

步骤2，采用2d-sl0重构算法解决最优化l0范数问题，重构出原始信号。

与现有技术相比，本发明的显著优点为：

本发明采用隐式的测量矩阵，同时提出改进的sl0重构算法——2d-sl0算法，来减小存储空间，减少成像时间，提高成像效率。

附图说明

图1为本发明雷达成像扫描方式示意图。

图2为本发明实施例中2d-sl0算法重构的点目标二维图像。

图3为本发明实施例中sl0算法重构的点目标二维图像。

具体实施方式

一种基于压缩感知重构的太赫兹雷达成像算法，包括以下步骤：

步骤1，根据回波信号的推导公式分离出两个隐藏式的测量矩阵，得到最优化l0范数问题；

步骤2，采用2d-sl0重构算法解决最优化l0范数问题，重构出原始信号。

进一步的，步骤1具体为：

如图1所示，雷达成像采用平面扫描方式，雷达发射天线在x′o′y′扫描平面发射单频thz信号，接收天线将目标反射的回波信号纪录下来，组成回波数据矩阵；

扫描平面内测量点(x′,y′)接收的回波信号为：

式中，k＝2πf/c，f为雷达信号频率，c为光速，σ为目标的散射系数，z0为扫描平面和成像目标之间的距离；(x,y)为目标平面成像点；

在满足远小于1这一条件时，2kr被表示为：

a为扫描平面长度，l为成像目标区域长度，λ为信号波长；

将式(2)代入式(1)可得：

测量目标是稀疏的点目标，其压缩感知算法是通过求解最优化l0范数问题来计算出σ：

min||σ||0s.t.φσ＝g(4)

式中，g表示测量信号，φ为测量矩阵；选取测量矩阵为：

由式(5)可知φ是可分离的，将测量矩阵分离为两个矩阵：

因此，式(4)的最优化l0范数问题表示为：

min||σ||0s.t.φ1σφ2＝g(8)

进一步的，最优化l0范数问题是一个np-hard问题，不能直接计算，需要寻找近似方法来求解。2d-sl0重构算法是用光滑的高斯函数来逼近离散的l0范数，从而将离散函数的最优化问题转化为连续函数的最优化问题，通过凸优化的方法对其求解，在迭代过程中采用最速下降法和梯度投影原理，经过多次迭代逐步逼近最优解。2d-sl0重构算法的具体步骤如下：

(1)初始化算法中的各个参数：

1)

2)选择外循环迭代值α以及确定算法终止的界限αmin；

(2)当α＞αmin时进行以下循环：

1)进行l次迭代，l＝1,2,…,l：

a，令

b，对重构信号进行更新s←s-μα²d，μ为常数；

c，根据梯度投影原理，得到

2)α←βα，其中β为递减因子；

(3)得到重构信号

下面结合具体实施例对本发明进行详细说明。

实施例

本实施例通过对点目标的太赫兹雷达成像仿真阐述基于二维sl0压缩感知重构算法的太赫兹雷达成像算法。

如图1所示，假设目标的参考坐标系为xyz，三个待成像点目标在xy平面内的坐标中的位置分别为(-0.02m,-0.02m)、(0m,0m)、(0.02m,0.02m)，雷达成像时扫描平面为x′y′，扫描平面与待测目标平面的距离z0为1.5m。雷达发射的单频信号频率为300ghz，设扫描平面x′y′的范围为100mm×100mm，现随机选取50mm×50mm的范围来接收回波数据，此时雷达回波信号可以表示为：

式中，k＝2πf/c，f表示雷达信号频率，c表示光速，σ表示目标的散射系数。

本实施例中扫描平面长度a为100mm，信号波长λ为1mm，假设成像目标区域长度l为0.1m，则满足了远小于1这一条件，那么2kr可以被表示为：

则g(x′,y′)可以表示为：

待成像目标是三个稀疏的点目标，因此基于压缩感知的雷达成像算法可归结为通过求解最优化l0范数问题计算σ：

min||σ||0s.t.φσ＝g(12)

式中，g表示测量信号，φ为测量矩阵。选取测量矩阵为：

由式(13)可以看出φ是可分离的，本发明中将测量矩阵分离为两个矩阵：

至此，式(12)的最优化l0范数问题又可以表示为：

min||σ||0s.t.φ1σφ2＝g(16)

下面通过2d-sl0算法求解上述最优化l0范数问题，具体步骤如下：

(1)初始化算法中的各个参数：

1)

2)选择一个合适的外循环迭代值α，以及确定算法终止的界限αmin；

(2)当α＞αmin时进行以下循环：

1)进行l次迭代，l＝1,2,…,l：

a.令

b.对重构信号进行更新s←s-μα²d，μ为常数，本实施例中取值为2；

c.根据梯度投影原理，得到

2)α←βα，其中β为递减因子，取值在0到1之间，本实施例中取0.8；

(3)得到重构信号

至此完成了由测量信号g重构目标散射系数σ的全部过程，接下来进行二维成像显示即完成全部算法过程。图2为2d-sl0算法重构的点目标二维图像。

为了验证本发明中2d-sl0算法在减少成像时间方面的优势，本实施例给出了sl0算法重构的点目标二维图像，如图3所示。本实例仿真时，2d-sl0重构算法耗时0.3588s，sl0重构算法耗时2.2152s，2d-sl0重构算法在减少成像时间方面的优势显著。