jrflower

木蟲 (正式寫手)

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

[求助] 求matlab求解一個(gè)方程組的問題

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

回復(fù)此樓

» 猜你喜歡

0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷已經(jīng)有8人回復(fù)
268求調(diào)劑已經(jīng)有8人回復(fù)
材料，紡織，生物（0856、0710），化學(xué)招生啦已經(jīng)有7人回復(fù)
301求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有7人回復(fù)
材料與化工求調(diào)劑已經(jīng)有6人回復(fù)
299求調(diào)劑已經(jīng)有4人回復(fù)
277調(diào)劑已經(jīng)有6人回復(fù)
341求調(diào)劑已經(jīng)有7人回復(fù)
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷已經(jīng)有8人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

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

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

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

dbb627

榮譽(yù)版主 (著名寫手)

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

【答案】應(yīng)助回帖

★
jrflower(金幣+20): 謝謝 2011-11-17 22:29:59
余澤成(金幣+1): 謝謝參與應(yīng)助！ 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人)

回復(fù)此樓

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

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

jrflower

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3767.3
散金: 50
紅花: 4
帖子: 597
在線: 230.8小時(shí)
蟲號(hào): 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 ...

十分感謝，這個(gè)結(jié)果肯定是對(duì)的。可是為什么我把你的式子復(fù)制到matlab里面就算不出來呢，得到下面的結(jié)果
>> [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 =

[]
請幫我看看這是怎么回事，另外我求解這個(gè)方程
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')也求解不出來，也請幫我看看是什么問題吧

贊一下

回復(fù)此樓

3樓2011-11-17 22:29:02

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

dbb627

榮譽(yù)版主 (著名寫手)

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

【答案】應(yīng)助回帖

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

回復(fù)此樓

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

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

jrflower

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3767.3
散金: 50
紅花: 4
帖子: 597
在線: 230.8小時(shí)
蟲號(hào): 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ǔ),b,c設(shè)置一下類型，比如syms a b c, real a b c,什么的，[ empty sym ]是不是當(dāng)成一個(gè)矩陣什么了，不知道怎么解決，另外看網(wǎng)上還有一個(gè)fsolve和fzero的命令，也是求解方程的，不知道是不是該用這兩個(gè)命令

贊一下

回復(fù)此樓

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

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

dbb627

榮譽(yù)版主 (著名寫手)

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

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

[ empty sym ]是空的符號(hào)矩陣

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)
解出來的結(jié)果也不對(duì)
估計(jì)因?yàn)槟隳欠匠叹蜎]有解

贊一下(1人)

回復(fù)此樓

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

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

jrflower

木蟲 (正式寫手)

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

引用回帖:

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

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) ...

第二個(gè)方程是我搞錯(cuò)了，里面少了括號(hào)，應(yīng)該是
>> 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
解出來和第一個(gè)方程組你幫我解出來的b值相等，這個(gè)式子就是我把前面的方程組合并了一下求的b

贊一下(1人)

回復(fù)此樓

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

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

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

[考研] 293求調(diào)劑 +7	zjl的號(hào) 2026-03-16	12/600	2026-03-17 18:22 by 重科小霸王
[考研] 289求調(diào)劑 +6	步川酷紫123 2026-03-11	6/300	2026-03-17 10:23 by Sammy2
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5	薛云鵬 2026-03-13	5/250	2026-03-17 10:15 by Sammy2
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識(shí)可樂 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 學(xué)碩285求調(diào)劑 +13	Wisjxn 2026-03-12	46/2300	2026-03-14 10:33 by JourneyLucky
[考研] 材料371求調(diào)劑 +9	鱷魚? 2026-03-11	11/550	2026-03-13 22:53 by JourneyLucky
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5	瑪卡巴卡啊哈 2026-03-11	5/250	2026-03-13 21:40 by JourneyLucky
[考研] 【考研調(diào)劑求收留】 +3	Ceciilia 2026-03-11	3/150	2026-03-13 20:18 by JourneyLucky
[考研] 材料工程調(diào)劑 +4	咪咪空空 2026-03-11	4/200	2026-03-13 19:57 by JourneyLucky
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 274求調(diào)劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz