clc;clear?all;close?all;
f0=1.501e9;
f1=1.502e9;
T=0.01;
C=3e8;?
REP=5;
for（rep=1:REP）
for（m=1:600）
???R=m;
???V=300;
???fd0=2*V*f0/C;
???fd1=2*V*f1/C;
???fs=10e5;Ts=1/fs;
???N=round（T/Ts）;
???t=linspace（-T/2T/2N）;
???f=linspace（-fs/2fs/2N）;

???echo=exp（j*2*pi*fd0*t-j*2*pi*f0*2*R/C）;
???echo=awgn（echo0‘measured‘）;
???echo_FFT=fftshift（fft（echo））;?
???position0=find（abs（echo_FFT）==max（abs（echo_FFT）））;
???A0（m）=angle（echo_FFT（position0））;

???echo1=exp（j*2*pi*fd1*t-j*2*pi*f1*2*R/C）;
???echo1=awgn（echo10‘measu