基于步进频波形的雷达定位算法

此算法针对利用步进频波形实现雷达目标定位，可在此仿真中进行改进，加入时延因子此算法可用于穿墙雷达目标定位及成像。

%计算并画出SFW的合成HRR像
%该函数使用长度为步长两倍的逆傅里叶变化IFFT
function [hl] = hrr_profile（nscatscat_rangescat_rcsndeltafprfvrnotewinid）

%nscat组成目标散射体的个数
%scat_range包含散射体距离的向量
%scat_rcs 包含散射体RCS的向量
%n步进数
%deltaf步进频率
%prf SFW的PRF
%v 目标速度
%rnote 像起始距离
%winid > 0为汉明窗  <0不加窗
%hl 距离像:是指目标散射中心在雷达视线上的投影，反映目标各散射中心间的径向位置关系，体现目标的部分结构特征，是雷达感知、识别目标的重要途径。

c = 3e8; %光速
num_pulses = n; %步进数
SNR_dB = 40;  
nfft = 256; %IFFT点数

carrier_freq = 9.5e9; %载波频率
freq_step = deltaf; %步进频率
V = v; %雷达速度
PRI = 1./prf; %脉冲重复间隔

if （nfft>2*num_pulses）
    num_pulses = nfft/2;
end

%初始化矩阵
Inphase = zeros（（2*num_pulses）1）;
Quadrature = zeros（（2*num_pulses）1）; 
Inphase_tgt = zeros（num_pulses1）;
Quadrature_tgt = zeros（num_pulses1）;
IQ_freq_domain = zeros（（2*num_pulses）1）;
Weighted_I_freq_domain = zeros（（num_pulses）1）;
Weighted_Q_freq_domain = zeros（（num_pulses）1）;
Weighted_IQ_freq_domain = zeros（（2*num_pulses）1）;
Weighted_IQ_time_domain = zeros（（2*num_pulses）1）;
abs_Weighted_IQ_time_domain = zeros（（2*num_pulses）1）;
dB_abs_Weighted_IQ_time_domain = zeros（（2*num_pulses）1）;
taur = 2.*rnote/c; %起始时间延迟

for jscat = 1:nscat
    ii = 0;
    for i = 1:num_pulses
        ii = ii + 1;
        rec_freq = （（i-1）*freq_step）; %fi
        %实数分量    利用欧拉公式
        Inphase_tgt（ii） = Inphase_tgt（ii） + sqrt（scat_rcs（jscat））*cos（-2*pi*rec_freq*（2.*scat_range（jscat）/c-2*（V/c）*（（i-1）*PRI + taur/2 + 2*scat_range（jscat）/c）））;
        %虚部分量
        Quadrature_tgt（ii） = Quadrature_tgt（ii） + sqrt（scat_rcs（jscat））*sin（-2*pi*rec_freq*（2.*scat_range（jscat）/c-2*（V/c）*（（i-1）*PRI + taur/2 + 2*scat_range（jscat）/c）））;
    end
end

if （winid >= 0） %加汉明窗
    window（1:num_pulses） = hamming（num_pulses）;
else
    window（1:num_pulses） = 1;
end

Inphase = Inphase_tgt;
Quadrature = Quadrature_tgt;
Weighted_I_freq_domain（1:num_pulses） = Inphase（1:num_pulses）.*window‘;
Weighted_Q_freq_domain（1:num_pulses） = Quadrature（1:num_pulses）.*window‘;
Weighted_IQ_freq_domain（1:num_pulses） = Weighted_I_freq_domain + Weighted_Q_freq_domain*j;   %虚部和实部的合并
Weighted_IQ_freq_domain（num_pulses:2*num_pulses） = 0.+0.i;
Weighted_IQ_time_domain = （ifft（Weighted_IQ_freq_domain））; %换算时域上
abs_Weighted_IQ_time_domain = abs（Weighted_IQ_time_domain）;
dB_abs_Weighted_IQ_time_domain = 20*log10（abs_Weighted_IQ_time_domain） + SNR_dB;

%计算未定义的范围窗口大小
Ru = c/2/deltaf;
hl = dB_abs_Weighted_IQ_time_domain;
numb = 2*num_pulses;
delx_meter = Ru/numb;
xmeter = 0:delx_meter:Ru-delx_meter;
plot（xmeterdB_abs_Weighted_IQ_time_domain‘k‘）;
xlabel（‘相对距离/m‘）;
ylabel（‘一维距离像‘）;
grid

        

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2831  2018-03-19 15:03  hrr_profile.m

     文件       1389  2018-03-22 19:56  HRR_profile_distance.m

----------- ---------  ---------- -----  ----

                 4220                    2