基于期望传播算法的逆合成孔径雷达成像方法技术

本发明专利技术公开了一种基于期望传播算法的逆合成孔径雷达成像方法，本发明专利技术的实现步骤是：(1)稀疏贝叶斯建模；(2)利用期望传播算法求解稀疏贝叶斯模型中每个距离单元的实数权向量；(3)重构逆合成孔径雷达ISAR图像散射系数矩阵；(4)转置并画出矩阵，得到回波缺损及低信噪比情况下的目标逆合成孔径雷达ISAR成像结果。本发明专利技术基于稀疏信号重构理论实现了逆合成孔径雷达ISAR成像，可用于在回波缺损及低信噪比情况下对空间与空中非合作目标的二维成像。

【技术实现步骤摘要】
基于期望传播算法的逆合成孔径雷达成像方法
本专利技术属于通信
，更进一步涉及雷达信号处理
中的一种基于期望传播算法的逆合成孔径雷达ISAR(InverseSyntheticApertureRadar)成像方法，可用于逆合成孔径雷达接收到目标回波信噪比低、存在缺损等情况下的空间目标二维成像。
技术介绍
由于具有全天时、全天候、高分辨率和远距离等特点，逆合成孔径雷达ISAR在航空与航天目标观测中发挥着重要作用。当受雷达工作模式限制而无法对目标进行连续观测时，会产生方位缺损回波，此时传统基于脉冲压缩的成像方法会产生高旁瓣。此外，当逆合成孔径雷达ISAR对远距离小目标进行探测时，回波信噪比较低，从而很难得到高质量成像结果。对于逆合成孔径雷达，目标强散射点数目远小于观测样本数，因此具有稀疏性。目前，充分利用其稀疏性，在回波缺损、低信噪比情况下实现空间目标高质量、高分辨成像是提高现有雷达对空间、空中目标探测与监视能力的关键技术。西安电子科技大学在其申请的专利技术专利文献“基于稀疏孔径的机动目标逆合成孔径雷达成像方法”(公开号：103901429A，申请号：201410140123.5)中公开了一种基于稀疏重构的成像方法。该方法具体步骤为：对回波数据进行距离压缩和运动补偿，利用正交匹配追踪算法得到稀疏孔径重构回波信号，接着对重构回波信号进行快速傅里叶变换，最终实现逆合成孔径雷达ISAR的距离-多普勒成像。该方法虽然能够在回波缺损情况下实现高分辨成像，但是，该方法仍然存在的不足之处是，正交匹配追踪算法稀疏表征权向量能力不足，在回波缺损及低信噪比情况下容易产生虚假点，无法获得聚焦良好的ISAR图像。吴称光，邓彬，苏伍各，王宏强，秦玉亮在其发表的论文“基于块稀疏贝叶斯模型的ISAR成像方法”(电子与信息学报2015，37(12)：2941-2947)中提出一种基于块稀疏贝叶斯模型的ISAR成像方法。该方法基于稀疏信号表示理论，将ISAR高分辨成像问题转化为稀疏信号表示问题，利用变分推断算法求解模型参数，最终实现目标逆合成孔径雷达ISAR成像。该方法虽然能够在回波缺损情况下实现二维ISAR成像，但该方法存在的不足之处是，在低信噪比情况下该方法参数估计误差较大，无法获得准确良好的ISAR像。
技术实现思路
本专利技术的目的是针对上述稀疏成像技术中数值优化方法稀疏表征能力不足，低信噪比条件下误差较高等局限性，提出一种基于期望传播算法的逆合成孔径雷达ISAR成像方法。该方法充分利用回波信号的稀疏性，通过回波信号的稀疏建模及稀疏系数求解最终实现目标回波缺损与低信噪比情况下的空间目标逆合成孔径雷达ISAR成像。为实现本专利技术的目的的具体思路是：基于稀疏信号重构理论，将ISAR成像问题转化为稀疏线性回归问题，利用伯努利高斯模型对目标回波建模，进而采用期望传播算法求解权值向量，最终实现回波缺损及低信噪比情况下的空间目标二维成像。为实现上述目的，本专利技术的主要步骤如下：(1)稀疏贝叶斯建模：(1a)计算距离向脉冲压缩后每个距离单元有效复数回波向量对应的实数向量；(1b)将待测目标周围环境中复数噪声向量对应的实数噪声向量的先验分布设置为高斯分布；(1c)将每个距离单元的实数权向量的稀疏先验分布设置为伯努利高斯分布；(2)利用期望传播算法计算每个距离单元的实数权向量：(2a)设置初始距离单元的序号q＝1；(2b)设置初始迭代次数k＝1，在第q个距离单元中，将先验分布R中所有高斯分布的方差值均初始化为p0vs，其中p0表示方位存在散射点的概率，初始化为5×10-3，vs表示散射点幅度变化的方差值，初始化为0.1，将距离单元实数权向量初始化为2Na×1维零向量，将噪声精度值初始化为4×10-4；(2c)利用迭代公式，依次计算当前迭代次数中的五个参数值，用当前迭代次数中的五个参数值分别替换上一次迭代次数中的五个参数值；(2d)将当前迭代与上次迭代所求得的距离单元实数权向量相减得到差向量，判断差向量中所有元素最大值是否小于10-5，若是，则停止迭代，得到距离单元对应的实数权向量，执行步骤(2e)，否则，在当前迭代次数上加1后执行步骤(2c)；(2e)更新距离单元序号q＝q+1，判断更新后的距离单元序号是否大于距离单元总数，若是，则停止对距离单元的搜索，执行步骤(3)；否则，执行步骤(2b)；(3)重构逆合成孔径雷达ISAR图像散射系数矩阵：(3a)产生一个维数为2Na×Nr的矩阵W，其中Na表示方位向采样总数，Nr表示距离向采样总数，矩阵W的第q列为第q个距离单元的实数权向量；(3b)利用公式W′＝W1+jW2，得到逆合成孔径雷达ISAR图像的散射系数矩阵W′，其中，W1表示矩阵W第1行到第Na行，第1列到第Nr列所有元素构成的矩阵，W2表示矩阵W第Na+1行到第2Na行，第1列到第Nr列所有元素构成的矩阵，j表示虚数单位符号；(4)转置并二维成像：对逆合成孔径雷达ISAR图像散射系数矩阵W′进行转置操作，将得到的Nr×Na维矩阵画出来，得到目标的二维逆合成孔径雷达ISAR图像。本专利技术与现有的技术相比，具有以下优点：第一，本专利技术在稀疏贝叶斯建模时，将待测目标周围环境中复数噪声向量对应的实数噪声向量的先验分布设置为高斯分布，以及将每个距离单元的实数权向量的稀疏先验分布设置为伯努利高斯分布，充分利用了噪声的先验信息和目标散射点分布的稀疏性，克服了现有技术在低信噪比条件下参数估计方法误差较大的问题，使得本专利技术在低信噪比条件下有效提高了成像质量。第二，本专利技术利用期望传播算法，计算每个距离单元的实数权向量，克服了现有技术采用数值优化方法与贝叶斯推断方法稀疏表征权向量能力不足、容易产生虚假点等问题，使得本专利技术在回波缺损等条件下，能够获得空间目标聚焦良好的逆合成孔径雷达ISAR图像。附图说明图1为本专利技术的流程图；图2为本专利技术的仿真图。具体实施方式下面结合附图对本专利技术做进一步的描述。参照附图1，本专利技术的具体实施步骤如下：步骤1，稀疏贝叶斯建模。A.按照下式，计算距离向脉冲压缩后每个距离单元有效复数回波向量对应的实数向量：yf＝Φwq+n其中，yf表示距离向脉冲压缩后第f个距离单元有效复数回波向量对应的实数向量，Φ表示与雷达方位向有效复数傅里叶字典对应的实数傅里叶字典矩阵，wq表示第q个距离单元与待重构逆合成孔径雷达ISAR图像复数散射系数向量对应的实数权向量，n表示与待测目标周围环境中复数噪声向量对应的实数噪声向量，f与q取值相等。所述距离向脉冲压缩后每个距离单元有效复数回波向量对应的实数向量由以下步骤得到：(1)逆合成孔径雷达发射的电磁波在传播过程中遇到目标后，所反射的回波在传播过程中受到遮挡、稀疏观测等影响，最终被雷达接收机接收并记录为目标缺损回波矩阵。(2)令逆合成孔径雷达到场景中心的距离为参考距离，令与逆合成孔径雷达发射信号载频、调频率相同，距离为参考距离的线性调频信号为参考信号Sref，将参考信号Sref取共轭后与有效回波矩阵Se相乘，得到解线频调后信号Srd：其中，为距离快时间，tm为方位慢时间，Sref(·)为参考信号，Srd为解线频调后的信号，*表示共轭操作。(3)对解线频调后的信号Srd沿快时间维作傅里叶变换，得到缺损回波距离脉冲压缩后的矩阵本文档来自技高网...

【技术保护点】
1.一种基于期望传播算法的逆合成孔径雷达成像方法，其特征在于，包括如下步骤：(1)稀疏贝叶斯建模：(1a)计算距离向脉冲压缩后每个距离单元有效复数回波向量对应的实数向量；(1b)将待测目标周围环境中复数噪声向量对应的实数噪声向量的先验分布设置为高斯分布；(1c)将每个距离单元的实数权向量的先验分布设置为伯努利高斯分布；(2)利用期望传播算法，计算每个距离单元的实数权向量：(2a)设置初始距离单元的序号q＝1；(2b)设置初始迭代次数k＝1，在第q个距离单元中，将先验分布R中所有高斯分布的方差值均初始化为p0vs，其中p0表示方位单元存在散射点的概率，初始化为5×10‑3，vs表示散射点幅度变化的方差值，初始化为0.1，将距离单元实数权向量初始化为2Na×1维零向量，将噪声精度值初始化为4×10‑4；(2c)利用迭代公式，依次计算当前迭代次数中的五个参数值，用当前迭代次数中的五个参数值分别替换上一次迭代次数中的五个参数值；(2d)将当前迭代与上次迭代所求得的距离单元实数权向量相减得到差向量，判断差向量中所有元素最大值是否小于10‑5，若是，则停止迭代，得到距离单元对应的实数权向量，执行步骤(2e)，否则，在当前迭代次数上加1后执行步骤(2c)；(2e)更新距离单元序号q＝q+1，判断更新后的距离单元序号是否大于距离单元总数，若是，则停止对距离单元的搜索，执行步骤(3)；否则，执行步骤(2b)；(3)重构逆合成孔径雷达ISAR图像散射系数矩阵：(3a)产生一个维数为2Na×Nr的矩阵W，其中Na表示方位向采样总数，Nr表示距离向采样总数，矩阵W的第q列为第q个距离单元的实数权向量；(3b)利用公式W′＝W1+jW2，得到逆合成孔径雷达ISAR图像的散射系数矩阵W′，其中，W1表示矩阵W第1行到第Na行，第1列到第Nr列所有元素构成的矩阵，W2表示矩阵W第Na+1行到第2Na行，第1列到第Nr列所有元素构成的矩阵，j表示虚数单位符号；(4)转置并二维成像：对逆合成孔径雷达ISAR图像散射系数矩阵W′进行转置操作，将得到的Nr×Na维矩阵画出来，得到目标的二维逆合成孔径雷达ISAR图像。...

【技术特征摘要】
1.一种基于期望传播算法的逆合成孔径雷达成像方法，其特征在于，包括如下步骤：(1)稀疏贝叶斯建模：(1a)计算距离向脉冲压缩后每个距离单元有效复数回波向量对应的实数向量；(1b)将待测目标周围环境中复数噪声向量对应的实数噪声向量的先验分布设置为高斯分布；(1c)将每个距离单元的实数权向量的先验分布设置为伯努利高斯分布；(2)利用期望传播算法，计算每个距离单元的实数权向量：(2a)设置初始距离单元的序号q＝1；(2b)设置初始迭代次数k＝1，在第q个距离单元中，将先验分布R中所有高斯分布的方差值均初始化为p0vs，其中p0表示方位单元存在散射点的概率，初始化为5×10-3，vs表示散射点幅度变化的方差值，初始化为0.1，将距离单元实数权向量初始化为2Na×1维零向量，将噪声精度值初始化为4×10-4；(2c)利用迭代公式，依次计算当前迭代次数中的五个参数值，用当前迭代次数中的五个参数值分别替换上一次迭代次数中的五个参数值；(2d)将当前迭代与上次迭代所求得的距离单元实数权向量相减得到差向量，判断差向量中所有元素最大值是否小于10-5，若是，则停止迭代，得到距离单元对应的实数权向量，执行步骤(2e)，否则，在当前迭代次数上加1后执行步骤(2c)；(2e)更新距离单元序号q＝q+1，判断更新后的距离单元序号是否大于距离单元总数，若是，则停止对距离单元的搜索，执行步骤(3)；否则，执行步骤(2b)；(3)重构逆合成孔径雷达ISAR图像散射系数矩阵：(3a)产生一个维数为2Na×Nr的矩阵W，其中Na表示方位向采样总数，Nr表示距离向采样总数，矩阵W的第q列为第q个距离单元的实数权向量；(3b)利用公式W′＝W1+jW2，得到逆合成孔径雷达ISAR图像的散射系数矩阵W′，其中，W1表示矩阵W第1行到第Na行，第1列到第Nr列所有元素构成的矩阵，W2表示矩阵W第Na+1行到第2Na行，第1列到第Nr列所有元素构成的矩阵，j表示虚数单位符号；(4)转置并二维成像：对逆合成孔径雷达ISAR图像散射系数矩阵W′进行转置操作，将得到的Nr×Na维矩阵画出来，得到目标的二维逆合成孔径雷达ISAR图像。2.根据权利要求1所述的基于期望传播算法的逆合成孔径雷达成像方法，其特征在于，步骤(1a)中所述的距离向脉冲压缩后每个距离单元有效复数回波向量对应的实数向量是按照下述公式计算得到的：yf＝Φwq+n其中，yf表示距离向脉冲压缩后第f个距离单元有效复数回波向量对应的实数向量，Φ表示与雷达方位向有效复数傅里叶字典对应的实数傅里叶字典矩阵，wq表示第q个距离单元与待重构逆合成孔径雷达ISAR图像复数散射系数向量对应的实数权向量，n表示与待测目标周围环境中复数噪声向量对应的实数噪声向量，f与q取值相等。3.根据权利要求2所述的基于期望传播算法的逆合成孔径雷达成像方法，其特征在于，所述的与雷达方位向有效复数傅里叶字典对应的实数傅里叶字典矩阵是由以下步骤得到的：第一步，按照下式，构造一个雷达方位向傅里叶字典：Φ′＝exp(...

【专利技术属性】
技术研发人员：白雪茹，王格，周峰，李小勇，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61