matlab仿真点目标 sar合成孔径雷达图像点目标仿真报告(附matlab代码).docx

sar合成孔径雷达图像点目标仿真报告(附matlab代码).docx

SAR图像点目标仿真报告徐一凡1SAR原理简介合成孔径雷达(SyntheticApertureRadar，简称SAR)是一种高分辨率成像雷达技术。它利用脉冲压缩技术获得高的距离向分辨率，利用合成孔径原理获得高的方位向分辨率，从而获得大面积高分辨率雷达图像。SAR回波信号经距离向脉冲压缩后，雷达的距离分辨率由雷达发射信号带宽决定：，式中表示雷达的距离分辨率，表示雷达发射信号带宽，表示光速。同2rrCBrrBC样，SAR回波信号经方位向合成孔径后，雷达的方位分辨率由雷达方位向的多谱勒带宽决定：，式中表示雷达的方位分辨率，表示雷达方位向多谱勒带宽，表示avaaav方位向SAR平台速度。在小斜视角的情况下，方位分辨率近似表示为，其中为2aD方位向合成孔径的长度。2SAR的几何关系雷达位置和波束在地面覆盖区域的简单几何模型如图1所示。此次仿真考虑的是正侧视的条带式仿真，也就是说倾斜角为零，SAR波束中心和SAR平台运动方向垂直的情况。图1雷达数据获取的几何关系建立坐标系XYZ如图2所示，其中XOY平面为地平面；SAR平台距地平面高H，以速度V沿X轴正向匀速飞行；P点为SAR平台的位置矢量，设其坐标为(x,y,z)；T点为目标的位置矢量，设其坐标为；由几何关系，目标与SAR平台的斜距为：(,)Txyz(1)222()()TTRyz由图可知：；令，其中为平台速度，s为慢时间变量(slow0,,0TyzHxvtime)，假设，其中表示SAR平台的x坐标为的时刻；再令，TxvsT2TrHy表示目标与SAR的垂直斜距，重写(1)式为：r(2)220(;)()PRrvs就表示任意时刻时，目标与雷达的斜距。一般情况下，，于是通过(;)Rsrs0vsr傅里叶技术展开，可将(2)式可近似写为：(3)222200(;)()()Rsrvsrs可见，斜距是的函数，不同的目标，也不一样，但当目标距SAR较远时，在观测和带内，可近似认为不变，即。rr图2：空间几何关系(a)正视图(b)侧视图图2(a)中，表示合成孔径长度，它和合成孔径时间的关系是。(b)LsarTsarLsarvT中，为雷达天线半功率点波束角，为波束轴线与Z轴的夹角，即波束视角，为近距点距离，为远距点距离，W为测绘带宽度，它们的关系为：minRmaxR(4)2in()aiHtg3SAR的回波信号模型SAR在运动中以一定的周期()发射和接收信号，具体过程如图3所示。发射1/PRF机以的时间发射啁啾脉冲，然后切换天线开关接收回波信号。l图3雷达发射脉冲串的时序当雷达不处于发射状态时，它接收3反射回波。发射和接收回波的时间序列如图4所示。在机载情况下，每个回波可以在脉冲发射间隔内直接接收到。但是在星载情况下，由于距离过大，某个脉冲的回波要经过6~10个脉冲间隔才能接收到。这里仿真为了方便，默认为机载情况。图4脉冲雷达的发射与接收周期假设为chirp信号持续时间，下标表示距离向；PRF为重复频率，PRT为重复周rTr期，等于。接收序列中，表示发射第个脉冲时，目标回波相对于1/PRF2*(;)nRsCi发射序列的延时。雷达的发射序列数学表达式为式(5):(5)2rjK()()trec)cnjfttspPTte式中，表示矩形信号，为距离向的chirp信号调频率，为载频。rect()Arcf雷达回波信号由发射信号波形，天线方向图，斜距，目标RCS，环境等因素共同决定，若不考虑环境因素，则单点目标雷达回波信号可写成式(6)所示：(6)()()rnnstwptPRT其中，表示点目标的雷达散射截面，表示点目标天线方向图双向幅度加权，表n示载机发射第n个脉冲时，电磁波再次回到载机时的延时，带入式(6)中得：2*(;)nRsrC(7)2(;)/()(exp[;/]4-(;)]exp[()]rnrrcntnPRTsCstwcjKtsjftT式(7)就是单点目标回波信号模型，其中，是2(;/]rKPRsrCchirp分量，它决定距离向分辨率；为多普勒分量，它决定方位向分辨4exp[-(;)]jRs率。对于任意一个脉冲，回波信号可表示为式(8)所示：(8)002(,)(2;)/()exp{4(;)/}exp{}rracstAwsrCwsjfRsrCjKR我们知道，由于随慢时间s的变化而变化，所以计算机记录到的回波数据存储(;)R形式如图5所示：图5目标照射时间内，单个点目标回波能量在信号处理器的二维存储器中的轨迹4距离徙动及校正根据图2可知，在倾斜角为零或很小的时候，目标与雷达的瞬时距离为，根据(;)Rsr几何关系可知，，根据泰勒级数展开可得：220(;)()Rsrvs(9)222200(;)()()vRsrvsrs由式(9)可知，不同慢时间对应着不同的，并且是一个双曲线形式或者近似为;R一个二次形式。如图5所示，同一目标的回波存储在计算机里不在同一直线上，存在距离徙动。从而定义距离徙动量：(10)220(,)()vRsrs为了进行方位向的压缩，方位向的回波数据必须在同一条直线上，也就是说必须校正距离徙动。由式(10)可知，不同的最近距离r对应着不同的，因此在(,)sr(,)Rsr时域处理距离徙动会非常麻烦。因此，对方位向进行傅里叶变换，对距离向不进行变换，得到新的域。由于方位向的频率即为多普勒频率，所以这个新的域也称为距离多普勒域。将斜距R写成多普勒fa的函数，即。众所周知，对最近距离为r的点目标(,)aRfrP，回波多普勒是倾斜角的函数，即，斜距，于是af2sinV(,)/cosaRf(11)2(,)/cos/1()8aarrfrfV所以距离多普勒域中的我距离徙动为=，可发现它不随慢时间变换，(,)Rr21(8arf同一最短距离r对应着相同大小的距离徙动。因此在距离多普勒域对一个距离徙动校正就是对一组具有相同最短距离的点目标的距离徙动校正，这样可以节省运算量。为了对距离徙动进行校正，需要得到距离徙动单元，即距离徙动体现在存储单元中的移动数值，距离徙动单元可以表示为，这个值通常为一个分数，由于存储单(,)/arf元都是离散的，所以不同通过在存储单元简单的移动得到准确的值。为了得到准确的徙动校正值，通常需要进行插值运算。本仿真采用了两种插值方法最近邻点插值和sinc插值，下面分别进行介绍。最近邻点插值法的优点是简单而快速，缺点是不够精确。，其中N为整数部(,)/=arRfn分，n为小数部分，整数部分徙动可以直接通过平移消除，对于小数部分则通过四舍五入的方法变为0或者1，这样就可以得到较为精确的插值。Sinc插值原理如下：在基带信号下，卷积核是sinc函数(12)sin()()i)xhxc插值信号为(13)()()sin)digxcxi即为所有输入样本的加权平均。可通过频域来理解，如图6所示，采样信号的频谱等于以采样率重复的()dgi()dGf信号频谱。为了重建信号，只需要一个周期频谱(