stft isar成像 matlab,基于时频分析的ISAR成像

1引言雷达目标的回波具有时变性,因此常用的频域或时域处理方法往往力不从心。解决该问题的主要工具联合时频技术应运而生。逆合成孔径雷达(ISAR)成像的基本方法为距离一多普勒法,距离一多普勒法采用DFT对回波作多普勒分析,假设观测期间其多普勒频率为常数,当目标做平稳运动时能够满足,但如果目标做机动运动,那么回波的多普勒频率将是时变的,这就需要对回波数据作时频分析,以得到各个时刻的瞬时多普勒频率分布。其中常用的时频分析有基于短时傅里叶变换(STFT)的成像方法,基于Wigner-Ville分布的瞬时成像方法和伪WignerVille分布。魏格纳-维尔分布(WVD)具有良好的时频集聚性,它可以获得高分辨的距离-瞬时多普勒像,但其主要问题是存在交叉项的干扰。目标任一距离单元内都有多个散射点,对于某一时刻的瞬时多普勒谱,交叉项混杂在其中,使其不能正确成像。在此对WVD变换进行改进,运用两个窗函数分别对WVD分布进行平滑操作称为平滑伪WignerVille分布,通过此方法可以获得最后抑制了交叉项干扰同时又保持较高时频分辨率的ISAR像。2ISAR成像的基本原理雷达目标成像的几何模型如图1所示,目标以O为旋转中心,X轴为方位轴,Y为雷达视线方向。图1ISAR成像原理设雷达距旋转中心的距离为0,目标绕O点旋转的角频率为,则目标上点P到雷达的距离为:r=[r20+r2a+2r0rasin(0+t)]12(1)当雷达与目标之间的距离远大于目标尺寸时,r的近似表达为:r=ra+x0sin(t)+y0cos(t)(2)假设雷达的发射信号为s(t)=ar(t)exp(j2(fct+12t2))(3)式中t为距离快时间,fc为发射信号载频,为发射信号的调频率。则雷达接收到P点的回波基频信号为:s(t,tm)=ar(t-2R(tm)c)aa(tm)exp(j(t-2R(tm)c)2)exp(-j4R(tm))(4)式中ar()和aa()分别为雷达线性调频信号的窗函数和方位窗函数,c为光速,为波长。经过距离匹配滤波之后,其变为s(t,tm)=Asinc[fr(t-2R(tm)c)]exp(-j4R(tm))(5)对由平动所引起的包络移动和越距离单元走动进行补偿。补偿之后为:s(t,tm)=Asinc[fr(t-2c(Ra+y0))]exp(-j4(Ra+y0))exp(j4x0tm)(6)然后对其进行方位向的FFT变换,就可得到目标的高分辨图像:|s(t,tm)|=|A'||sinc[fr(t-2c(Ra+y0))]||sinc[Ta2(fa+2x0)]|(7)3机动目标ISAR成像的时频分布成像算法中最主要的问题是如何平衡效率与分辨率。若想得到高的方位分辨率就只能牺牲效率。如果目标作机动运动,旋转速率不均匀或旋转轴变化时,运动过程就会变的很复杂,致使各个散射点子回波的多普勒频率发生变化,从而导致ISAR像变得模糊。如果能够得到各个距离单元里所有散射点子回波的时频分布,就可以从各个时刻的瞬时多普勒分布得到相应时刻的瞬时ISAR像。当目标作机动飞行时,在观测期间内目标回波的多普勒频率将是时变的,如果仍用传统的距离多普勒方法对其进行成像,将会使图像变的模糊不清,此时应当用时频分布对其作分析,在每一个距离单元里存在多个散射点,对该单元的回波数据进行时频分析,能够得到各散射点子回波的瞬时多普勒频率分布,从而构成目标在各个时刻的ISAR成像图,即距离-瞬时多普勒分布图。其流程框图如图2所示。图2距离-瞬时多普勒成像原理4ISAR成像的时频分析方法在ISAR的成像过程中,每一个距离单元都有多