步进频距离像拼接仿真

%% 同距离选大距离像拼接
clc;
close all;
clear all

%% 系统参数
tao = 1e-6/37.5;      % 时宽，对应带宽37.5MHz
PRT = 2e-6;         % 脉冲重复周期 2us
detf = 30e6;          % 步进频率间隔 30MHz
c = 3e8;              % 光速 
N = 50;               % 步进频率跳变次数 
B = N * detf;         % 雷达合成带宽 1.5GHz
Rmax = PRT * c / 2;   % 最大不模糊距离

f0 = 9.025e9;         % 载频起始频率9.025GHz
fs = 50e6;            % 回波采样频率50MHz       

nprt = ceil（（PRT-tao）*fs）;         % 1个prt内采样的点数 ceil 向上取整
ntao = ceil（tao*fs）;               % 1个脉宽内采样的点数

detR = c/（2*B）;       % 细化后的距离分辨率
Rfs = c/2/fs;         % 采样距离分辨率
Rtao = c*tao/2;       % 单脉冲分辨距离
Ri = c/2/detf;        % 单点不模糊距离
Ru = c/2/detf/N;
ts_num = Rfs/Ru;      %每个采样点代表的有效数据段，60%的有效，重叠20%，冗余20%
ttao_num = round（Rtao/Ru）;   %每个脉宽代表的有效数据段，80%的有效，冗余20%
Overlapping=ttao_num-ts_num;  % 每个脉冲重叠数据段

snr = 2;      %回波信号叠加信噪比大小
R0 = 30;      %目标所在距离
V0 =0;        %假设为静止目标回波
ntao = 3;     %回波一般都有展宽，从而ntao设为3，而不是2

%% 坐标轴参数设置
range_axis = 0 : Ru : （ts_num*nprt-1）*Ru;
range_axis = range_axis - 3;

%% 生成回波信号
x1 = zeros（Nnprt）;
for n = 0:N-1
    t = n*PRT;
    R（n+1） = R0-V0*t;
    DELAY = 2*R（n+1）/c;
    temp = abs（DELAY*fs）;
    if （ abs（temp-round（temp））        ntarget = round（temp）;
    else
        ntarget = ceil（temp）-1;   % ceil/floor 根据数据定 
    end
    %%%%%%%%%  高斯型的回波 %%%%%%%%%%
    x1（n+1:） = [zeros（1ntarget）exp（-1j*2*pi*（f0+n*detf）*DELAY）...
                 *Func_Gauss（ntao1）zeros（1nprt-ntao-ntarget）]; 
end    %%% 结束N个脉冲，产生回波

% % echo = x1;
echo = awgn（x1snr‘measured‘‘linear‘）;    %回波信号中加入噪声

%% 一维成像
echo_ifft = ifft（echoN）;      %IFFT处理
figure;
mesh（abs（echo_ifft））;
xlabel（‘采样点‘）;ylabel（‘距离分辨单元‘）;
zlabel（‘幅度‘）;
title（‘各频点IFFT相参脉压结果‘）