密立根油滴实验计算代码

  • 密立根油滴实验动态法计算公式和预设参数值
    • 计算代码(python实现):
    • 代码测试:

密立根油滴实验动态法计算公式和预设参数值

密立根油滴实验是一个著名的实验,其中有用非平衡法(动态法)测量油滴电荷量的步骤,在做实验时计算量繁琐,这里用一个python小程序提供辅助计算。
其中主要的计算公式和相关参数值如下:
注意:不同实验室里装置不同,相关参数会有所不同,如有不同请在程序开始位置进行修改。
求油滴电荷量的公式:q=18πη′32L322g(ρ−ρ′)⋅dU′(1tg+1te)(1tg)12油滴半径:r=9ηvg2(ρ−ρ′)gη的修正函数:η′=η1+brP{\begin{aligned} 求油滴电荷量的公式 :q &= \frac{ 18\pi\eta'^{\frac{3}{2}} L^{\frac{3}{2}} }{\sqrt{2g(\rho-\rho')} } \cdot\frac{d}{U'} \big( \frac{1}{t_g} + \frac{1}{t_e} \big)\big( \frac{1}{t_g} \big)^{\frac{1}{2}}\\ \\ 油滴半径:r&=\sqrt{\frac{9\eta v_g}{2(\rho-\rho')g}}\\ \\ \eta的修正函数:\eta'&=\frac{\eta}{1+\frac{b}{rP}}\\ \end{aligned} }\\ 求油滴电荷量的公式:q油滴半径:rη的修正函数:η′​=2g(ρ−ρ′)​18πη′23​L23​​⋅U′d​(tg​1​+te​1​)(tg​1​)21​=2(ρ−ρ′)g9ηvg​​​=1+rPb​η​​
相关参数的默认值:(有不同的需要在代码中进行简单修改)空气压强:P=76.0cmHg修正常数:b=6.17×10−6m⋅cmHg粘滞系数:η=1.83×10−5kg⋅m−1⋅s−1平行极板间距离:d=5.00×10−3m匀速上升/下降的距离:L=1.6mm重力加速度:g=9.794m/s2油的密度:ρ=981kg⋅m−3空气密度:ρ′=1.2928kg⋅m−3匀速上升的时间te:测量得到匀速下降的时间tg:测量得到油滴的上升电压U′:测量得到\begin{aligned} 相关参数的默认值&:(有不同的需要在代码中进行简单修改)\hspace{5cm}\\ 空气压强&:P=76.0cmHg\\ 修正常数&:b=6.17\times 10^{-6}m\cdot cmHg\\ 粘滞系数&:\eta = 1.83\times10^{-5}kg\cdot m^{-1} \cdot s^{-1}\\ 平行极板间距离&:d=5.00\times 10^{-3}m\\ 匀速上升/下降的距离&:L=1.6mm\\ 重力加速度&:g=9.794m/s^{2}\\ 油的密度&:\rho=981kg\cdot m^{-3}\\ 空气密度&:\rho'=1.2928kg\cdot m^{-3}\\ 匀速上升的时间t_e&:测量得到\\ 匀速下降的时间t_g&:测量得到\\ 油滴的上升电压U'&:测量得到\\ \end{aligned} 相关参数的默认值空气压强修正常数粘滞系数平行极板间距离匀速上升/下降的距离重力加速度油的密度空气密度匀速上升的时间te​匀速下降的时间tg​油滴的上升电压U′​:(有不同的需要在代码中进行简单修改):P=76.0cmHg:b=6.17×10−6m⋅cmHg:η=1.83×10−5kg⋅m−1⋅s−1:d=5.00×10−3m:L=1.6mm:g=9.794m/s2:ρ=981kg⋅m−3:ρ′=1.2928kg⋅m−3:测量得到:测量得到:测量得到​

计算代码(python实现):

使用时确保这些相关参数与自己实验的参数设置值相同,不同需要要对数值进行相应修改。

测量次数默认为5次,也可调整。

from decimal import *
import math
getcontext().prec = 24text_times = 5 # 设置每个油滴上下往返时间的测量次数,默认为5
# 相关参数值(不同的需要修改)
P = Decimal('76')  # 空气压强
b = Decimal("6.17e-6")  # 修正常数
n = Decimal("1.83e-5")  # 粘滞系数
d = Decimal('0.005')  # 平行极板间距离
L = Decimal('0.0016')  # 匀速上升/下降的距离
g = Decimal('9.794')  # 重力加速度
rou = Decimal(918)  # 油的密度
rou2 = Decimal('1.2928')  # 空气密度
pi = Decimal(str(math.pi))  # 圆周率π
e1 = Decimal("1.6021892e-19")  # 理论元电荷值data_tg = []  # 用于存储下降时间
data_te = []  # 用于存储下降时间
for i in range(1, text_times+ 1 ):tg = input("输入匀速下降时间:tg" + str(i) + "= ")data_tg.append(Decimal(tg))
tgp = sum(data_tg)/Decimal(text_times)
print("计算得到的tg平均值:"+str(sum(data_tg)/Decimal(text_times)) )
for i in range(1, text_times+ 1 ):te = input("输入匀速上升时间:te" + str(i) + "= ")data_te.append(Decimal(te))
tep = sum(data_te)/Decimal(text_times)
print("计算得到的te平均值:"+str(sum(data_te)/Decimal(text_times)) )
U = input("输入上升电压:U'= ")
U = Decimal(U)#  计算公式
r = (9*n*L / (2*(rou - rou2)* g * tgp) ) ** Decimal('0.5')
n2 = n/(1+b/(r*P))  # 修正后的粘滞系数
K = ( 18*pi*(n2**Decimal('1.5'))*(L**Decimal('1.5')) / ( (2*g*(rou-rou2))**Decimal('0.5') ) ) * d
q = (K/U) * (1/tgp + 1/tep) * ((1/tgp)**Decimal('0.5'))
N = q/e1
#  结果显示
print("油滴半径:r=",end='');print("%e" % r)
print("r=",end='');print(r)
print("修正后的粘滞系数:n2=",end='');print("%e" % n2)
print("n2=",end='');print(n2)
print("油滴所带的电荷量:q=",end='');print('%e'% q)
print("q=",end='');print(q)
print('q/e:N=',end='');print(N)
input("结束,按任意键继续")
# 这是一组测试数据
# 12.59
# 12.61
# 12.29
# 12.94
# 12.88
# 7.97
# 8.11
# 8.07
# 8.21
# 8.05
# 342

程序中为了保证数据计算的精确性,引入了decimal库。

代码测试:

用代码中的测试数据进行测试:

结果中有两个值,一个是进行了科学计数法表示,一个是原始数据。

有问题欢迎交流指正,如果觉得有用的话还求一个小小的赞(〃‘▽’〃)。

密立根油滴实验数据处理,油滴电荷量计算,简单复制即可用相关推荐

  1. 密立根油滴实验数据处理(平衡法和动态法)

    密立根油滴实验数据处理(平衡法和动态法) 平衡法 //平衡法 #include <iostream> #include <cmath> using namespace std; ...

  2. 基于python对密立根油滴实验数据处理改进

    基于python对密立根油滴实验数据处理改进 近似最大公约数 import math ac=1.34 def accuracy(x, y):return math.fabs(x - y) < a ...

  3. 科大奥锐密立根油滴实验数据_密立根油滴实验数据处理分析

    第 01 卷第 3期 1 9 9 4 年 9 月 黄 淮 H U A N G H U A I 学 刊 J O U R N A L V o l . 1 0 N o . 3 Se P t . 1 9 9 ...

  4. 密立根油滴实验数据处理(基于Python)

    import math as m import numpy as nu QQ_average = [] # 电荷量均值 NN_average = [] # 电子数量均值 QQNN = [] # 电子数 ...

  5. 科大奥瑞物理实验——密立根油滴实验

    实验名称:密立根油滴实验 1. 实验目的: 学习控制油滴在视场中的运动,并选择合适的油滴测量元电荷,要求测得6个不同的油滴8次以上. 2. 实验器材: 密立根油滴仪 显示器 油滴管 实验总体装置 3. ...

  6. 密里根油滴实验实战代码(静态法与动态法)

    学习控制油滴在视场中的运动,并选择合适的油滴测量元电荷. 1. 选择适当的油滴并测量油滴上所带电荷 要做好油滴实验,所选的油滴体积要适中,大的油滴虽然比较亮,但下降速度快,不容易测准确:太小则受布朗运 ...

  7. 密立根测油滴实验c语言测试,密立根油滴实验-平衡法-数据处理-C程序源代码

    <密立根油滴实验-平衡法-数据处理-C程序源代码>由会员分享,可在线阅读,更多相关<密立根油滴实验-平衡法-数据处理-C程序源代码(2页珍藏版)>请在人人文库网上搜索. 1.* ...

  8. 密立根测油滴实验c语言测试,密立根油滴实验

    一.实验目的 1.了解密立根油滴仪的结构,油滴实验测定电子电荷的设计思想和方法. 2.了解CCD图像传感器的原理和电视显微测量方法. 3.通过对带电油滴在重力场和静电场中运动的测量,来测量基本电荷 的 ...

  9. 计算机课件模拟实验,计算机模拟仿真密立根油滴实验课件.pdf

    计算机模拟仿真密立根油滴实验 倪敏 , (上海师范大学数理信息学院上海200234) 提要:根据教学实际,用计算机模拟仿真密立根油滴实验.使实验中的原理,仪器.操作方法,数据处理等形 象地呈现在计算机 ...

最新文章

  1. 第四次作业 结对编程-黄金点游戏
  2. 【SharePoint 2010】SharePoint 2010开发方面的课堂中整理有关问题
  3. C++开发中类的定制
  4. Shell命令-文件及内容处理之split、paste
  5. spring 同时配置hibernate and jdbc 事务
  6. 【Tika基础教程之一】Tika基础教程
  7. scroll-view——小程序横向滚动
  8. 分词消除歧义_角色标题消除歧义
  9. 面向物联网的可重构流式深度卷积神经网络加速器
  10. python爬取英雄联盟所有皮肤_python爬取lol英雄皮肤
  11. 《一个人的朝圣》:爱,才是救赎一切的良药
  12. 微信扫描自定义二维码关注公众号并登录
  13. RN 组件 Animated
  14. 懂车帝视频连接--JS逆向
  15. XXO客户端携带WUP包体进行协议测试
  16. mac m1上esc键失灵不能退出vi解决方法
  17. 增长率超 100%!东软数据可视化到底什么样?
  18. 联想小新pad安装Google Play教程
  19. proxychans的使用
  20. 端口号被占用的解决方法

热门文章

  1. 51系列—基于51单片机的校园教室打铃系统(代码+仿真+电路图超级全)
  2. Servlet生命周期及请求、响应
  3. MATLAB type文件名,Matlab产生IGES文件代码
  4. 刘未鹏逃出肖申克系列链接
  5. texturePacker 指令 .pvr.ccz 转 .png, 报错 : Not all sprites could be packed into the texture!
  6. 桌面word excel ppt图标变白怎么解决?
  7. JZOJ 3896. 【NOIP2014模拟10.26】战争游戏
  8. 潍坊医学院本科毕业论文答辩PPT模板
  9. 为什么只有义乌个体户可以无限制结汇美金,什么原因
  10. 最简单Visio导出无边框无白边pdf