当前位置:首页 > 通信技术 > 通信技术
[导读]现有的雷达成像超分辨算法是基于目标回波信号的二维正弦信号模型,所以模型误差,特别是距离走动误差,将使算法性能严重下降或失效.为此,本文采用距离走动误差下的一阶近似雷达成像二维信号模型,提出了一种基于非线

现有的雷达成像超分辨算法是基于目标回波信号的二维正弦信号模型,所以模型误差,特别是距离走动误差,将使算法性能严重下降或失效.为此,本文采用距离走动误差下的一阶近似雷达成像二维信号模型,提出了一种基于非线性最小二乘准则的参数化超分辨算法.在算法中,距离走动误差补偿与目标参量估计联合进行.文中同时给出了算法估计性能的Cramer-Rao界及仿真结果.

关键词:距离走动误差;补偿;超分辨;雷达成像

A Super Resolution Radar Imaging Algorithm Based on the 2-D Approximate Model

SUN Chang-yin,BAO Zheng

(Kay Laboratory for Radar Signal Processing,Xidian University,Xi’an 710071,China)

Abstract:The recently proposed super resolution radar imaging algorithms,which are based on the 2-D sinusoid signal model,often suffer from the motion through resolution cell error(MTRC) and failed completely.In this paper,an algorithm is proposed based on the 2-D approximate radar imaging model.By minimizing a nonlinear least-squares cost function,the algorithm combines the parameter estimation with the compensation of MTRC errors.The Cramer-Rao bounds are derived and simulation results are also presented to demonstrate the performance of the algorithm.

Key words:motion through resolution cell error;compensation;super resolution;radar imaging

一、引  言

雷达成像基于目标的散射点模型.雷达通常发射长时宽的线频调(chirp)信号,然后用参考信号对回波作解线频调(dechirp)处理,再将解线频调的回波作横向排列,则在一定条件下它可近似为二维正弦信号模型,通过二维傅里叶变换,可以重构目标的二维像;采用超分辨算法[1~3],还可得到更精细的二维目标像.

应当指出,上述二维模型是假设散射点在成像期间不发生超越分辨单元走动,近似认为散射点的移动只影响回波的相移,而子回波包络则固定不变.这种近似,只适用于小观察角时参考点附近有限小尺寸目标成像.

如果目标较大,特别是在离参考点较远处,越分辨单元移动(MTRC)便会发生,从而使得用简单二维模型获得的图像模糊.传统解决的方法是按目标转动用极坐标-直角坐标插值.插值不可避免地会有误差,而超分辨算法通常基于参数化估计,对误差较为敏感,这会影响成像质量.

本文介绍一种近似度较高的二维模型,并利用该模型通过超分辨算法成像,可获得较好的结果.

二、维回波模型

设目标有K个散射点,雷达以平面波自下向上照射目标(图1).目标以参考点为原点相对雷达射线转动,经过N次脉冲发射,散射点Pk点移至P′k点,移动中第n次脉冲时该散射点的垂直坐标为:

ykn=yk+Δykn=xksin(nδθ)+ykcos(nδθ),n=0,1,…,N-1 (1)

式中δθ为相邻脉冲的转角,总观测角Δθ=(N-1)δθ.考虑到雷达发射的是长时宽的线频调信号,以原点为参考作解线频调处理,并对信号以 的频率采样,得目标的回波信号(离散形式)为:

 

 (2)

 

式中Ak为第k个散射点子回波信号的复振幅;fc、γ分别是雷达载频和调频率,c为光速;e(m,n)为加性噪声.

 

 

图1 二维雷达目标几何图

由于观测角Δθ很小,取近似sin(nδθ)≈nδθ和cos(nδθ)≈1,则式(2)可近似写成:

 

 (3)

 

式中

 

式(3)指数项中的第三项是时频耦合项,它是线频调信号(其模糊函数为斜椭圆)所特有的,如果采用窄脉冲发射,则该项不存在.将该项忽略,则式(3)成为常用的回波二维正弦信号模型.

实际上,式(3)的第三项系“距离移动”项,它与散射点的横坐标xk成正比,目标区域大时必须考虑,而且这还远远不够,散射点的多普勒移动也必须考虑.为此,令sin(nδθ)≈nδθ和cos(nδθ)≈1-(nδθ)2/2,则式(2)较精确的近似式可写成:

 


 

式(4)与式(3)相比较,指数中增加了两项,其中前一项是“多普勒移动”项,纵坐标yk越大,影响也越大,这可以补充式(3)之不足;而后项是时频耦合的多普勒移动项,由于Mγ/Fs<

 

 (5)

 

需要指出,每个散射点的参数之间存在下述关系:ωk/μk=2γ/Fsfcδθ2和

k/vk=fcFs/γδθ.由于雷达参数(fc,γ,Fs)和运动参数(δθ)均已知,所以待估计的五个参数中只有三个是独立的.本文假设五个参数是独立的,而在成像计算中已考虑参数之间的关系.

 

设{ξk}Kk=1≡{αk,ωk,

k,μk,vk}Kk=1,现在我们要从y(m,n)中估计参量{ξk}Kk=1.

 

三、二维推广的RELAX算法

对于(5)式所示的信号模型,令:

Y=[y(m,n)]M×N

 (6)

 

式中

 

 

设ξk估计值为

,则ξk的估计问题可通过优化下述代价函数解决:

 

 

 (7)

 

式中‖.‖F表示矩阵的Frobenius范数,⊙表示矩阵的Hadamard积.

上式中C1的最优化是一个多维空间的寻优问题,十分复杂.本文将RELAX[3]算法推广以求解.为此,首先做以下准备工作,令:

 

 (8)

 

即假定{

i}i=1,2,…,K,i≠k已经求出,则式(7)C1的极小化等效于下式的极小化:

 

C2(ξk)=‖Yk-αk(aM(ωk)bTN(

k)Pk)⊙Dk(vk)‖2F (9)

 

令:  Zk=YkP-1k⊙Dk(-vk) (10)

由于Pk为酉矩阵,矩阵Dk的每个元素的模|Dk(m,n)|=1,显然矩阵Yk与Zk的F范数相同,故C2的极小化等效于下式的极小化:

C3=‖Zk-αkaM(ωk)bTN(

k)‖2F (11)

 

对上式关于αk求极小值就获得αk的估计值

k:

 

 

k=aHM(ωk)Zkb*N(

k)/(MN) (12)
 


 

2.SAR成像模拟

雷达参数为:中心频率f0=24.24GHz,调频率γ=33.357×1011Hz/s,带宽B=133.5MHz,脉冲宽度tp=40μs.四个点目标作正方形放置,间隔50米,左下角的点作为参考点.雷达与目标间隔1公里,观察角Δθ=3.15,数据长度为128×128.采用FFT成像方法时,其纵向和横向距离分辨率为ρr=ρa=1.123米,防止MTRC现象发生所需的目标最大范围为[4]:纵向尺寸Dr<4ρ2r/λ=40米,横向尺寸Da<4ρ2a/λ=40米.采用常规超分辨方法时,目标尺寸Dr=Da>10米则出现明显的性能下降.图2、图3分别给出了RELAX方法及本文推广的RELAX(Extended RELAX)算法的成像结果.可以看出,由于目标远离参考中心,已在横向和纵向出现距离走动,采用常规超分辨的RELAX算法产生图像模糊,对于本文算法,则得到基本正确的成像结果.图4和图5则比较了RELAX算法和推广的RELAX算法的散射点强度估计结果,可以看到,RELAX算法由于距离走动影响,散射点(除参考点以外)的强度降低.对于本文算法,散射点强度接近真实值.

 

          

 

图2 距离走动误差下的RELAX成像结果图3 距离走动误差下的

 

           

 

图4 RELAX方法估计的信号强度推广RELAX成像结果图5 推广RELAX方法估计的信号强度

五、结束语

现有的雷达成像超分辨算法是基于目标回波信号的二维正弦信号模型,所以仅适用于目标位于参考点附近很小区域时的情形.当目标远离参考点时,模型误差,特别是距离走动误差,将使算法性能严重下降或失效.为此,本文提出一种基于雷达成像近似二维模型的超分辨算法,从而扩大了超分辨算法的适用范围.本文进一步的工作包括SAR实测数据成像及ISAR机动目标成像,结果将另文报道.

附 录:参数估计的C-R界

下面我们给出式(5)所示的二维信号参量估计的C-R界表达式.同时假设式(5)中加性噪声为零均值高斯色噪声,其协方差矩阵未知.令:

y=vec(Y) (A.1)

e=vec(E) (A.2)

dk=vec(Dk) (A.3)

式中vec(X)=(xT1,xT2,…,xTN)T,向量xn(n=1,2,…,N)为矩阵X的列向量.我们将式(5)改写为如下向量形式:

 

 (A.4)

 

式中

表示Kronecker积,Ω=[{[P1bN(

1)]

aM(ω1)}⊙d1…{[PkbN(

K)]

aM(ωK)}⊙dK],α=(α1,α2,…,αK)T.

 

令Q=E(eeH)为e的协方差矩阵,则对于由式(A.4)所示的二维信号模型,其Fisher信息阵(FIM)的第ij个元素推广的Slepian-Bangs公式为[5,6]:

(FIM)ij=tr(Q-1Q′iQ-1Q′j)+2Re[(αHΩH)′iQ-1(Ωα)′j] (A.5)

式中X′i表示矩阵X对第i个参数求导,tr(X)为矩阵的迹,Re(X)为矩阵的实部.由于Q与Ωα中的参量无关,而Ωα亦与Q的元素无关,显然FIM为一块对角阵.所以待估计参量的C-R界矩阵由(A.5)式的第二项得到.

令:η=([Re(α)]T[Im(α)]TωT

TμTvT)T (A.6)

 

式中ω=(ω1,ω2,…,ωK)T,μ=(μ1,μ2,…,μK)T,

=(

1,

2,…,

K)T,v=(v1,v2,…,vK)T.

 

令:F=[Ω jΩ DωΘ D

Θ DμΘ DvΘ] (A.7)

 

式中矩阵Dω、D

、Dμ、Dv的第k列分别为:

[{[PkbN(

k)]

aM(ωk)}⊙dk]/

ωk、

[{[PkbN(

k)]

aM(ωk)}⊙dk]/

k、

[{[PkbN(

k)]

aM(ωk)}⊙dk]/

μk、

[{[PkbN(

k)]

aM(ωk)}⊙dk]/

vk,Θ=diag{α1 α2 … αK}.则关于参量向量η的CRB矩阵为

 

CRB(η)=[2Re(FHQ-1F)]-1 (A.8)

本站声明: 本文章由作者或相关机构授权发布,目的在于传递更多信息,并不代表本站赞同其观点,本站亦不保证或承诺内容真实性等。需要转载请联系该专栏作者,如若文章内容侵犯您的权益,请及时联系本站删除。
换一批
延伸阅读

9月2日消息,不造车的华为或将催生出更大的独角兽公司,随着阿维塔和赛力斯的入局,华为引望愈发显得引人瞩目。

关键字: 阿维塔 塞力斯 华为

加利福尼亚州圣克拉拉县2024年8月30日 /美通社/ -- 数字化转型技术解决方案公司Trianz今天宣布,该公司与Amazon Web Services (AWS)签订了...

关键字: AWS AN BSP 数字化

伦敦2024年8月29日 /美通社/ -- 英国汽车技术公司SODA.Auto推出其旗舰产品SODA V,这是全球首款涵盖汽车工程师从创意到认证的所有需求的工具,可用于创建软件定义汽车。 SODA V工具的开发耗时1.5...

关键字: 汽车 人工智能 智能驱动 BSP

北京2024年8月28日 /美通社/ -- 越来越多用户希望企业业务能7×24不间断运行,同时企业却面临越来越多业务中断的风险,如企业系统复杂性的增加,频繁的功能更新和发布等。如何确保业务连续性,提升韧性,成...

关键字: 亚马逊 解密 控制平面 BSP

8月30日消息,据媒体报道,腾讯和网易近期正在缩减他们对日本游戏市场的投资。

关键字: 腾讯 编码器 CPU

8月28日消息,今天上午,2024中国国际大数据产业博览会开幕式在贵阳举行,华为董事、质量流程IT总裁陶景文发表了演讲。

关键字: 华为 12nm EDA 半导体

8月28日消息,在2024中国国际大数据产业博览会上,华为常务董事、华为云CEO张平安发表演讲称,数字世界的话语权最终是由生态的繁荣决定的。

关键字: 华为 12nm 手机 卫星通信

要点: 有效应对环境变化,经营业绩稳中有升 落实提质增效举措,毛利润率延续升势 战略布局成效显著,战新业务引领增长 以科技创新为引领,提升企业核心竞争力 坚持高质量发展策略,塑强核心竞争优势...

关键字: 通信 BSP 电信运营商 数字经济

北京2024年8月27日 /美通社/ -- 8月21日,由中央广播电视总台与中国电影电视技术学会联合牵头组建的NVI技术创新联盟在BIRTV2024超高清全产业链发展研讨会上宣布正式成立。 活动现场 NVI技术创新联...

关键字: VI 传输协议 音频 BSP

北京2024年8月27日 /美通社/ -- 在8月23日举办的2024年长三角生态绿色一体化发展示范区联合招商会上,软通动力信息技术(集团)股份有限公司(以下简称"软通动力")与长三角投资(上海)有限...

关键字: BSP 信息技术
关闭
关闭