加窗都脉冲压缩的影响0.5 1 1.5 2 2.5 3 3.5x 10-4-60-50-40-30-20-100t - seconds 冲冲db冲 冲 冲 冲 冲 冲 冲 冲冲 冲 冲冲 hamming冲0.5 1 1.5 2 2.5 3 3.5x 10-4-60-50-40-30-20-100t - seconds 冲冲db冲 冲 冲 冲 冲 冲 冲 冲冲 冲 冲冲 hanning冲0.5 1 1.5 2 2.5 3 3.5x 10-4-60-50-40-30-20-100t - seconds 冲冲db冲 冲 冲 冲 冲 冲 冲 冲冲 冲 冲冲 kaiser冲0.5 1 1.5 2 2.5 3 3.5x 10-4-70-60-50-40-30-20-100t - seconds 冲冲db冲 冲 冲 冲 冲 冲 冲 冲冲 冲 冲冲 blackman冲加各种窗函数脉压的结果对比主瓣宽度 峰值 主副瓣比 主瓣损失不加窗 1.8us 22.77dB 10.83dB 0.00dB加 hamming 窗 2.0us 20.26dB 18.96dB 2.51dB加 hanning 窗 2.2us 19.99dB 20.67dB 2.78dB加 kaiser 窗 1.9us 21.78dB 15.66dB 0.99dB加 blackman 窗 3.1us 19.20dB 20.16dB 3.67dB可见加窗脉压对副瓣抑制比较好。
程序如下：close all; clear all;eps = 1e-10;%=====================================================================% 雷达参数 %%====================================================================C=3.0e8; %光速(m/s)Lambda=0.03;%雷达工作波长f0=C/Lambda;BandWidth=1.0e6; %发射信号带宽TimeWidth=200e-6; %发射信号时宽mu=BandWidth/TimeWidth %调频率Fs=2*BandWidth; %采样频率Ts=1/Fs;Ns=fix(Fs*TimeWidth);%计算一个脉冲周期的采样点数；N=1024; %FFT 点数t=0:Ts:TimeWidth-Ts;%====================================================================y=exp(j*pi*mu*t.^2);yfft = fft(y,N) ;h=zeros(1,Ns);for i=1:Nsh(i)=conj(y(Ns-i+1));endhfft= fft(h,N); % 匹配滤波器的频域响应ycomp =abs(ifft(yfft .*hfft)); %脉冲压缩 maxval = max (ycomp);ycomp = eps + ycomp ./ maxval; % 利用最大值归一化ycomp_db=20*log10(ycomp); %取对数%%%%%%%%%%%%%% 加窗处理 %%%%%%%win = hamming(Ns)';h_w=h.*win; % 加窗hfft_w=fft(h_w,N); % 加窗的匹配滤波器的频域响应ycomp_w = abs(ifft(yfft .*hfft_w)); %脉冲压缩 maxval1 = max(ycomp_w);val=ycomp_w ;ycomp_w = eps + ycomp_w ./ maxval; % 利用 ycomp 的最大值归一化ycomp_w1 = eps + val./ maxval1; % 利用 ycomp_w 的最大值归一化ycomp_w_db=20*log10(ycomp_w); %取对数ycomp_w1_db=20*log10(ycomp_w1); %取对数%%%%%%%%%%%%%%%%tt =0:Ts:2*TimeWidth-Ts;figure(1)plot (tt,ycomp_db(1:2*Ns),'b')axis([.2*TimeWidth 1.8*TimeWidth -60 0] )xlabel ('t - seconds ');ylabel('幅度 db')title('未加窗的脉冲压缩输出')grid onfigure(2)plot (tt,ycomp_w1_db(1:2*Ns),'r')axis([.2*TimeWidth 1.8*TimeWidth -60 0] )xlabel ('t - seconds ');ylabel(' 幅度 db')title('加窗的脉冲压缩输出')grid onfigure(3)plot (tt,ycomp_db(1:2*Ns),'b',tt,ycomp_w_db(1:2*Ns),'r')axis([.2*TimeWidth 1.8*TimeWidth -60 0] )xlabel ('t - seconds ');ylabel(' 幅度 db')legend('未加窗 ','加窗');title('脉冲压缩输出对比')grid on。