| 查看: 270 | 回复: 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 |
回复此楼» 猜你喜欢
求调剂
已经有3人回复
-
265求调剂
已经有4人回复
-
085700资源与环境308求调剂
已经有6人回复
-
一志愿吉林大学材料学硕321求调剂
已经有12人回复
-
286分人工智能专业请求调剂愿意跨考！
已经有3人回复
-
329求调剂
已经有5人回复
-
申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复
已经有4人回复
-
材料学硕318求调剂
已经有5人回复
-
一志愿中国海洋大学，生物学，301分，求调剂
已经有6人回复
-
081700化工学硕调剂
已经有3人回复
10