FMCW SAR 距离-多普勒算法的matlab程序

以下是FMCW SAR距离-多普勒算法的MATLAB程序：

%% FMCW SAR距离-多普勒算法
% 假设目标是静止的
% 参数设置
c = 3e8; % 光速
fc = 10e9; % 中心频率
B = 100e6; % 带宽
T = 1/B; % 脉冲宽度
Rmax = 250; % 最大距离
Vmax = 50; % 最大速度
Np = 256; % 脉冲个数
Nr = 512; % 接收点数
dR = Rmax/Nr; % 距离分辨率
dV = Vmax/(Np/2); % 速度分辨率
Fs = 2*B; % 采样率
Ts = 1/Fs; % 采样时间间隔
t = linspace(0, Np*T, Np*Fs); % 时间序列
f = linspace(-Fs/2, Fs/2, Nr); % 频率序列
v = linspace(-Vmax, Vmax, Np/2); % 速度序列

% 产生线性调频信号
s = exp(1j*pi*B*t.^2);

% 产生目标场景
Rt = 100; % 目标距离
Vt = 0; % 目标速度
Dr = 10; % 目标距离差
Dv = 5; % 目标速度差
target = exp(1j*2*pi*fc*(t - 2*Rt/c)) .* exp(1j*pi*B*(t - 2*Rt/c - Dr/c).^2) .* ...
    exp(1j*2*pi*B*Vt*(t - 2*Rt/c - Dr/c)) .* ...
    exp(1j*2*pi*fc*(t - 2*(Rt+Dr)/c)) .* exp(1j*pi*B*(t - 2*(Rt+Dr)/c - Dr/c).^2) .* ...
    exp(1j*2*pi*B*(Vt+Dv)*(t - 2*(Rt+Dr)/c - Dr/c)) .* ...
    exp(1j*2*pi*fc*(t - 2*(Rt+2*Dr)/c)) .* exp(1j*pi*B*(t - 2*(Rt+2*Dr)/c - Dr/c).^2) .* ...
    exp(1j*2*pi*B*(Vt+2*Dv)*(t - 2*(Rt+2*Dr)/c - Dr/c));

% 产生接收信号
Rx = zeros(Np, Nr);
for i = 1:Np
    for j = 1:Nr
        tau = 2*(Rt + (j-1)*dR)/c;
        fd = 2*B*Vt/c;
        Rx(i,j) = target(i) * exp(-1j*2*pi*fc*tau) * exp(-1j*pi*B*(tau-T).^2) * ...
            exp(-1j*2*pi*fd*(t(i)-tau)) * exp(-1j*pi*B*(t(i)-tau-T).^2);
    end
end

% 进行STFT处理
S = zeros(Np, Np/2);
for i = 1:Np
    x = squeeze(Rx(i,:));
    X = fftshift(fft(x, Nr));
    S(i,:) = abs(X(Nr/2+1:end)).^2;
end

% 显示距离-多普勒谱
figure;
imagesc(v, f, S);
axis xy;
xlabel('速度 (m/s)');
ylabel('频率 (Hz)');
title('距离-多普勒谱');
colorbar;