基于NSL0重构算法的THz-FMCW雷达成像算法制造技术

本发明专利技术提供基于NSL0重构算法的太赫兹调频连续波雷达成像算法，步骤1、设计一测量矩阵，该测量矩阵确定THz FMCW雷达在扫描平面内的测量位置点；步骤2、根据步骤1中设计的测量矩阵，在对应的空间位置进行THz FMCW雷达回波信号采集，得到测量信号；步骤3、对步骤2中得到的测量信号采用NSL0重构算法进行重构，得到变换域稀疏信号；步骤4、对步骤3中重构得到的变换域稀疏信号进行处理，实现原始信号的重构；步骤5、将重构得到的原始信号进行三维成像处理。本发明专利技术中采用的NSL0重构算法，可克服SL0重构算法存在的最优解索路径呈“锯齿状”、函数收敛速度降慢等缺点，具有重构速度快、成像时间短等优点。

【技术实现步骤摘要】
基于NSL0重构算法的THz-FMCW雷达成像算法
本专利技术涉及太赫兹调频连续波雷达成像算法领域，特别是涉及基于NSL0重构算法的THz-FMCW雷达成像算法。
技术介绍
太赫兹波(THz)是指频率在0.1THz—10THz之间的电磁波，是一个尚未被人类完全认知和利用的电磁波段。THz雷达具有大的工作带宽，能够实现高分辨率成像，同时可以穿透其他频段的电磁波不能穿透的物质，适用于无损检测、安检等一系列实际应用当中。压缩感知理论利用信号的稀疏性，通过较少的测量数据即可准确恢复出原始信号。在雷达成像领域，尽管THz成像技术可极大提高成像分辨率，但同时对采样率、方位采样间隔提出了更高要求，由此带来数据获取时间长、数据传输困难等问题。依据压缩感知理论，如果雷达场景稀疏或者可通过某种变换稀疏表示，则可通过采集远低于传统奈奎斯特采样方法要求的数据量进行场景重构，即利用较少的采样数据来实现对原始完整信号的近似逼近，获得较理想的恢复重建。由于压缩采样方法的采用，压缩感知雷达成像技术可以降低系统采样率、增大测绘带宽度、降低采样数据量。所以本专利技术中将压缩感知理论应用到THz雷达成像中以减少数据采集样本，进而提高成像速度。目前已经有了一系列基于压缩感知的THz雷达成像算法，比如“基于压缩感知的太赫兹调频连续波雷达成像算法研究”论文中研究了基于SL0压缩感知重构算法的THz-FMCW雷达成像算法。专利方面，专利号CN201710782421.8，专利名称基于压缩感知重构的太赫兹雷达成像算法则采用2D-SL0重构算法尽管目前已有不少基于压缩感知的雷达成像方法的专利，但对于太赫兹调频连续波成像雷达，其基于压缩感知的成像算法较少，比较相关的是本研究团队的“基于压缩感知重构的太赫兹雷达成像算法”、“基于压缩感知的太赫兹调频连续波雷达三维成像算法”等专利，及“基于压缩感知的太赫兹调频连续波雷达成像算法研究”等论文。与对比文献相比，本申请专利主要是采用了NSL0重构算法。由于SL0重构算法在迭代的过程中采用的是最速下降法，沿着最速下降方向搜索最优解，其搜索路径呈“锯齿状”，导致函数收敛速度降低。基于NSL0重构算法的太赫兹调频连续波雷达成像算法，可以克服SL0重构算法中函数收敛速度降低的问题，具有重构速度快、成像时间短等优点。
技术实现思路
为了解决以上问题，本专利技术提供基于NSL0重构算法的THz-FMCW雷达成像算法，本专利技术中采用的NSL0重构算法，可克服SL0重构算法存在的最优解索路径呈“锯齿状”、函数收敛速度降慢等缺点，具有重构速度快、成像时间短等优点，为达此目的，本专利技术提供基于NSL0重构算法的THz-FMCW雷达成像算法，具体步骤如下，其特征在于：步骤1、设计一测量矩阵，该测量矩阵确定THzFMCW雷达在扫描平面内的测量位置点；步骤2、根据步骤1中设计的测量矩阵，在对应的空间位置进行THzFMCW雷达回波信号采集，得到测量信号y；步骤3、对步骤2中得到的测量信号采用NSL0重构算法进行重构，得到变换域稀疏信号；步骤4、对步骤3中重构得到的变换域稀疏信号进行处理，实现原始信号的重构；步骤5、将重构得到的原始信号进行三维成像处理。本专利技术的进一步改进，步骤3中由较少的测量信号样本y重构原始信号g，然而g是不稀疏的，由y直接重构出g的过程不能成立，此时，将原始信号进行稀疏表示：g＝Ψx(1)式中，x为原始信号g在Ψ变换域的稀疏表示，Ψ表示稀疏变换基，本专利技术中采用的稀疏变换基为离散余弦变换基。考虑到测量测量公式y＝Φg，则测量信号y与变换域稀疏信号x之间存在着如下的关系：y＝ΦΨx＝Ax(2)式中，A＝ΦΨ为感知矩阵，表示测量矩阵和稀疏变换基的乘积；将y看作是稀疏信号x关于“测量矩阵”A的测量值，通过求解一个类似最优l0范数问题来重构稀疏信号x，即本专利技术的进一步改进，步骤3中基于NSL0的压缩感知重构算法总结如下；输入：(1)M×N的感知矩阵A，其中M＜＜N，A＝ΦΨ；(2)长度为M的测量值y；输出：长度为N的重构稀疏信号满足初始化各参数：(1)x0＝AT(AAT)-1y；(2)选择合适的σ递减序列{σ1,σ2,...,σj},σj＝βσj-1，β(0＜β＜1)为递减因子；算法步骤：(1)进行外循环：j＝1,2,3,...,JA、σ＝σj；B、xj＝xj-1；C、r0＝0；D、内部循环：n＝1,2,3,...,Na)计算修正后的牛顿方向b)对重构信号进行更新c)根据梯度投影原理，得到d)计算余量，e)判断余量的差值，若||r-r0||＜e则退出当前σ值下的内循环，否则r0＝r；E、(2)得到重构稀疏信号本专利技术的进一步改进，步骤4中将重构的稀疏信号乘上稀疏变换基Ψ进行稀疏反变换得到重构的原始信号本专利技术公开了一种基于NSL0重构算法的太赫兹调频连续波雷达成像算法，该算法首先设计好测量矩阵并完成THzFMCW雷达回波信号采集；采用NSL0重构算法对测量雷达回波信号进行重构，对得到的变换域稀疏信号进行处理实现原始信号重构；将重构信号进行三维成像处理实现目标的三维成像。本专利技术与SL0重构算法相比，具有收敛速度快，成像效率高等优点。附图说明图1为基于压缩感知的THz-FMCW雷达三维成像算法的流程图。图2为一具体实施例中稀疏率为20％时重构的剪刀图像。图3为一具体实施例中稀疏率为30％时重构的剪刀图像。图4为一具体实施例中稀疏率为50％时重构的剪刀图像。图5为不同稀疏率条件下，SL0算法与NSL0算法的成像重构时间对比曲线图。具体实施方式下面结合附图与具体实施方式对本专利技术作进一步详细描述：本专利技术提供基于NSL0重构算法的THz-FMCW雷达成像算法，本专利技术中采用的NSL0重构算法，可克服SL0重构算法存在的最优解索路径呈“锯齿状”、函数收敛速度降慢等缺点。本实施实例将通过THz-FMCW雷达成像系统的实际测量数据对本专利技术中提出的算法进行验证说明。为验证本专利技术中提出的基于压缩感知的THz-FMCW雷达三维成像算法，利用THz-FMCW雷达三维成像系统进行实测成像实验，成像目标为距离扫描平面0.25m处的剪刀。成像系统的工作频率范围为170GHz-220GHz，扫描平面的扫描范围为151mm*151mm。基于NSL0的THzFMCW雷达压缩感知成像重构算法，算法流程如图1所示，具体操作步骤如下：步骤1、设计一测量矩阵Φ，该测量矩阵确定THzFMCW雷达在扫描平面内的测量位置点；步骤2、根据步骤1中设计的测量矩阵，在对应的空间位置进行THzFMCW雷达回波信号采集，得到测量信号y；步骤3、由较少的测量信号样本y重构原始信号g，然而g是不稀疏的，由y直接重构出g的过程不能成立。此时，将原始信号进行稀疏表示：g＝Ψx(1)式中，x为原始信号g在Ψ变换域的稀疏表示，Ψ表示稀疏变换基，本专利技术中采用的稀疏变换基为离散余弦变换基。考虑到测量测量公式y＝Φg，则测量信号y与变换域稀疏信号x之间存在着如下的关系：y＝ΦΨx＝Ax(2)式中，A＝ΦΨ为感知矩阵，表示测量矩阵和稀疏变换基的乘积；将y看作是稀疏信号x关于“测量矩阵”A的测量值，通过求解一个类似最优l0范数问题来重构稀疏信号x，即基于NSL0的压缩感知重构算法总结如下。输入：(1)M×N的感知矩阵A，其中本文档来自技高网...

【技术保护点】
1.基于NSL0重构算法的THz‑FMCW雷达成像算法，具体步骤如下，其特征在于：步骤1、设计一测量矩阵，该测量矩阵确定THz FMCW雷达在扫描平面内的测量位置点；步骤2、根据步骤1中设计的测量矩阵，在对应的空间位置进行THz FMCW雷达回波信号采集，得到测量信号y；步骤3、对步骤2中得到的测量信号采用NSL0重构算法进行重构，得到变换域稀疏信号；步骤4、对步骤3中重构得到的变换域稀疏信号进行处理，实现原始信号的重构；步骤5、将重构得到的原始信号进行三维成像处理。

【技术特征摘要】
1.基于NSL0重构算法的THz-FMCW雷达成像算法，具体步骤如下，其特征在于：步骤1、设计一测量矩阵，该测量矩阵确定THzFMCW雷达在扫描平面内的测量位置点；步骤2、根据步骤1中设计的测量矩阵，在对应的空间位置进行THzFMCW雷达回波信号采集，得到测量信号y；步骤3、对步骤2中得到的测量信号采用NSL0重构算法进行重构，得到变换域稀疏信号；步骤4、对步骤3中重构得到的变换域稀疏信号进行处理，实现原始信号的重构；步骤5、将重构得到的原始信号进行三维成像处理。2.根据权利要求1所述的基于NSL0重构算法的THz-FMCW雷达成像算法，其特征在于：步骤3中由较少的测量信号样本y重构原始信号g，然而g是不稀疏的，由y直接重构出g的过程不能成立，此时，将原始信号进行稀疏表示：g＝Ψx(1)式中，x为原始信号g在Ψ变换域的稀疏表示，Ψ表示稀疏变换基，本发明中采用的稀疏变换基为离散余弦变换基。考虑到测量测量公式y＝Φg，则测量信号y与变换域稀疏信号x之间存在着如下的关系：y＝ΦΨx＝Ax(2)式中，A＝ΦΨ为感知矩阵，表示测量矩阵和稀疏变换基的乘积；将y看作是稀疏信号...

【专利技术属性】
技术研发人员：顾姗姗，杨忠，葛玲玉，司海飞，陈维娜，
申请(专利权)人：金陵科技学院，
类型：发明
国别省市：江苏,32