jrflower

木蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 3767.3
散金: 50
紅花: 4
帖子: 597
在線: 230.8小時
蟲號: 414707
注冊: 2007-06-28
性別: MM
專業(yè): 金屬材料的合金相、相變及

[求助] 求matlab求解一個方程組的問題

想用matlab求解一個包含3個式子的方程組，結果沒算出來，不知道是怎么回事，我對matlab了解非常少，想著應該是我寫的式子有問題，請有經(jīng)驗的朋友幫幫忙
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的值

回復此樓

» 猜你喜歡

一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有3人回復
290求調(diào)劑已經(jīng)有11人回復
材料與化工求調(diào)劑已經(jīng)有4人回復
293求調(diào)劑已經(jīng)有11人回復
工科材料085601 279求調(diào)劑已經(jīng)有3人回復
301求調(diào)劑已經(jīng)有3人回復
材料專碩306英一數(shù)二已經(jīng)有8人回復
材料與化工專碩調(diào)劑已經(jīng)有5人回復
085600材料與化工已經(jīng)有4人回復
085601材料工程專碩求調(diào)劑已經(jīng)有3人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

matlab數(shù)值求解邊界條件微分方程組已經(jīng)有7人回復
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復
【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題已經(jīng)有21人回復
【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(shù)解！已經(jīng)有18人回復
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復
【求助】用fortran求解大型線性方程組時出現(xiàn)的錯誤【已解決】已經(jīng)有11人回復
【求助】求MAtlab求解程序�。�！已經(jīng)有20人回復

1樓 2011-11-17 21:32:46

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
程序強帖: 2
應助: 289 (大學生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學
管轄: 計算模擬

【答案】應助回帖

★
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

贊一下(1人)

回復此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2樓2011-11-17 22:12:06

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

jrflower

木蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 3767.3
散金: 50
紅花: 4
帖子: 597
在線: 230.8小時
蟲號: 414707
注冊: 2007-06-28
性別: MM
專業(yè): 金屬材料的合金相、相變及

引用回帖:

2樓: Originally posted by dbb627 at 2011-11-17 15:12:06:
[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.083464305892 ...

十分感謝，這個結果肯定是對的�？墒菫槭裁次野涯愕氖阶訌椭频絤atlab里面就算不出來呢，得到下面的結果
>> [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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
程序強帖: 2
應助: 289 (大學生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學
管轄: 計算模擬

【答案】應助回帖

★
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')

贊一下(1人)

回復此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

4樓2011-11-17 22:39:41

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

jrflower

木蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 3767.3
散金: 50
紅花: 4
帖子: 597
在線: 230.8小時
蟲號: 414707
注冊: 2007-06-28
性別: MM
專業(yè): 金屬材料的合金相、相變及

引用回帖:

4樓: Originally posted by dbb627 at 2011-11-17 15:39:41:
可能是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 s ...

我的版本很老，是matlab7.0,難道需要在開始把a,b,c設置一下類型，比如syms a b c, real a b c,什么的，[ empty sym ]是不是當成一個矩陣什么了，不知道怎么解決，另外看網(wǎng)上還有一個fsolve和fzero的命令，也是求解方程的，不知道是不是該用這兩個命令

贊一下

回復此樓

5樓2011-11-17 23:48:43

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
程序強帖: 2
應助: 289 (大學生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學
管轄: 計算模擬

★
jjdg(金幣+1): 感謝參與 2011-11-18 01:06:49

[ empty sym ]是空的符號矩陣

fsolve和fzero的命令是求數(shù)值解的
你可以用help fsolve查看用法

f=@(x)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;
fsolve(f,0)
解出來的結果也不對
估計因為你那方程就沒有解

贊一下(1人)

回復此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

6樓2011-11-17 23:56:55

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

jrflower

木蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 3767.3
散金: 50
紅花: 4
帖子: 597
在線: 230.8小時
蟲號: 414707
注冊: 2007-06-28
性別: MM
專業(yè): 金屬材料的合金相、相變及

引用回帖:

6樓: Originally posted by dbb627 at 2011-11-17 16:56:55:
[ empty sym ]是空的符號矩陣

fsolve和fzero的命令是求數(shù)值解的
你可以用help fsolve查看用法

f=@(x)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) ...

第二個方程是我搞錯了，里面少了括號，應該是
>> 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

贊一下(1人)

回復此樓

7樓2011-11-18 00:12:30

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 材料與化工求調(diào)劑 +4	為學666 2026-03-16	4/200	2026-03-17 14:28 by houyaoxu
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +3	小李0854 2026-03-16	3/150	2026-03-17 13:40 by 熱情沙漠
[考研] 304求調(diào)劑 +7	小熊joy 2026-03-14	7/350	2026-03-17 08:53 by 霧散后相遇lc
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 0703化學調(diào)劑，求各位老師收留 +8	秋有木北 2026-03-14	8/400	2026-03-16 15:21 by 哦哦123
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 344求調(diào)劑 +3	knight344 2026-03-16	3/150	2026-03-16 09:42 by 無際的草原
[考研] 機械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 288求調(diào)劑 +4	奇點0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學習一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 材料與化工085600調(diào)劑求老師收留 +9	jiaanl 2026-03-11	9/450	2026-03-13 20:22 by JourneyLucky
[考研] 求調(diào)劑 +5	一定有學上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[論文投稿] 投稿問題 5+4	星光燦爛xt 2026-03-12	6/300	2026-03-13 14:17 by god_tian
[考研] 0856化工原理 +6	z2839474511 2026-03-10	6/300	2026-03-13 10:41 by houyaoxu
[考研] 求調(diào)劑資源與環(huán)境 285 +3	未名考生 2026-03-10	3/150	2026-03-13 10:31 by houyaoxu
[考研] 420求調(diào)劑 +4	莫向外求11 2026-03-10	6/300	2026-03-12 14:41 by ruiyingmiao
[考研] 大連大學化學專業(yè)研究生調(diào)劑 +3	琪久. 2026-03-10	8/400	2026-03-11 10:02 by 琪久.