| 查看: 880 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
8楼2010-06-16 12:19:31
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
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
这是错误吗？ 5楼2010-06-16 10:19:54
100