| 查看: 258 | 回复: 0 | |||
xb490284223金虫 (正式写手)
蘑菇头
|
[求助]
密封舱浓度足迹法 求解方程问题
|
|
各位大神，对于下图的足迹法，我编写了如下的程序，但计算结果不满意，求各位帮忙。程序是 二分法求根的。 Cequ=1.5; slope=-0.21/3600; Intercept=-1.75; hm=0.0025; L=0.0038; B=0.33333; a=pi/2+0.00000001; b=pi*3/2-0.000000001; m=1; while(abs(m)>10^-10) q1=(a+b)/2; Dm=-slope*L*L/(q1*q1); %%求Dm Bim=hm*L/Dm; %%求Bim K=q1/((q1*q1/Bim-B)*tan(q1)); %%求K A1=(K*B-q1*q1*K/Bim+1)*cos(q1)-(1+2*K/Bim)*q1*sin(q1); m=log(-2*(K*B+1)*sin(q1)/q1/A1)-Intercept; Cm0=Cequ*K+1/B; %%求Cm0 if(m>0) a=q1; else b=q1; end if(abs(m)<10^-10) q1end=q1; Kend=K; Dmend=Dm; Cm0end=Cm0; m=abs(m); break; end end disp('q1是') q1end disp('初始散发浓度C0是') Cm0end disp('分离系数K是') Kend disp('扩散系数是') Dmend disp('m是') m  密封舱浓度足迹法流程图.png |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复