close?all;?clear?all;?clc;?

T?=?0.26;
fs?=?20e3;
fc?=?5e3;
B?=?fc/3.2;
t?=?0:1/fs:T-1/fs;

type_index?=?3;?%?Choosing?signal:?1?for?pcw?2?for?lfm?3?for?hfm?4?for?hfm+pcw
type_name?=?{‘PCW‘‘LFM‘‘HFM‘‘HFM+PCW‘};

switch?type_index
????case?1
????%?1.?pcw
????x?=?sin（2*pi*fc*t）;
????
????case?2
????%?2.?lfm
????k?=?B/T;
????f0?=?fc-B/2;
????x?=?sin（2*pi*（f0*t+k/2*t.^2））;
????
????case?3
????%?3.?hfm
????f0?=?fc+B/2;
????beta?=?B/f0/（fc-B/2）/T;
%?????f0=f0*（fc-B/2）/fc;
????x?=?sin（2*pi/beta*log（1+beta*fc*t））;
????
????case?4
????%?4.?hfm?+?pcw
????f0?=?fc+B/2;
????beta?=?B/f0/（fc-B/2）/T;
????x?=?sin（2*pi/beta*log（1+beta*f0*t））?+?sin（2*pi*fc*t）;
end

re_fs?=?0.9*fs:10:1.1*fs;
alpha?=?re_fs/fs;???????%?Doppler?ratio?alpha?=?1-2*v/c
doppler?=?（1-alpha）*fc;	%?Doppler?=?2v/c*fc?=?（1-alpha）*fc

N_a?=?length（?resample（xfsmin（re_fs））?）;
N?=?N_a?+?length（x）-1;
afmag?=?zeros（length（alpha）N）;

hh?=?waitbar（0?‘Calculating?the?Wide-band?Ambiguity?Function?Waiting?please!?‘）;
for?i?=?1:length（alpha）
????x_alpha?=?resample（xfsre_fs（i））;
????x_temp?=?zeros（1length（x_alpha））;
????for?j?=?1:lengt