去斜脉冲压缩算法研究

本仿真对脉冲压缩进行改进，提出去斜处理。与脉冲压缩技术不同，此算法可用于带宽极大的雷达定位：穿墙雷达，探地雷达等。定位效果极佳，还可为后续压缩感知及后投影BP算法提供研究思路。

function [y] = stretch（nscattaupf0bscat_rangerrecscat_rcswinid）
eps = 1e-16;
htau = taup/2;
c = 3e8;
trec = 2.*rrec/c;
n = fix（2.*trec*b）;%fix向0取整
m = power_integer_2（n）;
nfft = 2.^m;
%初始化矩阵
x（nscat1:n） = 0;
y（1:n） = 0;
%确定适当的窗口
if （winid == 0）
    win（1:n） = 1;
    win = win‘;    %不加窗
else
    if（winid == 1）
    win = hamming（n）;   %加汉明
    else
        if（winid == 2）
            win = kaiser（npi）;  %参数为pi的kaiser窗
        else
            if（winid == 3）
                win = chebwin（n60）;   %旁瓣为-60db的切比雪夫窗
            end
        end
    end
end

deltar = c/2/b;
max_rrec = deltar*nfft/2;
maxr = max（scat_range）;

if（rrec>max_rrec|maxr>=rrec）
    ‘Error.Receive window is too large;or scatterers fall outside window‘;
    return
end

t = linspace（0taupn）;
for j = 1:1:nscat
    range = scat_range（j）;
    psi1 = 4.*pi*range*f0/c-4.*pi*b*range*range/c/c/taup;
    psi2 = （2*4.*pi*b*range/c/taup）.*t;
    x（j:） = scat_rcs（j）.*exp（i*psi1+i.*psi2）;
    y = y + x（j:）;
end

figure
plot（treal（y）‘k‘）;
xlabel（‘时间延迟‘）;
ylabel（‘未压缩的回波信号‘）;

ywin = y.*win‘;
yfft = fft（yn）./n;
out = fftshift（abs（yfft））;
figure
delinc = rrec/n;
dist = linspace（（-rrec/2）rrec/2n）;
plot（distout‘k‘）;
xlabel（‘相对距离/m‘）;
ylabel（‘压缩回波信号‘）;

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

     文件       1355  2018-03-21 21:20  stretch.m

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

                 1355                    1