雷达SAR成像仿真的应用（Matlab代码实现）

%% 
%    线性调频信号的脉冲压缩
% 
%    介绍：对线性调频信号进行仿真，输出其时频域的相关信息，并模拟回波信号，
%    对其进行脉冲压缩和加窗处理。
% 
%    实验记录：
%     1.线性调频信号时域包络、相位；实部、虚部
%     2.线性调频信号频谱幅频、相频特性；实部、虚部
%     3.两个目标回波的时域和频域波形
%     4.信号通过匹配滤波器的输出结果（脉冲压缩）。
%     5.用Hamming窗抑制脉冲压缩结果副瓣
%%  基本参数 
    clc;clear all;close all;

   T = 10e-6; % LFM周期/脉宽 10us
    B = 60e6; % LFM带宽 60Mhz
    fs = 100e6; % 采样率 100MHz
    K = B/T;
%%  模拟发射信号
    n = round(15*T*fs);
    t = linspace(-10*T, 10*T,n);
    lfmT = rectpuls(t,T).*exp(1j*pi*K*t.^2);
    lfmF = fftshift(fft(fftshift(lfmT)));
    f = linspace(-fs,fs,n);
    %% 时域绘图
        figure();
        plot(diff(phase(lfmT)));
        title('LFM信号的时间-频率变化趋势图');
        xlabel('时间');
        ylabel('频率');
        xlim([7200,7800])
    % 包络
        figure();
        subplot(2,2,1);
        plot(t,abs(lfmT));
        title('LFM信号时域包络');
        xlabel('t/s');
        ylabel('幅度');
        xlim([-1e-5,1e-5])
        ylim([-0.5,1.5])
    % 相位
        subplot(2,2,2);
        plot(t,phase(lfmT));
        title('LFM信号时域相位');
        xlabel('t/s');
        ylabel('相位');
        xlim([-5e-6,5e-6])
    % 实部
        subplot(2,2,3);
        plot(t,real(lfmT));
        title('LFM信号时域实部');
        xlabel('t/s');
        ylabel('幅度');
        xlim([-1.5e-6,1.5e-6]);
        ylim([-1,1]);
    % 虚部
        subplot(2,2,4);
        plot(t,imag(lfmT));
        title('LFM信号时域虚部');
        xlabel('t/s');
        ylabel('幅度');
        xlim([-1.5e-6,1.5e-6]);
        ylim([-1,1]);
    %% 频域绘图
        figure();
        subplot(2,2,1);
        plot(f,abs(lfmF));
        title('LFM信号幅频特性');
        xlabel('Hz');
        ylabel('幅度');
        subplot(2,2,2);
        plot(unwrap(angle(lfmF)));
        title('LFM信号相频特性');
        xlabel('Hz');
        ylabel('相位');
        subplot(2,2,3);
        plot(f,real(lfmF));
        title('LFM信号频谱实部');
        xlabel('Hz');
        ylabel('幅度');    
        xlim([-3e7,3e7]);
        subplot(2,2,4);
        plot(f,imag(lfmF));
        title('LFM信号频谱虚部');
        xlabel('Hz');
        ylabel('幅度');    
        xlim([-3e7,3e7]);