| 查看: 872 | 回复: 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变化情况
已经有15人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
-
filecode
已经有19人回复
-
基金中了
已经有11人回复
-
听说今天filecode变了
已经有44人回复
-
固定端突然变了，今天
已经有8人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
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
50