| 24小時(shí)熱門(mén)版塊排行榜 |
| 8 | 1/1 | 返回列表 |
| 查看: 1288 | 回復(fù): 7 | ||
xiaohsu2006木蟲(chóng) (著名寫(xiě)手)
|
[求助]
求1stopt擬合動(dòng)力學(xué)參數(shù)
|
|
小弟沒(méi)有高版本的,求高手幫忙,感激不盡! 反應(yīng)動(dò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 分別擬合出三個(gè)溫度下的反應(yīng)速率常數(shù)k1,k2,k3. |
木蟲(chóng) (著名寫(xiě)手)
鐵桿木蟲(chóng) (職業(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)計(jì)(F-Statistic): 0.302254888086507 參數(shù) 最佳估算 -------------------- ------------- k1 0.000571090248522173 k2 0.00237024445720707 k3 0.00354220354969495 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) cc 計(jì)算 cc 目標(biāo) cd 計(jì)算 cd 目標(biāo) ce 計(jì)算 ce 目標(biāo) cf 計(jì)算 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 |
木蟲(chóng) (著名寫(xiě)手)
鐵桿木蟲(chóng) (職業(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)計(jì)(F-Statistic): 0.553612983479765 參數(shù) 最佳估算 -------------------- ------------- k1 0.000745071719863344 k2 0.00228132894455383 k3 0.00353913908564415 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) cc 計(jì)算 cc 目標(biāo) cd 計(jì)算 cd 目標(biāo) ce 計(jì)算 ce 目標(biāo) cf 計(jì)算 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)計(jì)(F-Statistic): 2.70389954875217 參數(shù) 最佳估算 -------------------- ------------- k1 0.000850324466622732 k2 0.00300797213925827 k3 0.00339594133977107 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) cc 計(jì)算 cc 目標(biāo) cd 計(jì)算 cd 目標(biāo) ce 計(jì)算 ce 目標(biāo) cf 計(jì)算 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 |
木蟲(chóng) (著名寫(xiě)手)
|
你好,能不能幫我再跑個(gè)代碼,寫(xiě)好了!麻煩你了,先謝過(guò)! 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 |
鐵桿木蟲(chóng) (職業(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)計(jì)(F-Statistic): 0.302815847593867 參數(shù) 最佳估算 -------------------- ------------- k1 0.000571141146818135 k2 0.00237051259796401 k3 0.00354131349016907 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) cc 計(jì)算 cc 目標(biāo) cd 計(jì)算 cd 目標(biāo) ce 計(jì)算 ce 目標(biāo) cf 計(jì)算 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 |
木蟲(chóng) (著名寫(xiě)手)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703 化學(xué) 求調(diào)劑,一志愿山東大學(xué) 342 分 +5 | Shern—- 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +5 | BrightLL 2026-03-29 | 5/250 |
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
[考研] 297求調(diào)劑 +11 | 田洪有 2026-03-26 | 11/550 |
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +5 | SYA! 2026-03-23 | 5/250 |
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
[考研] 打過(guò)很多競(jìng)賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚(yú)567 2026-03-24 | 4/200 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|