《合成孔径雷达成像算法与实现》Figure5.17

clc
clear
close all

距离向参数
R_eta_c = 20e3;             % 景中心斜距
Tr = 25e-6;                 % 发射脉冲时宽
Kr = 0.25e12;               % 距离向调频率
Fr = 7.5e6;                 % 距离向采样率
Nrg = 256;                  % 距离线采样点数
Bw = abs(Kr*Tr);            % 距离信号带宽

方位向参数
Delta_f_dop = 80;           % 多普勒带宽
Fa = 104;                   % 方位向采样率
Naz = 256;                  % 距离线数；
theta_sq_c = [0,22.8];      % 波束斜视角
eta_c = [0,-51.7];          % 波束中心偏移时间
f_eta_c = [0,2055];         % 多普勒中心频率
Vr = 150;                   % 等效速度

其他参数
c = 3e8;
f0 = 5.3e9;
lambda = c/f0;
R0 = R_eta_c*cosd(theta_sq_c(1));   % 最短斜距
La = 0.886*(2*Vr*cosd(theta_sq_c(1)))/Delta_f_dop;  % 真实孔径
theta_bw = 0.886*lambda/La;         % 方位向波束宽度
Trr = Nrg/Fr;                       % 发射脉冲时宽
Taa = Naz/Fa;                       % 目标照射时间
Ka = 2*Vr^2/lambda/R0;              % 方位向调频率
d_t_tau = 1/Fr;                     % 距离采样时间间隔
d_t_eta = 1/Fa;                     % 方位采样时间间隔
d_f_tau = Fa/Nrg;                   % 距离采样频率间隔
d_f_eta = Fa/Naz;                   % 方位采样频率间隔

变量设置
tau = (-Trr/2:d_t_tau:Trr/2-d_t_tau)+2*R_eta_c/c;     % 距离时间变量
eta = (-Taa/2:d_t_eta:Taa/2-d_t_eta)+eta_c(1);        % 方位时间变量

坐标设置
[tauX,etaY] = meshgrid(tau,eta);        % 二维网格坐标   

信号设置
R_eta = R0+Vr^2*etaY.^2/2/R0;           % 瞬时斜距
A0 = 1;                                 % 后向散射系数幅度
wr = (abs(tauX-2*R_eta/c) <= Tr/2);     % 距离向包络
wa = sinc(0.886*atan(Vr*(etaY-eta_c(1))/R0)/theta_bw).^2;   % 方位包络
srt = A0*wr.*wa.*exp(-1j*4*pi*R0/lambda)...
               .*exp(-1j*pi*Ka*etaY.^2)...
               .*exp(1j*pi*Kr*(tauX-2*R_eta/c).^2);
srt_z = A0*wr.*wa.*exp(-1j*4*pi*R0/lambda)...
               .*exp(-1j*pi*Ka*etaY.^2)...
               .*exp(1j*pi*Kr*(tauX-2*R_eta/c).^2); % 正扫频
srt_f = A0*wr.*wa.*exp(-1j*4*pi*R0/lambda)...
               .*exp(-1j*pi*Ka*etaY.^2)...
               .*exp(-1j*pi*Kr*(tauX-2*R_eta/c).^2); % 负扫频
Srf_rd = fft(srt);
Srf_2d = fft2(srt);

绘图
figure,subplot(221)
imagesc(abs(Srf_rd))
% shading interp; % 使用插值着色
% colorbar
subplot(222)
imagesc(angle(Srf_rd))
subplot(223)
imagesc(abs(Srf_2d))
subplot(224)
imagesc(angle(Srf_2d))