24小时热门版块排行榜    

查看: 926  |  回复: 14
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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 ]
回复此楼
已阅   关注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楼2009-09-20 19:58:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

monitor2885

至尊木虫 (知名作家)

队长


fishfly80(金币+2):谢谢参与
上帝啊,这是什么啊?解析解?
Retirement
2楼2009-09-20 13:43:24
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoshans

铁杆木虫 (小有名气)

★ ★
fishfly80(金币+2):谢谢参与
fishfly80(金币+1,VIP+0): 9-20 13:53
看不懂,楼主再说说吧
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode +8 documentary 2026-08-10 8/400 2026-08-10 19:20 by loufangrui
[基金申请] 帮忙看看fileCode +3 wwncly 2026-08-10 5/250 2026-08-10 19:06 by 2000zf36392
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 11/550 2026-08-10 17:09 by 六两废铜
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +14 Lanmanbaby 2026-08-09 22/1100 2026-08-10 14:57 by wwncly
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 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
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
请填处理意见