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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】这个解PR方程的程序错哪里了？【已完成】
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 549  |  回复: 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.

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

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

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的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本科211，293分请求调剂 +4 莲菜就是藕吧 2026-04-03 4/200 2026-04-04 05:18 by gswylq
[考研] 265求调剂 +17 林深温澜 2026-04-01 20/1000 2026-04-04 01:09 by userper
[考研] 282求调剂 +20 ycy1201 2026-04-01 22/1100 2026-04-04 00:42 by userper
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 322求调剂 +6 FZAC123 2026-04-03 6/300 2026-04-03 22:23 by 科研小专家
[考研] 281求调剂 +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
[考研] 085601，一志愿厦大334复试被刷求调剂 +10 曾仰之 2026-04-03 12/600 2026-04-03 21:41 by zhq0425
[考研] 295求调剂 +6 FZAC123 2026-04-03 6/300 2026-04-03 21:01 by zhq0425
[考研] 材料调剂 +15 一样YWY 2026-04-02 15/750 2026-04-03 15:39 by 晴空210210
[考博] 申博求助 +3 Reee1Llll 2026-04-01 3/150 2026-04-02 22:29 by 这是一个无聊的�
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
[考研] 379求调剂 +3 ?苦瓜不苦 2026-04-01 3/150 2026-04-01 20:09 by 暮云清寒
[考研] 353求调剂 +4 拉钩不许变 2026-04-01 4/200 2026-04-01 18:10 by 记事本2026
[考研] 材料专业调剂 +5 啦啦啦哭 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[考研] 352分-085602-一志愿985 +6 海纳百川Ly 2026-03-29 6/300 2026-03-31 21:06 by yuq
[考研] 材料科学与工程求调剂 +13 深V宿舍吧 2026-03-29 13/650 2026-03-31 19:50 by Dyhoer
[考研] 江苏苏北高校诚邀调剂同学 +3 zzll406 2026-03-31 3/150 2026-03-31 16:54 by 及时行乐fan
[考研] 0856 335分 +9 cccchenso 2026-03-29 9/450 2026-03-31 16:37 by lishahe
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-03-31 4/200 2026-03-31 11:56 by jp9609
[考研] 本科211总分289，08工学真心求调剂 +3 utopiaE 2026-03-30 3/150 2026-03-30 23:42 by ms629
信息提示
关闭
请填处理意见
关闭