声速的测量数据处理代码

#说明，我的距离都是负数，所以你们如果是正数需要修改for循环的内容
l1=[100.26,95.46,90.63,85.47,80.91,76.15,71.59,66.88,61.98,57.27,52.61,47.83,43.06,38.24,33.56,28.71]  #共振干涉法
l2=[100.14,95.56,90.37,86.08,81.26,76.44,71.8,67.03,62.12,57.56,52.76,48.07,43.28,38.44,33.79,28.87]  #相位比较法
l3=[34.71,28.37,22.08,16.92]  #时差法
t=[354,372,388,404]
f=[36378,36378,36378,36379,36378]  #谐振频率
T=23.4  #室温######以下内容根据数据修改##########
sum1=0
sum2=0
dL1=[]
dL2=[]
f1=sum(f)/5
vr=331.45+0.59*T
print('Li+8-Li分别为')
for i in range(8):dL1.append(l1[i]-l1[i+8])sum1+=l1[i]-l1[i+8]print(l1[i]-l1[i+8],end='  ')
print('\nLi+8-Li分别为')
for i in range(8):dL2.append(l2[i]-l2[i+8])sum2+=l2[i]-l2[i+8]print(l2[i]-l2[i+8],end='  ')
L1=sum1/8
L2=sum2/8
print('\n△L分别为 '+f'{L1}  {L2}')
r1=sum1/32
r2=sum2/32
print('λ分别为 '+f'{r1}  {r2}')
v1=f1*r1*0.001
v2=f1*r2*0.001
print('测量声速为 '+f'{v1}  {v2}')
E1=(v1-vr)/vr*100
E2=(v2-vr)/vr*100
print("E分别为 "+f'{round(E1,2)}%  {round(E2,2)}%')ufa=(((f[0]-f1)**2+(f[1]-f1)**2+(f[2]-f1)**2+(f[3]-f1)**2+(f[4]-f1)**2)/20)**0.5
ufb=1/(3)**0.5
ufc=((ufa)**2+(ufb)**2)**0.5
ufr=ufc/f1*100
print('频率不确定度a,b,c,r分别为 '+f'{ufa}  {ufb}  {ufc}  {ufr}%')uLa=(((dL1[0]-L1)**2+(dL1[1]-L1)**2+(dL1[2]-L1)**2+(dL1[3]-L1)**2+(dL1[4]-L1)**2+(dL1[5]-L1)**2+(dL1[6]-L1)**2+(dL1[7]-L1)**2)/56)**0.5
uLb=0.02/(3)**0.5
uLc=((uLa)**2+(uLb)**2)**0.5
uLr=uLc/L1*100
print('共振干涉法△L不确定度a,b,c,r分别为 '+f'{uLa}  {uLb}  {uLc}  {uLr}%')
uvr1=((ufc/f1)**2+(uLc/L1)**2)**0.5*100
uvc1=uvr1*v1/100
uLa=(((dL2[0]-L2)**2+(dL2[1]-L2)**2+(dL2[2]-L2)**2+(dL2[3]-L2)**2+(dL2[4]-L2)**2+(dL2[5]-L2)**2+(dL2[6]-L2)**2+(dL2[7]-L2)**2)/56)**0.5
uLb=0.02/(3)**0.5
uLc=((uLa)**2+(uLb)**2)**0.5
uLr=uLc/L1*100
print('相位比较法△L不确定度a,b,c,r分别为 '+f'{uLa}  {uLb}  {uLc}  {uLr}%')uvr2=((ufc/f1)**2+(uLc/L2)**2)**0.5*100
uvc2=uvr2*v2/100
print('共振干涉法v不确定度r,c分别为 '+f'{uvr1}%  {uvc1}')
print('相位比较法v不确定度r,c分别为 '+f'{uvr2}%  {uvc2}')v3=[]
for i in range(3):v3.append((l3[i]-l3[i+1])/(t[i+1]-t[i])*1000)
print('时差法测量声速为 '+f'{v3}')
print('平均声速为 '+f'{sum(v3)/3}')
print('相对百分误差为 '+f'{(sum(v3)/3-vr)/vr*100}%')

效果图如下

声速的测量数据处理代码相关推荐

大学物理实验_超声声速的测量_数据处理软件_源代码
大学物理实验_超声声速的测量_数据处理软件_源代码1 前言程序适用于大学物理实验_超声声速的测量的有关数据处理.本程序以NUAA大学物理实验中心所提供的的器材以及实验报告模板为基础编写. 代码(基于 ...
科大奥瑞物理实验——声速的测量
实验名称:声速的测量 1. 实验目的: (1)了解超声波的发射和接收方法. (2)加深对振动合成.波动干涉等理论知识的理解. (3)掌握用驻波法和相位法测声速. 2. 实验器材: 示波器声速仪信号 ...
声速的测量的实验原理和应用_大学物理实验报告声速测量
实验报告声速的测量 [ 实验目的] 1.学会用共振干涉法.相位比较法以及时差法测量介质中的声速 2.学会用逐差法进行数据处理; 3.了解声速与介质参数的关系. [ 实验原理] 由于超声波具有波 ...
pandas数据处理代码_使用Pandas方法链接提高代码可读性
pandas数据处理代码 We have been talking about using the Pandas pipe function to improve code readability. ...
声速的测量的实验原理和应用_声速的测定实验报告心得体会
测量声速的实验报告 1. 提出问题如何测出声音的速度? 2. 猜想与假设如果在一定距离内听到声音要多少时间? 3. 实验步骤步骤应该就是实施实验,第三是实验器材的话,就是要秒表. 4. 实施实验 ...
python测试代码与模块_测量Python代码运行时间
Python 社区有句俗语: "python自己带着电池" ,别自己写计时框架. Python 2.3 具备一个叫做 timeit 的完美计时工具可以测量python代码的运行时间 ...
测量数据的计算机处理实验报告,测量数据处理实验报告..doc
测量数据处理实验报告. 院系____________________ 班级____________________ 姓名____________________ 学号_________________ ...
canbus数据处理代码分析
canbus数据处理代码分析下面以刹车控制为例进行分析在**_controller.cc中, brake_command_102_ = dynamic_cast<Brakecommand10 ...
中南大学实用测量数据处理实验
目录实验一:总体最小二乘估计实验二:线性回归分析实验三:非线性回归分析实验四:平滑与滤波中南大学<实用测量数据处理>实验(自用) 实验一:总体最小二乘估计为了确定直线方程 y ...
声速的测量的实验原理和应用_示波器的原理和使用声速测量实验报告.docx
示波器的原理和使用声速测量实验报告.docx 示波器的原理和使用.声速测量一实验目的1了解示波器的基本结构及其工作原理,学习并掌握示波器的基本使用方法2学习电信号有关参数的基本概念及其测量3了解声波在 ...

声速的测量数据处理代码

声速的测量数据处理代码相关推荐

最新文章

热门文章