时域脉冲压缩与频域脉冲压缩比较附matlab源码.docVIP

首先在脉冲压缩之前要搞清楚以下两点，而不是上来就自己想着设计一个匹配滤波器进行脉冲压缩，这样往往会导致错误的结果! 1. LFM波形是up-chirp的，还是down-chirp的？ 2. LFM的带宽B，脉冲持续时间Te（不是Tr噢），mu = B/Te和TS(时域采样间隔) 其次，要清楚频域脉冲压缩的实质是快速计算时域脉冲压缩，这两个是一回事！ 1 时域脉冲压缩，采用xcorr(等长)，或conv(不等长) nscat = 3;rrec = 200;b = 100e6;scat_range = [10, 30, 100];scat_rcs = [1 1 2];taup = 0.005e-3; f0 = 5.6e9; % carrier freq ,由于exp(-i*2*pi*f0*2*range/c)与t无关，此项是个固定的相移，不起作用，建模中可以忽略c = 3.e8; % speed of lightfs = 2*b; % sampling freqsampling_interval = 1/fs;?n = fix(taup/sampling_interval);freqlimit = 0.5*fs;freq = linspace(-freqlimit,freqlimit,n); % notice freq_sampling_interval = fs/n = 1/taup;t = linspace(-taup/2,taup/2,n);% initialize input, output and replica vectorsx(nscat,1:n) = 0.; % x is a nscat-by-n matrixy(1:n) = 0.;replica(1:n) = 0.;% baseband lfm signalreplica = exp(i * pi * (b/taup) .* t.^2); % 我以前不知道原来接收的一个脉冲中包含了多个目标的回波for j = 1:1:nscat??? range = scat_range(j) ;%对于函数Y来说，平移都是加．．%其实，平移都是针对函数Y来说的，而不是x．．．%因为x是y的自变量,分析的时候将对函数Y的加，等效的变成了自变量的减！%MATLAB针对自变量X的移动 (t - tau) 必须变成针对函数Y本身的移动 (t + tau)% t +(2*range/c) means target locate at tau = 2*range/c%? t +(2*range/c) MUST greater than t, 不是课本上自变量的移动 t - (2*range/c) ！% f(t-t0)是自变量t的范围不动,比如观测范围永远是t=0:100ms，用于不会有t0=200ms的图像出现在f(t-t0)的图像中% 所以f(t-t0)是通过调整自变量实现曲线的移动的，而且自变量t的范围固定，这样移动曲线y不具有物理意义。% 要采用直接移动y的方式移动曲线，要t的范围不固定，比如将t=0:100ms的回波波形，移动到t=100:200，移动后要画图的话，t的范围要变大。% 也就是说t凡是书上针对自变量 y = f(x - x0)移动曲线的写法在Matlab中建模雷达回波时要写成 y = f(x + x0)，才有物理意义，注意x的范围已经变大 ???? x(j,:) = scat_rcs(j) .*exp(-i*2*pi*f0*2*range/c).* exp(i * pi * (b/taup) .* (t +(2*range/c)).^2) ;???? y = x(j,:)? + y;?end out =xcorr(replica, y);?% 注意：xcorr(replica, y)中 replica和y的次序不能颠倒！out = out ./ n; % 除以n是有偏自相关估计的要求s = taup * c /2; %最大不模糊距离Npoints = ceil(rrec * n /s); %?最大不模糊距离s上采了n个点，问在接收窗rrec大小的距离上采了多少点？是个比例式 % n是LFM信号taup上的采样点数，按说 Ru = c*taup/2 = s上面应该对应 n/2点才对啊，为什么脉冲压缩后这里s对应了n点呢？ % 后面做了个实验（见最后）验证了下，人家这里是对的我自己的理解是错的。先死记住这个结论：xcorr脉冲压缩后的n点对应最大不模糊距离(0~Te/2)！ %?原因是经过脉冲压缩后原来长度为n的序列变成了长度为2n-1的序列，前边的n-1点加上0点这n个点对应着(-Te/2~0)这部