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 根据需求进行修改