24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

慢三儿

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 4256.4
散金: 21
红花: 1
帖子: 234
在线: 91.3小时
虫号: 433712
注册: 2007-08-18
专业: 化学反应工程

[交流] 【求助】这个解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 ]

回复此楼

» 猜你喜欢

1楼 2009-09-28 20:05:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

慢三儿

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 4256.4
散金: 21
红花: 1
帖子: 234
在线: 91.3小时
虫号: 433712
注册: 2007-08-18
专业: 化学反应工程

引用回帖:

Originally posted by 慢三儿 at 2009-9-29 09:49:
还有为啥分别取液态和气态的初值，计算的结果是一样滴？？

气态初值：V0=R*T/P
液态初值： V0=2*b

难道是超临界，气液不分了，就出现这样的情况？

赞一下

回复此楼

5楼2009-09-29 10:29:47

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 7 个回答

nono2009

超级版主 (文学泰斗)

No gains, no pains.

专家经验: +21105
仿真EPI: 1
应助: 28684 (院士)
贵宾: 513.911
金币: 2555220
散金: 27828
红花: 2147
沙发: 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

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

赞一下(13人)

回复此楼

2楼2009-09-28 20:26:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖