f2812做的超声波流量计 matlab仿真部分 有图

%产生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