加窗都脉冲压缩及影响.docVIP

加窗都脉冲压缩的影响 加各种窗函数脉压的结果对比 主瓣宽度 峰值 主副瓣比 主瓣损失 不加窗 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:Ns h(i)=conj(y(Ns-i+1)); end hfft= 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 on figure(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 on figure(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