2017.3.31
%%%%%%%??利用頻域處理方法進行脈沖壓縮??%%%%%%%
clear?all
clc
clf
eps?=?1e-10;
B=100e-6;???????%信號帶寬
Fm=1e6;???????????%調頻
k=Fm/B;??????????%調頻斜率
Ts=1/（5*Fm）;????????%采樣周期
Ns=fix（B/Ts）;????????%采樣點數
Nf=1024;???????????????%?fft點數
t=0:Ts:B-Ts;?
y=exp（j*pi*k*t.^2）;???%脈沖壓縮前的線形調頻信號
yfft?=?fft（yNf）?;
h=zeros（1Ns）;
%%?%%%%%%%%%%%%%%%%%%%%%%%%%Hamming窗%%%%%%%%%%%%%%%%%%%%%%%%%%%?
for?i=1:Ns
????h（i）=conj（y（Ns-i+1））;
end
hfft=?fft（hNf）;?????%?匹配濾波器的頻域響應
lfm?=abs（ifft（yfft?.*hfft））;?%脈沖壓縮????
maxval?=?max?（lfm）;
lfm?=?eps?+?lfm?./?maxval;????%?利用最大值歸一化
lfm_db=20*log10（lfm）;???%取對數
%%%%%%%%%%%%%%?加窗處理?%%%%%%%
win?=?hamming（Ns）‘;
h_w=h.*win;???????%?加窗
hfft_w=fft（h_wNf）;?????%?加窗的匹配濾波器的頻域響應
lfm_w?=?abs（ifft（yfft?.*hfft_w））;?%脈沖壓縮?
maxval1?=?max（lfm_w）;
val=lfm_w?;
lfm_w?=?eps?+?lfm_w?./?maxval;????%?利用lfm的最大值歸一化
lfm_w1?=?eps?+?val./?maxval1;????%?利用lfm_w的最大值歸一化
lfm_w_db=20*log10（lfm_w）;???%取對數
lfm_w1_db=20*log10（lfm_w1）;???%取對數
%%%%%%%%%%%%%%%%
tt?=0:Ts:2*B-Ts;
figure（1）
plot?（ttlfm_db（1:2*Ns）‘b‘）
axis（[.2*B?1.8*B?-60?0]?）
xlabel?（‘t?-?seconds?‘）;
ylabel（‘?db‘）
title（‘沒有加Hamming窗的脈沖壓縮輸出‘）
grid?on
figure（2）
plot?（ttlfm_w1_db（1:2*Ns）‘r‘）
axis（[.2*B?1.8*B?-60?0]?）
xlabel?（‘t?-?seconds?‘）;
ylabel（‘?db‘）
title（‘加Hamming窗的脈沖壓縮輸出‘）
grid?on
figure（3）
plot?（ttlfm_db（1:2*Ns）‘b‘ttlfm_w_db（1:2*Ns）‘r‘）
axis（[.7*B??1.3*B?-60?0]?）
xlabel?（‘t?-?seconds?‘）;
ylabel（‘?db‘）
legend（‘未加Hamming窗‘‘加Hamming窗‘）;
title（‘脈沖壓縮輸出對比‘）
grid?on

%%?%%%%%%%%%%%%%%%%%%%%%%%%%Hanning窗%%%%%%%%%%%%%%%%%%%%%%%%%%%?
for?i=1:Ns
????h（i）=conj（y（Ns-i+1））;
end
hfft=?fft（hNf）;?????%?匹配濾波器的頻域響應
lfm?=abs（ifft（yfft?.*hfft））;?%脈沖壓縮????
maxval?=?max?（lfm）;
lfm?=?eps?+?lfm?./?maxval;????%?利用最大值歸一化
lfm_db=20*log10（lfm）;???%取對數
%%%%%%%%%%%%%%?加窗處理?%%%%%%%
win?=?hanning（Ns）‘;
h_w=h.*win;???????%?加窗
hfft_w=fft（h_wNf）;?????%?加窗的匹配濾波器的頻域響應
lfm_w?=?abs（ifft（yfft?.*hfft_w）