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

clc
clear
close all%% 参数设置
Ta = 64;                                                                    % 脉冲持续时间
Ka = -1.56e-2;                                                              % 方位向调频率
Delta_f_dop = abs(Ta*Ka);                                                   % 方位向带宽
alpha_OS = 1;                                                               % 过采样率
Fa = alpha_OS*Delta_f_dop;                                                  % 采样率
N = 2*ceil(Fa*Ta/2);                                                        % 方位向采样点数
dt = Ta/N;                                                                  % 采样时间间隔
df = Fa/N;                                                                  % 采样频率间隔%% 变量设置
t1 = -Ta/2:dt:Ta/2-dt;                                                      % 时间轴
f1 = -Fa/2:df:Fa/2-df;                                                      % 频率轴
t2 = -3/4*Ta:dt:1/4*Ta-dt;                                              % 时间轴
f2 = -3/4*Fa:df:1/4*Fa-df;                                              % 频率轴
n = 0:N-1;                                                                  % 采样轴%% 信号表达
st1 = exp(1j*pi*Ka*t1.^2);                                                  % chirp信号
st2 = exp(1j*pi*Ka*t2.^2);                                                  % chirp信号
omega_a = sinc(1.5*t1/Ta).^2;                                               % 双程天线波束方向图
omega_a_nor = abs(omega_a)/max(abs(omega_a));                               % 归一化
omega_a_nor_log = db(omega_a_nor);                                          % 对数化
st1 = omega_a.*st1;                                                         % 天线方向图调制的chirp信号
st2 = omega_a.*st2;                                                         % 天线方向图调制的chirp信号% 频域
Sf1 = fft(st1);                                                             % 频域表达式
Sf2 = fft(st2);                                                             % 频域表达式%% 绘图
subplot(221),plot(n,real(st1)),xlabel('(a)信号实部（斜视角为0）')
subplot(222),plot(f1,abs(Sf1)),xlabel('(b)信号频谱（斜视角为0）')
subplot(223),plot(n,real(st2)),xlabel('(a)信号实部（斜视角非0）')
subplot(224),plot(f2,abs(Sf2)),xlabel('(d)信号频谱（斜视角非0）')

注意：斜视即等价于时间轴平移，即现在时刻可得到之前未来时刻的结果