| 5 | 1/1 | 返回列表 |
| 查看: 1850 | 回復(fù): 10 | |||
| 【懸賞金幣】回答本帖問題,作者ZXR319將贈送您 40 個(gè)金幣 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
ZXR319新蟲 (小有名氣)
|
[求助]
請會1stopt 的大佬幫幫我 已有2人參與
|
||
|
我在做 lh 機(jī)理動力學(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)測量值的差(比如殘差平方和)最小。 @月只藍(lán) 發(fā)自小木蟲Android客戶端 |
新蟲 (小有名氣)
送紅花一朵 |
老師,能不能麻煩您給個(gè)源代碼?我自己弄的有問題,擬合出來反應(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客戶端 |
版主 (知名作家)

新蟲 (小有名氣)
|
老師您的意思是數(shù)據(jù)里每一個(gè)時(shí)間下應(yīng)該寫出來對應(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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩290求調(diào)劑 +3 | 杰尼龜aaa 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿鄭州大學(xué),學(xué)碩,物理化學(xué), 333求調(diào)劑 +4 | 李魔女斗篷 2026-03-04 | 4/200 |
|
|
[考博] 26申博-目前4篇SCI一作 +3 | chen_2024 2026-03-02 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | G212 2026-03-03 | 5/250 |
|
|
[考研] 環(huán)境調(diào)劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
[基金申請] 請問大家,研究風(fēng)險(xiǎn)與應(yīng)對措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 266材料化工求調(diào)劑 +3 | 哇塞王帥 2026-03-03 | 3/150 |
|
|
[考研] 一志愿東北大學(xué)化學(xué)314分求調(diào)劑 +5 | lr1212.. 2026-03-02 | 5/250 |
|
|
[考研] 一志愿中科大能動297求調(diào)劑,本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
[考研]
材料工程專碩283求調(diào)劑
5+8
|
,? 2026-03-02 | 10/500 |
|
|
[考研] 324求調(diào)劑 +4 | wxz2 2026-03-03 | 5/250 |
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
[考研] 調(diào)劑材料學(xué)碩 +4 | 詞凝Y 2026-03-02 | 4/200 |
|
|
[考研] 288求調(diào)劑 +3 | 少71.8 2026-03-02 | 5/250 |
|
|
[考博] 博士自薦 +4 | kkluvs 2026-02-28 | 5/250 |
|
|
[基金申請] 此成果不能導(dǎo)入原因:元數(shù)據(jù)必填信息不完整,可 進(jìn)行補(bǔ)充。 +4 | Kittylucky 2026-03-02 | 5/250 |
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|
|
[考研] 275求調(diào)劑 +3 | L-xin? 2026-03-01 | 6/300 |
|
|
[考研] 295復(fù)試調(diào)劑 +3 | 簡木ChuFront 2026-03-01 | 3/150 |
|
|
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6 | Soliloquy_Q 2026-02-28 | 11/550 |
|