xwl196

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2985.5
帖子: 109
在線: 52.2小時
蟲號: 1805601
注冊: 2012-05-09
性別: MM
專業(yè): 能源化工

[求助] matlab求解含中間變量的非線性方程組已有1人參與

方程組如附件所示，其中：
T=328.2;
ρmc=1.0254;
Eii=0.432*195.2;
C=1.23;
Eis=-1250;
ρads為中間變量
求大神指點：1.如何求解想x和y；2.或者如何從方程(1)中把ρads表示成關(guān)于x的表達式。

1.jpg

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有3人回復(fù)
265求調(diào)劑已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有6人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有12人回復(fù)
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有3人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 已經(jīng)有4人回復(fù)
材料學(xué)碩318求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有6人回復(fù)
081700化工學(xué)碩調(diào)劑已經(jīng)有3人回復(fù)

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

matlab 如何求解多元二階微分方程組已經(jīng)有7人回復(fù)
求助如何用matlab求非線性方程組的解已經(jīng)有8人回復(fù)
求助matlab高手指點：大型非線性方程組（牛頓拉佛森法）已經(jīng)有7人回復(fù)
MATLAB解非線性方程組已經(jīng)有9人回復(fù)
MATLAB求解五元非線性方程組，在線等已經(jīng)有7人回復(fù)
求用牛頓迭代法解非線性方程組的matlab程序已經(jīng)有4人回復(fù)
matlab非線性最小二乘法求解已經(jīng)有12人回復(fù)
非線性二階微分方程組求解，matlab 已經(jīng)有12人回復(fù)
matlab 非線性微分方程求解已經(jīng)有3人回復(fù)
matlab非線性方程組該如何編程。求助已經(jīng)有10人回復(fù)
用matlab求解一個非線性方程組的解已經(jīng)有4人回復(fù)
請教一個matlab求解非線性方程組的問題已經(jīng)有9人回復(fù)
非常簡單��！MATLAB求解非線性一元方程已經(jīng)有5人回復(fù)
用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？已經(jīng)有11人回復(fù)
matlab如何求解一個非線性微分方程組已經(jīng)有8人回復(fù)
求教matlab線性優(yōu)化求最小值的問題已經(jīng)有7人回復(fù)
MATLAB求解非線性方程組已經(jīng)有5人回復(fù)
求Matlab解非線性方程代碼已經(jīng)有12人回復(fù)
matlab求解非線性方程組，錯誤提示怎么解決已經(jīng)有5人回復(fù)
求助matlab---fsolve解非線性方程組已經(jīng)有6人回復(fù)
matlab求解非線性方程組已經(jīng)有16人回復(fù)
非線性方程組matlab擬合已經(jīng)有9人回復(fù)
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(shù)解！已經(jīng)有18人回復(fù)
【求助】求教matlab解非線性方程組已經(jīng)有9人回復(fù)

1樓 2015-04-09 10:04:17

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

dingd

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

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

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

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

nex有數(shù)值嗎，多少？

贊一下

回復(fù)此樓

2樓2015-04-09 14:36:36

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

xwl196

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2985.5
帖子: 109
在線: 52.2小時
蟲號: 1805601
注冊: 2012-05-09
性別: MM
專業(yè): 能源化工

引用回帖:

2樓: Originally posted by dingd at 2015-04-09 14:36:36
nex有數(shù)值嗎，多少？

ρb=x=[[0.0116
0.0240
0.0502
0.0804
0.1170
0.1608
0.2152
0.2972
0.4123
0.5257
0.6100];
所以問題轉(zhuǎn)化為在方程(1)中用matlab的solve求解ρads，然后將ρads帶入方程(2)中求解nex就好了！
matlab程序：
>> T=328.2;
rhomc=1.0254;
Eii=0.432*195.2;
C=1.23;
Eis=-1250;
x=[0.0116
0.0240
0.0502
0.0804
0.1170
0.1608
0.2152
0.2972
0.4123
0.5257
0.6100];
>> syms rhoads
>> [rhoads]=solve(log(rhoads*(rhomc-x)/(x*(rhomc-rhoads)))+(7*rhoads-8*x)/rhomc*Eii/T+Eis/T==0)
運行后出現(xiàn)的提示是：
??? Error using ==> mupadmex
Error in MuPAD command: Array sizes must match.

Error in ==> sym.plus at 14
X = mupadmex('symobj::zip',A.s,B.s,'_plus');
沒找到問題出現(xiàn)在哪？

贊一下

回復(fù)此樓

3樓2015-04-09 15:50:39

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

dingd

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

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

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xwl196: 金幣+10, ★★★很有幫助 2015-04-09 16:40:44

Solve命令是求符號解，不一定有的�？梢允菙�(shù)值解。
下面是用1stOpt求解，很簡單：

CODE:

Constant T=328.2,rhomc=1.0254,Eii=0.432*195.2,C=1.23,Eis=-1250;

LoopConstant x=[0.0116,0.0240,0.0502,0.0804,0.1170,0.1608,0.2152,0.2972,0.4123,0.5257,0.6100];

PassParameter nex=2*C*(rhoads-x)/rhomc;

Function

ln(rhoads*(rhomc-x)/(x*(rhomc-rhoads)))+(7*rhoads-8*x)/rhomc*Eii/T+Eis/T;

Ñ-»·3￡êy x rhoads nex
0.0116 0.257851015773962 0.590771892728639
0.024 0.378121213925542 0.849559378054255
0.0502 0.520370658449439 1.1279693970993
0.0804 0.618855300376744 1.29178860827656
0.117 0.699976313181814 1.3985973575456
0.1608 0.769244989810605 1.45969833716997
0.2152 0.831497991791605 1.47853818978677
0.2972 0.896153280504329 1.43692712116311
0.4123 0.952359132344292 1.2956363034591
0.5257 0.984696573619486 1.10116205490924
0.61 0.999729281534294 0.934985403329788

贊一下

回復(fù)此樓

4樓2015-04-09 16:05:54

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

xwl196

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2985.5
帖子: 109
在線: 52.2小時
蟲號: 1805601
注冊: 2012-05-09
性別: MM
專業(yè): 能源化工

引用回帖:

4樓: Originally posted by dingd at 2015-04-09 16:05:54
Solve命令是求符號解，不一定有的�？梢允菙�(shù)值解。
下面是用1stOpt求解，很簡單：

Constant T=328.2,rhomc=1.0254,Eii=0.432*195.2,C=1.23,Eis=-1250;
LoopConstant x=;
PassParameter nex=2*C*(rhoads-x)/rh ...

你用的1stopt是哪個版本的，我用1stopt1.5運行的結(jié)果如下：
"model nex"

====== 結(jié)果 ======

迭代數(shù): 16
計算用時(時:分:秒:毫秒): 00:00:00:50
計算中止原因: 達到收斂判定標準
優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法
函數(shù)表達式 ln(rhoads*(1.0254-x)/(x*(1.0254-rhoads)))+(7*rhoads-8*x)/1.0254*84.3264/328.2+(-1250)/328.2
目標函數(shù)值: 0
rhoads: 0.485072231445505
x: 0.0421283782771957

====== 計算結(jié)束 ======

沒有得到你運行的結(jié)果，是不是我用的1stopt版本太低的原因？

贊一下

回復(fù)此樓

5樓2015-04-09 16:29:27

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

dingd

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

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

引用回帖:

5樓: Originally posted by xwl196 at 2015-04-09 16:29:27
你用的1stopt是哪個版本的，我用1stopt1.5運行的結(jié)果如下：
"model nex"

====== 結(jié)果 ======

迭代數(shù): 16
計算用時(時:分:秒:毫秒): 00:00:00:50
計算中止原因: 達到收斂判定標準
優(yōu)化算法: 麥 ...

版本問題。

回復(fù)此樓

6樓2015-04-09 16:33:16

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

xwl196

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2985.5
帖子: 109
在線: 52.2小時
蟲號: 1805601
注冊: 2012-05-09
性別: MM
專業(yè): 能源化工

引用回帖:

6樓: Originally posted by dingd at 2015-04-09 16:33:16
版本問題。...

高版本的1stopt找不到，只能回到matlab中繼續(xù)奮斗了

感謝！

贊一下

回復(fù)此樓

7樓2015-04-09 16:40:33

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

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

[考研] 329求調(diào)劑 +5	想上學(xué)吖吖 2026-03-19	5/250	2026-03-19 23:56 by 23Postgrad
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +18	rare12345 2026-03-18	18/900	2026-03-19 22:36 by 學(xué)員8dgXkO
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	10/500	2026-03-19 16:05 by 余麟余
[考研] 266求調(diào)劑 +5	陽陽哇塞 2026-03-14	10/500	2026-03-19 15:08 by 陽陽哇塞
[考研] 一志愿天大材料與化工（085600）總分338 +5	蔡大美女 2026-03-13	5/250	2026-03-19 10:44 by 是小劉呀～
[考研] 332求調(diào)劑 +3	ydfyh 2026-03-17	3/150	2026-03-19 10:14 by 功夫瘋狂
[考研] 304求調(diào)劑 +6	司空. 2026-03-18	6/300	2026-03-18 23:03 by 星空星月
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 085600材料與化工 +5	安全上岸！ 2026-03-16	5/250	2026-03-18 15:33 by cmz0325
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考博] 環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名 +3	QGZDSYS 2026-03-13	5/250	2026-03-18 11:13 by QGZDSYS
[考研] 0703化學(xué)求調(diào)劑總分331 +3	ZY-05 2026-03-13	3/150	2026-03-18 10:58 by macy2011
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 288求調(diào)劑 +4	奇點0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky