大学物理实验——密立根油滴实验数据计算Python代码
先手准备
import numpy as np
编写函数
def f(t,U):# 基本物理量准备rou = 981.0g = 9.80l = 2.00*10**-3b = 8.23*10**-3p = 1.01*10**5d = 5.00*10**-3n = 1.83*10**-5r = np.sqrt((9*n*l)/(2*rou*g*t))e = 1.60*10**-19# 计算准备k = (18*np.pi)/np.sqrt(2*rou*g)x = (n*l)/(t*(1+b/(p*r)))y = d/Uprint("r",l/(1+b/(p*r)))# 开始计算q = k*x**(3/2)*ynum = (q/e)e1 = q/numprint("电荷带电量为：",q)print("可能有这些个电荷：",num)print("我算出来的e：",e1)print("=="*10+"我是分割线"+"=="*10)return l
成果展示
time = np.array([23.9,24.8,23.9,24.1,24.6,18.5,19.0,18.7,18.2,19.0,18.7,19.1,19.5,19.1,19.6,27.9,28.4,28.0,28.6,28.7,23.6,23.9,23.9,24.1,23.9])time.shape = (5,5)U = np.array([239,214,223,262,189])for i in range(5):l = f(time[i].mean(),U[i])
