雷达成像近似二维模型及其超分辨算法

2011-04-18 23:51:34来源: 互联网

现有的雷达成像超分辨算法是基于目标回波信号的二维正弦信号模型,所以模型误差,特别是距离走动误差,将使算法性能严重下降或失效.为此,本文采用距离走动误差下的一阶近似雷达成像二维信号模型,提出了一种基于非线性最小二乘准则的参数化超分辨算法.在算法中,距离走动误差补偿与目标参量估计联合进行.文中同时给出了算法估计性能的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)δθ.考虑到雷达发射的是长时宽的线频调信号,以原点为参考作解线频调处理,并对信号以 的频率采样,得目标的回波信号(离散形式)为:

g84.gif (2442 bytes) (2)

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

t84-1.gif (1497 bytes)

图1 二维雷达目标几何图

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

g85-1.gif (1711 bytes) (3)

式中g85-2.gif (622 bytes)
  式(3)指数项中的第三项是时频耦合项,它是线频调信号(其模糊函数为斜椭圆)所特有的,如果采用窄脉冲发射,则该项不存在.将该项忽略,则式(3)成为常用的回波二维正弦信号模型.
  实际上,式(3)的第三项系“距离移动”项,它与散射点的横坐标xk成正比,目标区域大时必须考虑,而且这还远远不够,散射点的多普勒移动也必须考虑.为此,令sin(nδθ)≈nδθ和cos(nδθ)≈1-(nδθ)2/2,则式(2)较精确的近似式可写成:

g85-3.gif (2315 bytes) (4)

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

g85-4.gif (1774 bytes) (5)

  需要指出,每个散射点的参数之间存在下述关系:ωk/μk=2γ/Fsfcδθ2和ts85-1.gif (92 bytes)k/vk=fcFs/γδθ.由于雷达参数(fc,γ,Fs)和运动参数(δθ)均已知,所以待估计的五个参数中只有三个是独立的.本文假设五个参数是独立的,而在成像计算中已考虑参数之间的关系.
  设{ξk}Kk=1≡{αk,ωk,ts85-1.gif (92 bytes)k,μk,vk}Kk=1,现在我们要从y(m,n)中估计参量{ξk}Kk=1.

三、二维推广的RELAX算法
  对于(5)式所示的信号模型,令:

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

g85-5.gif (1215 bytes) (6)
式中

g85-6.gif (3131 bytes)

  设ξk估计值为g85-7.gif (561 bytes),则ξk的估计问题可通过优化下述代价函数解决:

g85-8.gif (1588 bytes) (7)

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

g85-9.gif (1243 bytes) (8)

  即假定{ts85-6.gif (102 bytes)i}i=1,2,…,K,i≠k已经求出,则式(7)C1的极小化等效于下式的极小化:

C2(ξk)=‖Yk-αk(aM(ωk)bTN(ts85-1.gif (92 bytes)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(ts85-1.gif (92 bytes)k)‖2F (11)

  对上式关于αk求极小值就获得αk的估计值ts82.gif (98 bytes)k:

ts82.gif (98 bytes)k=aHM(ωk)Zkb*N(ts85-1.gif (92 bytes)k)/(MN) (12)

  从式(12)可以看出:ts82.gif (98 bytes)是Zk归一化的二维离散傅里叶变换在{ωk,ts85-1.gif (92 bytes)k}处的值,所以只要得到估计值{ts85-2.gif (93 bytes)k,ts85-3.gif (102 bytes)k,ts85-4.gif (105 bytes)k,ts85-5.gif (93 bytes)k},即可通过2D-FFT获得ts82.gif (98 bytes)k.
  将估计值ts82.gif (98 bytes)k代入式(11)后,估计值{ts85-2.gif (93 bytes)k,ts85-3.gif (102 bytes)k,ts85-4.gif (105 bytes)k,ts85-5.gif (93 bytes)k}可由下式寻优得到:

g85-10.gif (1590 bytes) (13)

  由上式可见,对于固定的{μk,vk}取值,估计值{ts85-2.gif (93 bytes)k,ts85-3.gif (102 bytes)k}为归一化的周期图|aHM(ωk)Zkb*N(ts85-1.gif (92 bytes)k)|2/(MN)主峰处的二维频率值.这样,式(13)的优化问题归结为:在(μk,vk)平面上可能的取值范围内寻找一点{ts85-4.gif (105 bytes)k,ts85-5.gif (93 bytes)k},在该点处周期图|aHM(ωk)Zkb*N(ts85-1.gif (92 bytes)k)|2/(MN)的主峰值比其余各点处的主峰值都大.所以,我们通过上述二维寻优获得{μk,vk}的估计值{ts85-4.gif (105 bytes)k,ts85-5.gif (93 bytes)k},再由式(13)得到{ωk,ts85-1.gif (92 bytes)k}的估计值{ts85-2.gif (93 bytes)k,ts85-3.gif (102 bytes)k}.
  实际中,为了加快运算速度,二维(μk,vk)平面的寻优可以用Matlab中的函数Fmin()实现.
  在做了以上的准备工作以后,基于推广的RELAX算法的参量估计步骤如下:
  第一步:假设信号数K=1,分别利用式(13)和式(12)计算ts85-6.gif (102 bytes)1.
  第二步(2):假设信号数K=2,首先将第一步计算所得到的ts85-6.gif (102 bytes)1代入式(8)求出Y2,再利用式(13)和式(12)计算ts85-6.gif (102 bytes)2;将计算的ts85-6.gif (102 bytes)2代入式(8)求出Y1,然后利用式(13)和式(12)重新计算ts85-6.gif (102 bytes)1,这个过程反复叠代,直至收敛.
  第三步:假设信号数K=3,首先将第二步计算所得到的ts85-6.gif (102 bytes)1和ts85-6.gif (102 bytes)2代入式(8)求出Y3,再利用式(13)和式(12)计算ts85-6.gif (102 bytes)3;将计算的ts85-6.gif (102 bytes)3和ts85-6.gif (102 bytes)2代入式(8)求出Y1,然后利用式(13)和式(12)重新计算ts85-6.gif (102 bytes)1;将计算的ts85-6.gif (102 bytes)1和ts85-6.gif (102 bytes)3代入式(8)求出Y2,然后利用式(13)和式(12)重新计算ts85-6.gif (102 bytes)2,这个过程反复叠代,直至收敛.
  剩余步骤:令K=K+1,上述步骤持续进行,直到K等于待估计信号数.
  上述过程中的收敛判据与RELAX算法的收敛判据相同,即比较代价函数C1在两次叠代过程中的变化值,如果这个变换值小于某个值,如ε=10-3,则认为过程收敛.

四、数值模拟
  1.算法参数估计性能模拟
  模拟数据由式(5)产生,M=10,N=10,信号数K=2.信号参数和实验条件如表1所示,为复高斯白噪声.注意两信号的频率差小于FFT的分辨率Δf=Δω/(2π)=0.1.表1给出了信号参数估计均方根误差的统计结果及相应情形时的C-R界,可见,估计均方根误差与CR界十分接近.另外表中还给出了估计均值,与真实值也非常接近.

表1 二维信号的参数估计、CRB及与均方根差的比较

  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算法由于距离走动影响,散射点(除参考点以外)的强度降低.对于本文算法,散射点强度接近真实值.

     t84-2.gif (2232 bytes)          t86-1.gif (2338 bytes)

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

图3 距离走动误差下的

     t86-2.gif (3010 bytes)           t86-3.gif (2871 bytes)

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

图5 推广RELAX方法估计的信号强度

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

[1] [2]

关键字:分辨  算法

编辑:神话 引用地址:http://www.eeworld.com.cn/mndz/2011/0418/article_7238.html
本网站转载的所有的文章、图片、音频视频文件等资料的版权归版权所有人所有,本站采用的非本站原创文章及图片等内容无法一一联系确认版权者。如果本网所选内容的文章作者及编辑认为其作品不宜公开自由传播,或不应无偿使用,请及时通过电子邮件或电话通知我们,以迅速采取适当措施,避免给双方造成不必要的经济损失。
论坛活动 E手掌握
微信扫一扫加关注
论坛活动 E手掌握
芯片资讯 锐利解读
微信扫一扫加关注
芯片资讯 锐利解读
推荐阅读
全部
分辨
算法

小广播

独家专题更多

富士通铁电随机存储器FRAM主题展馆
富士通铁电随机存储器FRAM主题展馆
馆内包含了 纵览FRAM、独立FRAM存储器专区、FRAM内置LSI专区三大部分内容。 
走,跟Molex一起去看《中国电子消费品趋势》!
走,跟Molex一起去看《中国电子消费品趋势》!
 
带你走进LED王国——Microchip LED应用专题
带你走进LED王国——Microchip LED应用专题
 
电子工程世界版权所有 京ICP证060456号 京ICP备10001474号 电信业务审批[2006]字第258号函 京公海网安备110108001534 Copyright © 2005-2016 EEWORLD.com.cn, Inc. All rights reserved