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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】急，还是matlab 解方程求助！【已完成】
151/2返回列表12下一页
查看: 970  |  回复: 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下一页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-31 00:21 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:30 by hZiFeudZyoGR
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:21 by hZiFeudZyoGR
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jCd0dEvKHShX 2026-08-29 3/150 2026-08-30 19:26 by hZiFeudZyoGR
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 19:01 by hZiFeudZyoGR
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:50 by hZiFeudZyoGR
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:40 by hZiFeudZyoGR
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 17:34 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:03 by hZiFeudZyoGR
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 18/900 2026-08-30 16:47 by 黄鸟于飞Chao
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:19 by l0VvVHGBGRLv
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
信息提示
关闭
请填处理意见
关闭