| 查看: 946 | 回复: 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 |
回复此楼» 猜你喜欢
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
4楼2010-06-16 09:57:27
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
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
这是错误吗？ 5楼2010-06-16 10:19:54
6楼2010-06-16 10:27:07
30