 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 583 | 回復(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é)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有239人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有72人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
人生已過半程，不惑之年漸近，回望來時路，滿心皆是感恩。
+5/1680
-
持續(xù)創(chuàng)業(yè)者，A8離異男征婚，事業(yè)遇到瓶頸，尋找創(chuàng)業(yè)生活伴侶
+1/461
-
上海真誠脫單，這是一篇很認真的相親貼，找一個靈魂伴侶
+1/178
-
香港中文大學(xué)（深圳）靳羽華教授交叉實驗室高薪誠聘博士后（光致變色方向）
+2/138
-
電力全國重點實驗室雙一流A類長江學(xué)者團隊招2026年全日制博士1名（4月3日截止）
+2/88
-
天津理工大學(xué)國家杰青王鐵課題組招收2026年博士研究生
+1/82
-
文獻分享，歡迎引用
+1/82
-
浙江大學(xué)侯陽教授團隊誠聘博士后---高性能固態(tài)電解質(zhì)與固態(tài)電池器件設(shè)計
+1/80
-
哈工大鄭州高等研究院劉紹琴與查正寶教授課題組誠招合成生物學(xué)背景相關(guān)的優(yōu)秀人才加盟
+1/30
-
研究生 調(diào)劑 看過來- 機械、控制工程、力學(xué)類（流體力學(xué)、一般力學(xué)等） 211 雙一流
+1/28
-
下周見一個合作者，應(yīng)該問些什么
+1/19
-
招聘科研助理——儲能方向（儲能、流體傳熱、CFD模擬）
+1/16
-
浙江中醫(yī)藥大學(xué)副校長團隊招聘藥學(xué)/生物信息學(xué)等方向高校教師7名，報名通道已開啟
+1/15
-
臨沂大學(xué)數(shù)學(xué)與統(tǒng)計學(xué)院招2026研究生
+1/11
-
寧波東方理工大學(xué)余鵬課題組2026年上海交大聯(lián)培博士生、博士后、科研助理招聘
+1/9
-
天津工業(yè)大學(xué)雙一流高校招博士、碩士、博士后
+1/7
-
大連海事大學(xué)船舶潔凈能源研究中心2026年博士研究生招生啟事
+1/6
-
香港大學(xué)David Srolovitz教授課題組招聘「香江學(xué)者」博士后
+1/5
-
2026考研的同學(xué)，分數(shù)不夠？想調(diào)劑上岸的同學(xué)，來看
+1/2
-
歡迎加入贛南師范大學(xué)先進光電子研究中心
+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ā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +6 | 生物工程調(diào)劑 2026-03-17 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +5 | Zhao anqi 2026-03-22 | 5/250 |
|
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
|
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗） +3 | 1911236844 2026-03-17 | 3/150 |
|
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研] 生物學(xué)調(diào)劑招人�。。� +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|