| 5 | 1/1 | 返回列表 |
| 查看: 4055 | 回復(fù): 17 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
SAYRG金蟲 (小有名氣)
|
[求助]
使用DSLF擬合吸附等溫線時(shí),怎麼樣才能讓得到的參數(shù)為identifiable(參數(shù)估計(jì)值唯一)
|
||
|
想向大夥兒請(qǐng)教的是使用DSLF模型擬合等溫吸附曲線時(shí),有無其他條件需注意? 才會(huì)使得經(jīng)由擬合所得到的參數(shù)具有再現(xiàn)性、信度、一致性以及identifiable(參數(shù)估計(jì)值唯一)的性質(zhì)。 (即只得到一組參數(shù),而非有多組參數(shù)都可將該曲線擬合得很好的情況。) 目前只能確定限制條件之一就是所有的參數(shù)的數(shù)值都必須大於0。 在閱讀一些文獻(xiàn)後發(fā)現(xiàn),有透過DSLF模型擬合得到參數(shù)的文獻(xiàn)裡,其報(bào)導(dǎo)參數(shù)的方式大致上有兩種。以使用DSLF模型的為例。 DSLF Q = (q1*b1*x^n1)/(1 + b1*x^n1) + (q2*b2*x^n2)/(1 + b2*x^n2) 或 Q = (q1*b1*x^(1/n1))/(1 + b1*x^(1/n1)) + (q2*b2*x^(1/n2))/(1 + b2*x^(1/n2)) 其一,直接報(bào)導(dǎo)使用DSLF模型擬合不同溫度時(shí)該吸附材的某氣體等溫吸附曲線所得到的參數(shù)。 溫度1 q1 b1 q2 b2 n1 n2 溫度2 q1' b1' q2' b2' n1' n2' 溫度3 q1" b1" q2" b2" n1" n2" 其二,不像上頭有多組參數(shù),只會(huì)有一組參數(shù)。 q1 b10 E1 n1 q2 b20 E2 n2 b1=b10 exp(E1/RT), ln(b1) = ln(b10)+E/RT b2=b20 exp(E2/RT), ln(b2) = ln(b20)+E/RT 這兩種報(bào)導(dǎo)方式讓我想到兩個(gè)問題! 第一,為何會(huì)有這兩種報(bào)導(dǎo)方式。而第二種方式又是怎麼擬合與分析的呢? 第二,bi肯定會(huì)隨著溫度的變化而有不同,但是q1、q2就其物理意義而言, 應(yīng)該不會(huì)依溫度的高低而有所變動(dòng)。這兩種報(bào)導(dǎo)方式似乎表示著其看待參數(shù) 物理意義的角度有所不同。 而在使用DSLF擬合文獻(xiàn)*1(可在本樓附件一下載)中的化合物1於298K時(shí)吸附CO2的等溫吸附曲線時(shí),無論怎麼擬合,都無法得到其*1之Table S2.所報(bào)導(dǎo)的參數(shù)值。而我用?磉M(jìn)行擬合的軟體為Origin 8.6版。(化合物1之CO2以及N2的單組成吸附資料在本樓附件二或是可取用附件一的Table S3.,而我自己擬合的結(jié)果在本樓末) 另外有想過以下事情。 Q = (q1*b1*x^n1)/(1 + b1*x^n1) + (q2*b2*x^n2)/(1 + b2*x^n2) 或 Q = (q1*b1*x^(1/n1))/(1 + b1*x^(1/n1)) + (q2*b2*x^(1/n2))/(1 + b2*x^(1/n2)) 對(duì)於同一個(gè)吸附劑在不同溫度下吸附某氣體時(shí),會(huì)隨溫度有變化的參數(shù)應(yīng)該就只有the affinity coefficients of site :bi (因?yàn)?bi = bi0 exp(Ei/RT) ; ln(bi) = ln(bi0)+Ei/RT )。因此將q1 q2 (或是q1 q2 n1 n2)看成共享參數(shù),使用DSLF模型以O(shè)rigin中的Global fitting功能進(jìn)行擬合,依舊得不到Table S2.中化合物1吸附CO2時(shí)的DSLF參數(shù)。 因此覺得在擬合時(shí)除了所有參數(shù)都大於0之外, 或許還有其他限制條件需要考慮進(jìn)?懟2艜?huì)蕼厦由擬合得到的參數(shù)具有再現(xiàn)性、一致性以及identifiable(參數(shù)估計(jì)值唯一)的性質(zhì)。 (即只得到一組參數(shù),而非有多組參數(shù)都可將該曲線擬合得很好的情況。) 在使用DSLF模型擬合曲線時(shí)需注意的事情與整個(gè)詳細(xì)過程想和大家請(qǐng)教! *1 Enhanced CO2 Binding Affinity of a High-Uptake rht-Type Metal−Organic Framework Decorated with Acylamide Groups Supporting Information中的Figure S10與Table S2 http://pubs.acs.org/doi/suppl/10.1021/ja110042b |
興趣 |
金蟲 (小有名氣)
|
?L???Origin ?? Global fitting????????У??????F?o?????(?????)??????? ???????J??Global fitting????????]?e?????????Щ????????????? ??????????*1??????o?????c?c???h?????x??? *1 http://www.gaoyang168.com/bbs/viewthread.php?tid=5201919&page=1#pid7 |
榮譽(yù)版主 (著名寫手)
彩色掛圖
![]() |
專家經(jīng)驗(yàn): +2 |

|
本帖內(nèi)容被屏蔽 |
金蟲 (小有名氣)
送鮮花一朵|
不知道能不能再向你請(qǐng)教一些問題? 1. " 擬合的時(shí)候不同溫度的q1,q2應(yīng)該設(shè)為常數(shù) " 我個(gè)人非常同意這個(gè)看法,q1,q2應(yīng)該不會(huì)隨著溫度而有變化。 a. q1, q2是先以其他方式求得的嗎?那方式又是什麼呢? b. 還是在擬合過程當(dāng)中將q1, q2視作共享參數(shù),進(jìn)而透過Global fitting得到? c. 我想了解怎麼得到q1, q2的。 2. " 比較好的方法是通過分子模擬估計(jì)吸附熱等參數(shù),確定數(shù)據(jù)可能確定范圍再代入擬合。 " a. 目前有辦法透過維里(Virial)方程式計(jì)算零覆蓋率(zero coverage)與不同覆蓋率時(shí)的等量吸附熱(Qst, Isosteric heat or absorption of heat)。但是不知道如何推估E1和E2?煞窈臀艺f說? |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 26申博 +4 | 加油沖。 2026-03-26 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 349求調(diào)劑 +4 | zwjjjjjj 2026-03-31 | 4/200 |
|
|
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +11 | 秋二十二 2026-03-30 | 11/550 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +18 | 荷包蛋hyj 2026-03-26 | 18/900 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +20 | Mr. Z 2026-03-25 | 21/1050 |
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +6 | 星沉uu 2026-03-26 | 7/350 |
|
|
[考研] 293分求調(diào)劑,外語(yǔ)為俄語(yǔ) +5 | 加一一九 2026-03-31 | 5/250 |
|
|
[考研] 335求調(diào)劑 +3 | 321* 2026-03-31 | 3/150 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 壓國(guó)家一區(qū)線,求導(dǎo)師收留,有恩必謝! +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|