SAR方位向成像matlab代码及注释

SAR方位向成像matlab代码及注释及详细注释 对应有博客的讲解

        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        %        CROSS-RANGE IMAGING         %
        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


colormap（gray（256））
cj=sqrt（-1）;
pi2=2*pi;
%
c=3e8;              % propagation speed
fc=200e6;           % frequency
lambda=c/fc;        % Wavelength
k=pi2/lambda;       % Wavenumber
Xc=1.e3;            % Range distance to center of target area
%
% Case 1:
L=400;              % synthetic aperture is 2*L
Y0=100;             % target area in cross-range is within [Yc-Y0Yc+Y0]
Yc=0;               % Cross-range distance to center of target area
%
% Case 2: 这时候需要进行zero-padding
% L=150;
% Y0=200;
% Yc=.5e3;
%
theta_c=atan（Yc/Xc）;  % squint angle to center of target area
Rc=sqrt（Xc^2+Yc^2）;   % squint range to center of target area
kus=2*k*sin（theta_c）; % Doppler frequency shift in ku domain due to squint
                      % 正侧视时候的正对着的目标斜视角为零 多普勒偏移也就为零
                      % 但正侧视情况下每一个目标都有被正视的情况 这个事情暂时是没有考虑的
                      % 毕竟这里只研究了方位向成像 且只是研究了一个合成孔径内部发生的事情
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   Program performs slow-time compression to save PRF   %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
Xcc=Xc/（cos（theta_c）^2）;       % redefine Xc by Xcc for squint processing
%
du=（Xcc*lambda）/（4*（Y0+L））;    % sample spacing in aperture domain 回波信号的理论的采样间隔
duc=（Xcc*lambda）/（4*Y0）;       % sample spacing in aperture domain 回波信号的压缩形式的采样间隔
                               % for compressed signal 这个和书上的公式不太一样？
DY=（Xcc*lambda）/（4*L）;         % Cross-range resolution
%
L_min=max（Y0L）;               % Zero-padded aperture
%这样为什么相当于补零呢？如果Y0较大，那么后面计算压缩信号的采样个数的时候利用较大的Y0
%　而不是Ｌ，但是信号其实在Ｌ长度之外是没有值的，也就是为零

% u domain parameters and arrays for compressed signal
%
mc=2*ceil（L_min/duc）;             % number of samples on aperture 
                                  % 用duc去采样回波信号 会aliased
                                  % 所以才要进行压缩处理
uc=duc*（-mc/2:mc/2-1）;            % synthetic aperture array 压缩信号的下标索引
dkuc=pi2/（mc*duc）;                % sample spacing in ku domain
kuc=dkuc*（-mc/2:mc/2-1）;          % kuc array 压缩信号的频域序列索引
%
dku=dkuc;                         % sample spacing in ku domain
%不确定这么理解对不对：频域的补零对应时域的内插，所以时域的信号长度是不变的，那么频域的间隔也就不变
%
% u domain parameters and arrays for Synthetic aperture signal  下面的参数对应经压缩-升采样-解压缩后的回波信号
%
m=2*ceil（pi/（du*dku））;            % number of samples on aperture pi/du 对应频域的范围 范围除以间隔得到采样个数 
%这里的意思就是说回波信号（未压缩）的理论采样间隔决定了其频域的范围 这个是确定的
%然后，回波信号按照duc采样后混叠了，解决办法是压缩 上采样 再解压缩

du=pi2/（m*dku）;%频域的而采样间隔乘以采样的个数是频域的范围 他的倒数对应时域的采样间隔

u=du*（-m/2:m/2-1）;                % synthetic aperture array 信号的下标索引
                                  % 这个是回波信号（经压缩-升采样-解压缩后的）的索引 注意不是直接采样得到的
ku=dku*（-m/2:m/2-1）;              % ku array 
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%          SIMULATION      目标的基本信息 %%%%%%%          
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
ntarget=4;             % Number of targets
%
% Targets‘ coordinates and reflectivity

yn（1）=0;             fn（1）=1;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       9130  2016-05-07 12:47  crange.m

----------- ---------  ---------- -----  ----

                 9130                    1