 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 587 | 回復(fù): 1 | |||
[交流]
簡單的多元化學(xué)平衡模擬，請各位進來幫忙
|
|
一個簡單的matlab程序問題，麻煩各位進來指點 各位專家，請教關(guān)于擬合方面的問題： 有一個方程組： (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組關(guān)于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ù)測值， 下面用化學(xué)語言重新表達一遍上述問題。 y是吸附劑，x1,x2,x3,x4是4個被吸附物質(zhì)，每個物質(zhì)的吸附都符合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同時放入反應(yīng)器，簡單計算后各物質(zhì)平衡濃度會同時滿足以下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是各物質(zhì)的初始濃度，x1,x2,x3,x4吸附平衡后各物質(zhì)濃度，a,b,c,d,e是吸附常數(shù)，現(xiàn)有26組實驗，已測各物質(zhì)的初始濃度和吸附平衡之后的濃度，現(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)預(yù)測值， 數(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 非常感謝各位專家?guī)兔? |
回復(fù)此樓» 猜你喜歡
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有89人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有29人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有26人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
UJN物理學(xué)專業(yè)調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復(fù)
-
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名
已經(jīng)有0人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
陜西理工大學(xué)資源與環(huán)境、材料化工部分調(diào)劑名額，理工兼收
+1/184
-
廣西師范大學(xué)-王艷課題組2026年招收藥學(xué)學(xué)碩 、生物與醫(yī)藥專碩調(diào)劑生
+1/89
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
+2/48
-
導(dǎo)師招生
+1/38
-
三峽大學(xué)國家青年人才課題組招2026年化學(xué)/材料/化工類碩士調(diào)劑生
+1/38
-
天津理工大學(xué)功能晶體研究院（晶體材料全國重點實驗室）杰青團隊招收2026年碩士研究生
+1/38
-
湖北汽車工業(yè)學(xué)院接受碩士調(diào)劑！
+1/37
-
2026年武漢紡織大學(xué)材料學(xué)院鄧波教授團隊接收調(diào)劑（碩士生4名）
+1/36
-
天津理工大學(xué)功能晶體研究院（晶體材料全國重點實驗室）杰青團隊招收2026年博士研究生
+1/18
-
天津科技大學(xué)-生物質(zhì)資源利用與制漿技術(shù)團隊2026年碩士招生
+1/9
-
中國農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院 - 招聘杰出人才（教授-專技四級）、青年研究員、優(yōu)秀人才
+1/8
-
浙江寧波 生物與醫(yī)藥 接受生物學(xué)、園藝學(xué)、作物學(xué)、農(nóng)業(yè)調(diào)劑
+1/7
-
江西理工大學(xué) 稀土學(xué)院 材料科學(xué)與工程 招收2026屆調(diào)劑碩士研究生
+1/6
-
陜西理工大學(xué)材料學(xué)院新能源團隊招調(diào)劑生�。�！
+1/6
-
北京石油化工學(xué)院機器人應(yīng)用技術(shù)研究所2026年招收調(diào)劑研究生！
+1/5
-
哈工程煙臺院材化專業(yè)2026調(diào)劑信息
+1/4
-
蘇州科技大學(xué)胡老師課題組招2026年碩士生若干名
+1/4
-
廣東工業(yè)大學(xué)生態(tài)環(huán)境與海洋學(xué)院 楊盼教授團隊招收博士研究生
+1/2
-
2026年北京石油化工學(xué)院環(huán)境學(xué)科-水污染控制工程方向課題組歡迎您
+1/1
-
2026年江西科技師范大學(xué)Tide課題組招收化學(xué)、生物、藥學(xué)、材料相關(guān)學(xué)碩或?qū)４T調(diào)劑生
+1/1
★
lilei119(金幣+2):謝謝參與
lilei119(金幣+8): 謝謝，太感謝你了。記得你上次回過我的鐵，能交個朋友嗎？ 我同濟大學(xué)環(huán)境學(xué)院的研究生。 2011-04-11 06:17:39
lilei119(金幣+2):謝謝參與
lilei119(金幣+8): 謝謝，太感謝你了。記得你上次回過我的鐵，能交個朋友嗎？ 我同濟大學(xué)環(huán)境學(xué)院的研究生。 2011-04-11 06:17:39
|
本帖內(nèi)容被屏蔽 |
2樓2011-04-10 21:20:14
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 lilei119 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600，材料與化工321分調(diào)劑 +8 | 大饞小子 2026-03-27 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 【求調(diào)劑】新能源材料本科，一志愿211，初試321 +6 | 求調(diào)劑學(xué)校， 2026-04-02 | 6/300 |
|
|
|
[考研] 化學(xué)工程專碩324分，一志愿中國礦業(yè)大學(xué)求調(diào)劑 +7 | 耿耿1314 2026-04-01 | 7/350 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +14 | 一定上岸aaa 2026-03-27 | 18/900 |
|
|
[考研] 303分 0807學(xué)碩求調(diào)劑 +3 | TYC3632 2026-04-01 | 3/150 |
|
|
|
[考研] 350求調(diào)劑 +7 | 阿佳～ 2026-03-31 | 7/350 |
|
|
|
[考研] 一志愿鄭大085600，310分求調(diào)劑 +6 | 李瀟可 2026-03-26 | 6/300 |
|
|
|
[考研] 330分求調(diào)劑 +11 | qzenlc 2026-03-29 | 11/550 |
|
|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行 +25 | lllllos 2026-03-30 | 26/1300 |
|
|
|
[考研] 一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑 +5 | 雙馬尾痞老板2 2026-03-31 | 5/250 |
|
|
|
[考研] 江蘇蘇北高校誠邀調(diào)劑同學(xué) +3 | zzll406 2026-03-31 | 3/150 |
|
|
|
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級。 +10 | XX小鄧 2026-03-29 | 10/500 |
|
|
|
[考研] 求收留 +8 | 1943443204 2026-03-28 | 8/400 |
|
|
|
[考研] 286求調(diào)劑 +6 | Faune 2026-03-30 | 6/300 |
|
|
|
[考研] 0703化學(xué) +20 | 妮妮ninicgb 2026-03-27 | 20/1000 |
|
|
|
[考研] 調(diào)劑求院校招收 +7 | 鶴鯨鴿 2026-03-28 | 7/350 |
|
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +6 | 星沉uu 2026-03-26 | 7/350 |
|
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
|
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|