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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】急，还是matlab 解方程求助！【已完成】
151/2返回列表12下一页
查看: 930  |  回复: 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下一页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 F4zcb7xRCEJA 2026-08-11 3/150 2026-08-11 21:51 by pMlf9MbdMHyc
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 yPv6ndDtiRxQ 2026-08-11 3/150 2026-08-11 21:09 by thgx1vqMqrLw
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 yPv6ndDtiRxQ 2026-08-11 4/200 2026-08-11 21:03 by thgx1vqMqrLw
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 yPv6ndDtiRxQ 2026-08-11 3/150 2026-08-11 20:51 by thgx1vqMqrLw
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 GpvwGNVshCOp 2026-08-10 7/350 2026-08-11 19:39 by thgx1vqMqrLw
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] FileCode能看出啥？ +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
关闭
请填处理意见
关闭