%產生CW脈沖信號模型
clc
clf
f=4*10^6;??????%信號頻率
T=1/f;?????????%信號周期
Fs=10*f;???????%10倍采樣
Ts=1/Fs;
P=30;
time=（0:P*T*Fs）/Fs;
N=length（time）;
var=sqrt（1/exp（1.0））;
noise=var*rand（1N）;
signal=2*sin（2*pi*f*time）;
%signal=2*signal.*hanning（N）‘;
subplot（331）
plot（timesignal）;grid
axis（[0?10*10^-6?-3?3]）

%理想時系統參數
v=0.003;
c=1500;
d=0.15;
tans_time1=d/（c+v）??%順流時間延遲
tans_time2=d/（c-v）??%逆流時間延遲

%順流接收信號

signal_1=2*sin（2*pi*f*（time-tans_time1））+noise;
[ba]=fir1（32[0.19?0.21]）;
signal_1=filter（basignal_1）;
subplot（332）
plot（timesignal_1）;grid
axis（[0?10*10^-6?-3?3]）

%逆流接收信號

signal_2=2*sin（2*pi*f*（time+tans_time2））+noise;
[ba]=fir1（32[0.19?0.21]）;
signal_2=filter（basignal_2）;

subplot（333）
plot（timesignal_2）;grid
axis（[0?10*10^-6?-3?3]）

%插值求相關
B=100;?????????????????????????%250倍插值
tt=（0:B*P*T*Fs）/（B*Fs）;???????%插枝點


xx=interp1（ti