wanyshenger

金蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 981.4
散金: 998
紅花: 7
帖子: 478
在線: 64.3小時(shí)
蟲號(hào): 2557404
注冊(cè): 2013-07-21
性別: GG
專業(yè): 高分子合成化學(xué)

[求助] 請(qǐng)教大神：matlab怎么求解有符號(hào)變量的方程組已有2人參與

MATLAB剛剛?cè)腴T，不太會(huì)用，老板讓我計(jì)算絕熱火焰溫度，現(xiàn)在遇到的問題是燃燒生成的CO2,H2O會(huì)分解，但化學(xué)平衡的方程組不會(huì)解，簡(jiǎn)單來說就是方程組里有符號(hào)變量如何求解
例如：
x*[c+0.5*(x+y)]^0.5/(a-x)=K1;
y*[c+0.5*(x+y)]^0.5/(b-y)=K2;
其中a,b,c,K1,K2 均為過程中賦值的變量，求解x，y，請(qǐng)問這種情況應(yīng)該如何用MATLAB求解。請(qǐng)各位大神不吝指教，感謝感謝！

回復(fù)此樓

» 猜你喜歡

085700資源與環(huán)境308求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)調(diào)劑招人�。�！已經(jīng)有4人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有17人回復(fù)
一志愿武漢理工材料工程專碩調(diào)劑已經(jīng)有5人回復(fù)
0856調(diào)劑，是學(xué)校就去已經(jīng)有7人回復(fù)
申博26年已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有10人回復(fù)
0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷已經(jīng)有12人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有13人回復(fù)
288求調(diào)劑已經(jīng)有14人回復(fù)

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

救助關(guān)于幾率分布、圖形擬合的問題已經(jīng)有1人回復(fù)

USAwaitforme

1樓 2015-07-30 04:02:14

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

斷劍殘雪

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 25 (小學(xué)生)
貴賓: 0.057
金幣: 8521.4
散金: 9306
紅花: 191
沙發(fā): 5
帖子: 3566
在線: 1232.8小時(shí)
蟲號(hào): 2187840
注冊(cè): 2012-12-15
專業(yè): 代數(shù)學(xué)

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

感謝參與，應(yīng)助指數(shù) +1

syms x y a b c K1 K2; %定義變量
[x,y]=solve(x*sqrt(c+0.5*(x+y))/(a-x)==K1 , y*sqrt(c+0.5*(x+y))/(b-y)==K2,x,y) %給定方程組，并返回x，y

贊一下

回復(fù)此樓

靜若處子

2樓2015-07-30 08:14:25

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

斷劍殘雪

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 25 (小學(xué)生)
貴賓: 0.057
金幣: 8521.4
散金: 9306
紅花: 191
沙發(fā): 5
帖子: 3566
在線: 1232.8小時(shí)
蟲號(hào): 2187840
注冊(cè): 2012-12-15
專業(yè): 代數(shù)學(xué)

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

output：

x =

(K1*a + K2*b + 2.0*K2*c)/(K1 - 1.0*K2) - (2.0*z1^3)/(K1 - 1.0*K2) + (2.0*c*z1)/(K1 - 1.0*K2) - (2.0*K2*z1^2)/(K1 - 1.0*K2)

y =

(2.0*z1^3)/(K1 - 1.0*K2) + (2.0*K1*z1^2)/(K1 - 1.0*K2) - (2.0*c*z1)/(K1 - 1.0*K2) - (1.0*K1*a + 1.0*K2*b + 2.0*K1*c)/(K1 - 1.0*K2)

贊一下

回復(fù)此樓

靜若處子

3樓2015-07-30 08:15:49

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

斷劍殘雪

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 25 (小學(xué)生)
貴賓: 0.057
金幣: 8521.4
散金: 9306
紅花: 191
沙發(fā): 5
帖子: 3566
在線: 1232.8小時(shí)
蟲號(hào): 2187840
注冊(cè): 2012-12-15
專業(yè): 代數(shù)學(xué)

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

z1 = RootOf(z^4 + z^3*(K1 + K2) - z^2*(c - K1*K2) - z*((K1*a)/2 + (K2*b)/2 + K1*c + K2*c) - K1*K2*c - (K1*K2*b)/2 - (K1*K2*a)/2, z)

贊一下

回復(fù)此樓

靜若處子

4樓2015-07-30 08:17:42

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

斷劍殘雪

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 25 (小學(xué)生)
貴賓: 0.057
金幣: 8521.4
散金: 9306
紅花: 191
沙發(fā): 5
帖子: 3566
在線: 1232.8小時(shí)
蟲號(hào): 2187840
注冊(cè): 2012-12-15
專業(yè): 代數(shù)學(xué)

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

★ ★ ★ ★ ★ ★ ★ ★
wanyshenger: 金幣+8, ★★★很有幫助 2015-07-30 13:00:56

最終版：
syms z a b c K1 K2 x y;
a=1 ;b=1 ;c=1 ;K1=1 ;K2=1; %比如已知它們的值，要求x，y。那么把前面的注釋去掉再運(yùn)行即可，它們的值也可以自己設(shè)定。
[x,y]=solve(x*sqrt(c+0.5*(x+y))/(a-x)==K1 , y*sqrt(c+0.5*(x+y))/(b-y)==K2,x,y); %給定方程組，并返回x，y的符號(hào)表達(dá)式;
z1 = solve(z^4 + z^3*(K1 + K2) - z^2*(c - K1*K2) - z*((K1*a)/2 + (K2*b)/2 + K1*c + K2*c) - K1*K2*c - (K1*K2*b)/2 - (K1*K2*a)/2==0, z)
x=subs(x,'z1',z1) %返回x的值
y=subs(y,'z1',z1) %返回y的值

贊一下

回復(fù)此樓

靜若處子

5樓2015-07-30 08:34:25

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

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
wanyshenger: 金幣+2 2015-07-30 13:01:34

不存在解析解的。
如果給出a,b,c,K1,K2 的具體數(shù)值，可以求出數(shù)值解。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

6樓2015-07-30 12:20:01

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

wanyshenger

金蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 981.4
散金: 998
紅花: 7
帖子: 478
在線: 64.3小時(shí)
蟲號(hào): 2557404
注冊(cè): 2013-07-21
性別: GG
專業(yè): 高分子合成化學(xué)

引用回帖:

6樓: Originally posted by 月只藍(lán) at 2015-07-30 12:20:01
不存在解析解的。
如果給出a,b,c,K1,K2 的具體數(shù)值，可以求出數(shù)值解。

請(qǐng)問：我會(huì)在過程中給a b c K1 K2賦值，不過因?yàn)橛袀€(gè)循環(huán)，所以這五個(gè)的值是會(huì)變化的，所以我想問的就是怎么在某一個(gè)循環(huán)下（就是a b c K1 K2被賦值），求解這步循環(huán)所應(yīng)該得到的x和y。謝謝！

贊一下

回復(fù)此樓

USAwaitforme

7樓2015-07-30 12:46:50

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

wanyshenger

金蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 981.4
散金: 998
紅花: 7
帖子: 478
在線: 64.3小時(shí)
蟲號(hào): 2557404
注冊(cè): 2013-07-21
性別: GG
專業(yè): 高分子合成化學(xué)

引用回帖:

5樓: Originally posted by 斷劍殘雪 at 2015-07-30 08:34:25
最終版：
syms z a b c K1 K2 x y;
a=1 ;b=1 ;c=1 ;K1=1 ;K2=1; %比如已知它們的值，要求x，y。那么把前面的注釋去掉再運(yùn)行即可，它們的值也可以自己設(shè)定。
=solve(x*sqrt(c+0.5*(x+y))/(a-x)==K1 , y*sqrt(c ...

謝謝你的幫忙，還想請(qǐng)問一下，如果我這個(gè)方程更加復(fù)雜，比如
x*[c+0.5*(x+y)]^0.5/(a-x)/(a+b+c+0.5*(x+y))^0.5=K1;
y*[c+0.5*(x+y)]^0.5/(b-y)/(a+b+c+0.5*(x+y))^0.5=K2;
這樣的話也可以這樣做嗎？這樣會(huì)不會(huì)效率不高呀，因?yàn)槲疫€有循環(huán)（每一步循環(huán)中a b c K1 K2會(huì)變），這樣會(huì)不會(huì)計(jì)算很慢呀，謝謝

贊一下

回復(fù)此樓

USAwaitforme

8樓2015-07-30 13:00:43

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

斷劍殘雪

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 25 (小學(xué)生)
貴賓: 0.057
金幣: 8521.4
散金: 9306
紅花: 191
沙發(fā): 5
帖子: 3566
在線: 1232.8小時(shí)
蟲號(hào): 2187840
注冊(cè): 2012-12-15
專業(yè): 代數(shù)學(xué)

引用回帖:

8樓: Originally posted by wanyshenger at 2015-07-30 13:00:43
謝謝你的幫忙，還想請(qǐng)問一下，如果我這個(gè)方程更加復(fù)雜，比如
x*^0.5/(a-x)/(a+b+c+0.5*(x+y))^0.5=K1;
y*^0.5/(b-y)/(a+b+c+0.5*(x+y))^0.5=K2;
這樣的話也可以這樣做嗎？這樣會(huì)不會(huì)效率不高呀，因?yàn)槲疫€有循環(huán) ...

給定矩陣，每一次讓它自己提取即可。

贊一下

回復(fù)此樓

靜若處子

9樓2015-07-30 13:08:44

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

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

引用回帖:

7樓: Originally posted by wanyshenger at 2015-07-30 12:46:50
請(qǐng)問：我會(huì)在過程中給a b c K1 K2賦值，不過因?yàn)橛袀€(gè)循環(huán)，所以這五個(gè)的值是會(huì)變化的，所以我想問的就是怎么在某一個(gè)循環(huán)下（就是a b c K1 K2被賦值），求解這步循環(huán)所應(yīng)該得到的x和y。謝謝！...

方程求解的程序放進(jìn)循環(huán)中，a b c K1 K2變化一次，求解一次方程。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

10樓2015-07-30 13:16:29

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

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

[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4	脫穎而出 2026-03-16	12/600	2026-03-19 16:17 by 脫穎而出
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 266求調(diào)劑 +5	陽陽哇塞 2026-03-14	10/500	2026-03-19 15:08 by 陽陽哇塞
[考研] 317求調(diào)劑 +3	申子申申 2026-03-19	6/300	2026-03-19 14:16 by 申子申申
[考研] 求調(diào)劑 +3	Mqqqqqq 2026-03-19	3/150	2026-03-19 14:11 by peike
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無懈可擊111
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 299求調(diào)劑 +5	△小透明* 2026-03-17	5/250	2026-03-18 11:49 by 盡舜堯1
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu