《合成孔径雷达成像算法与实现》Figure3.12

clc
clear
close all

% 参数设置
TBP = 724;                              % 时间带宽积
T = 42e-6;                              % 脉冲持续时间
t_0 = 1e-6;                             % 脉冲回波时延
Nfft = 2^11;                            % fft长度

% 参数计算
B = TBP/T;                              % 信号带宽
K = B/T;                                % 线性调频频率
alpha_os = 1.07;                        % 过采样率
F = alpha_os*B;                         % 采样频率
N = 2*ceil(F*T/2);                      % 采样点数
dt = T/N;                               % 采样时间间隔
df = F/N;                               % 采样频率间隔

% 变量设置
t = -T/2:dt:T/2-dt;                     % 时间变量
f = -F/2:df:F/2-df;                     % 频率变量

% 信号表达
st = exp(1j*pi*K*t.^2);
Sf = fft(st);

% 窗函数
window = kaiser(N,2.5)';                %时域窗函数
Window = fftshift(window);              %频域窗函数

% 信号变换
Hf_3 = Window.*exp(1j*pi*f.^2/K);
Hf_phi = pi*f.^2/K;
Hf_fre = 2*pi*f/K;

%绘图
subplot(311),plot(abs(Hf_3)),axis([-15 791,0 1.2])
subplot(312),plot(1:388,Hf_fre(388:775)*1e+4),hold on
plot(389:776,Hf_fre(1:388)*1e+4)
subplot(313),plot(fftshift(Hf_phi))