| 查看: 276 | 回复: 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 |
回复此楼» 猜你喜欢
08工科求调剂290分
已经有12人回复
-
286求调剂
已经有16人回复
-
材料类284调剂
已经有24人回复
-
考研调剂-材料类-284
已经有14人回复
-
327求调剂
已经有11人回复
-
一志愿985初试354分生物调剂
已经有3人回复
-
一志愿2110，化学学硕310分，本科重点双非求调剂
已经有13人回复
-
353求调剂
已经有8人回复
-
308求调剂
已经有16人回复
-
334求调剂
已经有14人回复
