申请/专利权人：电子科技大学

申请日：2023-08-31

公开（公告）日：2023-12-05

公开（公告）号：CN117169881A

主分类号：G01S13/89

分类号：G01S13/89;G06F17/15;G06F17/16

优先权：

专利状态码：在审-实质审查的生效

法律状态：2023.12.22#实质审查的生效;2023.12.05#公开

摘要：本发明公开了一种基于全变差约束的扫描雷达递归快速超分辨成像方法，首先将方位回波建模为目标散射分布与天线辐射函数的卷积，然后在正则化框架下，选择全变差范数作为约束导出优化问题，最后将矩阵求逆运算转化为矩阵乘法的迭代运算，利用实时的扫描回波有效地递归更新重建结果，实现前视区域目标实时递归更新。本发明的方法与现有批处理方法相比，在不损失成像性能的同时显著降低计算复杂度，节省运算时间，适合扫描雷达的快速成像，与现有TV正则化方法相比，则可以随着扫描回波的不断更新而加强整体的重建结果，提高重建结果的鲁棒性，同时迭代更新过程中的各个矩阵维度是固定不变的，仅需根据当前波束回波去更新，提供了恒定的计算和存储成本。

主权项：1.一种基于全变差约束的扫描雷达递归快速超分辨成像方法，具体步骤如下：步骤一、机载扫描雷达方位回波建模；基于机载扫描雷达的运动几何模型，扫描雷达发射线性调频信号，经过下变频处理，接收到的回波信号为： 其中，σx,y表示场景Ω中点目标的散射系数，x,y表示场景中某一点的坐标，wt表示天线方向图函数调制，rect·表示矩形窗函数，τ表示距离向时间采样向量，t表示方位向时间采样向量，Tp表示发射信号的脉冲时宽，λ表示载频波长，c表示电磁波传播速度，k表示线性调频率，nτ,t表示加性高斯白噪声，Rt表示目标的距离历史，且R0表示目标起始距离，v表示机载平台运动速度，θ0表示目标空间方位角；经过脉冲压缩和距离走动校正后，方位向回波视为目标散射系数与天线方向图构成的测量矩阵的类卷积，具体如下：y＝Hx+n2其中，y表示回接收到的方位向回波向量，H表示由天线方向图构成的卷积测量矩阵，x表示目标散射系数，n表示噪声向量；步骤二、构建正则化目标函数；在正则化框架下，选择全变差范数构建目标函数，其表达式为： 其中，表示恢复的目标散射系数分布，表示数据表征项，表示向量二范数的平方：η表示正则化参数；||Dx||1表示全变差范数，D表示梯度矩阵；选择交替方向乘子法ADMM对目标函数进行求解，引入新变量z将式3中的目标函数重写为如下形式： 对于式4，它的增广拉格朗日形式可写为： 其中，·T表示转置操作；p表示增广拉格朗日参数；x表示目标散射系数，u表示对偶变量；将式5表示的优化问题分解为三个子问题，通过交替更新变量xi+1、zi+1和ui+1来求解，具体更新过程如下： 其中，表示门限，表示收缩算子，i表示变量的迭代次数；步骤三、矩阵求逆转化为递归快速迭代运算；对方位向回波y和天线测量矩阵H进行分块处理，则yn＝[y0,y1,…,yn]T∈Cn×1，Hn＝[h1,h2,…,hn]T∈Cn×N；其中，yn表示n时刻方位向接收到的总的回波信号，yn表示n时刻总回波信号yn的第n个元素，即方位上接收到的第n个回波值；C表示复数域，Cn×1表示复数域上维度为n×1大小的矩阵，Hn表示将天线测量矩阵H进行行分块处理得到的矩阵，hn∈C1×N表示H的第n行行向量，N表示天线测量矩阵H的列数；则式6中关于x的更新公式可转换为如下形式： 令可以得到pn+1和pn的关系如下： 利用矩阵反演恒等式可得： 将上式带入到式7可得： 观察式10可以发现，矩阵反演被转换成迭代处理；步骤四、波束回波递归更新重建结果；基于上述推导，在每次回波到来导致回波矩阵更新时，利用下述公式进行递归更新： 其中，选取作为初始值，关于变量z和u的初始值如下： 经过I次迭代后，直至相邻两次的超分辨结果的误差不大于设定的误差值χ，循环结束；循环终止条件为： 每次回波数据更新，重复步骤四，输出实时超分辨处理结果，直至所有回波更新完毕。

全文数据：

权利要求：

百度查询： 电子科技大学 一种基于全变差约束的扫描雷达递归快速超分辨成像方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD