| 查看: 653 | 回复: 6 | |||||
| 当前主题已经存档。 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
慢三儿木虫 (小有名气)
|
[交流]
【求助】这个解PR方程的程序错哪里了？【已完成】
|
||||
|
程序如下： Tc=304.19;Pc=7.381*10^6;w=0.225;R=8.3145; T=140+273.15; P=1440/145*10^6; Tr=T/Tc; alfa=(1+(0.37464+1.54226*w-0.26992*w^2)*(1-Tr^0.5))^2; ac=0.457235*(R*Tc)^2/Pc; b=0.077796*(R*Tc)/Pc; a=ac*alfa; V0=R*T/P; V=fzero(@PReq,V0,P,T,a,b,R) function f=PReq(V,P,T,a,b,R) f=P-R*T/(V-b)+a/(V*(V+b)+b*(V-b)); 运行时老提示 Error in ==> SPR at 10 V=fzero(@PReq,V0,P,T,a,b,R); why？？？ [ Last edited by nono2009 on 2009-10-5 at 08:11 ] |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有25人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有8人回复
-
分享一下我之前已中青C的计划书的filecode
已经有4人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有3人回复
-
为什么网上很多人说本周 12号出结果
已经有6人回复
-
听说今天filecode变了
已经有49人回复
-
帮忙看看fileCode
已经有12人回复
-
FileCode能看出啥？
已经有17人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有6人回复
-
filecode
已经有9人回复
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 仿真EPI: 1
- 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
|
估计是fzero的使用问题，再仔细读一下Matlab online help吧。 另外，我对比了fsolve和fzero的计算结果（如下），说明fsolve比fzero对初值的稳定性较高。 >> Tc=304.19;Pc=7.381e+6;w=0.225;R=8.3145; T=140+273.15; P=1440/145*1e6; Tr=T/Tc; alfa=(1+(0.37464+1.54226*w-0.26992*w^2)*(1-Tr^0.5))^2; ac=0.457235*(R*Tc)^2/Pc; b=0.077796*(R*Tc)/Pc; a=ac*alfa; V0=R*T/P; f=@(V)P-R*T/(V-b)+a/(V*(V+b)+b*(V-b)); V=fzero(f,V0); V V = 2.8961e-004 结果正确！ >> Tc=304.19;Pc=7.381e+6;w=0.225;R=8.3145; T=140+273.15; P=1440/145*1e6; Tr=T/Tc; alfa=(1+(0.37464+1.54226*w-0.26992*w^2)*(1-Tr^0.5))^2; ac=0.457235*(R*Tc)^2/Pc; b=0.077796*(R*Tc)/Pc; a=ac*alfa; V0=R*T/P; f=@(V)P-R*T/(V-b)+a/(V*(V+b)+b*(V-b)); V=fsolve(f,V0); V Optimization terminated: norm of relative change in X is less than max(options.TolX^2,eps) and sum-of-squares of function values is less than sqrt(options.TolFun). V = 2.8961e-004 结果正确！ >> Tc=304.19;Pc=7.381e+6;w=0.225;R=8.3145; T=140+273.15; P=1440/145*1e6; Tr=T/Tc; alfa=(1+(0.37464+1.54226*w-0.26992*w^2)*(1-Tr^0.5))^2; ac=0.457235*(R*Tc)^2/Pc; b=0.077796*(R*Tc)/Pc; a=ac*alfa; V0=2*b; f=@(V)P-R*T/(V-b)+a/(V*(V+b)+b*(V-b)); V=fzero(f,V0); V V = 2.6658e-005 结果错误！ >> f(V) ans = -1.1779e+019 残差很大！ >> Tc=304.19;Pc=7.381e+6;w=0.225;R=8.3145; T=140+273.15; P=1440/145*1e6; Tr=T/Tc; alfa=(1+(0.37464+1.54226*w-0.26992*w^2)*(1-Tr^0.5))^2; ac=0.457235*(R*Tc)^2/Pc; b=0.077796*(R*Tc)/Pc; a=ac*alfa; V0=2*b; f=@(V)P-R*T/(V-b)+a/(V*(V+b)+b*(V-b)); V=fsolve(f,V0); V Optimization terminated: norm of relative change in X is less than max(options.TolX^2,eps) and sum-of-squares of function values is less than sqrt(options.TolFun). V = 2.8961e-004 结果正确！ [ Last edited by nono2009 on 2009-9-30 at 00:23 ] |
7楼2009-09-30 00:21:54
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 仿真EPI: 1
- 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
试试这个
★ ★ ★ ★ ★ ★ ★
kuhailangyu(金币+2,VIP+0):感谢double no的积极参与，o(∩_∩)o... 9-28 22:12
慢三儿(金币+5,VIP+0):谢谢！偶还想问问偶那段程序哪里错了呢？偶刚学，小菜一个...望高手指教(*^__^*) 9-29 09:14
kuhailangyu(金币+2,VIP+0):感谢double no的积极参与，o(∩_∩)o... 9-28 22:12
慢三儿(金币+5,VIP+0):谢谢！偶还想问问偶那段程序哪里错了呢？偶刚学，小菜一个...望高手指教(*^__^*) 9-29 09:14
|
Tc=304.19;Pc=7.381e+6;w=0.225;R=8.3145; T=140+273.15; P=1440/145*1e6; Tr=T/Tc; alfa=(1+(0.37464+1.54226*w-0.26992*w^2)*(1-Tr^0.5))^2; ac=0.457235*(R*Tc)^2/Pc; b=0.077796*(R*Tc)/Pc; a=ac*alfa; V0=R*T/P; f=@(V)P-R*T/(V-b)+a/(V*(V+b)+b*(V-b)); V=fsolve(f,V0); V |
2楼2009-09-28 20:26:50
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 仿真EPI: 1
- 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
3楼2009-09-28 22:00:46
慢三儿
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4256.4
- 散金: 21
- 红花: 1
- 帖子: 234
- 在线: 91.3小时
- 虫号: 433712
- 注册: 2007-08-18
- 专业: 化学反应工程
4楼2009-09-29 09:49:00
