应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】这个解PR方程的程序错哪里了？【已完成】
北京石油化工学院2026年研究生招生接收调剂公告
71/1返回列表
查看: 545  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

慢三儿

木虫 (小有名气)

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

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主

试试这个

★ ★ ★ ★ ★ ★ ★
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的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主
★
kuhailangyu(金币+1,VIP+0):欢迎补充 9-28 22:12
补充一点：fzero对初值比较敏感，不推荐。
赞一下(10人) 回复此楼
3楼2009-09-28 22:00:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

慢三儿

木虫 (小有名气)
nono2009(金币+0,VIP+0):可能还是计算的初值选得有问题。将计算条件贴出来看看。 9-29 10:23
还有 为啥分别取液态和气态的初值，计算的结果是一样滴？？
赞一下(6人) 回复此楼
4楼2009-09-29 09:49:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

慢三儿

木虫 (小有名气)
引用回帖:
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的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主
★ ★ ★ ★ ★
慢三儿(金币+5,VIP+0):华丽丽的金币全送给no版主吧,呵呵~~版主能给诉我我那段程序错误的地方吗？ 9-29 16:34
压力和温度给了，自然只有一个比容值（两相区除外）。

对超临界，一般只称超临界流体，没有气液之分。
赞一下(1人) 回复此楼
6楼2009-09-29 11:44:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主
★
小木虫(金币+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 ]
赞一下(1人) 回复此楼
7楼2009-09-30 00:21:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 慢三儿 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 281求调剂 +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
[考研] 11408,335分，本科211，求调剂，可转专业 +3 鳄梨大鳄鱼 2026-04-03 3/150 2026-04-03 21:18 by zhq0425
[考研] 310求调剂 +18 争取九点睡 2026-03-30 18/900 2026-04-03 18:35 by ls刘帅
[考研] 301求调剂 +15 骆驼男人 2026-04-02 15/750 2026-04-03 18:26 by ls刘帅
[考研] 工科 267求调剂 +5 wanwan00 2026-04-02 7/350 2026-04-03 14:14 by zhangdingwa
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +5 蓝笺片 2026-04-02 6/300 2026-04-02 21:37 by dongzh2009
[考研] 0856材料与化工调剂，339 +14 10213207 2026-03-31 14/700 2026-04-02 21:01 by 1104338198
[考研] 260求调剂 +6 朱芷琳 2026-04-02 6/300 2026-04-02 20:27 by 6781022
[考研] 。 +4 雾与海 2026-04-02 5/250 2026-04-02 19:16 by 土木硕士招生
[考研] 农学考研求调剂 +3 dkdkxm 2026-04-01 3/150 2026-04-02 16:04 by wangjagri
[考研] 367求调剂 +8 芋泥啵啵… 2026-03-28 8/400 2026-04-02 09:20 by olim
[考研] 材料求调剂 +8 呢呢妮妮 2026-04-01 8/400 2026-04-02 07:13 by yjolah
[考研] 086000生物与医药298调剂求助 +4 元元青青 2026-03-31 6/300 2026-04-01 11:13 by syh9288
[考研] 化学0703 调剂 306分 一志愿211 +12 26要上岸 2026-03-28 12/600 2026-04-01 11:10 by chemdavid
[考研] 一志愿武理材料工程302调剂环化或化工 +15 Doleres 2026-03-31 16/800 2026-04-01 09:49 by lfj11
[考研] 【调剂】一志愿厦大生物与医药调剂 +3 Echo虾米 2026-03-31 3/150 2026-04-01 08:40 by JourneyLucky
[考研] 调剂申请 +8 张张张张zy 2026-03-31 9/450 2026-04-01 08:29 by zjbkx
[考研] 080200学硕，机械工程专业277分，求带走！ +4 瓶子PZ 2026-03-31 4/200 2026-03-31 20:16 by vgtyfty
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
信息提示
关闭
请填处理意见
关闭