求声速

// 相位法和驻波法只是记录数据的方法不同处理数据是一样的
L=input('请输入所有的L组成的数组（单位是mm）L=');
n=length(L)/2;%确定数组长度
Ln=L(n+1:2*n);%分成后半和前半数组，方便使用逐差法处理数据
Lk=L(1:n);
lambda_2=(Ln-Lk)./n;%用数组减法对应元素相互减，点除得波长一半
lambda_2_=sum(lambda_2)/n%求平均值
f=input('请输入频率（单位是kHz）f=');
v1=lambda_2_*2*f%测量出的声速
t=input('请输入温度（单位是℃）t=');
T0=273.15;%K
v0=331.45;%m/s
v2=v0*sqrt((t+T0)/T0)%当前温度下理论声速
epsilon=abs(v1-v2)/v2%abs是求绝对值函数

如果有更好的思路或者一些建议，请私信我，我是在自行摸索，难免会有不规范或者有差错的地方。