基于线性调频信号的超表面天线三维成像方法与流程

本发明属于雷达
技术领域：
，更进一步一种超表面天线三维成像方法，可用于超表面孔径雷达的成像重构。
背景技术：
：近年来，超表面孔径雷达成像受到研究学者以及工程人员越来越多的关注。不同于传统超表面孔径成像雷达基于矢量网络分析仪，采用步进频信号，基于线性调频信号的超表面成像系统采用线性调频的发射机和自混频的接收机，可显著增强对波形的控制能力。随着对超表面天线三维成像方法研究的深入，很多成像方法被提出来，但是一些超表面三维成像方法存在一定的问题，如：大多数压缩感知算法同时对距离、方位和俯仰进行处理，字典维度高，计算量大。yurduseven,o.和gollub,j.n.在其发表的论文"softwarecalibrationofafrequency-diverse,multistatic,computationalimagingsystem"(ieeeaccess,2016,4,pp.2488–2497)中使用矢量网络分析仪与天线进行相连，进行逐频点的激励和回波数据的采集，但是该方法不适用于连续扫频的应用，对波形控制力比较低。zhenhuawu和leizhang,在其发表的论文"rangedecouplingalgorithmforacceleratingmetamaterialapertures-basedcomputationalimaging[j]"(ieeesensorsjournal,2018,18(9):3619-3631)中提出了一种距离解耦合并行处理算法，由于同时对距离、方位和俯仰进行处理，所以计算复杂度较高。且该算法是对场景空间进行并行分块处理，如果不同目标在距离上比较接近，则每一个距离段的点数很少，对应的频点数也很少，由于超材料天线的横向自由度是根据频点个数决定的，因而会减少横向可分辨目标的个数，降低目标分辨率。技术实现要素：本发明的目的在于针对现有技术存在的缺陷，提出一种基于线性调频信号的超表面天线三维成像方法，以增强波形控制能力，减小计算复杂度，提高目标分辨能力。为实现上述目的，本发明的技术方案包括如下步骤：一种基于成像系统的超表面天线三维成像方法，所述成像系统包括超表面天线、喇叭天线和混频器，其特征在于，成像步骤包括如下：(1)超表面天线发射线性调频信号xt对目标进行探测，喇叭天线接收该线性调频信号信号xr并与混频器进行混频，得到目标回波信号yr：其中，f2为目标回波信号的幅度，fc为回波信号的中心频率，μ＝b/t是调频率，b是有效带宽，t是有效时宽，τ为回波延时，τ＝2r/c，r是目标与雷达间的距离，c是光速，为发射信号xt的共轭信号；(2)对目标回波信号进行距离估计和检测，得到频域回波信号：(2a)根据(1)中目标回波信号的相位部分exp(-j4πμtτ)得到频率导向矩阵am×n：其中，表示矩阵中第m行第k列元素，m＝1,2,...,m，m表示信号的长度，k＝1,2,...,n，n表示搜索频点的个数，t表示时间，j表示虚数，τk表示第k个目标的时延；(2b)对目标回波信号yr进行m点均匀采样，得到长度为m的采样序列y：y＝[y1,y2,...ym,...,ym]t其中，ym表示第m个采样值，t表示转置；(2c)根据频率导向矩阵am×k，计算频率导向矩阵的协方差矩阵：其中，为频率导向矩阵的共轭转置，为目标信号初始功率；(2d)计算在τk处的信号其中，k表示搜索频点的个数，为协方差矩阵的逆矩阵，am(τk)为频率导向矩阵的第k列，为频率导向矩阵第k列的共轭转置；(2e)根据(2c)中所求信号求解在τk处的功率：(2f)重复(2c)到(2e)，得到目标信号功率以及目标信号其中，表示τk处目标的功率，表示τk处的目标信号；(2g)设置检测门限值为pf，将目标信号功率与该检测门限值pf进行比较：如果则认为在τk处存在目标，反之，则在τk处不存在目标；(2h)根据(2g)的检测结果，在距离窗内，根据存在目标的τk重构频域回波信号s：其中，sk是在τk处重构的信号，sk＝f2exp(-j2π(fcτk+2μtτk-μτk2)，l表示在距离窗内，存在的目标个数；(3)根据(2)中重构的频域回波信号s，结合超表面天线频率敏感方向图，完成方位、俯仰的二维成像。本发明与现有的技术相比具有以下优点：第一，由于本发明设计的实际成像系统采用的是时宽较大的宽频带信号，即线性调频信号，相比现有的成像系统采用步进频的形式，提高了波形控制能力。第二，本发明由于先对目标回波进行距离估计，再进行方位、俯仰向的二维成像，实现了降维成像，相比于现有的大多数算法，比如三维压缩感知、距离解耦等算法，减小了计算复杂度，提高了目标分辨能力。附图说明图1为本发明的实现流程图；图2为本发明中所用的成像系统原理框图；图3为用本发明对五个标准角反射器的仿真成像图。具体实施方式以下结合附图对本发明的实施例及效果做进一步详细描述本实例是基于成像系统实现，该成像系统如图2所示，其主要由超表面发射天线、喇叭接收天线、混频器、上位机、锁相环、压控振荡器、功放、低噪声放大器和模数转换器等组成。该成像系统完成信号的发射和接收，其中：发射过程，是由上位机控制产生随时间线性增长的的电压信号，经由锁相环pll、压控振荡器vco调制锁定，产生高频调制信号，再经功率放大器pa进行功率放大，最后通过超表面天线发射出去，同时产生一参考信号进入混频器。接收过程，是由目标回波信号经由喇叭天线接收，然后通过低噪声放大器lna进行功率放大，在混频器内与参考信号混频后，经模数转换器a/d采样后，传输到上位机。该成像系统的参数如表1。表1工作带宽33-37ghz天线面板尺寸250mm*250mm谐振单元数目125*125频率采样间隔10mhz距离分辨率4cm方位角分辨率0.4°俯仰向分辨率0.4°参照图1，本实例的实现步骤如下：步骤1，成像系统的超表面天线发射线性调频信号xt对目标进行探测，喇叭天线接收该线性调频信号xr并与混频器进行混频，得到目标回波信号yr：(1.1)超表面天线发射线性调频信号xt，喇叭天线接收线性调频信号xr：其中，f为线性调频信号的幅度，fc为信号的中心频率，μ＝b/t是调频率，b是有效带宽，t是有效时宽，τ为回波延时，τ＝2r/c，r是目标与雷达间的距离，c是光速，j表示虚数；(1.2)对喇叭天线接收到线性调频信号xr作解线频调处理，即与超表面发射天线发射信号xt的共轭信号xt*作混频处理，得到目标回波信号yr：步骤2，对目标回波信号yr作距离估计和检测，得到目标信号功率和目标信号(2.1)将目标回波信号yr展开，进一步表示为：yr＝f2exp(-j2πfcτ)·exp(-j4πμtτ)·exp(j2πμτ2)<4>由于目标距离比较近，由τ＝2r/c可知，exp(j2πμτ2)可忽略不计，故可将式<4>中exp(-j4πμtτ)作为目标回波信号yr的相位部分，f2exp(-j2πfcτ)作为目标回波信号yr的幅度部分；(2.2)根据目标回波信号yr的相位部分exp(-j4πμtτ)，可得到频率导向矩阵am×n：其中，表示矩阵中第m行第k列元素，m＝1,2,...,m，m表示信号的长度，k＝1,2,...,n，n表示搜索频点的个数，t表示时间，j表示虚数，τk表示第k个目标的时延；(2.3)对式<4>的目标回波信号yr进行均匀采样，得到长度为m的采样序列y：y＝[y1,y2,...ym,...,ym]t<6>其中，ym表示第m个采样值，t表示转置；(2.4)根据频率导向矩阵am×n，计算频率导向矩阵的协方差矩阵rm×m：其中，为频率导向矩阵的共轭转置，为目标信号初始功率，表示如下：其中，为第k个目标信号初始功率，τk表示第k个目标的时延,y(n)为长度为m的回波信号中的第n个值，am(τk)为频率导向矩阵的第k列，为频率导向矩阵第k列的共轭转置，∑|·|2表示对模的平方求和；(2.5)计算在τk处的信号其中，为协方差矩阵的逆矩阵，am(τk)为频率导向矩阵am×n的第k列，为频率导向矩阵am×n的第k列的共轭转置；(2.6)根据(2.5)中所求信号求解在τk处的功率：(2.7)重复(2.3)到(2.5)，得到目标信号功率以及目标信号其中，表示τk处目标的功率。步骤3，根据得到目标信号功率以及目标信号得到重构的频域回避。(3.1)设置检测门限值为pf，将目标信号功率与该检测门限值pf进行比较：如果则认为在τk处存在目标，反之，则在τk处不存在目标；(3.2)根据(3.1)的检测结果，在距离窗内，根据存在目标的τk重构频域回波信号s：其中，sk是在τk处重构的信号，sk＝f2exp(-j2π(fcτk+2μtτk-μτk2)，l表示在距离窗内，存在的目标个数；步骤4，根据重构的频域回波信号s，结合超表面天线频率敏感方向图，进行方位、俯仰二维成像。(4.1)测试超表面天线频率敏感方向图矩阵r，并与重构频域信号s一起构成如下代价函数x(λ)：其中，表示求解l2范数的平方，||·||1表示求解l1范数，λ为正则化参数，x是待重构的目标信号；(4.2)求解(4.1)代价函数的最小值：(4.2.1)从t＝1开始迭代，设x的初始稀疏迭代估计设x的初始残差z0＝s；(4.2.2)计算x的第t次非稀疏噪声估计其中，rh为超表面天线频率敏感方向图矩阵r的共轭转置，是x的第t-1次稀疏迭代估计；x的第t次稀疏迭代估计为：其中，是软阈值函数，表示的相位角，zt-1表示x的第t-1次迭代残差；x的第t次残差迭代为：其中，ηr和ηi分别是软阈值函数的实部和虚部，是ηr关于xr的偏导，是ηi关于xi的偏导，为欠采样因子，<·>表示求解平均值；(4.2.3)设误差阈值tol＝10-5，令t＝t+1，重复(4.2.2)直到输出x的第t次稀疏迭代估计其中，||·||2表示求解l2范数。本发明的优点可通过以下仿真实验进一步说明：1.实验条件本发明的仿真实验的硬件平台是：超表面天线雷达成像系统，matlabr2017a。2.实验内容将五个标准角反射器以间距4cm放置于成像系统正前方1.2m处，利用超表面成像系统对其进行探测，再用本发明方法对其进行成像，结果如图3所示，其中图3a表示实验场景图，图3b表示目标的距离估计像，图3c表示重构的频域回波，图3d表示最终的成像结果。由图3d可见，本发明能将全部五个目标均重构出来，即在超表面天线距离分辨率为4cm、方位角和俯仰角的分辨率均为0.4°的情况下，实现了超分辨，提高了目标分辨能力。当前第1页12