《合成孔径雷达成像——算法与实现》之【10】仿真图3.13

频域匹配滤波有三种方式：

• 方式1的匹配滤波器：时间反褶后取复共轭，计算Nfft点补零DFT
  H1(f) = fft(conj(s(-t)),Nfft)

• 方式2的匹配滤波器：补零后计算Nfft点DFT，对结果取复共轭
  H2(f) = conj(fft(s(t),Nfft))

• 方式3的匹配滤波器：直接在频域生成匹配滤波器
  H3(f) = exp(1i*pi*f.^2/K)列表内容

滤波器的形式不对，滤波后脉冲峰值的位置也不同，方式1的脉冲在信号尾部，方式2的脉冲在信号前部，方式3的脉冲与信号零频位置对准。

% SAR_Figure_3_13% 2016/08/31close all;clear all;clc;T = 10e-6;          % 脉冲持续时间B = 15e6;           % 脉冲带宽K = B/T;            % 调频率ratio = 5;          % 过采样率Fs = ratio*B;       % 采样频率dt = 1/Fs;          % 采样间隔Nr = ceil(T/dt);    % 采样点数t0 = ((0:Nr-1)-Nr/2)/Nr*T;          % 基本时间轴st0 = exp(1i*pi*K*(t0-T/5).^2);     % 基本信号space1 = zeros(1,round(Nr/5));      % 生成空信号space2 = zeros(1,Nr);               % 生成空信号st = [space1,st0,space2,st0,space2,st0,space1];     % 实际信号N = length(st);             % 实际信号长度n = 0:N-1;                  % 样本轴f = ((0:N-1)-N/2)/N*Fs;     % 基本频率轴Sf = fftshift(fft(st));                         % 实际信号的傅里叶变换Hf1 = fftshift(fft(conj(fliplr(st0)),N));       % 方式1的匹配滤波器：时间反褶后取复共轭，计算N点补零DFTHf2 = fftshift(conj(fft(st0,N)));               % 方式2的匹配滤波器：补零后计算DFT，对结果取复共轭Hf3 = exp(1i*pi*f.^2/K);                        % 方式3的匹配滤波器：直接在频域生成匹配滤波器out1 = ifft(ifftshift(Sf.*Hf1));out2 = ifft(ifftshift(Sf.*Hf2));out3 = ifft(ifftshift(Sf.*Hf3));figure,set(gcf,'Color','w');subplot(4,1,1),plot(n,real(st));axis tighttitle('(a)输入阵列信号的实部');ylabel('幅度');subplot(4,1,2),plot(n,abs(out1));axis tighttitle('(b)方式1的匹配滤波输出');ylabel('幅度');subplot(4,1,3),plot(n,abs(out2));axis tighttitle('(c)方式2的匹配滤波输出');ylabel('幅度');subplot(4,1,4),plot(n,abs(out3));axis tighttitle('(d)方式3的匹配滤波输出');xlabel('时间(采样点)');ylabel('幅度');