 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 仿真模擬 » MATLAB » 【求助】多元方程組求解及擬合
| 5 | 1/1 | 返回列表 |
| 查看: 2084 | 回復: 9 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
lilei119新蟲 (正式寫手)
|
[交流]
【求助】多元方程組求解及擬合 已有2人參與
|
||
|
多元非線性方程組： (x10-x1)/y=abx1/(1+bx1+cx2+dx3+ex4) (x20-x2)y=acx2/(1+bx1+cx2+dx3+ex4) (x30-x3)/y=adx3/(1+bx1+cx2+dx3+ex4) (x40-x4)/y=aex4/(1+bx1+cx2+dx3+ex4) 其中x10,x20,x30,x40 都是已知的常數(shù)，abcde是待擬合的常數(shù)，現(xiàn)已知一系列的x1,x2,x3,x4和y值，要擬合出上述abcde5個參數(shù)，上面的方程組看樣子是無法手工消元然后用一元非線性擬合了，類似的文獻說用牛頓-拉弗遜迭代法，各位專家?guī)涂纯催@個問題解決起來的難度？并告知具體或者大概步驟。不甚感激 |
回復此樓» 猜你喜歡
291求調(diào)劑
已經(jīng)有7人回復
-
材料與化工272求調(diào)劑
已經(jīng)有13人回復
-
283求調(diào)劑（080500）
已經(jīng)有5人回復
-
材料科學與工程求調(diào)劑
已經(jīng)有5人回復
-
一志愿北京化工大學材料與化工（085600）296求調(diào)劑
已經(jīng)有17人回復
-
考研調(diào)劑
已經(jīng)有7人回復
-
343求調(diào)劑085601
已經(jīng)有3人回復
-
08工科求調(diào)劑286
已經(jīng)有3人回復
-
337求調(diào)劑
已經(jīng)有3人回復
-
329求調(diào)劑，一志愿西北工業(yè)大學，材料工程（085601）
已經(jīng)有10人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 大家看看我的excell擬合方程，怎么能合理的解釋相關性呢？
已經(jīng)有8人回復
- matlab擬合方程參數(shù)時初值的選擇
已經(jīng)有15人回復
- lsqnonlin函數(shù)擬合微分方程組參數(shù)擬合問題
已經(jīng)有10人回復
- 非線性方程組matlab擬合
已經(jīng)有9人回復
- 請教Matlab高手解多元非線性方程組中的參數(shù)
已經(jīng)有6人回復
- 【求助】擬合動力學方程求助
已經(jīng)有13人回復
- 請教怎樣在Design-Expert 7.0 求解擬合方程？
已經(jīng)有9人回復
- 【求助】求教 多元指數(shù)方程如何解 有軟件否？
已經(jīng)有6人回復
- 【求助】 求MAtlab求解程序�。。�
已經(jīng)有20人回復
lilei119
新蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 1734.1
- 散金: 1168
- 紅花: 5
- 帖子: 751
- 在線: 531.3小時
- 蟲號: 569135
- 注冊: 2008-06-04
- 專業(yè): 環(huán)境工程
★
xiegangmai(金幣+1): 鼓勵討論交流！ 2011-04-10 21:42:25
xiegangmai(金幣+1): 鼓勵討論交流！ 2011-04-10 21:42:25
|
謝謝回復， 問題是這樣的， 有一個方程組： (x10-x1)/y=a*b*x1/(1+b*x1+c*x2+d*x3+e*x4) (x20-x2)y=a*c*x2/(1+b*x1+c*x2+d*x3+e*x4) (x30-x3)/y=a*d*x3/(1+b*x1+c*x2+d*x3+e*x4) (x40-x4)/y=a*e*x4/(1+b*x1+c*x2+d*x3+e*x4) x10,x20,x30,x40,y是自變量，x1,x2,x3,x4是因變量，a,b,c,d,e是常數(shù)，現(xiàn)已給出26組關于x1,x2,x3,x4x10,x20,x30,x40,y數(shù)據(jù)，要擬合a,b,c,d,e 這5個 參數(shù)。 擬合時的約束條件是讓SSD這個值最�。� SSD = Σ[ (x1 – x1 (model))2 / (x10)2+ (x2- x2(model))2 / (x20)2+ (x3- x3(model))2 /(x30)2 + (x4- x4(model))2 / (x40) 2] 這里：x1(model)預測值， 下面用化學語言重新表達一遍上述問題。 y是吸附劑，x1,x2,x3,x4是4個被吸附物質，每個物質的吸附都符合Langmuir 模型， 即單獨吸附x1的時候，滿足公式： (x10-x1)/y=a*b*x1/(1+b*x1) ，x10是x1的初始濃度，x1是x1吸附平衡時的濃度，y是吸附劑的投量。a，b是吸附常數(shù)。 現(xiàn)把x1,x2,x3,x4和y同時放入反應器，簡單計算后各物質平衡濃度會同時滿足以下4個方程： (x10-x1)/y=a*b*x1/(1+b*x1+c*x2+d*x3+e*x4) (x20-x2)y=a*c*x2/(1+b*x1+c*x2+d*x3+e*x4) (x30-x3)/y=a*d*x3/(1+b*x1+c*x2+d*x3+e*x4) (x40-x4)/y=a*e*x4/(1+b*x1+c*x2+d*x3+e*x4) 即：x10,x20,x30,x40是各物質的初始濃度，x1,x2,x3,x4吸附平衡后各物質濃度，a,b,c,d,e是吸附常數(shù)，現(xiàn)有26組實驗，已測各物質的初始濃度和吸附平衡之后的濃度，現(xiàn)要擬合a,b,c,d,e 這5個 參數(shù)。 擬合時的約束條件是讓SSD這個值最�。� SSD = Σ[ (x1 – x1 (model))2 / (x10)2+ (x2- x2(model))2 / (x20)2+ (x3- x3(model))2 /(x30)2 + (x4- x4(model))2 / (x40) 2] 這里：x1(model)預測值， 數(shù)據(jù)如下 Y X10 X1 X20 X2 X30 X3 X40 X4 0.098709898 0 0 0 0 0.001071429 1.2854E-05 0 0 0.034516927 0 0 0 0 0.001071429 2.27537E-05 0 0 0.023548457 0 0 0 0 0.001071429 4.57405E-05 0 0 0.009339838 0 0 0 0 0.001071429 0.000158461 0 0 0.001373626 0 0 0 0 0.001071429 0.000714286 0 0 0.000446429 0 0 0 0 0.001071429 0.000928571 0 0 2.67857E-05 3.74667E-05 0.0000324 0.000103226 9.03226E-05 0 0 0 0 5.35714E-05 3.81333E-05 2.86667E-05 0.000103226 7.41935E-05 0 0 0 0 8.92857E-05 3.69333E-05 0.000022 0.0001 5.51613E-05 0 0 0 0 0.000142857 0.0000364 0.0000136 9.83871E-05 2.99355E-05 0 0 0 0 0.00018125 0.0000368 8.66667E-06 0.000102903 1.67742E-05 0 0 0 0 0.000196429 0.000036 5.73333E-06 9.6129E-05 9.6129E-06 0 0 0 0 0.000221429 3.54667E-05 3.06667E-06 9.45161E-05 4.80645E-06 0 0 0 0 0.000267857 0.000034 2.53333E-07 9.35484E-05 2.77419E-06 0 0 0 0 0.000175 6.66667E-06 3.08889E-07 9.03226E-05 2.49336E-06 0 0 0 0 0.0001125 6.66667E-06 2.21111E-06 9.03226E-05 2.32827E-05 0 0 0 0 3.39286E-05 6.66667E-06 4.85556E-06 9.03226E-05 6.57495E-05 0 0 0 0 7.14286E-05 6.66667E-06 3.67778E-06 9.03226E-05 4.78178E-05 0 0 0 0 0.00011 6.66667E-06 2.28889E-06 9.03226E-05 2.57306E-05 0 0 0 0 0.000151786 6.66667E-06 0.0000012 9.03226E-05 1.05882E-05 0 0 0 0 0.000210714 6.66667E-06 0.00000011 9.03226E-05 1.0019E-06 0 0 0 0 8.03571E-05 2.26667E-06 1.5569E-06 9.67742E-05 4.96451E-05 0.001066071 0.001066071 0.0000044 4.30976E-06 0.000139286 0.0000034 1.13333E-06 9.67742E-05 1.22203E-05 0.00105 0.00105 3.26667E-06 3.06667E-06 0.000214286 4.06667E-06 4.46554E-07 9.67742E-05 1.06928E-06 0.001028571 0.001028571 0.0000026 2.22222E-06 0.000321429 4.53333E-06 1.51089E-07 9.67742E-05 2.6058E-07 0.000985714 0.000985714 2.13333E-06 1.3468E-06 0.000535714 4.99333E-06 4.66667E-08 9.67742E-05 9.67742E-08 0.000942857 0.000942857 1.67333E-06 6.06667E-07 |
3樓2011-04-10 08:18:58
★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
xiegangmai(金幣+2): 鼓勵討論交流！ 2011-03-12 20:20:51
zzuwangshilei(金幣+1): 鼓勵交流 2011-03-12 20:26:58
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
xiegangmai(金幣+2): 鼓勵討論交流！ 2011-03-12 20:20:51
zzuwangshilei(金幣+1): 鼓勵交流 2011-03-12 20:26:58
|
本帖內(nèi)容被屏蔽 |
2樓2011-03-12 20:19:36
★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
xiegangmai(金幣+2): 謝謝應助！ 2011-04-10 21:42:38
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
xiegangmai(金幣+2): 謝謝應助！ 2011-04-10 21:42:38
|
本帖內(nèi)容被屏蔽 |
4樓2011-04-10 21:30:21
lilei119
新蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 1734.1
- 散金: 1168
- 紅花: 5
- 帖子: 751
- 在線: 531.3小時
- 蟲號: 569135
- 注冊: 2008-06-04
- 專業(yè): 環(huán)境工程
5樓2011-04-11 06:27:46
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +9 | 張zz111 2026-03-27 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 309求調(diào)劑 +7 | 誰不是少年 2026-03-29 | 7/350 |
|
|
[考研] 一志愿鄭州大學，080500學碩，總分317分求調(diào)劑 +8 | 舉個栗子oi 2026-03-24 | 9/450 |
|
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +9 | 芋泥小鈴鐺 2026-03-28 | 9/450 |
|
|
|
[考研] 材料學碩333求調(diào)劑 +11 | 北道巷 2026-03-24 | 11/550 |
|
|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
|
[考研] 生物學學碩，一志愿湖南大學，初試成績338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
|
[考研] 一志愿華北電力大學能動專碩，293，求調(diào)劑 +3 | 15537177284 2026-03-23 | 5/250 |
|
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
|
[考研] 0703本科鄭州大學求調(diào)劑 +3 | nhj_ 2026-03-25 | 3/150 |
|
|
|
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經(jīng)歷，求調(diào)劑！ +4 | Zuhui0306 2026-03-25 | 4/200 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +4 | 來財aa 2026-03-25 | 4/200 |
|
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分，前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯條 2026-03-23 | 5/250 |
|