| 查看: 312 | 回复: 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 |
回复此楼» 猜你喜欢
关于Filecode分析方法
已经有6人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有9人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
FileCode能看出啥？
已经有21人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有35人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
