| 8 | 1/1 | 返回列表 |
| 查看: 1292 | 回復(fù): 7 | ||
xiaohsu2006木蟲 (著名寫手)
|
[求助]
求1stopt擬合動力學(xué)參數(shù)
|
|
小弟沒有高版本的,求高手幫忙,感激不盡! 反應(yīng)動力學(xué): A+B - 2.5C+3.5D+0.5F+G A+D - C+F 3C+D - E+A r1=k1*CB r2=k2*CD r3=k3*CD dCAdt=-r1-r2+r3 dCBdt=-r1 dCCdt=2.5*r1+r2-3*r3 dCDdt=3.5*r1-r2-r3 dCEdt=r3 dCFdt=0.5r1+r2 dCG=r1 初始條件 CA=4.0799, CB=0.0707,CC=CD=CE=CF=CG=0 1 t/s CC CD CE CF 30 2.35E-04 3.36E-03 1.15E-04 4.53E-03 60 1.78E-03 8.62E-03 9.35E-04 6.27E-03 120 2.92E-03 1.20E-02 2.27E-03 6.74E-03 180 4.43E-03 1.54E-02 2.93E-03 6.99E-03 2 t/s CC CD CE CF 30 7.33E-04 5.87E-03 4.50E-04 5.01E-03 60 3.07E-03 1.36E-02 2.67E-03 6.25E-03 120 5.07E-03 1.66E-02 4.42E-03 6.62E-03 180 5.86E-03 1.71E-02 5.30E-03 7.29E-03 3 t/s CC CD CE CF 30 1.67E-03 9.37E-03 1.08E-03 4.77E-03 60 5.55E-03 1.68E-02 4.33E-03 5.91E-03 120 8.29E-03 1.77E-02 6.02E-03 7.50E-03 180 1.00E-02 1.68E-02 6.28E-03 8.12E-03 分別擬合出三個溫度下的反應(yīng)速率常數(shù)k1,k2,k3. |
木蟲 (著名寫手)
鐵桿木蟲 (職業(yè)作家)
|
代碼及第一組數(shù)據(jù)結(jié)果: 均方差(RMSE):0.00206339294816877 殘差平方和(SSE):6.81214473368416E-5 相關(guān)系數(shù)(R): 0.857496435844223 相關(guān)系數(shù)之平方(R^2): 0.735300137485547 確定系數(shù)(DC): -2.62619250526438 F統(tǒng)計(F-Statistic): 0.302254888086507 參數(shù) 最佳估算 -------------------- ------------- k1 0.000571090248522173 k2 0.00237024445720707 k3 0.00354220354969495 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標 cc 計算 cc 目標 cd 計算 cd 目標 ce 計算 ce 目標 cf 計算 cf 1 30 0.000235 0.00250994279834437 0.00336 0.00385071637240616 0.000115 0.000211282929925092 0.00453 0.000741861291302684 2 60 0.00178 0.0041032640107171 0.00862 0.00701015349688254 0.000935 0.000793940277216935 0.00627 0.00172202530354545 3 120 0.00292 0.00514480375607497 0.012 0.0116906247177385 0.00227 0.00281539943089081 0.00674 0.00422532664962184 4 180 0.00443 0.00410080160481611 0.0154 0.0147450856032192 0.00293 0.00564846729296284 0.00699 0.00723295135707574 |
木蟲 (著名寫手)
鐵桿木蟲 (職業(yè)作家)
|
第二組: 均方差(RMSE):0.00231720632138809 殘差平方和(SSE):8.59111221740947E-5 相關(guān)系數(shù)(R): 0.868731725547655 相關(guān)系數(shù)之平方(R^2): 0.754694810973007 確定系數(shù)(DC): -2.87359587747822 F統(tǒng)計(F-Statistic): 0.553612983479765 參數(shù) 最佳估算 -------------------- ------------- k1 0.000745071719863344 k2 0.00228132894455383 k3 0.00353913908564415 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標 cc 計算 cc 目標 cd 計算 cd 目標 ce 計算 ce 目標 cf 計算 cf 1 30 0.000733 0.00325877285433546 0.00587 0.00501713253581565 0.00045 0.000275173337049803 0.00501 0.000958759952701426 2 60 0.00307 0.00529425321594842 0.0136 0.00911952654539591 0.00267 0.00103303889142659 0.00625 0.00221139144787796 3 120 0.00507 0.00650713396046731 0.0166 0.0151522454277829 0.00442 0.00365541671324483 0.00662 0.00537970201979195 4 180 0.00586 0.00495183891595699 0.0171 0.0190254293617908 0.0053 0.00731587838390978 0.00729 0.00915254679655734 第三組: 均方差(RMSE):0.00290108831059639 殘差平方和(SSE):0.000134661014174064 相關(guān)系數(shù)(R): 0.87612822607638 相關(guān)系數(shù)之平方(R^2): 0.767600668527745 確定系數(shù)(DC): -1.04816539387532 F統(tǒng)計(F-Statistic): 2.70389954875217 參數(shù) 最佳估算 -------------------- ------------- k1 0.000850324466622732 k2 0.00300797213925827 k3 0.00339594133977107 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標 cc 計算 cc 目標 cd 計算 cd 目標 ce 計算 ce 目標 cf 計算 cf 1 30 0.00167 0.00381899279623351 0.00937 0.00566811232455994 0.00108 0.000299316725528757 0.00477 0.00115548563185028 2 60 0.00555 0.00643102513322945 0.0168 0.0102027330692526 0.00433 0.00111646934774825 0.00591 0.00274722131198251 3 120 0.00829 0.00889289413916752 0.0177 0.0166430893826131 0.00602 0.00390344884962953 0.0075 0.0068866477146232 4 180 0.01 0.00874715868387847 0.0168 0.0205465263667531 0.00628 0.00772723662236323 0.00812 0.0118613247042733 |
木蟲 (著名寫手)
|
你好,能不能幫我再跑個代碼,寫好了!麻煩你了,先謝過! ConstStr r1=k1*CB,r2=k2*CD,r3=k3*CD; Variable t,CC,CD,CE,CF; InitialODEValue t=0, CA=4079.9, CB=70.7,CC=0, CD=0, CE=0, CF=0, CG=0; ODEFunction CA'=-r1-r2+r3; CB'=-r1; CC'=2.5*r1+r2-3*r3; CD'=3.5*r1-r2-r3; CE'=r3; CF'=0.5*r1+r2; CG'=r1; Data; //t/s CC CD CE CF 30 0.235008523 3.363239986 0.115121201 4.531327956 60 1.780010877 8.616751392 0.934615204 6.268066726 120 2.923654401 12.0012815 2.266319353 6.737997984 180 4.433302192 15.40323797 2.929517334 6.994715676 |
鐵桿木蟲 (職業(yè)作家)
|
均方差(RMSE):2.06284664311107 殘差平方和(SSE):68.0853803679137 相關(guān)系數(shù)(R): 0.858287409589619 相關(guān)系數(shù)之平方(R^2): 0.736657277460058 確定系數(shù)(DC): -2.62363520042484 F統(tǒng)計(F-Statistic): 0.302815847593867 參數(shù) 最佳估算 -------------------- ------------- k1 0.000571141146818135 k2 0.00237051259796401 k3 0.00354131349016907 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標 cc 計算 cc 目標 cd 計算 cd 目標 ce 計算 ce 目標 cf 計算 cf 1 30 0.235008523 2.51033673349595 3.363239986 3.85109149173353 0.115121201 0.211249848393764 4.531327956 0.74194373907573 2 60 1.780010877 4.10425865484917 8.616751392 7.01089022534024 0.934615204 0.793820390023803 6.268066726 1.72224301149238 3 120 2.923654401 5.14738702699264 12.0012815 11.6919857722601 2.266319353 2.81501817652468 6.737997984 4.22595972144955 4 180 4.433302192 4.10534301010109 15.40323797 14.746973106494 2.929517334 5.64773674450559 6.994715676 7.23413221207872 |
木蟲 (著名寫手)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 346求調(diào)劑 一志愿070303有機化學(xué) +5 | 蘿卜燉青菜 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 329求調(diào)劑 +8 | 星野? 2026-03-26 | 8/400 |
|
|
[考研] 0703化學(xué) +14 | 妮妮ninicgb 2026-03-27 | 14/700 |
|
|
[碩博家園] 求調(diào)劑 有機化學(xué)考研356分 +10 | Nadiums 2026-03-25 | 11/550 |
|
|
[考研] 310求調(diào)劑 +9 | 爭取九點睡 2026-03-30 | 9/450 |
|
|
[考研] 0703 化學(xué) 求調(diào)劑,一志愿山東大學(xué) 342 分 +6 | Shern—- 2026-03-28 | 6/300 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +7 | 傾____蕭 2026-03-27 | 8/400 |
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
[考研] 壓國家一區(qū)線,求導(dǎo)師收留,有恩必謝! +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
[考研] 一志愿上海理工能源動力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 26考研-291分-廈門大學(xué)(085601)-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|