| 5 | 1/1 | 返回列表 |
| 查看: 1851 | 回復(fù): 10 | ||
| 【懸賞金幣】回答本帖問題,作者ZXR319將贈(zèng)送您 40 個(gè)金幣 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
ZXR319新蟲 (小有名氣)
|
[求助]
請(qǐng)會(huì)1stopt 的大佬幫幫我 已有2人參與
|
|
|
我在做 lh 機(jī)理動(dòng)力學(xué),已經(jīng)有一部分?jǐn)?shù)據(jù)了,目前已經(jīng)用 matlab 處理了,擬合初值是隨便給的,我現(xiàn)在想用1stopt 軟件擬合一下,擬合的結(jié)果作為 matlab 的參數(shù)初值,以便更準(zhǔn)確。求大佬幫我擬合一下,最好有源代碼,金幣我可以盡量給,有多少都可以給,好商量。 數(shù)據(jù)如下:t ca cb 0 0.1679 0 15 0.1511 0.0164 30 0.1322 0.0376 45 0.115 0.0512 60 0.1037 0.06 90 0.1003 0.069 120 0.09486 0.07368 還需要提供什么可隨時(shí)聯(lián)系我,拜謝各位! 圖2是我用 matlab 擬合的思路,目標(biāo)函數(shù)就是計(jì)算值和實(shí)驗(yàn)測(cè)量值的差(比如殘差平方和)最小。 @月只藍(lán) 發(fā)自小木蟲Android客戶端 |
版主 (知名作家)
|
假設(shè) dcA/dt=-r 模型公式: ca' = -(k*ka*ca*sqrt(kh*3/kbc))/((1+ka*ca+kb*cb+sqrt(kh*3/kbc))^2) 殘差平方和(SSE): 0.000192586899315775 均方差(RMSE): 0.005665493495948 相關(guān)系數(shù)(R): 0.959232883604991 相關(guān)系數(shù)之平方(R^2): 0.920127724989146 修正R平方(Adj. R^2): 0.866879541648577 確定系數(shù)(DC): 0.91769457654023 F統(tǒng)計(jì)(F-Statistic): 2.54023664039648 參數(shù) 最佳估算 ---- -------- k 2837450.21546848 ka 2.67743713648235E-7 kh 67.0303631498266 kbc 0.100651239217611 kb 687.010099643109 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) ca 計(jì)算 ca 1 15 0.1511 0.143509362937322 2 30 0.1322 0.129868469396974 3 45 0.115 0.120093597136612 4 60 0.1037 0.112222162940757 5 90 0.1003 0.0997387665702967 6 120 0.09486 0.0893241064493534 |

版主 (知名作家)

新蟲 (小有名氣)
|
老師您的意思是數(shù)據(jù)里每一個(gè)時(shí)間下應(yīng)該寫出來對(duì)應(yīng)的反應(yīng)速率嗎?我通過實(shí)驗(yàn)只能得到不同時(shí)間的反應(yīng)物濃度和生成物濃度,那在1stopt 中可以通過代碼得到不同時(shí)刻的反應(yīng)速率嗎?我在 Matlab 中是通過 ODE45函數(shù)求反應(yīng)速率的 發(fā)自小木蟲Android客戶端 |
新蟲 (小有名氣)
|
Function r = (k * KA * cA * sqrt(KH * PH2 / Kbc)) / ((1 + KA * cA + KB * cB + sqrt(KH * PH2 / Kbc))^2); 這是模型方程的代碼形式 發(fā)自小木蟲Android客戶端 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)求調(diào)劑,總分320分,一志愿中南民族大學(xué) +3 | 19286712440 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工(高分子)290分求調(diào)劑 +4 | wengyujian 2026-03-04 | 4/200 |
|
|
[考研] 化學(xué)工程求調(diào)劑 +6 | 化工人999 2026-03-04 | 6/300 |
|
|
[考研] 292求調(diào)劑 +9 | yhk_819 2026-02-28 | 9/450 |
|
|
[考研] 材料328求調(diào)劑 +7 | 一個(gè)蘿卜02 2026-03-03 | 7/350 |
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
[考研] 325求調(diào)劑 +5 | 學(xué)家科 2026-03-04 | 5/250 |
|
|
[論文投稿]
EST拒稿重投
5+3
|
15102603076 2026-03-02 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 280分 +10 | yyqqhh 2026-03-03 | 10/500 |
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
[考研] 一志愿中科大能動(dòng)297求調(diào)劑,本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
[論文投稿]
通訊作者寫誰,問題是你意想不到的問題
15+3
|
阿爾法啊 2026-03-01 | 3/150 |
|
|
[考研] 338求調(diào)劑 +5 | 18162027187 2026-03-02 | 6/300 |
|
|
[考研] 0854復(fù)試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
[考研] 321求調(diào)劑一志愿東北林業(yè)大學(xué)材料與化工英二數(shù)二 +5 | 蟲蟲蟲蟲蟲7 2026-03-01 | 9/450 |
|
|
[考研] 272求調(diào)劑 +9 | 材紫有化 2026-02-28 | 9/450 |
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|
|
[考研] 328求調(diào)劑 +3 | aaadim 2026-03-01 | 5/250 |
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 3/150 |
|