%%本振信號?
%------------------------------------------------------------------------

T=1e-4;?????????????????????????????????%調頻周期
B=30e6;???????????????????????????????????%調頻帶寬?30MHz?
K=B/T;??????????????????????????????????????%調頻斜率?
Fs=2.5*B;Ts=1/Fs;??????????????????????%采樣頻率?
N=T/Ts;
t0=linspace（-T/2T/2N）;
St0=exp（j*pi*K*t0.^2）;
f0=linspace（-Fs/2Fs/2N）;
%plot（t0*1e4real（St0））;?
%xlabel（‘Time?in?u?sec‘）;?
%title（‘Real?part?of?chirp?signal‘）;?
%grid?on;axis?tight;
%figure（2）;
%plot（f0abs（fftshift（fft（St0length（f0）））））

%-------------------------------------------------------------------------

%%目標設定
%-------------------------------------------------------------------------

c=3e8;????????????%光速
mm=[306090];?????%三個目標的距離
tao=2.*mm./c;??????%時間延遲
rmax=100;
BSt=zeros（1N）;
for?i=1:3
????temp=[];
????k0=ceil（tao（i）/Ts）;
????temp（k0:k0+N-1）=exp（j*pi*K*t0.^2）;
????temp1=temp（1:N）;
????BSt=temp1+BSt;
end
%plot（t0*1e4real（BSt）‘r‘）;

%%加入雜波干擾項
noise=randn（1N）;????%設信噪比0.2
BSt=BSt+noise;
%%混頻濾波獲取差拍信號
%-------------------------------------------------------------------------
u=real（BSt）.*real（St0）;
%f=linspace（-Fs/2Fs/2N）;
plo