lose 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