小木蟲zb

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時
蟲號: 1227944
注冊: 2011-03-10
性別: GG
專業(yè): 分離過程

[求助] 求助matlab解一個非線性方程組

四個方程如圖，其中qm1=1.13788,qm2=1.12991,
k1=0.58955,k2=0.54607,c0=0.5,1,2,4,6,8,10,15,
20,25,30這幾個點，分別求這幾個點上的c1,
q1,c2,q2,然后作圖(橫坐標為c,縱坐標為q),我
寫的matlab程序如上傳文件binaryequlibrum.m中所
示，但提示有錯，求大神指教
方程組.jpg

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : binaryequlibrum.m

2013-01-04 11:12:37, 736 bytes

» 猜你喜歡

中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有154人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有2人回復(fù)

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

求Matlab高手解決線性方程組的迭代求解問題已經(jīng)有17人回復(fù)
matlab求解非線性方程組，錯誤提示怎么解決已經(jīng)有5人回復(fù)
求助 MATLAB解方程組-fslove 已經(jīng)有7人回復(fù)
求助matlab---fsolve解非線性方程組已經(jīng)有6人回復(fù)
求助！matlab用fsolve函數(shù)求解非線性方程組的問題！已經(jīng)有19人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
matlab求解非線性方程組已經(jīng)有16人回復(fù)
高金求助matlab解微分方程組已經(jīng)有12人回復(fù)
matlab解微分方程組已經(jīng)有15人回復(fù)
求matlab求解一個方程組的問題已經(jīng)有6人回復(fù)
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題已經(jīng)有21人回復(fù)
【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(shù)解！已經(jīng)有18人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】如何求解下列偏微分方程已經(jīng)有13人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)
【求助】求MAtlab求解程序�。。� 已經(jīng)有20人回復(fù)
【求助】100金幣求助高人一個程序解方程已經(jīng)有15人回復(fù)
【求助】求教matlab解非線性方程組已經(jīng)有9人回復(fù)
【求助】一個時變微分方程組的matlab實現(xiàn) 已經(jīng)有8人回復(fù)

1樓 2013-01-04 11:12:46

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

小木蟲zb

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時
蟲號: 1227944
注冊: 2011-03-10
性別: GG
專業(yè): 分離過程

A=importdata 部分可以刪除，那是我用來和其他數(shù)據(jù)比較

贊一下

回復(fù)此樓

2樓2013-01-04 11:18:15

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

dingd

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

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
小木蟲zb: 金幣+10, ★有幫助 2013-01-04 12:00:01
csgt0: 金幣+2, 謝謝 2013-01-05 10:01:47

用1stOpt跑了下，供參考：

CODE:

LoopConstant c0=[0.5,1,2,4,6,8,10,15,20,25,30];

Constant qm=[1.13788,1.12991],

         k=[0.58955,0.54607];

PlotLoopData c1[x],q1;

Function

q1=qm[1]*k[1]*c1/(1+k[1]*c1+k[2]*c2);

q2=qm[2]*k[2]*c2/(1+k[1]*c1+k[2]*c2);

q1=(c0-c1)*10/0.1;

q2=(c0-c2)*10/0.1;

循環(huán)常數(shù) c0 q1 c1 c2 q2
0.5 0.213344987124234 0.497866550128758 0.498037063635551 0.196293636444935
1 0.31363598030533 0.996863640196947 0.997114571953084 0.288542804691647
2 0.409862629704777 1.99590137370295 1.99622962277441 0.377037722559132
4 0.484017906040076 3.9951598209396 3.99554775905092 0.445224094907642
6 0.515050983766089 5.99484949016234 5.99526243482094 0.473756517905549
8 0.532101170278375 7.99467898829722 7.99510567894046 0.489432105953705
10 0.542881105828988 9.99457118894171 9.99500657297151 0.499342702848826
15 0.557948647226742 14.9944205135277 14.9948680519249 0.513194807506159
20 0.565798599049311 19.9943420140095 19.9947958861554 0.520411384458806
25 0.57061485553382 24.9942938514447 24.9947516101101 0.524838988988562
30 0.573871227314218 29.9942612877269 29.9947216743709 0.527832562908426

v1.JPG

贊一下(1人)

回復(fù)此樓

3樓2013-01-04 11:53:42

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

小木蟲zb

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時
蟲號: 1227944
注冊: 2011-03-10
性別: GG
專業(yè): 分離過程

引用回帖:

3樓: Originally posted by dingd at 2013-01-04 11:53:42
用1stOpt跑了下，供參考：

LoopConstant c0=;
Constant qm=,
k=;
PlotLoopData c1,q1;
Function
q1=qm*k*c1/(1+k*c1+k*c2);
q2=qm*k*c2/(1+k*c1+k*c2);
q1=(c0-c1)*10/0.1;
q2=(c0-c2)*10/0.1; ...

這個軟件沒用過，能給個matlab的編程嗎

贊一下

回復(fù)此樓

4樓2013-01-04 12:00:39

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
小木蟲zb: 金幣+10, ★★★★★最佳答案 2013-01-05 10:45:41

CODE:

function binaryequlibrum

clear all;clc

syms q1 q2 c1 c2;

qm1=1.13788;

qm2=1.12991;%DBT,IDL單組份擬合Langmuir曲線結(jié)果

k1=0.58955;

k2=0.54607;%DBT,IDL單組份擬合Langmuir曲線結(jié)果

c0=[0.5,1,2,4,6,8,10,15,20,25,30];

[c1,c2,q1,q2]=solve('q1=qm1*k1*c1/(1+k1*c1+k2*c2)',...

'q2=qm2*k2*c2/(1+k1*c1+k2*c2)','q1=(c0-c1)*10/0.1',...

'q2=(c0-c2)*10/0.1','c1','c2','q1','q2');

c11=subs(c1(3));

c22=subs(c2(3));

q11=subs(q1(3));

q22=subs(q2(3));

plot(c11,q11,'r-',c22,q22,'b-');

xlabel('c/mol/m3'),ylabel('q/mol/kg'),

legend('DBT-ExtendLangmuir','IDL-ExtendLangmuir','Location','best') ,title('雙組份DBT+IDL平衡曲線')

grid

贊一下(1人)

回復(fù)此樓

showmethemoney

5樓2013-01-05 10:33:07

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

小木蟲zb

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時
蟲號: 1227944
注冊: 2011-03-10
性別: GG
專業(yè): 分離過程

引用回帖:

5樓: Originally posted by csgt0 at 2013-01-05 10:33:07
function binaryequlibrum
clear all;clc
syms q1 q2 c1 c2;
qm1=1.13788;
qm2=1.12991;%DBT,IDL單組份擬合Langmuir曲線結(jié)果
k1=0.58955;
k2=0.54607;%DBT,IDL單組份擬合Langmuir曲線結(jié)果
c0=;
=solve('q ...

非常感謝，問題解決了，但我有點疑問，為什么是c11=subs(c1(3))？c1(3)指的是什么？

贊一下

回復(fù)此樓

6樓2013-01-05 10:47:22

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

小木蟲zb

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時
蟲號: 1227944
注冊: 2011-03-10
性別: GG
專業(yè): 分離過程

引用回帖:

非常感謝，問題解決了，但我有點疑問，為什么是c11=subs(c1(3))？c1(3)指的是什么？

贊一下

回復(fù)此樓

7樓2013-01-05 10:47:58

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

引用回帖:

7樓: Originally posted by 小木蟲zb at 2013-01-05 10:47:58
非常感謝，問題解決了，但我有點疑問，為什么是c11=subs(c1(3))？c1(3)指的是什么？...

有3組解，前兩組是負的，我沒要。你可以把3組分別畫出來看看

贊一下

回復(fù)此樓

showmethemoney

8樓2013-01-05 10:49:58

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

小木蟲zb

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時
蟲號: 1227944
注冊: 2011-03-10
性別: GG
專業(yè): 分離過程

引用回帖:

8樓: Originally posted by csgt0 at 2013-01-05 10:49:58
有3組解，前兩組是負的，我沒要。你可以把3組分別畫出來看看...

哦，我明白了，謝謝大神

贊一下

回復(fù)此樓

9樓2013-01-05 11:14:38

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

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

[考研] 材料專碩 085600求調(diào)劑 +7	BBQ233 2026-03-30	7/350	2026-03-30 17:44 by oooqiao
[考研] 279求調(diào)劑 +8	j的立方 2026-03-29	8/400	2026-03-30 17:21 by 無際的草原
[考研] 303求調(diào)劑 +4	DLkz1314. 2026-03-30	4/200	2026-03-30 16:18 by JourneyLucky
[考研] 求調(diào)劑 +7	青春裁為三截 2026-03-29	7/350	2026-03-30 13:14 by laoshidan
[考研] 282求調(diào)劑 +4	wcq131415 2026-03-24	4/200	2026-03-30 10:39 by Delta2012
[考研] 求調(diào)劑 +10	張zz111 2026-03-27	11/550	2026-03-30 09:17 by 無際的草原
[考研] 290求調(diào)劑 +3	dfffsar 2026-03-29	3/150	2026-03-29 22:38 by 毛毛毛阿莫2
[考研] 材料科學(xué)與工程求調(diào)劑 +6	深V宿舍吧 2026-03-29	6/300	2026-03-29 20:51 by 唐沐兒
[考研] 【求調(diào)劑】085601材料工程專碩 \| 總分272 \| +7	腳滑的守法公民 2026-03-27	7/350	2026-03-29 20:21 by dophin1985
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +4	12.27 2026-03-27	4/200	2026-03-29 14:40 by zhshch
[考研] 22408 359分調(diào)劑 +4	Qshers 2026-03-27	6/300	2026-03-29 14:16 by Qshers
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4	小張zxy張 2026-03-26	8/400	2026-03-28 19:25 by lbsjt
[考研] 311（085601）求調(diào)劑 +4	liziyeyeye 2026-03-28	4/200	2026-03-28 18:50 by 535743368
[考研] 一志愿北化085600材料專碩275\|有文章專利｜求調(diào)劑 +7	Micky11223 2026-03-25	7/350	2026-03-28 18:34 by 無際的草原
[考研] 085701環(huán)境工程，267求調(diào)劑 +16	minht 2026-03-26	16/800	2026-03-28 12:16 by zllcz
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3	@taotao 2026-03-28	3/150	2026-03-28 10:26 by JourneyLucky
[考研] 08開頭275求調(diào)劑 +4	拉誰不重要 2026-03-26	4/200	2026-03-27 14:12 by Delta2012
[考研] 304材料求調(diào)劑 +4	鐘llll 2026-03-26	4/200	2026-03-27 03:42 by wxiongid
[考研] 求b區(qū)院校調(diào)劑 +4	周56 2026-03-24	5/250	2026-03-25 17:12 by yishunmin
[考研] 各位老師您好：本人初試372分 +5	jj涌77 2026-03-25	6/300	2026-03-25 14:15 by mapenggao