版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

夜雨tju

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 852.3
散金: 40
帖子: 77
在線: 52.3小時(shí)
蟲號(hào): 3085703
注冊(cè): 2014-03-24
性別: GG
專業(yè): 分離過程

[求助] matlab求解非線性方程組已有2人參與

方程組有36個(gè)，未知數(shù)為x1-x34,xa,xb，除了未知數(shù)，其他都已知，其中兩個(gè)方程組是定積分，兩個(gè)上限為未知量，分別為,xa,xb，下限皆為0.初值為
x0=[0.036 0.055 0.04 0.032 0.028 0.031 0.06 0.05 0.05 0.03 360 0.04 0.06 0.05 0.034 0.028 0.03 0.062 0.05 0.055 0.04 365 0.91 0.87 0.89 0.86 0.9 0.88 0.85 0.89 0.88 0.9 5900 20]；
xa0=0.01;
xb0=0.01;
開始做了不合理的簡(jiǎn)化，只解了前f(35),f(36),f(1)-f(11)，解釋合理的，現(xiàn)在不做簡(jiǎn)化，四個(gè)組分一起列式子，方成數(shù)變多，方程組總是不能得到合理的解，請(qǐng)高手指點(diǎn)，如何解這個(gè)方程組，小女子不勝感激，由于積分顯示不出來，我上傳了附件和圖片。謝謝各位大神
fff(35)=( K/k1)*∫_0^xa▒〖CA0/(K*CA0*CB0-CC0*CD0-(K*(CA0^2+CA0*CB0)+(CA0*CC0+CA0*CD0)).*x+(K-1)*CA0^2.*x.^2)dx〗 -Lc1/U
fff(36)= =( K/k1)*∫_0^xb▒〖x(11)/(K*x(11)*x(22)-x(33)*x(34)-(K*(x(11)^2+x(11)*x(22))+(x(11)*x(33)+x(11)*x(34))).*x+(K-1)*x(11)^2.*x.^2)dx〗 –Lc2/U

fff(1)=xa-(xa1-x(8))/xa1;
fff(2)=h1/HETP1-(x(4)-x(6))/(K1*x(7)-x(6));
fff(3)=L*(S1+S3+S5)/V*S-(x(4)-x(6))/(xa1-x(7));
fff(4)=x(2)-(x(7)*(S1+S3+S5)+x(8)*(S2+S4))/S;
fff(5)=h2/HETP2-mi*(x(6)-x(5))/(K1*x(11)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4)-mi*x(5));
fff(6)=L/V-mi*(x(6)-x(5))/(mi*x(2)-x(11)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4));
fff(7)=h3/HETP3-(HETP4*x(5)-h4*(K1*x(1)-ya1)-ya1*HETP4)/(K1*HETP4*x(9)-h4*(K1*x(1)-ya1)-ya1*HETP4);
fff(8)=L*(S2+S4)/V*S-(HETP4*x(5)-h4*(K1*x(1)-ya1)-ya1*HETP4)*mi/(HETP4*(xa1*M1+xb1*M2+xc1*M3+xd1*M4-mi*x(9)));
fff(9)=xb-(x(11)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4)-mi*x(10))/(x(11)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4));
fff(10)=x(3)-(x(10)*(S1+S3+S5)+x(9)*(S2+S4))/S;
fff(11)=L/V-(h4/HETP4)*(K1*x(1)-ya1)/(x(3)-x(1));
fff(12)=xa-(xb1-x(19))/xa1;
fff(13)=h1/HETP1-(x(15)-x(17))/(K2*x(18)-x(17));
fff(14)=L*(S1+S3+S5)/V*S-(x(15)-x(17))/(xb1-x(18));
fff(15)=x(13)-(x(18)*(S1+S3+S5)+x(19)*(S2+S4))/S;
fff(16)=h2/HETP2-mi*(x(17)-x(16))/(K2*x(22)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4)-mi*x(16));
fff(17)=L/V-mi*(x(17)-x(16))/(mi*x(13)-x(22)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4));
fff(18)=h3/HETP3-(HETP4*x(16)-h4*(K2*x(12)-yb1)-yb1*HETP4)/(K2*HETP4*x(20)-h4*(K2*x(12)-yb1)-yb1*HETP4);
fff(19)=L*(S2+S4)/V*S-(HETP4*x(16)-h4*(K2*x(12)-yb1)-yb1*HETP4)*mi/(HETP4*(xa1*M1+xb1*M2+xc1*M3+xd1*M4-mi*x(20)));
fff(20)=xb-(x(22)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4)-mi*x(21))/(x(11)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4));
fff(21)=x(14)-(x(21)*(S1+S3+S5)+x(20)*(S2+S4))/S;
fff(22)=L/V-(h4/HETP4)*(K2*x(12)-yb1)/(x(14)-x(12));
fff(23)=xa-(x(30)-xc1)/xa1;
fff(24)=h1/HETP1-(x(26)-x(28))/(K3*x(29)-x(28));
fff(25)=L*(S1+S3+S5)/V*S-(x(26)-x(28))/(xc1-x(29));
fff(26)=x(24)-(x(29)*(S1+S3+S5)+x(30)*(S2+S4))/S;
fff(27)=h2/HETP2-mi*(x(28)-x(27))/(K3*x(33)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4)-mi*x(27));
fff(28)=L/V-mi*(x(28)-x(27))/(mi*x(24)-x(33)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4));
fff(29)=h3/HETP3-(HETP4*x(27)-h4*(K3*x(23)-yc1)-yc1*HETP4)/(K3*HETP4*x(31)-h4*(K3*x(23)-yc1)-yc1*HETP4);
fff(30)=L*(S2+S4)/V*S-(HETP4*x(27)-h4*(K3*x(23)-yc1)-yc1*HETP4)*mi/(HETP4*(xa1*M1+xb1*M2+xc1*M3+xd1*M4-mi*x(31)));
fff(31)=xb-(mi*x(32)-x(33)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4))/(x(11)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4));
fff(32)=x(25)-(x(32)*(S1+S3+S5)+x(31)*(S2+S4))/S;
fff(33)=L/V-(h4/HETP4)*(K3*x(23)-yc1)/(x(25)-x(23));
fff(34)=(mi-(x(11)-x(22)-x(33))*(xa1*M1+xb1*M2+xc1*M3+xd1*M4))/(xa1*M1+xb1*M2+xc1*M3+xd1*M4)-x(34);

回復(fù)此樓

» 本帖附件資源列表

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

2015-08-16 11:38:06, 14.64 K

1樓 2015-08-16 11:47:05

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

hsdwlxhx

金蟲 (正式寫手)

應(yīng)助: 20 (小學(xué)生)
金幣: 2986.8
散金: 204
紅花: 6
帖子: 430
在線: 318.5小時(shí)
蟲號(hào): 1537165
注冊(cè): 2011-12-14
性別: GG
專業(yè): 配位化學(xué)

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

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

你是這是我遇到的用matlab解的最復(fù)雜問題，還是想想怎么簡(jiǎn)化吧，或者找個(gè)新的算法。fsolve這個(gè)函數(shù)估計(jì)是不夠用了。

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

BringLightonChip

2樓2015-08-16 14:45:21

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

cooooldog

鐵桿木蟲 (著名寫手)

ส็็็

應(yīng)助: 237 (大學(xué)生)
金幣: 6101.9
散金: 1114
紅花: 39
帖子: 1380
在線: 553.8小時(shí)
蟲號(hào): 506699
注冊(cè): 2008-02-18
專業(yè): 模式識(shí)別

隨機(jī)優(yōu)化中選擇合適的，
差分演化，粒子群
再結(jié)合fmincon或fsolve里的方法進(jìn)一步迭代

可能要反復(fù)嘗試才能找到最合適的參數(shù)。

贊一下

回復(fù)此樓

ส็็็็็็็็็็็็็็็็็็็็

3樓2015-08-16 16:10:17

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

夜雨tju

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 852.3
散金: 40
帖子: 77
在線: 52.3小時(shí)
蟲號(hào): 3085703
注冊(cè): 2014-03-24
性別: GG
專業(yè): 分離過程

引用回帖:

3樓: Originally posted by cooooldog at 2015-08-16 16:10:17
隨機(jī)優(yōu)化中選擇合適的，
差分演化，粒子群
再結(jié)合fmincon或fsolve里的方法進(jìn)一步迭代

可能要反復(fù)嘗試才能找到最合適的參數(shù)。

大神，我是化工的，對(duì)于matlab學(xué)的不是很精，你能幫我大體寫一下程序么，我自己再摸索一下。這個(gè)程序已經(jīng)困擾我好幾個(gè)月了，謝謝，謝謝

贊一下

回復(fù)此樓

4樓2015-08-16 16:45:40

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

cooooldog

鐵桿木蟲 (著名寫手)

ส็็็

應(yīng)助: 237 (大學(xué)生)
金幣: 6101.9
散金: 1114
紅花: 39
帖子: 1380
在線: 553.8小時(shí)
蟲號(hào): 506699
注冊(cè): 2008-02-18
專業(yè): 模式識(shí)別

引用回帖:

4樓: Originally posted by 夜雨tju at 2015-08-16 16:45:40
大神，我是化工的，對(duì)于matlab學(xué)的不是很精，你能幫我大體寫一下程序么，我自己再摸索一下。這個(gè)程序已經(jīng)困擾我好幾個(gè)月了，謝謝，謝謝...

這個(gè)不是不幫你，而是工作量太大了。

建議先在版上找擅長(zhǎng)用1stOpt的同學(xué)試試算算看，也是differential evolution之類的算法
如果效果不錯(cuò)就將就用

如果必須用matlab之類，
也還是找現(xiàn)成的開源代碼比較好。
這類算法現(xiàn)在有不少開源代碼放出來，針對(duì)不同的問題，效率和效果差別很大
要做不少嘗試。

贊一下

回復(fù)此樓

ส็็็็็็็็็็็็็็็็็็็็

5樓2015-08-16 17:06:06

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

夜雨tju

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 852.3
散金: 40
帖子: 77
在線: 52.3小時(shí)
蟲號(hào): 3085703
注冊(cè): 2014-03-24
性別: GG
專業(yè): 分離過程

內(nèi)容已刪除

回復(fù)此樓

6樓2015-08-16 20:09:45

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

夜雨tju

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 852.3
散金: 40
帖子: 77
在線: 52.3小時(shí)
蟲號(hào): 3085703
注冊(cè): 2014-03-24
性別: GG
專業(yè): 分離過程

我是化工類的，不太擅長(zhǎng)這個(gè)，大神可不可以指點(diǎn)細(xì)一些，我自己去查資料也可以，大神覺得什么函數(shù)可以，我之前也看了智能算法的書。

贊一下

回復(fù)此樓

7樓2015-08-16 20:11:03

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

dingd

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

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

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

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

其它一直常數(shù)都一并給出來啊，比如：s1, s2..., M1, M2,..., Lc1, Lc2, U等等。

贊一下

回復(fù)此樓

8樓2015-08-16 22:38:05

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

dingd

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

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

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

抱歉，應(yīng)該是“已知常數(shù)”

贊一下

回復(fù)此樓

9樓2015-08-16 22:38:40

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

cooooldog

鐵桿木蟲 (著名寫手)

ส็็็

應(yīng)助: 237 (大學(xué)生)
金幣: 6101.9
散金: 1114
紅花: 39
帖子: 1380
在線: 553.8小時(shí)
蟲號(hào): 506699
注冊(cè): 2008-02-18
專業(yè): 模式識(shí)別

引用回帖:

7樓: Originally posted by 夜雨tju at 2015-08-16 20:11:03
我是化工類的，不太擅長(zhǎng)這個(gè)，大神可不可以指點(diǎn)細(xì)一些，我自己去查資料也可以，大神覺得什么函數(shù)可以，我之前也看了智能算法的書。

@dingd 是高手，你把方程和詳細(xì)情況發(fā)給他，用1stOpt幫你算算。
你要把方程寫清楚，常數(shù)的數(shù)值都寫出來，所有常數(shù)的物理意義和取值范圍最好都有

贊一下

回復(fù)此樓

ส็็็็็็็็็็็็็็็็็็็็

10樓2015-08-17 06:56:36

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主夜雨tju 的主題更新

返回列表

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

[考研] 0856調(diào)劑，是學(xué)校就去 +4	sllhht 2026-03-19	5/250	2026-03-19 17:45 by 星空星月
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷 +11	曦熙兮 2026-03-15	11/550	2026-03-19 17:13 by yywzz
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 281求調(diào)劑（0805） +9	煙汐憶海 2026-03-16	19/950	2026-03-19 11:42 by laoshidan
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +10	SUICHILD 2026-03-12	10/500	2026-03-19 09:51 by Delta2012
[考研] 304求調(diào)劑 +6	司空. 2026-03-18	6/300	2026-03-18 23:03 by 星空星月
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7	步川酷紫123 2026-03-13	7/350	2026-03-18 17:12 by 盡舜堯1
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 268求調(diào)劑 +6	簡(jiǎn)單點(diǎn)0 2026-03-17	6/300	2026-03-18 09:04 by 無際的草原
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky

24小時(shí)熱門版塊排行榜

[求助] matlab求解非線性方程組 已有2人參與

» 本帖附件資源列表

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

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

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

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

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

[求助] matlab求解非線性方程組已有2人參與

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