一种基于复块贝叶斯算法的雷达成像方法技术

本发明专利技术涉及探地雷达信号处理技术领域，具体涉及一种基于复块贝叶斯算法的雷达成像方法，包括建立探地雷达反演的真实模型；基于真实模型进行正演运算，获得观测雷达数据，生成观测数据矩阵；将成像区域划分为若干网格，并计算不同阵元到每个网格的时延，将时延作为字典矩阵；基于观测数据矩阵和字典矩阵利用复块稀疏贝叶斯算法计算系数向量；对系数向量进行高分辨成像，通过观测数据矩阵和字典矩阵利用复块稀疏贝叶斯算法计算系数向量完成对观测数据矩阵的目标反射系数的重构，降低了计算的复杂度，提高了成像的精度，解决了现有的雷达成像方法通过利用回波信息直接计算的计算量大，从而减缓了计算速度的问题。从而减缓了计算速度的问题。从而减缓了计算速度的问题。

【技术实现步骤摘要】
一种基于复块贝叶斯算法的雷达成像方法


[0001]本专利技术涉及探地雷达信号处理
，尤其涉及一种基于复块贝叶斯算法的雷达成像方法。

技术介绍

[0002]探地雷达可以无损探测地下浅层目标，地下目标成像具有广泛应用前景。
[0003]目前传统地下目标成像算法是通过回波时延计算出各道回波的幅值，进而对幅值相干叠加完成成像。
[0004]但是这些雷达成像技术利用回波信息直接计算，虽然技术实现简单，但是在这种情况下，需要大量计算回波的幅值，导致计算量大、计算速度缓慢。

技术实现思路

[0005]本专利技术的目的在于提供一种基于复块贝叶斯算法的雷达成像方法，旨在解决现有的雷达成像方法通过利用回波信息直接计算的计算量大，从而减缓了计算速度的问题。
[0006]为实现上述目的，本专利技术提供了一种基于复块贝叶斯算法的雷达成像方法，包括以下步骤：
[0007]S1建立探地雷达反演的真实模型；
[0008]S2基于所述真实模型进行正演运算，获得观测雷达数据，生成观测数据矩阵；
[0009]S3将成像区域划分为若干网格，并计算不同阵元到每个网格的时延，将时延作为字典矩阵；
[0010]S4基于所述观测数据矩阵和所述字典矩阵利用复块稀疏贝叶斯算法计算系数向量；
[0011]S5对所述系数向量进行高分辨成像。
[0012]其中，所述基于所述真实模型进行正演运算，获得观测雷达数据，生成观测数据矩阵的具体方式为：
[0013]S21获取所述真实模型的若干个测量位置；
[0014]S22设定发射天线的发射频率，得到设定发射天线；
[0015]S23使用所述设定发射天线分别在若干个所述测量位置发射电磁波，并使用接收天线分别在若干个所述测量位置对经过目标反射后所述电磁波进行接收，得到若干个目标回波信号；
[0016]S24将若干个所述目标回波信号进行整合，得到观测雷达数据；
[0017]S25基于所述观测雷达数据生成观测数据矩阵。
[0018]其中，所述将成像区域划分为若干网格，并计算不同阵元到每个网格的时延，将时延作为字典矩阵的具体方式为：
[0019]S31将成像区域划分为若干网格；
[0020]S32使用所述设定发射天线在每个所述网格发射电磁波，使用所述接收天线对经
每个所述网格反射后的所述电磁波进行接收，得到若干网格回波信号；
[0021]S33基于若干所述网格回波信号计算不同阵元到每个网格的时延；
[0022]S34将所述时延作为字典矩阵。
[0023]其中，所述基于所述观测数据矩阵和所述字典矩阵利用复块稀疏贝叶斯算法计算系数向量的具体方式为：
[0024]S41基于所述观测数据矩阵和所述字典矩阵构建回波信号模型；
[0025]S42输入所述回波信号模型和所述字典矩阵，得到输入公式；
[0026]S43对所述回波信号模型的信号的方差及噪声的方差进行初始化，得到初始化值；
[0027]S44利用块稀疏贝叶斯基权向量公式迭代更新所述初始化值，每迭代一次，信号的均值、相关性、协方差矩阵及噪声方差更改一次。
[0028]S45判断所述当前基权向量中的信号均值与上次迭代得到的均值次数差值的最大值是否小于设置的终止阈值，若是，执行步骤S47，否则，执行步骤S46；
[0029]S46判断所述当前迭代次数是否大于设置的迭代次数，若是，执行步骤S47，否则，将所述当前迭代次数加1后执行步骤S44；
[0030]S47所述信号均值即为需要求得的系数向量。
[0031]本专利技术的一种基于复块贝叶斯算法的雷达成像方法，通过建立探地雷达反演的真实模型；基于所述真实模型进行正演运算，获得观测雷达数据，生成观测数据矩阵；将成像区域划分为若干网格，并计算不同阵元到每个网格的时延，将时延作为字典矩阵；基于所述观测数据矩阵和所述字典矩阵利用复块稀疏贝叶斯算法计算系数向量；对所述系数向量进行高分辨成像，通过所述观测数据矩阵和所述字典矩阵利用复块稀疏贝叶斯算法计算系数向量完成对观测数据矩阵的目标反射系数的重构，降低了计算的复杂度，提高了成像的精度，解决了现有的雷达成像方法通过利用回波信息直接计算的计算量大，从而减缓了计算速度的问题。，解决了现有的雷达成像方法通过利用回波信息直接计算的计算量大，从而减缓了计算速度的问题。
附图说明
[0032]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0033]图1是本专利技术提供的一种基于复块贝叶斯算法的雷达成像方法的流程图。
[0034]图2是基于所述真实模型进行正演运算，获得观测雷达数据，生成观测数据矩阵的流程图。
[0035]图3是将成像区域划分为若干网格，并计算不同阵元到每个网格的时延，将时延作为字典矩阵的流程图。
[0036]图4是基于所述观测数据矩阵和所述字典矩阵利用复块稀疏贝叶斯算法计算系数向量的流程图。
[0037]图5是单目标模型场景示意图。
[0038]图6是复块稀疏贝叶斯恢复单目标模型仿真结果图.
[0039]图7是多目标模型场景示意图。
[0040]图8是复块稀疏贝叶斯恢复多目标模型仿真结果图。
[0041]图9是不同信噪比下的均方误差图。
[0042]图10为探地雷达成像模型。
具体实施方式
[0043]下面详细描述本专利技术的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本专利技术，而不能理解为对本专利技术的限制。
[0044]请参阅图1至图5，本专利技术提供一种基于复块贝叶斯算法的雷达成像方法，包括以下步骤：
[0045]S1建立探地雷达反演的真实模型；
[0046]具体的，将所述真实模型设置为两层介质，一层为空气，设置介电常数ε1，另一层为土壤，设置介电常数ε2。
[0047]S2基于所述真实模型进行正演运算，获得观测雷达数据，生成观测数据矩阵；
[0048]具体方式为：
[0049]S21获取所述真实模型的若干个测量位置；
[0050]具体的，假设在整个测量过程中共有M个测量位置。
[0051]S22设定发射天线的发射频率，得到设定发射天线；
[0052]具体的，设置发射天线发射频率从f0到f
H
的电磁波。
[0053]S23使用所述设定发射天线分别在若干个所述测量位置发射电磁波，并使用接收天线分别在若干个所述测量位置对经过目标反射后所述电磁波进行接收，得到若干个目标回波信号；
[本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于复块贝叶斯算法的雷达成像方法，其特征在于，包括以下步骤：S1建立探地雷达反演的真实模型；S2基于所述真实模型进行正演运算，获得观测雷达数据，生成观测数据矩阵；S3将成像区域划分为若干网格，并计算不同阵元到每个网格的时延，将时延作为字典矩阵；S4基于所述观测数据矩阵和所述字典矩阵利用复块稀疏贝叶斯算法计算系数向量；S5对所述系数向量进行高分辨成像。2.如权利要求1所述的基于复块贝叶斯算法的雷达成像方法，其特征在于，所述基于所述真实模型进行正演运算，获得观测雷达数据，生成观测数据矩阵的具体方式为：S21获取所述真实模型的若干个测量位置；S22设定发射天线的发射频率，得到设定发射天线；S23使用所述设定发射天线分别在若干个所述测量位置发射电磁波，并使用接收天线分别在若干个所述测量位置对经过目标反射后所述电磁波进行接收，得到若干个目标回波信号；S24将若干个所述目标回波信号进行整合，得到观测雷达数据；S25基于所述观测雷达数据生成观测数据矩阵。3.如权利要求2所述的基于复块贝叶斯算法的雷达成像方法，其特征在于，所述将成像区域划分为若干网格，并计算不同阵元到每个网格的时延，将时延作为字典矩阵的具体方式为：S3...

【专利技术属性】
技术研发人员：刘庆华，杜文静，廖可非，谢宁波，
申请(专利权)人：桂林电子科技大学，
类型：发明
国别省市：