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

1、精品文档SAR图像点目标仿真报告徐一凡1 SAR原理简介合成孔径雷达(SyntheticApertureRadar ，简称 SAR)是一种高分辨率成像雷达技术。它利用脉冲压缩技术获得高的距离向分辨率，利用合成孔径原理获得高的方位向分辨率，从而获得大面积高分辨率雷达图像。SAR 回波信号经距离向脉冲压缩后，雷达的距离分辨率由雷达发射信号带宽决定：Cr 表示雷达的距离分辨率，Br 表示雷达发射信号带宽，C 表示光速。同r，式中2Br样，SAR回波信号经方位向合成孔径后， 雷达的方位分辨率由雷达方位向的多谱勒带宽决定：ava ，式中a 表示雷达的方位分辨率，Ba 表示雷达方位向多谱勒带宽，va 表示

2、方位向BaSAR平台速度。在小斜视角的情况下，方位分辨率近似表示为成孔径的长度。Da，其中 D 为方位向合22 SAR的几何关系雷达位置和波束在地面覆盖区域的简单几何模型如图1 所示。此次仿真考虑的是正侧视的条带式仿真，也就是说倾斜角为零，SAR波束中心和SAR平台运动方向垂直的情况。图 1 雷达数据获取的几何关系建立坐标系XYZ如图 2 所示，其中XOY平面为地平面；SAR平台距地平面高H，以速度V 沿 X 轴正向匀速飞行；P 点为 SAR平台的位置矢量，设其坐标为(x,y,z) ； T 点为目标的位置矢量，设其坐标为( xT , yT , zT ) ；由几何关系，目标与SAR平台的斜距为：

3、uuur( x xT )2( y yT )2zT )2RPT( z(1)由图可知： y 0, zH , zT0 ；令 xv s ，其中 v 为平台速度，s 为慢时间变量（ slow。1欢迎下载精品文档time ），假设 xTvs，其中 s 表示 SAR平台的 x 坐标为 xT 的时刻； 再令 rH 2yT2， r表示目标与 SAR的垂直斜距，重写 (1)式为：uuurr 2v2 (ss0 )2PTR(s ; r )(2)R( s ; r ) 就表示任意时刻s时，目标与雷达的斜距。一般情况下，v s s0r ，于是通过傅里叶技术展开，可将(2)式可近似写为：R(s ; r )r 2v2(ss0

4、)2rv2(s s0 )2(3)2r可见，斜距是 s和r 的函数，不同的目标，r 也不一样，但当目标距SAR较远时，在观测带内，可近似认为r 不变，即 r R0 。图 2：空间几何关系(a)正视图 (b) 侧视图图 2(a) 中， Lsar 表示合成孔径长度，它和合成孔径时间Tsar 的关系是 LsarvTsar。(b)中，为雷达天线半功率点波束角，为波束轴线与Z 轴的夹角，即波束视角，R min 为近距点距离， Rmax为远距点距离， W为测绘带宽度，它们的关系为：R minH tg (2)R maxH tg (2 )(4)WR maxR min3 SAR的回波信号模型SAR在运动中以一定的

5、周期( 1/ PRF ) 发射和接收信号，具体过程如图3 所示。发射机。2欢迎下载精品文档以 l 的时间发射啁啾脉冲，然后切换天线开关接收回波信号。图 3 雷达发射脉冲串的时序当雷达不处于发射状态时，它接收3 反射回波。发射和接收回波的时间序列如图4 所示。在机载情况下，每个回波可以在脉冲发射间隔内直接接收到。但是在星载情况下，由于距离过大，某个脉冲的回波要经过610 个脉冲间隔才能接收到。这里仿真为了方便，默认为机载情况。图 4脉冲雷达的发射与接收周期假设 Tr 为 chirp 信号持续时间，下标r 表示距离向； PRF为重复频率， PRT为重复周期，等于 1/ PRF 。接收序列中，n2*

6、R( s; r )表示发射第 i 个脉冲时，目标回波相对于C发射序列的延时。雷达的发射序列数学表达式为式(5):s(t )p(tn * PRT )n（5）rect( tp(t)e j K r t 2e j 2 f ctTr式中， rect(g) 表示矩形信号，K r为距离向的 chirp信号调频率， fc 为载频。雷达回波信号由发射信号波形，天线方向图，斜距，目标RCS，环境等因素共同决定，若不考虑环境因素，则单点目标雷达回波信号可写成式(6)所示：sr (t )nwp(tnPRTn )(6)其中，表示点目标的雷达散射截面，w 表示点目标天线方向图双向幅度加权，n 表示载机发射第n 个脉冲时，

7、电磁波再次回到载机时的延时2* R(s; r )n，带入式 (6)中C得：。3欢迎下载精品文档s (t )w rect( tn PRT2R( s; r ) / C )rnTrexp jK r (t nPRT 2R( s; r ) / C)2 (7)exp- j 4R(s; r )exp j 2f c (tnPRTn )式 (7)就是单点目标回波信号模型，其中，exp jK r ( tnPRT2R( s; r ) / C ) 2 是chirp 分量，它决定距离向分辨率；exp- j4R(s; r ) 为多普勒分量，它决定方位向分辨率。对于任意一个脉冲，回波信号可表示为式(8) 所示：sr (t,

8、 s)A0 wr (2R( s; r ) / C)wa (s sc ) expj 4f0 R(s; r ) / Cexp j K r (2R(s; r ) / C) 2（8）我们知道，由于R(s; r ) 随慢时间 s 的变化而变化，所以计算机记录到的回波数据存储形式如图5 所示：图 5 目标照射时间内，单个点目标回波能量在信号处理器的二维存储器中的轨迹4 距离徙动及校正根据图 2 可知，在倾斜角为零或很小的时候，目标与雷达的瞬时距离为R( s; r ) ，根据几何关系可知，R(s ;r )r 2v2 (ss0 )2 ，根据泰勒级数展开可得：。4欢迎下载精品文档R(s ; r )r 2v2 (

9、ss0 )2rv2(s s0 )2(9)2r由式 (9) 可知，不同慢时间对应着不同的R(s ; r ) ，并且是一个双曲线形式或者近似为一个二次形式。如图 5 所示，同一目标的回波存储在计算机里不在同一直线上，存在距离徙动。从而定义距离徙动量：v2(s s0 )2（ 10）R( s, r )2r为了进行方位向的压缩，方位向的回波数据必须在同一条直线上，也就是说必须校正距离徙动R(s,r ) 。由式（ 10）可知，不同的最近距离r 对应着不同的R(s, r ) ，因此在时域处理距离徙动会非常麻烦。因此，对方位向进行傅里叶变换，对距离向不进行变换，得到新的域。由于方位向的频率即为多普勒频率，所以

10、这个新的域也称为距离多普勒域。将斜距 R 写成多普勒 fa 的函数，即 R( f a , r ) 。众所周知，对最近距离为r 的点目标P，回波多普勒 f a 是倾斜角的函数，即 f a2V sin，斜距 R( fa , r ) r / cos ，于是R( f a , r )r / cosr /1sin 2r /1(fa)2(11)2Vr1() 2 rf a28 V所以距离多普勒域中的我距离徙动为(f a,r)122 ，可发现它不随慢时间变换，R= ()rf a8 V同一最短距离 r 对应着相同大小的距离徙动。因此在距离多普勒域对一个距离徙动校正就是对一组具有相同最短距离的点目标的距离徙动校正，

11、这样可以节省运算量。为了对距离徙动进行校正，需要得到距离徙动单元，即距离徙动体现在存储单元中的移动数值，距离徙动单元可以表示为R( f a , r ) /r ，这个值通常为一个分数，由于存储单元都是离散的，所以不同通过在存储单元简单的移动得到准确的值。为了得到准确的徙动校正值，通常需要进行插值运算。本仿真采用了两种插值方法最近邻点插值和sinc 插值，下面分别进行介绍。最近邻点插值法的优点是简单而快速，缺点是不够精确。R( fa ,r ) /r =Nn ，其中 N 为整数部分， n 为小数部分，整数部分徙动可以直接通过平移消除，对于小数部分则通过四舍五入的方法变为 0 或者 1，这样就可以得到

12、较为精确的插值。Sinc 插值原理如下：在基带信号下，卷积核是sinc 函数sin(x)h( x)sin c(x)（12）x插值信号为。5欢迎下载精品文档g( x)gd (i )sin c( xi )（13）i即为所有输入样本的加权平均。可通过频域来理解，如图6 所示，采样信号 gd (i ) 的频谱 Gd ( f ) 等于以采样率重复的信号频谱。为了重建信号g (x) ，只需要一个周期频谱（如基带周期），因此需要理想矩形低通滤波器在频域中提取基带频谱（如图6）所示。已知该理想滤波器在时域中是sinc 函数。由于频域相乘相当于时域卷积，故插值可以通过与sinc 核的卷积来实现。图 6 理想低通

13、滤波器怎样对采样信号进行插值5 点目标成像 matlab 仿真5.1 距离多普勒算法距离多普勒算法 （ RDA）是在 1976 年至 1978 年为民用星载SAR提出的，它兼顾了成熟、简单、 高效和精确等因素，至今仍是使用最广泛的成像算法。它通过距离和方位上的频域操作，到达了高效的模块化处理要求，同时又具有了一维操作的简便性。图 7 示意了 RDA的处理流程。这里主要讨论小倾斜角及短孔径下的基本RDA处理框图。1. 当数据处在方位时域时， 可通过快速卷积进行距离压缩。 也就是说， 距离 FFT 后随即进行距离向匹配滤波，再利用距离IFFT 完成距离压缩。回波信号为：s0 (t , s)A0 wr t2R(s) / cwa (ssc )exp- j 4f0 R(s) / cexpj K r (t - 2R( s) / c)2 （ 14）距离向压缩后的信号为：src (t , s)IFFTt S0 ( f t, s) H ( f t )（ 15）A0r t2R(s) / cwa (ssc )ex