| 查看: 875 | 回复: 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
已经有6人回复
-
关于代码变化问题，想知道的进来
已经有8人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有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
3楼2010-06-16 09:25:44
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
6楼2010-06-16 10:27:07
7楼2010-06-16 10:27:28
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
- 专业: 化学反应工程
★
zzuwangshilei(金币+1):多谢给出建议 2010-06-16 16:18:57
zzuwangshilei(金币+1):多谢给出建议 2010-06-16 16:18:57
|
不懂就不要装懂，安心坐下来看书就是一种罪过吗？！！ 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; options=optimset('display','off'); x1=fsolve(@(x)r*x^3-s*x^2+t*x+u-(3*r*x^2-2*s*x+t),0.3,options); p1=l*y1/x1; p2=l*(1-y1)/(1-x1); [ Last edited by change0618 on 2010-6-16 at 13:34 ] |
9楼2010-06-16 13:31:36
10楼2010-06-16 14:35:28
100