一种分布式阵列雷达迭代自适应高分辨成像方法技术

技术编号:22329242 阅读:19 留言:0更新日期:2019-10-19 12:06
本发明专利技术提供了一种分布式阵列雷达迭代自适应高分辨成像方法,属于雷法成像技术领域。本发明专利技术针对传统分布式阵列雷达的构型依赖问题,首先利用分布式阵列雷达稀疏布阵的空间关系,导出分布式雷达阵列空间谱分布;其次,利用分布式雷达空间谱分布和系统参数,构造分布式阵列雷达系统观测矩阵;最后,为降低成像处理过程的构型依赖性,采用迭代自适应成像处理方法,实现分布式阵列雷达高分辨成像。本发明专利技术通过目标求解过程中的迭代自适应收敛,降低了分布式阵列雷达对构型的依赖性。

【技术实现步骤摘要】
一种分布式阵列雷达迭代自适应高分辨成像方法
本专利技术属于雷达成像
,特别涉及一种分布式阵列雷达迭代自适应高分辨成像方法。
技术介绍
雷达利用反射电磁波中携带的目标信息,获得雷达图像,在灾情监测、目标追踪等场景中得到广泛应用。但单个雷达平台由于其天线孔径的限制,难以获得高分辨图像。分布式阵列雷达利用多个平台在空间上分集关系形成大的观测孔径,为雷达高分辨成像提供条件。而分布式阵列雷达的空间分布关系,对目标观测信息的获取产生着至关重要的影响。现有技术中,通过建立一种分布式多通道雷达成像观测模型,采用非线性傅里叶变化方法获得目标成像结果,但是该方法对雷达构型分布的要求较为严格。另外。可基于MUSIC方法的分布式雷达成像方法,但是场景中目标散射点的数量需要通过盖尔圆盘准则进行估计,该方法对目标数量和雷达构型分布均较为敏感。综上可知,上述方法能够用于分布式阵列雷达成像应用中,但是上述方法对阵列雷达空间构型依赖性强,工程应用难度大。
技术实现思路
本专利技术的目的在于针对分布式雷达在成像过程中对空间构型的依赖性,提出一种分布式阵列雷达迭代自适应高分辨成像方法,通过在目标求解过程中的迭代自适应收敛,降低分布式阵列雷达成像的空间构型依赖性。一种分布式阵列雷达迭代自适应高分辨成像方法,包括以下步骤:S1、根据发射信号和回波信号,构造滤波因子,得到频域回波信号;S2、构建分布式阵列雷达测量矩阵,融合多成像通道的频域回波信号,将所述频域回波信号转化为矩阵形式;S3、构建自相关矩阵,迭代自适应目标散射系数反演;S4、输出目标成像结果。进一步地,所述步骤S1包括:所述分布式阵列雷达包括M个发射机和N个接收机,在球坐标中,第m个发射机的位置为第n个接收机的位置为目标散射点的位置为其中,及分别表示发射机、接收机及点目标的斜距、俯仰角及方位角;在二维成像场景中,发射带宽为B的窄带信号,第m个发射机的发射信号为其中,um(t)、fm∈[-B/2,B/2]及φm分别表示第m个发射机的幅度包络、发射信号载频及初始相位,不同发射信号的幅度包络相互正交;第n个接收机接收到的第m个发射机的回波信号为其中,Ω表示二维成像区域,σ(r)表示在极坐标系内场景中点的目标散射系数,τmn(r)表示发射信号的传播延迟,nmn(t)表示加性高斯白噪声;将所述回波信号gmn(t)和发射信号sm(t)处理后构造滤波因子,滤去与积分变量r无关的参数,得到频域回波信号Ymn(f)其中,σxy表示直角坐标系内r=(x,y)位置处的目标散射系数,Nmn(f)表示加性高斯白噪声的频谱,和表示空间谱,表示为其中,c表示电磁波的传播速度。进一步地,所述步骤S2包括:所述发射信号带宽B内的Ws个频率采样点分别为将成像场景划分为Kx行Ky列共K=Kx×Ky个目标点,第k(k=1,2,...,K)个目标点的目标散射系数为σk,将所有目标散射系数依次排列为列向量形式,即σ=(σ1,σ2,...,σK)T其中,(·)T表示转置运算;构建分布式阵列雷达测量矩阵H,融合M个发射机、目标点、N个接收机构成的个成像通道的回波信号,所述回波频域回波信号转化为矩阵形式,即Y=Hσ+N其中,N表示噪声矩阵,Y=(Y11,Y12,...,Y1N,Y21,...,YMN)T表示回波矢量矩阵,H=(H11,H12,...,H1N,H21,...,HMN)T表示分布式阵列雷达测量矩阵,表示为其中,(ws=1,...,WS;k=1,...,K)。进一步地,所述步骤S3包括:构建自相关矩阵R其中,hk表示矩阵H的第k列数据,矩阵P=diag(p1,p2,...,pK),pk=|σk|2;定义代价函数为其中,Qk=R-pkhk(hk)H表示相关矩阵;初始化矩阵P为单位矩阵,对于k=1,2,...,K,采用最小二乘法得到目标点散射系数为每经过一次迭代,所述自相关矩阵R更新,目标散射系数迭代自适应收敛。本专利技术的有益效果:本专利技术提供了一种分布式阵列雷达迭代自适应高分辨成像方法,首先利用分布式阵列雷达稀疏布阵的空间关系,导出分布式雷达阵列空间谱分布;其次,利用分布式雷达空间谱分布和系统参数,构造分布式阵列雷达系统观测矩阵;最后,为降低成像处理过程的构型依赖性,采用迭代自适应成像处理方法,实现分布式阵列雷达高分辨成像。本专利技术通过目标求解过程中的迭代自适应收敛,降低了分布式阵列雷达对构型的依赖性。附图说明图1为本专利技术实施例的流程图。图2为本专利技术实施例中的分布式雷达成像系统集合构型图。图3为本专利技术实施例中的目标原始场景图。图4为本专利技术实施例中几何构型A空间谱分布图。图5为本专利技术实施例中采用MUSIC方法在几何构型A下的成像结果图。图6为本专利技术实施例中采用本专利技术方法在几何构型A下的成像结果图。图7为本专利技术实施例中几何构型B空间谱分布图。图8为本专利技术实施例中采用MUSIC方法在几何构型B下的成像结果图。图9为本专利技术实施例中采用本专利技术方法在几何构型B下的成像结果图。具体实施方式下面结合附图对本专利技术的实施例做进一步的说明。请参阅图1,本专利技术提出了一种分布式阵列雷达迭代自适应高分辨成像方法,通过以下步骤实现:S1、根据发射信号和回波信号,构造滤波因子,得到频域回波信号和空间谱分布。本实施例中,分布式雷达成像系统几何关系如图2所示,原始成像场景中的5个目标点,其布设方式如图3所示。下表1和表2分别为雷达系统参数和对比参数表。表1分布式雷达系统参数表表2系统比较参数表本实施例中,采用表2中两组参数进行对比。本实施例中,分布式阵列雷达包括M个发射机和N个接收机,在球坐标中,第m个发射机的位置为第n个接收机的位置为目标散射点的位置为其中,及分别表示发射机、接收机及点目标的斜距、俯仰角及方位角。本实施例在二维成像场景下进行,有M=4个发射机均匀处在工作距离100km、俯仰角为0°、方位角为0°~5°的位置上,发射带宽为B=30MHz、载频为fm=8GHz的窄带信号,第m个发射信号为:其中,um(t)、fm∈[-B/2,B/2]及φm分别表示第m个发射机的幅度包络、发射信号载频及初始相位,不同发射信号的幅度包络相互正交,幅度包络满足有N=4个接收机均匀处在工作距离10km、俯仰角为0°的位置上,几何构型A的方位角处在90°~95°,几何构型B的方位角处在135°~140°。第m个发射信号经过目标点反射后,被第n个接收机接收,接收信号表示为:其中,Ω表示二维成像区域,σ(r)表示在极坐标系内场景中点的目标散射系数,表示发射信号经目标点反射后的时间延迟,dmn表示信号传播距离,nmn(t)表示加性高斯白噪声。将回波信号gmn(t)和发射信号sm(t)处理后构造滤波因子,滤去与积分变量r无关的部分,得到频域回波信号Ymn(f)其中,σxy表示直角坐标系内r=(x,y)位置处的目标散射系数,Nmn(f)表示加性高斯白噪声的频谱,和表示分布式雷达二维空间谱,表示为其中,c表示电磁波的传播速度。S2、构建分布式阵列雷达测量矩阵,融合多成像通道的频域回波信号,将频域回波信号转化为矩阵形式。本实施例中,矩阵H为分布式阵列雷达测量矩阵,由发射机、接收机和成像场景的相对位置以及发射信号的形式决定。本实施例中,将成像场景划分为K=Kx×Ky=40×40个目标点,本文档来自技高网...

【技术保护点】
1.一种分布式阵列雷达迭代自适应高分辨成像方法,其特征在于,包括以下步骤:S1、根据发射信号和回波信号,构造滤波因子,得到频域回波信号;S2、构建分布式阵列雷达测量矩阵,融合多成像通道的频域回波信号,将所述频域回波信号转化为矩阵形式;S3、构建自相关矩阵,迭代自适应目标散射系数反演;S4、输出目标成像结果。

【技术特征摘要】
1.一种分布式阵列雷达迭代自适应高分辨成像方法,其特征在于,包括以下步骤:S1、根据发射信号和回波信号,构造滤波因子,得到频域回波信号;S2、构建分布式阵列雷达测量矩阵,融合多成像通道的频域回波信号,将所述频域回波信号转化为矩阵形式;S3、构建自相关矩阵,迭代自适应目标散射系数反演;S4、输出目标成像结果。2.如权利要求1所述的分布式阵列雷达迭代自适应高分辨成像方法,其特征在于,所述步骤S1包括:所述分布式阵列雷达包括M个发射机和N个接收机,在球坐标中,第m个发射机的位置为第n个接收机的位置为目标散射点的位置为其中,及分别表示发射机、接收机及点目标的斜距、俯仰角及方位角;在二维成像场景中,发射带宽为B的窄带信号,第m个发射机的发射信号为其中,um(t)、fm∈[-B/2,B/2]及φm分别表示第m个发射机的幅度包络、发射信号载频及初始相位,不同发射信号的幅度包络相互正交;第n个接收机接收到的第m个发射机的回波信号为其中,Ω表示二维成像区域,σ(r)表示在极坐标系内场景中点的目标散射系数,τmn(r)表示发射信号的传播延迟,nmn(t)表示加性高斯白噪声;将所述回波信号gmn(t)和发射信号sm(t)处理后构造滤波因子,滤去与积分变量r无关的参数,得到频域回波信号Ymn(f)其中,σxy表示直角坐标系内r=(x,y)位置处的目标散射系数,Nmn(f)表示加性高斯白噪...

【专利技术属性】
技术研发人员:张永超徐帆云任叶涵毛德庆黄钰林裴季方张寅郭德明杨建宇
申请(专利权)人:电子科技大学
类型:发明
国别省市:四川,51

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1