基于最大似然估计的前视超分辨成像方法技术

本发明专利技术公开了一种基于最大似然估计的前视超分辨成像方法。其包括获取正交通道接收的雷达回波及距离向脉冲压缩，获取天线方向图信息，计算方位向任意点幅度的边缘概率密度函数，构造目标函数并计算目标函数梯度，初始化系统参数并利用目标函数梯度进行方位迭代处理对目标函数逼近最优解，判断距离压缩后的回波数据是否处理完毕。本发明专利技术首先根据数据采集过程将回波分解为I/Q双通道进行噪声建模，计算回波联合概率密度函数，从而构建准确的目标函数，再通过最大化目标函数获得相对精确的超分辨率结果；本发明专利技术通过采用更准确的回波统计模型，可以在低信噪比情况下更准确地恢复目标。

【技术实现步骤摘要】
基于最大似然估计的前视超分辨成像方法
本专利技术属于雷达成像
，尤其涉及一种基于最大似然估计的前视超分辨成像方法。
技术介绍
雷达成像技术凭借其全天时、全天候和高成像质量的优点，已成为当今探测领域不可取代的技术手段，在民用、军事领域都发挥着越来越重要的作用。由于常规的单站SAR和DBS技术无法实现前视成像，实波束扫描雷达超分辨技术成为实现前视成像的一种重要技术手段。文献“Weightedleastsquaresmethodforforward-lookingimagingofscanningradar”(ZhangY,HuangY,ZhaY,GeoscienceandRemoteSensingSymposium.IEEE,2014:714-716.)采用了一种非线性限制迭代逆滤波技术对实波束进行锐化，但此方法没有利用噪声和目标的先验信息，在低信噪比情况下性能较差；文献“Superresolutionofcoherentsourcesinreal-beamdata”(UttamS,GoodmanN.AerospaceandElectronicSystems,IEEETransactionson,2010,46(3):1557-1566.)利用幅度相位信息进行锐化，但此方法由于需要大量快拍数对自相关矩阵进行准确估计，故在单次或快拍数较少的情况下性能较差；文献“AnimprovedRichardson-Lucyalgorithmforradarangularsuper-resolution”(ZhaY,HuangY,YangJ,etal.RadarConference,2014IEEE,2014：0406-0410)在最大似然框架下，分别将噪声和目标建模为高斯分布和泊松分布，进而利用最大似然方法对目标散射系数进行估计，但由于采用的噪声模型并未考虑雷达正交双通道回波的情况，导致虚假目标现象仍然比较突出。
技术实现思路
本专利技术的专利技术目的是：为了解决现有技术中存在的缺陷，本专利技术提出了一种基于最大似然估计的前视超分辨成像方法，以期解决传统最大似然方法成像结果不准，特别是低信噪比情况下虚假目标较多的问题，从而实现较准确的超分辨成像。本专利技术的技术方案是：一种基于最大似然估计的前视超分辨成像方法，包括以下步骤：A、获取正交通道接收的雷达回波数据，对距离向进行脉冲压缩，得到距离压缩后的回波数据；B、获取天线方向图，构造天线方向图卷积矩阵；C、根据步骤B中天线方向图卷积矩阵计算方位向任意点幅度的边缘概率密度函数；D、根据步骤C中方位向任意点幅度的边缘概率密度函数计算回波联合概率密度函数，构造目标函数并计算目标函数梯度；E、初始化系统参数，利用步骤D中的目标函数梯度进行方位迭代处理对目标函数逼近最优解；F、判断步骤A中距离压缩后的回波数据是否处理完毕；若是，则输出成像结果；若否，则返回步骤E。进一步地，所述步骤B获取天线方向图，构造天线方向图卷积矩阵具体为：获取天线方向图信息，对天线方向图进行离散化，得到离散化的天线方向图，表示为h＝[h1,...,hm]T，其中m为天线方向图长度，h1,...,hm表示天线方向图的各个采样点；构造天线方向图卷积矩阵，表示为进一步地，所述步骤C根据步骤A中距离压缩后的回波数据及步骤B中天线方向图卷积矩阵计算方位向任意点幅度的边缘概率密度函数具体包括以下分步骤：C1、构造I/Q通道解调回波模型，表示为其中为方位向回波，x为待估计的目标，为多普勒相位向量，nc和ns分别为I、Q通道高斯噪声向量，⊙为点乘，j为虚数单位；C2、利用步骤C1中I/Q通道解调回波模型对回波采样点进行处理，得到I/Q通道解调模型，表示为其中si和φi分别为方位向第i采样点的随机幅度和随机相位，为多普勒相位，和分别为I、Q通道高斯噪声；C3、对I、Q通道高斯噪声根据二维随机变量的概率密度函数进行变换，得到联合概率密度函数，表示为其中σ为高斯噪声分布标准差，为和联合概率密度函数；C4、根据联合概率密度函数对方位向任意采样点的随机相位进行积分，得到方位向任意点幅度的边缘概率密度函数，表示为其中J0(·)为零阶贝叶斯函数。进一步地，所述步骤D中回波联合概率密度函数表示为其中，L为方位采样点数。进一步地，所述步骤D中目标函数表示为进一步地，所述步骤D中目标函数梯度表示为其中J1(·)是一阶贝塞尔函数。进一步地，所述步骤E初始化系统参数，利用步骤D中的目标函数梯度进行方位迭代处理对目标函数逼近最优解具体包括以下分步骤：E1、设置迭代步长α，对回波矩阵s第q行数据sq取幅值，并将其设置为迭代初始值x0＝sq，设置加速迭代过程的预测初始值为y0＝sq，设置控制预测点迭代速度的参数初值为t1＝1；E2、根据迭代公式对步骤D中目标函数梯度沿梯度方向进行方位迭代，所述迭代公式表示为其中，xk表示第k步迭代结果；E3、根据迭代初始值x0和迭代步长α计算第一步迭代结果x1；E4、根据控制预测点迭代速度的参数初值t1计算参数值t2，计算公式表示为E5、计算加速迭代过程的预测值y1，计算公式表示为E6、根据加速迭代过程的预测值y1和迭代步长α计算第二步迭代结果x2；E7、判断第一步迭代结果x1和第二步迭代结果x2是否满足迭代终止条件||x2-x1||2＜δ，其中δ为设定的阈值；若是，则进行下一步骤；若否，则将x2、x1、y1和t2分别赋值给x1、x0、y0和t1，返回步骤E2。本专利技术的有益效果是：本专利技术首先根据数据采集过程将回波分解为I/Q双通道进行噪声建模，计算回波联合概率密度函数，从而构建准确的目标函数，再通过最大化目标函数获得相对精确的超分辨率结果；本专利技术通过采用更准确的回波统计模型，可以在低信噪比情况下更准确地恢复目标。附图说明图1是本专利技术的基于最大似然估计的前视超分辨成像方法的流程示意图。图2是本专利技术实施例中前视扫描雷达系统结构示意图。图3是本专利技术实施例中目标场景设置一维剖面示意图。图4是本专利技术实施例中归一化天线方向示意图。图5是本专利技术实施例中场景实波束回波结果剖面示意图。图6是本专利技术实施例中对目标成像结果的剖面示意图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本专利技术，并不用于限定本专利技术。如图1所示，为本专利技术的基于最大似然估计的前视超分辨成像方法的流程示意图。一种基于最大似然估计的前视超分辨成像方法，包括以下步骤：A、获取正交通道接收的雷达回波数据，对距离向进行脉冲压缩，得到距离压缩后的回波数据；B、获取天线方向图，构造天线方向图卷积矩阵；C、根据步骤B中天线方向图卷积矩阵计算方位向任意点幅度的边缘概率密度函数；D、根据步骤C中方位向任意点幅度的边缘概率密度函数计算回波联合概率密度函数，构造目标函数并计算目标函数梯度；E、初始化系统参数，利用步骤D中的目标函数梯度进行方位迭代处理对目标函数逼近最优解；F、判断步骤A中距离压缩后的回波数据是否处理完毕；若是，则输出成像结果；若否，则返回步骤E。在步骤A中，本专利技术针对获取的二维回波数据S(τ,t)，其中τ表示快时间，t表示慢时间，利用雷达发射信号参数，构造匹配滤波器，对距离向进本文档来自技高网...

【技术保护点】
一种基于最大似然估计的前视超分辨成像方法，其特征在于，包括以下步骤：A、获取正交通道接收的雷达回波数据，对距离向进行脉冲压缩，得到距离压缩后的回波数据；B、获取天线方向图，构造天线方向图卷积矩阵；C、根据步骤B中天线方向图卷积矩阵计算方位向任意点幅度的边缘概率密度函数；D、根据步骤C中方位向任意点幅度的边缘概率密度函数计算回波联合概率密度函数，构造目标函数并计算目标函数梯度；E、初始化系统参数，利用步骤D中的目标函数梯度进行方位迭代处理对目标函数逼近最优解；F、判断步骤A中距离压缩后的回波数据是否处理完毕；若是，则输出成像结果；若否，则返回步骤E。

【技术特征摘要】
1.一种基于最大似然估计的前视超分辨成像方法，其特征在于，包括以下步骤：A、获取正交通道接收的雷达回波数据，对距离向进行脉冲压缩，得到距离压缩后的回波数据；B、获取天线方向图，构造天线方向图卷积矩阵；C、根据步骤B中天线方向图卷积矩阵计算方位向任意点幅度的边缘概率密度函数；D、根据步骤C中方位向任意点幅度的边缘概率密度函数计算回波联合概率密度函数，构造目标函数并计算目标函数梯度；E、初始化系统参数，利用步骤D中的目标函数梯度进行方位迭代处理对目标函数逼近最优解；F、判断步骤A中距离压缩后的回波数据是否处理完毕；若是，则输出成像结果；若否，则返回步骤E。2.如权利要求1所述的基于最大似然估计的前视超分辨成像方法，其特征在于，所述步骤B获取天线方向图，构造天线方向图卷积矩阵具体为：获取天线方向图信息，对天线方向图进行离散化，得到离散化的天线方向图，表示为h＝[h1,...,hm]T，其中m为天线方向图长度，h1,...,hm表示天线方向图的各个采样点；构造天线方向图卷积矩阵，表示为3.如权利要求2所述的基于最大似然估计的前视超分辨成像方法，其特征在于，所述步骤C根据步骤B中天线方向图卷积矩阵计算方位向任意点幅度的边缘概率密度函数具体包括以下分步骤：C1、构造I/Q通道解调回波模型，表示为其中为方位向回波数据，x为待估计的目标，为多普勒相位向量，nc和ns分别为I、Q通道高斯噪声向量，为点乘，j为虚数单位；C2、利用步骤C1中I/Q通道解调回波模型对回波采样点进行处理，得到I/Q通道解调模型，表示为其中si和φi分别为方位向第i采样点的随机幅度和随机相位，为多普勒相位，和分别为I、Q通道高斯噪声；C3、对I、Q通道高斯噪声根据二维随机变量的概率密度函数进行变换，得到联合概率密度函数，表示为其中σ为高斯噪声分布标准差；C4、根据联合概率密度函数对方位向任意采样点的随机相位进行积分，得到方位向任意点幅度的边缘概率密度函数，表示为其中J0(·)为零阶贝叶斯函数。4.如权利要求3所述的基于最大似然估...

【专利技术属性】
技术研发人员：李文超，贺凤祥，谭珂，张倩，张永超，毛德庆，张启平，徐帆云，黄钰林，杨建宇，武俊杰，张寅，杨海光，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51