| 查看: 2081 | 回复: 6 | |||||
jrflower木虫 (正式写手)
|
[求助]
求matlab求解一个方程组的问题
|
|
想用matlab求解一个包含3个式子的方程组，结果没算出来，不知道是怎么回事，我对matlab了解非常少，想着应该是我写的式子有问题，请有经验的朋友帮帮忙 solve('1/4596.75-8.31446/716700*log(a)=1/4435.46-8.31446/837700*log(b)', '1/4435.46-8.31446/837700*log(b)=1/4097.94-8.31446/820000*log(c)', 'a+b+c=1') 要求出a,b,c的值 |
回复此楼» 猜你喜欢
A期刊撤稿
已经有3人回复
-
职称评审没过，求安慰
已经有34人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
回收溶剂求助
已经有6人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab数值求解边界条件微分方程组
已经有7人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】用fortran求解大型线性方程组时出现的错误【已解决】
已经有11人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 程序强帖: 2
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★
jrflower(金币+20): 谢谢 2011-11-17 22:29:59
余泽成(金币+1): 谢谢参与应助！ 2011-11-19 14:27:53
jrflower(金币+20): 谢谢 2011-11-17 22:29:59
余泽成(金币+1): 谢谢参与应助！ 2011-11-19 14:27:53
|
[a,b,c]=solve('1/4596.75-8.31446/716700*log(a)=1/4435.46-8.31446/837700*log(b)', '1/4435.46-8.31446/837700*log(b)=1/4097.94-8.31446/820000*log(c)', 'a+b+c=1','a','b','c') a = 0.08346430589202496974797491804609 b = 0.12177617993836521593616873649657 c = 0.79475951416960981431585634545734 |
2楼2011-11-17 22:12:06
jrflower
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3767.3
- 散金: 50
- 红花: 4
- 帖子: 597
- 在线: 230.8小时
- 虫号: 414707
- 注册: 2007-06-28
- 性别: MM
- 专业: 金属材料的合金相、相变及
|
十分感谢，这个结果肯定是对的。可是为什么我把你的式子复制到matlab里面就算不出来呢，得到下面的结果 >> [a,b,c]=solve('1/4596.75-8.31446/716700*log(a)=1/4435.46-8.31446/837700*log(b)', '1/4435.46-8.31446/837700*log(b)=1/4097.94-8.31446/820000*log(c)', 'a+b+c=1','a','b','c') Warning: Explicit solution could not be found. > In solve at 140 a = [ empty sym ] b = [] c = [] 请帮我看看这是怎么回事，另外我求解这个方程 x=solve('exp(1/8.31446*(1/4596.75-1/4435.46)+1/837700*log(x))*716700+exp(1/8.31446*(1/4097.94-1/4435.46)+1/837700*log(x))*820000+x-1')也求解不出来，也请帮我看看是什么问题吧 |
3楼2011-11-17 22:29:02
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 程序强帖: 2
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★
jjdg(金币+1): 感谢参与 2011-11-18 01:06:58
jjdg(金币+1): 感谢参与 2011-11-18 01:06:58
|
可能是matlab版本的问题 我的是2010b x=solve('exp(1/8.31446*(1/4596.75-1/4435.46)+1/837700*log(x))*716700+exp(1/8.31446*(1/4097.94-1/4435.46)+1/837700*log(x))*820000+x-1','x') Warning: Explicit solution could not be found. > In solve at 81 x = [ empty sym ] 画图可知 没有解满足条件 ezplot('exp(1/8.31446*(1/4596.75-1/4435.46)+1/837700*log(x))*716700+exp(1/8.31446*(1/4097.94-1/4435.46)+1/837700*log(x))*820000+x-1') |
4楼2011-11-17 22:39:41
jrflower
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3767.3
- 散金: 50
- 红花: 4
- 帖子: 597
- 在线: 230.8小时
- 虫号: 414707
- 注册: 2007-06-28
- 性别: MM
- 专业: 金属材料的合金相、相变及
5楼2011-11-17 23:48:43
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 程序强帖: 2
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
6楼2011-11-17 23:56:55
jrflower
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3767.3
- 散金: 50
- 红花: 4
- 帖子: 597
- 在线: 230.8小时
- 虫号: 414707
- 注册: 2007-06-28
- 性别: MM
- 专业: 金属材料的合金相、相变及
|
第二个方程是我搞错了，里面少了括号，应该是 >> x=fsolve('exp((1/8.31446*(1/4596.75-1/4435.46)+1/837700*log(x))*716700)+exp((1/8.31446*(1/4097.94-1/4435.46)+1/837700*log(x))*820000)+x-1',0.13) Optimization terminated: first-order optimality is less than options.TolFun. x = 0.1218 解出来和第一个方程组你帮我解出来的b值相等，这个式子就是我把前面的方程组合并了一下求的b |
7楼2011-11-18 00:12:30