传统SAR算法研究报告文献.docx

传统SAR算法研究报告 丁晟 2015年2月 概述 不管是在模拟器领域（回波和雷达），还是在RCS测量、机载雷达，抑或是在地面监控雷达领域，合成孔径雷达将会成为公司的一个重点发展方向。事实上，基于多个收发通道的合成雷达可以视为更加广义的合成孔径雷达。因此，从这个意义上说，不管是MIMO、DBF、数字相控阵、组网雷达等体制都可以看成是广义合成孔径雷达，其处理的最终目的是在多个维度上进行“去高次”相位处理，然后通过傅里叶变化得到具有sinc形状的压缩波形，以提高雷达的分辨率（在很多场合，被视为获得极窄的波束）。因此，对合成孔径雷达的深入研究是极有必要的。 本文第二章节将回顾与SAR相关的信号处理内容，特别是提出一些我个人对于脉冲压缩理论的理解，这一章节的内容将会贯穿全文；第三章节将给出SAR回波信号的数学表达式；第四章节将给出距离-多普勒算法（RDA）和ω-k算法（ωkA）解???及说明；附件将给出所有的MATLAB源代码。 信号处理 脉冲压缩的理解 连续时间的傅里叶变换公式如下： 式  SEQ 式 \* ARABIC 1 式  SEQ 式 \* ARABIC 2 详细的傅里叶变换说明请参考相关信号处理教科书。本文采用傅里叶变换来重新解释脉冲压缩的基本理论。 当一个矩形脉冲进行傅里叶变换后，将会得到一个类似sinc的波形，如下图所示 图 SEQ 图 \* ARABIC 1 矩形脉冲与sinc函数的傅里叶变换关系 矩形脉冲与sinc函数的这一关系似乎可以用来进行脉冲压缩，但事实上，由于时域的尺度和频域的尺度具有一一对应的关系，因此在归一化以后，实质上sinc函数的主瓣宽度并不比矩形脉冲宽度窄多少。并且，由于有限时域对应无限频域，因此矩形脉冲在经过傅里叶变换后，在频域上变为无限延展的波形，副瓣的存在会降低压缩的意义。 但是，如果我们把矩形脉冲变化为正弦波点频脉冲，则会有很有意思的现象，下图为MATLAB仿真图。 图 SEQ 图 \* ARABIC 2 点频的傅里叶变换 由 REF _Ref412816081 \h 图2我们可以看到，当对一个正弦波点频脉冲进行傅里叶变换后，在频谱上会得到一个sinc函数。而这个sinc函数具有一个令人振奋的特点，即其主瓣宽度与脉冲点频正弦波的脉冲宽度成反比！比例关系为： 当脉冲宽度为5us时，其傅里叶变换的主瓣宽度为0.177MHz；而当脉冲宽度为50us时，其傅里叶变换的主瓣宽度为0.0177MHz。这一特点意味着通过傅里叶变换我们可以将一个很宽的脉冲压缩成一个极窄的脉冲，并且脉冲宽度越宽，则压缩后的窄脉冲宽度越窄。 值得注意的是，上述的变换是按照“时域-频域”这一顺序进行的。但傅里叶变换是一个数学变换，没有任何约束条件要求变换必须按照上述顺序执行。傅里叶变换完全可以在“频域-时域”这样的顺序执行。换言之，不管是在时域还是在频域，只要我们在一个域（domain）中获得了宽脉冲点频信号，那么通过一次傅里叶变换，我们就可以得到一个窄脉冲，从而完成处理工作！ 因此，雷达信号处理工作的目的就变得明晰起来： 我们要得到窄脉冲，以获得高分辨率； 要得到窄脉冲，就必须在某一个域内（可以是时域，频域，波数域，空间域等等）得到点频宽脉冲信号； 点频宽脉冲信号意味着两个方面：矩形包络、线性相位。 下面我们来考察一个线性调频脉冲回波信号： 式  SEQ 式 \* ARABIC 3 为目标延迟，为脉冲周期，为脉冲宽度，为调频斜率。为矩形脉冲框。 对上式进行傅里叶变换，可得：（通过相位驻定POSP原理推导） 式  SEQ 式 \* ARABIC 4 我们构造一个参考函数：，将此函数与相乘，则得： 式  SEQ 式 \* ARABIC 5 仔细观察这个函数，我们发现完全满足脉冲压缩的两个条件： 具有幅度一致矩形脉冲包络，脉冲宽度为； 相位对于变量是线性的。 因此只要再对进行一次FFT，即可以形成一个sinc函数，其主瓣宽度为。另外考察的相位具有一个负号（由一次傅里叶变换带来），因此进行一次逆傅里叶变换，可把负号矫正为正号。 我将上述过程用MATLAB实现一遍如下，结果如下： 图 SEQ 图 \* ARABIC 3 线性调频脉冲信号的脉冲压缩过程 再次总结一下脉冲压缩的处理过程： 将时域信号通过傅里叶变换转换到频域； 在频域中，将信号包络调节为矩形脉冲包络，将信号相位调节为对于频率成线性关系； 将频域信号通过傅里叶变换转换到时域，完成压缩。 这一脉冲压缩处理过程贯穿所有SAR算法，是SAR算法的核心。SAR算法的实质是通过对信号的二维谱进行处理，使得信号在距离向和方位向都满足脉冲压缩的两个条件（矩形包络、线性相位），然后通过逆傅里叶变换，完成距离向和方位向