密里根油滴实验数据matlab实现

clc,clear;

Ut=; %动态法中提升电压

Up=; %静态法中平衡电压

tr=; %为上升时间；

tf=; %为下降时间；

vf=0.002/tf;

vr=0.002/tr;

g=9.794;

rou1=981;
rou2=1.293;

eta=1.83*10^-5;

d2=0.002；

b=8.22*10^-3;

p=1.013*10^5;

d1=0.005;

q1=1.602*10^-19;

r=(9*eta*vf/2/g/(rou1-rou2))^0.5

k1=9*2^0.5*pi*d1*((eta*d2)^3/(rou1-rou2)/g)^0.5/Ut*(1/(1+b/p/r))^1.5;

disp('动态法带电量');

qd=k1*(1/tf+1/tr)*(1/tf)^0.5 %动态法带电量

k2=9*2^0.5*pi*d1*((eta*d2)^3/(rou1-rou2)/g)^0.5/Up*(1/(1+b/p/r))^1.5;

disp('静态法带电量');

qj=k2*(1/tf)^1.5 %静态法带电量

n1=round(qd/q1); %动态法电荷数

n2=round(qj/q1); %静态法电荷数

disp('动态法基本电荷带电量');

jqd=qd/n1 %动态法基本电荷带电量

disp('静态法基本电荷带电量');

jqj=qj/n2 %静态法基本电荷带电量

disp('动态法相对误差');

ud=abs(((jqd-q1)/q1)*100) %动态法相对误差

disp('静态法相对误差');

uj=abs(((jqj-q1)/q1)*100) %静态法相对误差

提供平衡电压，上升电压，下降时间，上升时间即可求出。

其余参数g rou1 rou2 eta p d 根据需求进行修改