大学物理实验——密立根油滴实验数据计算Python代码

1
2
import numpy as np

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
def f(t,U):
    # 基本物理量准备
    rou = 981.0
    g = 9.80
    l = 2.00*10**-3
    b = 8.23*10**-3
    p = 1.01*10**5
    d = 5.00*10**-3
    n = 1.83*10**-5
    r = 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/U
    print("r",l/(1+b/(p*r)))
    # 开始计算
    q = k*x**(3/2)*y
    num = (q/e)
    e1 = q/num
    print("电荷带电量为：",q)
    print("可能有这些个电荷：",num)
    print("我算出来的e：",e1)
    print("=="*10+"我是分割线"+"=="*10)
    return l

1
2
3
4
5
6
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])