當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

夜雨tju

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 852.3
散金: 40
帖子: 77
在線: 52.3小時(shí)
蟲號: 3085703
注冊: 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;
開始做了不合理的簡化，只解了前f(35),f(36),f(1)-f(11)，解釋合理的，現(xiàn)在不做簡化，四個(gè)組分一起列式子，方成數(shù)變多，方程組總是不能得到合理的解，請高手指點(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)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(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的回帖

dingd

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

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

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

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

贊一下

回復(fù)此樓

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

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

查看全部 16 個(gè)回答

hsdwlxhx

金蟲 (正式寫手)

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

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

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

你是這是我遇到的用matlab解的最復(fù)雜問題，還是想想怎么簡化吧，或者找個(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í)
蟲號: 506699
注冊: 2008-02-18
專業(yè): 模式識別

隨機(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í)
蟲號: 3085703
注冊: 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ù)。

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

贊一下

回復(fù)此樓

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

已閱回復(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
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 294求調(diào)劑材料與化工專碩 +14	陌の森林 2026-03-18	14/700	2026-03-19 22:38 by 學(xué)員8dgXkO
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +5	Doleres 2026-03-19	5/250	2026-03-19 20:14 by 制度的
[考研] 能源材料化學(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
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3	@taotao 2026-03-19	3/150	2026-03-19 14:07 by peike
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽w
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs