一种存在阵元幅相误差时的微波关联成像方法技术

本发明专利技术提供一种存在阵元幅相误差时的微波关联成像方法。技术方案包括以下步骤：第一步：接收雷达回波，设微波关联成像系统包括N个发射阵元和1个接收阵元，发射阵元同时发射一组相互独立的发射信号；第二步：划分网格，对成像区域进行均匀网格划分，网格大小由成像分辨率决定；第三步：推演辐射场，利用迭代方法计算辐射场参考信号；第四步：目标重构，上述迭代结计算目标散射系数矢量；第五步：估计幅相误差；第六步：如果迭代结果满足要求，利用目标散射系数矢量进行成像。本发明专利技术在存在幅相误差时可以对目标进行高精度成像。

A method of microwave correlation imaging in the presence of phase and amplitude errors of array elements

The invention provides a microwave correlation imaging method when the amplitude and phase errors of an array element exist. The technical scheme comprises the following steps: the first step: receiving a microwave radar echo, correlation imaging system consists of N emission array and 1 receiving array, signal transmitting array element and transmitting a group of independent; the second step: mesh, uniform mesh on the imaging area, the grid size is determined by imaging resolution; the third step: deduction of radiation field, radiation field calculation reference signal using iterative method; the fourth step: target reconstruction, the iteration calculation of scattering coefficient vector; the fifth step: the estimation of amplitude and phase errors; the sixth step: if the iterative results meet the requirements of imaging by scattering coefficient vector. The invention can carry out high-precision imaging of the target in the presence of amplitude phase errors.

【技术实现步骤摘要】
一种存在阵元幅相误差时的微波关联成像方法
本专利技术涉及微波成像
，特别是微波关联成像方法，更为具体地涉及提高存在阵元幅相误差时的微波关联成像质量。
技术介绍
微波关联成像(DongzeLi,XiangLi,YuliangQin,YongqiangCheng,andHongqiangWang,“RadarCoincidenceImaging:AnInstantaneousImagingTechniqueWithStochasticSignals,”IEEETransactionsonGeoscienceandRemoteSensing,vol.52,no.4,Apr.2014)作为一种新的凝视成像方法，不依赖于雷达与目标的相对运动，具有高分辨、抗截获、抗干扰等优势，在静止或准静止平台凝视成像、灾情监测、海洋监视、高分辨对地观测等应用领域具有广泛的应用前景。微波关联成像借鉴经典的光学关联成像的原理，通过对发射信号波前的调制，构造在时间和空间上随机分布的二维随机辐射场来模拟具有随机涨落的光场分布。然后将目标散射回波与二维随机辐射场进行关联处理，完成对波束内目标信息的提取与解耦，实现微波关联成像。相比于现有的微波成像方法，微波关联成像不依赖于目标的多普勒信息进行分辨，大大降低了对目标相对运动的要求，同时微波关联成像能够利用雷达系统的一次发射对目标成像，极短的成像时间大大降低了目标的非合作运动对成像质量的影响。微波关联成像在处理相对静止目标和非合作运动目标方面具有极大的发展潜力。在微波关联成像中，将雷达作用的区域称为成像区域，成像时，成像区域首先被划分为若干个网格，并假设目标由若干散射中心组成且位于网格中心。当不存在阵元幅相误差等模型误差时，利用参数化的微波关联成像方法能够获得高质量的目标图像。但是在实际微波凝视成像雷达系统中，各阵元之间的幅相误差广泛存在，成像过程被引入误差，降低了目标成像精度。因此，研究存在阵列幅相误差时的微波关联成像方法对于微波关联成像在实际场景中的应用具有重要的意义和价值。针对幅相误差，人们提出了多种解决方法，例如测角中经常用到的基于特征值分解的校正方法(AifeiLiu,GuishengLiao,CaoZeng,ZhiweiYang,andQing,“AnEigenstructureMethodforEstimatingDOAandSensorGain-PhaseErrors,”IEEETransactionsonSignalProcessing,vol.59,no.12,Dec.2011)，这类方法对幅相误差不敏感，但是在低信噪比、观测数据有限、目标距离较近时性能较差。
技术实现思路
本专利技术针对实际微波关联成像系统中普遍存在的阵列幅相误差导致目标成像精度降低的问题，提出了一种存在阵元幅相误差时的自聚焦微波关联成像方法。本方法原理简单，计算量小，对噪声不敏感，可以有效解决存在阵元幅相误差下的成像问题，直接应用于实际的微波关联成像场景。本专利技术的基本思路是，将成像过程分为两个步骤：(1)目标重构，估计目标散射中心的位置和散射强度；(2)失配误差估计，采用牛顿法精确估计阵元幅相误差，并进行补偿，从而提高成像精度和质量。本专利技术的技术方案是：一种存在阵元幅相误差时的微波关联成像方法。具体包括以下步骤：第一步：接收雷达回波设微波关联成像系统包括N个发射阵元和1个接收阵元，发射阵元在t时刻同时发射一组相互独立的发射信号Sn(t)，n＝1,2,…,N。每个发射信号Sn(t)在时域上独立分布，接收阵元接收的目标散射回波y(t)。从y(t)中均匀抽取M个数据获得离散序列得到目标散射回波矢量其中为时间采样点，M的取值根据实际情况确定。第二步：划分网格对成像区域进行均匀网格划分，设在方位向和距离向上的网格数分别为Kx、Ky，总的网格数为K＝Kx×Ky，网格大小由成像分辨率决定。各网格中心的位置矢量构成集合初始化，令迭代次数i＝1，幅度误差初始值a0＝[1,…,1]T，相位误差初始值第三步：推演辐射场设第i-1次迭代计算的幅度误差矢量为第i-1次迭代计算的相位误差矢量为利用公式一计算第i次迭代第k个网格中心tm时刻对应的辐射场参考信号其中，和分别为第i-1次迭代得到的第n个发射阵元与接收阵元之间的幅度误差和相位误差，和分别为第n个发射阵元和接收阵元的位置矢量，表示两个矢量之间的欧氏距离，c为电磁波的传播速度。则第i次迭代的参考信号矩阵即：第四步：目标重构利用下式计算第i次迭代的目标散射系数矢量上式中，x表示待求解的目标散射系数矢量，||·||2表示2范数，||·||1表示1范数；λ是正则化参数，根据实际情况确定。第五步：估计幅相误差利用下式计算第i次迭代的代价函数fi：利用公式五和公式六分别求解幅度误差和相位误差(即牛顿法)和分别表示第i次迭代的代价函数对幅度误差求解的梯度和Hessian矩阵，和分别表示第i次迭代的代价函数对相位误差求解的梯度和Hessian矩阵。第六步：判断收敛如果小于设定的阈值η，或迭代次数i达到最大迭代次数，则迭代结束，利用第i次迭代的目标散射系数矢量进行成像，得到成像结果。否则令i＝i+1，返回第三步。本专利技术的有益效果是：利用目标散射中心空域分布的稀疏性，将存在阵元幅相误差时的成像建模为目标散射系数、幅度误差和相位误差的联合优化问题。首先将与目标散射系数相对应的代价函数最小化来估计目标散射系数(即公式三)，其次通过牛顿法来估计幅度误差和相位误差(即公式五与公式六)，并更新参考信号矩阵。该方法通过交替迭代估计目标散射系数、幅度误差和相位误差，直到方法收敛。其中第四步所述计算目标散射系数矢量的方法是一个不需要信号稀疏度和噪声水平的自适应算法，且不需用收敛准则来终止迭代条件，在大尺度数据中明显比SBL算法快，并在低信噪比的条件下获得了优于SBL算法的精度。而第五步所述的牛顿法无须对模型进行近似，可以较准确地估计幅相误差。所以本文所述的方法在存在幅相误差时可以对目标进行高精度成像。附图说明图1是本专利技术所述的微波关联成像雷达成像几何示意图；图2是本专利技术所述的存在阵元幅相误差时的微波关联成像方法流程图；图3至图6是本专利技术所述的方法的成像结果；图7至图10是本专利技术所述的方法对幅相误差的估计效果和收敛性能；图11是本专利技术所述的方法在不同信噪比下的相对成像误差；图12至图13是本专利技术所述的方法在不同信噪比下对幅相误差的估计误差。具体实施方式下面结合附图对本专利技术所述的存在阵元幅相误差时的微波关联成像方法进行详细说明。图1是本专利技术所述的微波关联成像雷达成像几何示意图。图中雷达阵列包括N个发射阵元和1个接收阵元，发射阵元发射一组相互独立的信号Sn(t)，n＝1,2,…,N，Sn(t)是在时域上呈独立分布的随机信号。各发射阵元的发射信号在成像区域处形成在时间和空间上随机分布的二维随机辐射场，经过目标调制的随机辐射场被接收阵元接收得到回波y(t)。成像时，成像区域在方位向和距离向上被划分为若干个相同大小的网格单元，目标散射中心均位于网格中心，图中的黑点所示。图2是本专利技术所述的存在阵元幅相误差时的微波关联成像方法流程图。雷达首先发射随机信号，在成像区域形成随机辐射场，同时接收目标散射回波。在进行成像前本文档来自技高网...

【技术保护点】
一种存在阵元幅相误差时的微波关联成像方法，其特征在于，包括以下步骤：第一步：接收雷达回波：设微波关联成像系统包括N个发射阵元和1个接收阵元，发射阵元在t时刻同时发射一组相互独立的发射信号S

【技术特征摘要】
1.一种存在阵元幅相误差时的微波关联成像方法，其特征在于，包括以下步骤：第一步：接收雷达回波：设微波关联成像系统包括N个发射阵元和1个接收阵元，发射阵元在t时刻同时发射一组相互独立的发射信号Sn(t)，n＝1,2,…,N；每个发射信号Sn(t)在时域上独立分布，接收阵元接收的目标散射回波y(t)；从y(t)中均匀抽取M个数据获得离散序列得到目标散射回波矢量其中为时间采样点，N和M的取值根据实际情况确定；第二步：划分网格：对成像区域进行均匀网格划分，设在方位向和距离向上的网格数分别为Kx、Ky，总的网格数为K＝Kx×Ky，网格大小由成像分辨率决定；各网格中心的位置矢量构成集合初始化，即令迭代次数i＝1，幅度误差初始值a0＝[1,…,1]T，相位误差初始值第三步：推演辐射场：设第i-1次迭代计算的幅度误差矢量为第i-1次迭代计算的相位误差矢量为利用公式一计算第i次迭代第k个网格中心tm时刻对应的辐射场参考信号其中，和分别为第i-1...

【专利技术属性】
技术研发人员：秦玉亮，程永强，周小利，王宏强，范波，李彦鹏，曹凯程，
申请(专利权)人：中国人民解放军国防科学技术大学，
类型：发明
国别省市：湖南,43