 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1878 | 回復(fù): 10 | ||
| 【懸賞金幣】回答本帖問題，作者ZXR319將贈送您 40 個金幣 | ||
ZXR319新蟲 (小有名氣)
|
[求助]
請會1stopt 的大佬幫幫我 已有2人參與
|
|
|
我在做 lh 機理動力學(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 還需要提供什么可隨時聯(lián)系我，拜謝各位！ 圖2是我用 matlab 擬合的思路，目標(biāo)函數(shù)就是計算值和實驗測量值的差（比如殘差平方和）最小。 @月只藍(lán) 發(fā)自小木蟲Android客戶端 |
回復(fù)此樓» 猜你喜歡
幫我的英語口語老師找學(xué)生
已經(jīng)有4人回復(fù)
-
計算三階力常數(shù)用thirdorder擴胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有61人回復(fù)
-
2026年循環(huán)經(jīng)濟(jì)功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有1人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
獨孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機械動力學(xué)
- 管轄: 計算模擬
2樓2024-09-09 13:31:42
ZXR319
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 74.6
- 帖子: 72
- 在線: 7.5小時
- 蟲號: 33303510
- 注冊: 2023-03-13
- 專業(yè): 化學(xué)工程及工業(yè)化學(xué)
|
老師您的意思是數(shù)據(jù)里每一個時間下應(yīng)該寫出來對應(yīng)的反應(yīng)速率嗎？我通過實驗只能得到不同時間的反應(yīng)物濃度和生成物濃度，那在1stopt 中可以通過代碼得到不同時刻的反應(yīng)速率嗎？我在 Matlab 中是通過 ODE45函數(shù)求反應(yīng)速率的 發(fā)自小木蟲Android客戶端 |
3樓2024-09-09 14:57:27
ZXR319
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 74.6
- 帖子: 72
- 在線: 7.5小時
- 蟲號: 33303510
- 注冊: 2023-03-13
- 專業(yè): 化學(xué)工程及工業(yè)化學(xué)
|
Function r = (k * KA * cA * sqrt(KH * PH2 / Kbc)) / ((1 + KA * cA + KB * cB + sqrt(KH * PH2 / Kbc))^2); 這是模型方程的代碼形式 發(fā)自小木蟲Android客戶端 |
4樓2024-09-09 15:28:45
hzlhm
至尊木蟲 (著名寫手)
- 應(yīng)助: 387 (碩士)
- 金幣: 17808.5
- 紅花: 53
- 帖子: 2879
- 在線: 606.3小時
- 蟲號: 1517335
- 注冊: 2011-11-30
- 性別: GG
- 專業(yè): 常微分方程與動力系統(tǒng)
5樓2024-09-09 22:28:55
ZXR319
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 74.6
- 帖子: 72
- 在線: 7.5小時
- 蟲號: 33303510
- 注冊: 2023-03-13
- 專業(yè): 化學(xué)工程及工業(yè)化學(xué)
送紅花一朵 |
謝謝老師我做的是加氫反應(yīng)，1molA 和3mol 氫氣生成1molB，A 和 B 摩爾比是1:1，所以是不是-dcA/dt=r；dcB/dt=r。不知道這樣可以計算了嗎？ 發(fā)自小木蟲Android客戶端 |
6樓2024-09-10 08:18:34
獨孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機械動力學(xué)
- 管轄: 計算模擬
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
假設(shè) dcB/dt=r 模型公式: cb' = (k*ka*ca*sqrt(kh*3/kbc))/((1+ka*ca+kb*cb+sqrt(kh*3/kbc))^2) 殘差平方和(SSE): 0.000100907681851596 均方差(RMSE): 0.00410096902881087 相關(guān)系數(shù)(R): 0.978424525403226 相關(guān)系數(shù)之平方(R^2): 0.957314551910527 修正R平方(Adj. R^2): 0.928857586517545 確定系數(shù)(DC): 0.956042109592344 F統(tǒng)計(F-Statistic): 5.07803520426272 參數(shù) 最佳估算 ---- -------- k 3947916360160.45 ka 4.19389648022836E15 kh 0.406855254346946 kbc 8.49921244458393 kb 1.28426131584115E16 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) cb 計算 cb 1 15 0.0164 0.023012428915109 2 30 0.0376 0.0370171975730106 3 45 0.0512 0.0473206201899215 4 60 0.06 0.05546914388147 5 90 0.069 0.068244316554898 6 120 0.07368 0.0782291213083908 |
7樓2024-09-11 21:54:51
獨孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機械動力學(xué)
- 管轄: 計算模擬
【答案】應(yīng)助回帖
|
假設(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)計(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 計算 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 |
8樓2024-09-11 21:58:39
ZXR319
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 74.6
- 帖子: 72
- 在線: 7.5小時
- 蟲號: 33303510
- 注冊: 2023-03-13
- 專業(yè): 化學(xué)工程及工業(yè)化學(xué)
9樓2024-09-13 09:34:45
ZXR319
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 74.6
- 帖子: 72
- 在線: 7.5小時
- 蟲號: 33303510
- 注冊: 2023-03-13
- 專業(yè): 化學(xué)工程及工業(yè)化學(xué)
送紅花一朵 |
老師，能不能麻煩您給個源代碼？我自己弄的有問題，擬合出來反應(yīng)速率和濃度都是0。Title "動力學(xué)模型2"; Parameters k,KA,KB,KH,Kbc; Variable t,CA,CB; Constant PH2=3; Function r = (k * K_A * c_A * sqrt(K_H * P_H2 / K_bc)) / ((1 + K_A * c_A + K_B * c_B + sqrt(K_H * P_H2 / K_bc))^2); Data; t CA CB 0 0.1679 0 30 0.1314 0.0367 60 0.0964 0.0719 90 0.0759 0.096 120 0.0575 0.1113 STARTPROGRAM var i: integer; begin for i = 1 to length(t)-1 do r = -(c_A - c_A[i-1]) / (t - t[i-1]); end ENDPROGRAM FIT 這是Ai 給出的代碼，我不知道怎么改 發(fā)自小木蟲Android客戶端 |
10樓2024-09-15 09:34:32
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ZXR319 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程,326分，求調(diào)劑 +5 | KRSLSR 2026-03-10 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 341求調(diào)劑 +3 | 搗蛋豬豬 2026-03-11 | 3/150 |
|
|
[考研] 材料371求調(diào)劑 +7 | 鱷魚? 2026-03-11 | 9/450 |
|
|
|
[考研] 293求調(diào)劑，一志愿陜師大生物學(xué) +3 | ??????.?.??? 2026-03-09 | 3/150 |
|
|
|
[考研] 考研材料與化工，求調(diào)劑 +7 | 戲精丹丹丹 2026-03-09 | 7/350 |
|
|
|
[考研] 材料調(diào)劑，307分 +12 | 張泳銘1 2026-03-09 | 13/650 |
|
|
|
[考研] 312求調(diào)劑 +5 | 陌宸希 2026-03-10 | 5/250 |
|
|
|
[考研] 標(biāo)題:撿漏預(yù)警|08工科/09農(nóng)學(xué)調(diào)劑!英語要求低，過線即有機會! +7 | 馬超放煙花 2026-03-07 | 12/600 |
|
|
|
[考研] 化工0817調(diào)劑 +8 | 燦若星晨 2026-03-10 | 8/400 |
|
|
|
[考研] 327求調(diào)劑 +3 | Ffff03 2026-03-10 | 3/150 |
|
|
|
[考研] 材料專碩調(diào)劑 +7 | 慕辰123 2026-03-05 | 10/500 |
|
|
|
[考研] 求調(diào)劑 一志愿蘇州大學(xué)，0856化工323分 | 本科應(yīng)化 | 有專利/競賽/科研助手經(jīng)歷 | +7 | 橙子cyx 2026-03-06 | 9/450 |
|
|
|
[考博] 求材料讀博院校 +9 | yanglei131 2026-03-08 | 9/450 |
|
|
|
[考研] 環(huán)境調(diào)劑 +5 | 曉看天暮看云 2026-03-09 | 5/250 |
|
|
|
[考研] 家人們 調(diào)劑不迷路 看這里 +8 | likeihood 2026-03-09 | 13/650 |
|
|
|
[碩博家園] 2026級碩士研究生招生/調(diào)劑 +3 | 知足常樂的樂 2026-03-06 | 5/250 |
|
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
|
[考研] 第一志愿上海大學(xué)，專業(yè)化學(xué)工程與技術(shù)，總分288，求調(diào)劑 +3 | 1829197082 2026-03-07 | 3/150 |
|
|
|
[考研] 2026調(diào)劑】考試A區(qū)0703化學(xué)類323分 誠求接收 +3 | 卷柏卷柏 2026-03-05 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +4 | 呼呼？~+123456 2026-03-06 | 4/200 |
|