步进频信号一维像仿真Matlab

%% 步进频率信号一维成像仿真程序，模拟矢网直连线缆成像，发射频率17-17.5GHz，频点数10001
%单点目标
%2018.4.20
clear all
close all

% [hl] = hrr_profile （1     3.03       100    10001 0.05e6[][] []    1）; 
% [hl] = hrr_profile （nscat scat_range scat_rcs n deltaf prf v rnotewinid）

%% Parameter--constant常数定义
c                  =  3.00e8;                           %  speed of light （m/s）电磁波在真空中传播速度
cj                 =  sqrt（-1）;
K                 = 10001;                             %  脉冲数（跳频点数目），频点数
Nfft             = 20*2^nextpow2（K）;        %  逆傅里叶变换点数

%% Parameter--point targets点目标参数定义
nscat           = 1;                                     %  点目标数目矢网为port2端口反射
scat_range  = 3.034;                              %  线长3.03m，等效点目标位置在scat_range/2处
%scat_range2   = 3.034;                        %  线长3.034m
scat_rcs       = 100;                                 % 复散射系数ρ（t）此时仿真幅度无衰减，实际测量时，线有接近5dB的衰减

winid           = 1;                                    % 加窗为1不加窗为-1
%% ------------------第一步：步进频发射信号参数--------------------------

minF           = 17e9;                               % 初始频率17GHz
deltaf          = 0.05e6;                            % △ffi=f0+i△f
f                  = minF+（0:K-1）*deltaf;

% 初始化变量
Inphase               = zeros（Nfft1）;          % 发射I路信号
Quadrature         = zeros（Nfft1）;          % 发射Q路信号
Inphase_tgt         = zeros（K1）;              % 回波I路信号
Quadrature_tgt   = zeros（K1）;              % 回波Q路信号
IQ_freq_domain  = zeros（（Nfft）1）;        % 频域IQ混频后复信号

% Inphase_tgt    = exp（-cj*2*pi*f（:）*scat_range1/c）;
% Quadrature_tgt = exp（-cj*2*pi*f（:）*scat_range2/c）;

%taur = 2. * rnote / c;
% ii = 0;
%    for i = 1:num_pulses
%       ii = ii+1;
%       rec_freq = （（i-1）*deltaf）;%回波频率增量对应式中第二个exp项S（i） = A*exp（-j*2*pi*f0）*exp（-j*2*pi*deltaf*（i-1））;
%       %回波信号的I路以及Q路信号下式为滤掉高频信号的低频信号，矢网为直连计算每个频点时延τ（i） =  scat_range/c
%       Inphase_tgt（ii） = Inphase_tgt（ii） + sqrt（scat_rcs）*cos（2*pi*rec_freq*（scat_range/c））;
%       Quadrature_tgt（ii） = Quadrature_tgt（ii） + sqrt（scat_rcs）*sin（2*pi*rec_freq*（scat_range/c））;
%    end

%% ------------------第三步：回波信号IQ频域加窗混频--------------------------
%% 初始化加窗混频