资源描述：
《时域脉冲压缩和频域脉冲压缩比较附matlab源码》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、首先在脉冲压缩之前要搞清楚以下两点，而不是上来就自己想着设计一个匹配滤波器进行脉冲压缩，这样往往会导致错误的结果!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=[112
2、];taup=0.005e-3;f0=5.6e9;%carrierfreq,由于exp(-i*2*pi*f0*2*range/c)与t无关，此项是个固定的相移，不起作用，建模中可以忽略c=3.e8;%speedoflightfs=2*b;%samplingfreqsampling_interval=1/fs; n=fix(taup/sampling_interval);freqlimit=0.5*fs;freq=linspace(-freqlimit,freqlimit,n);%noticefreq_sampli
3、ng_interval=fs/n=1/taup;t=linspace(-taup/2,taup/2,n);%initializeinput,outputandreplicavectorsx(nscat,1:n)=0.;%xisanscat-by-nmatrixy(1:n)=0.;replica(1:n)=0.;%basebandlfmsignalreplica=exp(i*pi*(b/taup).*t.^2);%我以前不知道原来接收的一个脉冲中包含了多个目标的回波forj=1:1:nscat range=sc
4、at_range(j);%对于函数Y来说，平移都是加．．%其实，平移都是针对函数Y来说的，而不是x．．．%因为x是y的自变量,分析的时候将对函数Y的加，等效的变成了自变量的减！%MATLAB针对自变量X的移动(t-tau)必须变成针对函数Y本身的移动(t+tau)%t+(2*range/c)meanstargetlocateattau=2*range/c% t+(2*range/c)MUSTgreaterthant,不是课本上自变量的移动t-(2*range/c)！%f(t-t0)是自变量t的范围不动,比如观测范
5、围永远是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,:)=
6、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; endout=xcorr(replica,y); %注意：xcorr(replica,y)中replica和y的次序不能颠倒！out=out./n;%除以n是有偏自相关估计的要求s=taup*c/2;%最大不模糊距离Npoints=ceil(rrec*n/s);% 最大不模糊距离s上采了n个点，问在接收窗rrec大小的距离上
7、采了多少点？是个比例式%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)这部分结果，也是个最大不模糊距离，但是此部分由于对应负的延迟没有实际意义。而脉冲宽度为Te
8、的LFM真正有用的是(0~Te/2)这部分结果，该部分进行xcorr或者卷积的匹配滤波后点数从n/2变为n。阿弥陀佛dist=linspace(0,rrec,Npoints);figure(3)plot(dist,abs(out(n:n+Npoints-1)),'k');%用xcorr匹配滤波有n-1点的延迟，看看1:n-1点被跳过xlabel('Targetrelati
