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

SAR 图像点目标仿真报告徐一凡1 SAR 原理简介合成孔径雷达(Synthetic Aperture Radar ，简称SAR)是一种高分辨率成像雷达技术。

它利用脉冲压缩技术获得高的距离向分辨率，利用合成孔径原理获得高的方位向分辨率，从而获得大面积高分辨率雷达图像。

SAR 回波信号经距离向脉冲压缩后，雷达的距离分辨率由雷达发射信号带宽决定：2r rC B ρ=，式中r ρ表示雷达的距离分辨率，r B 表示雷达发射信号带宽，C 表示光速。

同样，SAR 回波信号经方位向合成孔径后，雷达的方位分辨率由雷达方位向的多谱勒带宽决定：a a av B ρ=，式中a ρ表示雷达的方位分辨率，a B 表示雷达方位向多谱勒带宽，a v 表示方位向SAR 平台速度。

在小斜视角的情况下，方位分辨率近似表示为2a D ρ=，其中D 为方位向合成孔径的长度。

2 SAR 的几何关系雷达位置和波束在地面覆盖区域的简单几何模型如图1所示。

此次仿真考虑的是正侧视的条带式仿真，也就是说倾斜角为零，SAR 波束中心和SAR 平台运动方向垂直的情况。

图1 雷达数据获取的几何关系建立坐标系XYZ 如图2所示，其中XOY 平面为地平面；SAR 平台距地平面高H ，以速度V 沿X 轴正向匀速飞行；P 点为SAR 平台的位置矢量，设其坐标为(x,y,z)； T 点为目标的位置矢量，设其坐标为(,,)T T T x y z ；由几何关系，目标与SAR 平台的斜距为：(R PT x ==(1)由图可知：0,,0T y z H z ===；令x v s =⋅，其中v 为平台速度，s 为慢时间变量（slowtime ），假设T x vs =，其中s 表示SAR 平台的x 坐标为T x的时刻；再令r =，r 表示目标与SAR 的垂直斜距，重写(1)式为：(;)PT R s r = = (2) (;)R s r 就表示任意时刻s 时，目标与雷达的斜距。

一般情况下，0vs s r -<<，于是通过傅里叶技术展开，可将(2)式可近似写为：220(;)()2v R s r r s s r =≈+- (3) 可见，斜距是s r 和的函数，不同的目标，r 也不一样，但当目标距SAR 较远时，在观测带内，可近似认为r 不变，即0r R =。

图2：空间几何关系 (a)正视图 (b)侧视图图2(a)中，Lsar 表示合成孔径长度，它和合成孔径时间Tsar 的关系是Lsar vTsar =。

(b)中，θ∆为雷达天线半功率点波束角，θ为波束轴线与Z 轴的夹角，即波束视角，min R 为近距点距离，max R 为远距点距离，W 为测绘带宽度，它们的关系为：22min (max ()max minR H tg R H tg W R R θθθθ∆∆=⋅-=⋅+=- (4)3 SAR 的回波信号模型SAR 在运动中以一定的周期(1/PRF )发射和接收信号，具体过程如图3所示。

发射机以l τ的时间发射啁啾脉冲，然后切换天线开关接收回波信号。

图3 雷达发射脉冲串的时序当雷达不处于发射状态时，它接收3反射回波。

发射和接收回波的时间序列如图4所示。

在机载情况下，每个回波可以在脉冲发射间隔内直接接收到。

但是在星载情况下，由于距离过大，某个脉冲的回波要经过6~10个脉冲间隔才能接收到。

这里仿真为了方便，默认为机载情况。

图4 脉冲雷达的发射与接收周期假设r T 为chirp 信号持续时间，下标r 表示距离向；PRF 为重复频率，PRT 为重复周期，等于1/PRF 。

接收序列中，2*(;)n R s r Cτ=表示发射第i 个脉冲时，目标回波相对于发射序列的延时。

雷达的发射序列数学表达式为式(5):2r 2j K ()(*)(t)rect()e c n j f tt rs t p t n PRT t p e T ππ∞=-∞=-=∑ （5） 式中，rect()表示矩形信号，r K 为距离向的chirp 信号调频率，c f 为载频。

雷达回波信号由发射信号波形，天线方向图，斜距，目标RCS ，环境等因素共同决定，若不考虑环境因素，则单点目标雷达回波信号可写成式(6)所示：()()r n n s t wp t n PRT στ∞=-∞=-⋅-∑ (6) 其中，σ表示点目标的雷达散射截面，w 表示点目标天线方向图双向幅度加权，n τ表示载机发射第n 个脉冲时，电磁波再次回到载机时的延时2*(;)n R s r Cτ=，带入式(6)中得：22(;)/()() exp[(2(;)/)]4 exp[-(;)]exp[2()]r n r r c n t n PRT R s r C s t w rect T j K t n PRT R s r C jR s r j f t n PRT σπππτλ∞=-∞-⋅-=⋅⋅⋅-⋅-⋅⋅-⋅-∑ (7) 式(7)就是单点目标回波信号模型，其中，2exp[(2(;)/)]r j K t n PRT R s r C π-⋅-是chirp 分量，它决定距离向分辨率；4exp[-(;)]j R s r πλ为多普勒分量，它决定方位向分辨率。

对于任意一个脉冲，回波信号可表示为式(8)所示：002(,)(2(;)/)()exp{4(;)/} exp{(2(;)/)}r r a c r s t s A w R s r C w s s j f R s r C j K R s r C τππτ=--⨯-⋅-（8）我们知道，由于(;)R s r 随慢时间s 的变化而变化，所以计算机记录到的回波数据存储形式如图5所示：图5 目标照射时间内，单个点目标回波能量在信号处理器的二维存储器中的轨迹 4 距离徙动及校正根据图2可知，在倾斜角为零或很小的时候，目标与雷达的瞬时距离为(;)R s r ，根据几何关系可知，(;)R s r =，根据泰勒级数展开可得：22(;)()2vR s r r s sr=≈+-(9)由式(9)可知，不同慢时间对应着不同的(;)R s r ，并且是一个双曲线形式或者近似为一个二次形式。

如图5所示，同一目标的回波存储在计算机里不在同一直线上，存在距离徙动。

从而定义距离徙动量：22(,)()2vR s r s sr∆=-（10）为了进行方位向的压缩，方位向的回波数据必须在同一条直线上，也就是说必须校正距离徙动(,)R s r∆。

由式（10）可知，不同的最近距离r对应着不同的(,)R s r∆，因此在时域处理距离徙动会非常麻烦。

因此，对方位向进行傅里叶变换，对距离向不进行变换，得到新的域。

由于方位向的频率即为多普勒频率，所以这个新的域也称为距离多普勒域。

将斜距R写成多普勒fa的函数，即(,)aR f r。

众所周知，对最近距离为r的点目标P，回波多普勒af是倾斜角θ的函数，即2sinaVfθλ=，斜距(,)/cosaR f r rθ=，于是22(,)/cos//1()8aaR f r r rrr rfVθλ===≈+(11)所以距离多普勒域中的我距离徙动为(,)aR f r∆=221()8arfVλ，可发现它不随慢时间变换，同一最短距离r对应着相同大小的距离徙动。

因此在距离多普勒域对一个距离徙动校正就是对一组具有相同最短距离的点目标的距离徙动校正，这样可以节省运算量。

为了对距离徙动进行校正，需要得到距离徙动单元，即距离徙动体现在存储单元中的移动数值，距离徙动单元可以表示为(,)/a rR f rρ∆，这个值通常为一个分数，由于存储单元都是离散的，所以不同通过在存储单元简单的移动得到准确的值。

为了得到准确的徙动校正值，通常需要进行插值运算。

本仿真采用了两种插值方法最近邻点插值和sinc插值，下面分别进行介绍。

最近邻点插值法的优点是简单而快速，缺点是不够精确。

(,)/=a rR f r N nρ∆+，其中N为整数部分，n为小数部分，整数部分徙动可以直接通过平移消除，对于小数部分则通过四舍五入的方法变为0或者1，这样就可以得到较为精确的插值。

Sinc插值原理如下：在基带信号下，卷积核是sinc函数sin()()sin()xh x c xxππ==（12）插值信号为()()sin ()d i g x g i c x i =-∑（13）即为所有输入样本的加权平均。

可通过频域来理解，如图6所示，采样信号()d g i 的频谱()d G f 等于以采样率重复的信号频谱。

为了重建信号()g x ，只需要一个周期频谱（如基带周期），因此需要理想矩形低通滤波器在频域中提取基带频谱（如图6）所示。

已知该理想滤波器在时域中是sinc 函数。

由于频域相乘相当于时域卷积，故插值可以通过与sinc 核的卷积来实现。

图6 理想低通滤波器怎样对采样信号进行插值5 点目标成像matlab 仿真5.1距离多普勒算法距离多普勒算法（RDA ）是在1976年至1978年为民用星载SAR 提出的，它兼顾了成熟、简单、高效和精确等因素，至今仍是使用最广泛的成像算法。

它通过距离和方位上的频域操作，到达了高效的模块化处理要求，同时又具有了一维操作的简便性。

图7示意了RDA 的处理流程。

这里主要讨论小倾斜角及短孔径下的基本RDA 处理框图。

1.当数据处在方位时域时，可通过快速卷积进行距离压缩。

也就是说，距离FFT 后随即进行距离向匹配滤波，再利用距离IFFT 完成距离压缩。

回波信号为：0020(,)[2()/]()exp{-4()/}exp{(-2()/)}r a c r s t s A w t R s c w s s j f R s c j K t R s c ππ=--⨯（14）距离向压缩后的信号为： 000(,){(,)()}[2()/]()exp{4()/}rc t t t r a c s t s IFFT S f s H f A t R s c w s s j f R s c ρπ==--- （15）20(){}exp{}exp{2}||t f f H f rect j j ft K T Kππ=-- （16） 2.通过方位FFT 将数据变换至距离多普勒域，多普勒中心频率估计以及大部分后续操作都在该域进行。

方位向傅里叶变换后信号为：10200(,){(,)}2() []()4 exp{-}exp{}c s s rc rd s r a s s s a S t f FFT s t s R f A p t W f f cf R f j j c K ππ==--⨯ （17） 3.在距离多普勒域进行随距离时间及方位频率变化的RCMC ，该域中同一距离上的一组目标轨迹相互重合。

RCMC 将距离徙动曲线拉直到与方位频率轴平行的方向。

这里可以采用最近邻点插值法或者sinc 插值法，具体插值方法见前面。