本发明专利技术公开了一种基于杂波的多输入多输出雷达阵列误差校正方法，主要解决现有技术在超表面成像系统波形控制差、成像处理过程计算量大和分辨率低的问题。其实现方案是：1)根据超表面成像系统获取三维线性调频回波信号；2)对回波信号进行迭代求解作距离估计，得到目标信号的功率；3)设置功率检测门限，对目标信号功率进行检测，重构频域回波信号；4)根据重构的频域回波信号，结合超表面频率敏感方向图矩阵，完成方位、俯仰向二维成像。本发明专利技术提高了波形控制能力，减小了计算复杂度，提高了目标分辨能力，可用于超表面孔径雷达的成像重构。
3D imaging method of hypersurface antenna based on LFM signal
【技术实现步骤摘要】
基于线性调频信号的超表面天线三维成像方法
本专利技术属于雷达
，更进一步一种超表面天线三维成像方法，可用于超表面孔径雷达的成像重构。
技术介绍
近年来，超表面孔径雷达成像受到研究学者以及工程人员越来越多的关注。不同于传统超表面孔径成像雷达基于矢量网络分析仪，采用步进频信号，基于线性调频信号的超表面成像系统采用线性调频的发射机和自混频的接收机，可显著增强对波形的控制能力。随着对超表面天线三维成像方法研究的深入，很多成像方法被提出来，但是一些超表面三维成像方法存在一定的问题，如：大多数压缩感知算法同时对距离、方位和俯仰进行处理，字典维度高，计算量大。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)对喇叭本文档来自技高网...
【技术保护点】
1.一种基于成像系统的超表面天线三维成像方法，所述成像系统包括超表面天线、喇叭天线和混频器，其特征在于，成像步骤包括如下：/n(1)超表面天线发射线性调频信号x
【技术特征摘要】
1.一种基于成像系统的超表面天线三维成像方法，所述成像系统包括超表面天线、喇叭天线和混频器，其特征在于，成像步骤包括如下：
(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的采样序列yM：
yM＝[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)中重构的频域回波...
【专利技术属性】
技术研发人员：戴奉周，冯高俊，刘宏伟，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。