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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】急，还是matlab 解方程求助！【已完成】
北京石油化工学院2026年研究生招生接收调剂公告
151/2返回列表12下一页
查看: 810  |  回复: 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下一页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +11 一样YWY 2026-04-04 11/550 2026-04-04 23:58 by 永字号
[考研] 285求调剂 +9 哦呦呼o 2026-04-04 9/450 2026-04-04 23:38 by lqwchd
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-04 20:37 by 蓝云思雨
[考研] 求调剂 +3 农业工程与信息� 2026-04-04 3/150 2026-04-04 12:19 by 舍而后得
[考研] 357求调剂 +13 1050389037 2026-04-03 13/650 2026-04-03 22:27 by 无际的草原
[考研] 材料专硕调剂 +18 椰椰。 2026-03-29 18/900 2026-04-03 16:45 by 玲玲0606
[考研] 机械专硕297 +3 Afksy 2026-04-03 3/150 2026-04-03 14:24 by 1753564080
[考研] 0705理学294求调剂 +3 成果成果cg5 2026-04-03 3/150 2026-04-03 14:04 by simons1972
[考研] 321求调剂 +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[考研] 315分 085602 求调剂 +15 26考研上岸版26 2026-04-02 15/750 2026-04-03 12:45 by xingguangj
[考研] 285求调剂 +7 AZMK 2026-04-02 9/450 2026-04-03 11:12 by wanwan00
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
[考研] 生物学308分求调剂（一志愿华东师大） +6 相信必会光芒万� 2026-03-31 7/350 2026-04-02 23:16 by JourneyLucky
[考研] 318求调剂 +3 笃行致远. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[考研] 377求调剂 +3 RASKIN 2026-04-02 3/150 2026-04-02 09:45 by zzchen2000
[考博] 26年申博 +3 staryer 2026-03-30 4/200 2026-04-01 23:21 by ai4pharm
[考研] 085600 一志愿9 总分351 求调剂学校 +7 czhcz 2026-03-31 9/450 2026-04-01 19:24 by 无际的草原
[考研] 350求调剂 +7 阿佳～ 2026-03-31 7/350 2026-04-01 16:12 by yanflower7133
[考研] 一志愿大连理工大学，机械工程学硕，341 +3 西瓜田的守望者 2026-03-30 3/150 2026-03-31 11:08 by asdfzly
[考研] 本科211总分289，08工学真心求调剂 +3 utopiaE 2026-03-30 3/150 2026-03-30 23:42 by ms629
信息提示
关闭
请填处理意见
关闭