夜雨tju

金蟲 (小有名氣)

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

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

方程組有36個，未知數(shù)為x1-x34,xa,xb，除了未知數(shù)，其他都已知，其中兩個方程組是定積分，兩個上限為未知量，分別為,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)在不做簡化，四個組分一起列式子，方成數(shù)變多，方程組總是不能得到合理的解，請高手指點，如何解這個方程組，小女子不勝感激，由于積分顯示不出來，我上傳了附件和圖片。謝謝各位大神
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);

回復此樓

» 本帖附件資源列表

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

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

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

matlab編程繪圖

» 猜你喜歡

307求調劑已經(jīng)有9人回復
288求調劑已經(jīng)有15人回復
085700資源與環(huán)境308求調劑已經(jīng)有5人回復
294求調劑材料與化工專碩已經(jīng)有14人回復
0817 化學工程 299分求調劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有18人回復
265求調劑已經(jīng)有3人回復
296求調劑已經(jīng)有6人回復
0703化學調劑已經(jīng)有13人回復
085600材料與化工求調劑已經(jīng)有15人回復
生物學調劑招人�。。� 已經(jīng)有4人回復

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

matlab求解非線性方程組已經(jīng)有5人回復
matlab求解非線性方程組的方法已經(jīng)有6人回復
Matlab 求解非線性方程組已經(jīng)有6人回復
用MATLAB求解非線性方程組已經(jīng)有0人回復
matlab求解非線性方程組報錯，請各位大神指點！已經(jīng)有8人回復
matlab求解非線性方程組已經(jīng)有3人回復
請教一個matlab求解非線性方程組的問題已經(jīng)有9人回復
matlab 求解非線性方程組未果已經(jīng)有10人回復
MATLAB求解非線性方程組已經(jīng)有5人回復
MATLAB求解非線性方程組時設置初始參數(shù)的問題已經(jīng)有1人回復
matlab求解非線性方程組的算法？？？？已經(jīng)有1人回復
matlab求解非線性方程組已經(jīng)有3人回復
matlab求解非線性方程組。我真的無能為力啊。求助已經(jīng)有1人回復
matlab求解非線性方程組，錯誤提示怎么解決已經(jīng)有5人回復
matlab 求解非線性方程組牛頓解法已經(jīng)有1人回復
用matlab求解非線性方程組問題已經(jīng)有2人回復
matlab求解非線性方程組，求助！已經(jīng)有6人回復
如何使用matlab求解非線性方程組的所有整數(shù)解？已經(jīng)有9人回復
matlab求解非線性方程組已經(jīng)有16人回復
matlab求解非線性方程組代碼已經(jīng)有2人回復
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復
matlab求解非線性方程組�。。� 已經(jīng)有2人回復
matlab求解非線性方程組已經(jīng)有4人回復

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

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

夜雨tju

金蟲 (小有名氣)

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

引用回帖:

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

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

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

贊一下

回復此樓

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

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

查看全部 16 個回答

hsdwlxhx

金蟲 (正式寫手)

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

【答案】應助回帖

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

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

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

贊一下

回復此樓

BringLightonChip

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

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

cooooldog

鐵桿木蟲 (著名寫手)

ส็็็

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

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

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

贊一下

回復此樓

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

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

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

cooooldog

鐵桿木蟲 (著名寫手)

ส็็็

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

引用回帖:

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

這個不是不幫你，而是工作量太大了。

建議先在版上找擅長用1stOpt的同學試試算算看，也是differential evolution之類的算法
如果效果不錯就將就用

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

贊一下

回復此樓

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

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

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

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

[考研] 085700資源與環(huán)境308求調劑 +5	墨墨漠 2026-03-18	5/250	2026-03-19 22:40 by 學員8dgXkO
[考研] 復試調劑 +4	z1z2z3879 2026-03-14	6/300	2026-03-19 17:18 by fei626-918
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 304求調劑 +3	曼殊2266 2026-03-18	3/150	2026-03-19 14:42 by peike
[考研] 求調劑 +3	Mqqqqqq 2026-03-19	3/150	2026-03-19 14:11 by peike
[考研] 一志愿天大材料與化工（085600）總分338 +5	蔡大美女 2026-03-13	5/250	2026-03-19 10:44 by 是小劉呀～
[教師之家] 焦慮 +9	水冰月月野兔 2026-03-13	13/650	2026-03-19 09:50 by otani
[考研] 328求調劑，英語六級551，有科研經(jīng)歷 +3	生物工程調劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 一志愿武理材料305分求調劑 +5	想上岸的鯉魚 2026-03-18	6/300	2026-03-18 17:53 by 無際的草原
[考研] 26調劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7	步川酷紫123 2026-03-13	7/350	2026-03-18 17:12 by 盡舜堯1
[考研] 312求調劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 303求調劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 301求調劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無際的草原
[考研] 278求調劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 考研求調劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 有沒有道鐵/土木的想調劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 0856專碩279求調劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 085601材料工程315分求調劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 復試調劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] ［0860］321分求調劑，ab區(qū)皆可 +4	寶貴熱 2026-03-13	4/200	2026-03-13 22:01 by 星空星月