《合成孔径雷达成像算法与实现》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）')注意：斜视即等价于时间轴平移，即现在时刻可得到之前未来时刻的结果