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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】急，还是matlab 解方程求助！【已完成】
北京石油化工学院2026年研究生招生接收调剂公告
151/2返回列表12下一页
查看: 816  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

fishfly80

铁虫 (小有名气)

[交流] 【求助】急，还是matlab 解方程求助！【已完成】

用了yeyunxiao提供的解一般方程 结果都对，可用来解pr方程，结果都是符号，这是为什么，数据也检验好几次都没错，本人刚刚学，望请高手不吝赐教，syms v;
f = 8.314* 77 / (v-0.00001649839676) - 0.025825/ (v^2+2*v*0.00001649839676 - 0.00001649839676^2) - 100000;
v = solve(f,v)

v =                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                8641733483570864607155275934720181524676991/(1914088313930278856981482166803696189440000000000*((242207590675292455113936516646304255883544074967378025331707206971586057602285690330878181207807615479640055200706696641271^(1/2)*10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000^(1/2))/10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000 + 12703585355945088984479032173729896672868785889114677239189359257/1324076147308127973621611831457867134646590601566632280064000000000000000)^(1/3)) + ((242207590675292455113936516646304255883544074967378025331707206971586057602285690330878181207807615479640055200706696641271^(1/2)*10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000^(1/2))/10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000 + 12703585355945088984479032173729896672868785889114677239189359257/1324076147308127973621611831457867134646590601566632280064000000000000000)^(1/3) + 2944691389338352349149/1383505805528216371200000
(3^(1/2)*i*(8641733483570864607155275934720181524676991/(1914088313930278856981482166803696189440000000000*((242207590675292455113936516646304255883544074967378025331707206971586057602285690330878181207807615479640055200706696641271^(1/2)*10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000^(1/2))/10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000 + 12703585355945088984479032173729896672868785889114677239189359257/1324076147308127973621611831457867134646590601566632280064000000000000000)^(1/3)) - (1/10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000*242207590675292455113936516646304255883544074967378025331707206971586057602285690330878181207807615479640055200706696641271^(1/2)*10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000^(1/2) + 12703585355945088984479032173729896672868785889114677239189359257/1324076147308127973621611831457867134646590601566632280064000000000000000)^(1/3)))/2 - ((242207590675292455113936516646304255883544074967378025331707206971586057602285690330878181207807615479640055200706696641271^(1/2)*10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000^(1/2))/10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000 + 12703585355945088984479032173729896672868785889114677239189359257/1324076147308127973621611831457867134646590601566632280064000000000000000)^(1/3)/2 - 8641733483570864607155275934720181524676991/(3828176627860557713962964333607392378880000000000*((242207590675292455113936516646304255883544074967378025331707206971586057602285690330878181207807615479640055200706696641271^(1/2)*10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000^(1/2))/10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000 + 12703585355945088984479032173729896672868785889114677239189359257/1324076147308127973621611831457867134646590601566632280064000000000000000)^(1/3)) + 2944691389338352349149/1383505805528216371200000
2944691389338352349149/1383505805528216371200000 - ((242207590675292455113936516646304255883544074967378025331707206971586057602285690330878181207807615479640055200706696641271^(1/2)*10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000^(1/2))/10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000 + 12703585355945088984479032173729896672868785889114677239189359257/1324076147308127973621611831457867134646590601566632280064000000000000000)^(1/3)/2 - (3^(1/2)*i*(8641733483570864607155275934720181524676991/(1914088313930278856981482166803696189440000000000*((242207590675292455113936516646304255883544074967378025331707206971586057602285690330878181207807615479640055200706696641271^(1/2)*10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000^(1/2))/10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000 + 12703585355945088984479032173729896672868785889114677239189359257/1324076147308127973621611831457867134646590601566632280064000000000000000)^(1/3)) - (1/10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000*242207590675292455113936516646304255883544074967378025331707206971586057602285690330878181207807615479640055200706696641271^(1/2)*10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000^(1/2) + 12703585355945088984479032173729896672868785889114677239189359257/1324076147308127973621611831457867134646590601566632280064000000000000000)^(1/3)))/2 - 8641733483570864607155275934720181524676991/(3828176627860557713962964333607392378880000000000*((242207590675292455113936516646304255883544074967378025331707206971586057602285690330878181207807615479640055200706696641271^(1/2)*10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000^(1/2))/10389200852564950573550137526405054421261259377954756445650428924947602420218991079181495985174733616228310114631680000000000000000000000000000 + 12703585355945088984479032173729896672868785889114677239189359257/1324076147308127973621611831457867134646590601566632280064000000000000000)^(1/3))

[ Last edited by nono2009 on 2009-9-21 at 19:28 ]
回复此楼

» 猜你喜欢
1楼 2009-09-19 23:31:35
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

monitor2885

至尊木虫 (职业作家)

队长
★
fishfly80(金币+2):谢谢参与
上帝啊，这是什么啊？解析解？
赞一下(1人) 回复此楼
Retirement
2楼2009-09-20 13:43:24
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoshans

铁杆木虫 (小有名气)
★ ★
fishfly80(金币+2):谢谢参与
fishfly80(金币+1,VIP+0): 9-20 13:53
看不懂，楼主再说说吧
赞一下(2人) 回复此楼
3楼2009-09-20 13:44:43
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

fishfly80

铁虫 (小有名气)

就是pr求解方程

p=rt/(v-b)-a/(v^2+2bv-b^2)这就是pr方程，是高压状态下的气体密度，已知p r t b a  求v ，方程应该有三个根，可打入后求出的根却是上面那个样子，这是为什么，请高手赐教啊。另外p r t b a 数据也没错。
赞一下 回复此楼
4楼2009-09-20 13:52:45
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主
★ ★ ★ ★
fishfly80(金币+2):谢谢参与
fishfly80(金币+2,VIP+0):太谢谢了 正是我要的结果。 9-20 19:15
kuhailangyu(金币+1,VIP+0):欢迎积极参与 9-20 22:16
solve给出的符号型的解。从上面的解来看，是一个实根，两个虚根。在高压状态下的气体密度，只可能是个实根。不妨直接求解如下

v = fzero(@(v)(8.314* 77 / (v-0.00001649839676) - 0.025825/ (v^2+2*v*0.00001649839676 - 0.00001649839676^2) - 100000),1e-6);
v
赞一下(11人) 回复此楼
5楼2009-09-20 15:09:44
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

fishfly80

铁虫 (小有名气)

版主，可是不管压强如何改，计算结果都一样的啊

★ ★ ★
nono2009(金币+3,VIP+0):谢谢指出问题！ 9-20 20:28
v = fzero(@(v)(8.314* 77 / (v-0.00001649839676) - 0.025827/ (v^2+2*v*0.00001649839676 - 0.00001649839676^2) - 10000),1e-6);
v

v =

  6.8339e-006

>> v = fzero(@(v)(8.314* 77 / (v-0.00001649839676) - 0.025827/ (v^2+2*v*0.00001649839676 - 0.00001649839676^2) - 1000000),1e-6);
v

v =

  6.8339e-006

>> v = fzero(@(v)(8.314* 77 / (v-0.00001649839676) - 0.025827/ (v^2+2*v*0.00001649839676 - 0.00001649839676^2) - 100),1e-6);
v

v =

  6.8339e-006
赞一下(6人) 回复此楼
6楼2009-09-20 19:58:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
kuhailangyu(金币+1,VIP+0):欢迎积极参与 9-20 22:16
是数值问题，将函数修改了一下，再试试：

v = fzero(@(v)(8.314* 77e+5 / (v-1.649839676) - 0.025827e+10/ (v-1.649839676)^2 - 10000),1.0);
v*1.0e-5
赞一下(11人) 回复此楼
7楼2009-09-20 20:27:18
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

fishfly80

铁虫 (小有名气)

结果还是一样的啊。

nono2009(金币+0,VIP+0):上面我搞错了，对PR方程不熟。Sorry. 9-20 22:01
v = fzero(@(v)(8.314* 77e+5 / (v-1.649839676) - 0.025827e+10/ (v^2+2*v*1.649839676 - 1.649839676^2) - 1000),1.0);
v*1.0e-5

ans =

  6.8339e-006
赞一下(6人) 回复此楼
8楼2009-09-20 20:44:19
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

fishfly80

铁虫 (小有名气)

斑竹，能否推荐关于如何推荐关于matlab解方程的知识啊

不胜感激啊
回复此楼
9楼2009-09-20 20:46:40
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
kuhailangyu(金币+2,VIP+0):欢迎积极参与 9-20 22:16
让我们回到PR方程：p=rt/(v-b)-a/(v^2+2bv-b^2)
你能告诉我你是如何取各参数的单位的吗？
例如：你的压力只取了10000，如果单位是Pa的话，那根本不是高压，而是真空；t应该是K，如果取77K，那是极低的温度。所以我想确认一下参数的单位。
赞一下(11人) 回复此楼
10楼2009-09-20 21:59:33
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fishfly80 的主题更新
151/2返回列表12下一页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:48 by 544594351
[考研] 一志愿华北电力大学（北京），材料科学与工程学硕265，求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[考研] 材料调剂 +11 吴棂颖！ 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 288求调剂 一志愿哈工大 材料与化工 +39 洛神哥哥 2026-03-31 41/2050 2026-04-03 21:51 by qlm5820
[考研] 考研调剂 +3 Draa 2026-04-03 3/150 2026-04-03 17:37 by hgwz7468
[考研] 一志愿生物与医药，296分，求调剂 +8 66鹿 2026-04-03 9/450 2026-04-03 14:22 by 化学化工硕士招�
[考研] 085501一志愿天工大，机械专硕求调剂，跨材料 +3 33上 2026-04-03 3/150 2026-04-03 14:08 by 1753564080
[考研] 311求调剂一志愿合肥工业大学 +15 秋二十二 2026-03-30 15/750 2026-04-03 10:19 by linyelide
[考研] 279求调剂 +6 qazplm0852 2026-04-02 6/300 2026-04-03 10:03 by 蓝云思雨
[考研] 一志愿安徽大学0817化学工程与技术，求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +5 蓝笺片 2026-04-02 6/300 2026-04-02 21:37 by dongzh2009
[考研] 0710生物学求调剂 +9 manman511 2026-04-01 9/450 2026-04-02 10:00 by zxl830724
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 13/650 2026-04-02 09:17 by olim
[考研] 307分求调剂 +14 (o~o) 2026-03-31 15/750 2026-04-01 20:43 by longlotian
[考研] 310分求调剂 +4 成功上岸wang 2026-04-01 4/200 2026-04-01 20:35 by liu823948201
[考研] 288资源与环境专硕求调剂，不限专业，有学上就行 +25 lllllos 2026-03-30 26/1300 2026-04-01 09:52 by 一只好果子?
[考研] 考研生物与医药调剂 +7 铁憨憨123425 2026-03-31 7/350 2026-04-01 08:45 by JourneyLucky
[考研] 复试调剂 +7 双马尾痞老板2 2026-03-31 7/350 2026-03-31 19:49 by Dyhoer
[考研] 一志愿华中师范化学332分求调剂 +3 Lyy930824@ 2026-03-29 3/150 2026-03-30 20:15 by DHUSHUAI
[考研] 0703化学321分求调剂 +10 三dd. 2026-03-30 11/550 2026-03-30 19:24 by markhwc
信息提示
关闭
请填处理意见
关闭