| 查看: 884 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】这个程序检查没有问题，但是运行显示的是BUSY? 已有3人参与
|
|||
|
%双线性LANGMUIR公式 %b1,b2是第一个LANGMUIR公式的常数，b3,b4是第二个LANGMUIR公式的常数 %得到只含有x1的三次方程，迭代求解 function [p1,p2,x1]=langmuir(b1,b2,b3,b4,l,y1) m=b1+b2; n=b1*b2; p=b3+b4; q=b3*b4; r=m*y1-p*y1+p; s=2*m*y1-n*p*y1^2-p*y1+p; t=m*y1-2*n*p*y1^2-q*p*(y1-1)^2; u=n*p*y1^2; x0=0.01; f=r*x0^3-s*x0^2+t*x0+u; g=3*r*x0^2-2*s*x0+t; f1=subs(sym(f),findsym(sym(f)),x0); f2=subs(sym(g),findsym(sym(g)),x0); x1=x0-(f1/f2); while((x1-x0)<1.0e-03) end p1=l*y1/x1; p2=l*(1-y1)/(1-x1); end |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
这是错误吗？ 5楼2010-06-16 10:19:54
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★
zyj8119(金币+5):多谢，程序这样写就更短了。。。 2010-06-16 09:26:31
robert2020(金币+2):多谢应助！辛苦了！ 2010-06-16 12:28:25
zyj8119(金币+5):多谢，程序这样写就更短了。。。 2010-06-16 09:26:31
robert2020(金币+2):多谢应助！辛苦了！ 2010-06-16 12:28:25
|
function [p1,p2,x1]=langmuir(b1,b2,b3,b4,l,y1) m=b1+b2; n=b1*b2; p=b3+b4; q=b3*b4; r=m*y1-p*y1+p; s=2*m*y1-n*p*y1^2-p*y1+p; t=m*y1-2*n*p*y1^2-q*p*(y1-1)^2; u=n*p*y1^2; x1=fsolve(@(x)r*x^3-s*x^2+t*x+u-(3*r*x^2-2*s*x+t),0.3); p1=l*y1/x1; p2=l*(1-y1)/(1-x1); |
2楼2010-06-16 08:49:07
4楼2010-06-16 09:57:27
6楼2010-06-16 10:27:07
100