| 查看: 2182 | 回復(fù): 11 | ||
蘇拉鐵蟲 (小有名氣)
|
[求助]
matlab擬合反應(yīng)動(dòng)力學(xué)參數(shù) 已有2人參與
|
|
|
各位大佬,可以幫我擬合一下動(dòng)力學(xué)參數(shù)嗎,可以有償。先謝過各位大佬。 t B C A 0 0.0024 0.0000 0.1894 6 0.0422 0.0000 0.1496 12 0.0838 0.0014 0.1066 18 0.1261 0.0061 0.0596 24 0.1470 0.0151 0.0297 30 0.1447 0.0361 0.0111 36 0.1279 0.0591 0.0047 42 0.1127 0.0767 0.0025 48 0.1009 0.0898 0.0011 54 0.0848 0.1070 0.0000 60 0.0721 0.1197 0.0000 ρ=0.0077 [H2]=4 r1 = k1*KA*[A]*KH2*[H2]/(1+KA*[A]+KB*[B]+KC*[C]+sqrt(KH2*[H2]))^3 r1 = k2*KB*[B]*KH2*[H2]/(1+KA*[A]+KB*[B]+KC*[C]+sqrt(KH2*[H2]))^3 1/ρ*d[A]/dt = -r1 1/ρ*d[B]/dt = r1-r2 1/ρ*d[C]/dt = r2 |
版主 (知名作家)
|
均方差(RMSE): 0.00213787148383087 殘差平方和(SSR): 0.000137114834441316 相關(guān)系數(shù)(R): 0.998778073301579 相關(guān)系數(shù)之平方(R^2): 0.997557639708014 修正R平方(Adj. R^2): 0.995605193318058 確定系數(shù)(DC): 0.996987006259065 F統(tǒng)計(jì)(F-Statistic): 161.269285477097 參數(shù) 最佳估算 -------------------- ------------- k1 13.5564432263352 k2 9.47763860792964 ka 28.7183814838917 kb 4.13770213410864 kc 3.5066823001231 kh2 2.55638641506165 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) b 計(jì)算 b 目標(biāo) c 計(jì)算 c 目標(biāo) a 計(jì)算 a 1 6 0.0422 0.0424108511588022 0 0.000565879792438171 0.1496 0.14882326904876 2 12 0.0838 0.0837866483079356 0.0014 0.00281912430154442 0.1066 0.10519422739052 3 18 0.1261 0.121738165846101 0.0061 0.00847370300303641 0.0596 0.0615881311508623 4 24 0.147 0.145287462165307 0.0151 0.0200631340157014 0.0297 0.0264494038189916 5 30 0.1447 0.145967706327343 0.0361 0.0378459854381184 0.0111 0.00798630823453884 6 36 0.1279 0.132131188851013 0.0591 0.0577272070532937 0.0047 0.00194160409569314 7 42 0.1127 0.115034522234822 0.0767 0.0763300773851729 0.0025 0.0004354003800049 8 48 0.1009 0.0989887391692515 0.0898 0.0927163359628322 0.0011 9.49248679162901E-5 9 54 0.0848 0.0848715515617263 0.107 0.106908019877197 0 2.04285610764949E-5 10 60 0.0721 0.0726606420421092 0.1197 0.119134999676661 0 4.35828122956682E-6 |

鐵蟲 (小有名氣)
|
你好,這里求出來的是全局最優(yōu)解還是局部最優(yōu)解,因?yàn)槲矣貌煌某踔禂M合也能收斂,只是每次求出來的參數(shù)可能都不一樣如(k1 = 17.932955894914, k2 = 801.597983738672, kA = 8.49463586176274, KB = 0.0324248271839244, kc = 2.11449075468859, kH2 = 0.537868687305816)。另外,可以加上帶5%誤差分布的擬合參數(shù)嗎,是用MATLAB計(jì)算的嗎,可以的話還想求一份計(jì)算的代碼,麻煩你了,我可以向你付費(fèi)咨詢的。 |
版主 (知名作家)
|
這是用1stopt擬合的,采用通用全局優(yōu)化算法,至于每次結(jié)果有些許差別那是正常的,只要總體趨勢滿足,擬合的統(tǒng)計(jì)參數(shù)符合要求就可以了。擬合只是尋找理想的最優(yōu)解,結(jié)果往往只有更優(yōu),沒有最優(yōu)。 發(fā)自小木蟲Android客戶端 |

鐵蟲 (小有名氣)
至尊木蟲 (知名作家)
|
https://ww2.mathworks.cn/help/matlab/ref/polyfit.html https://ww2.mathworks.cn/products/curvefitting/features.html https://zhuanlan.zhihu.com/p/21754490 |
鐵蟲 (小有名氣)
鐵蟲 (小有名氣)
鐵蟲 (小有名氣)
|
大哥,數(shù)據(jù)和模型發(fā)了,你抽空幫我用程序跑一下吧,非常感激。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +8 | 土土小蟲 2026-03-03 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 材料復(fù)試調(diào)劑 +7 | 學(xué)材料的點(diǎn) 2026-03-01 | 8/400 |
|
|
[考研] 080500材料科學(xué)與工程 +8 | 202114020319 2026-03-03 | 8/400 |
|
|
[考研] 282求調(diào)劑 +5 | 2103240126 2026-03-02 | 8/400 |
|
|
[考研] 0855機(jī)械工程324求調(diào)劑 +3 | 鄭SK 2026-03-02 | 7/350 |
|
|
[考研] 學(xué)碩材料275調(diào)劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +6 | xiiiia 2026-03-04 | 6/300 |
|
|
[考研] 一志愿中科大080500總分324求調(diào)劑 +3 | jorna 2026-03-03 | 6/300 |
|
|
[考研] 293求調(diào)劑 +4 | 是樂渝哇 2026-03-03 | 4/200 |
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
[考研] 085700資環(huán)求調(diào)劑,初始279,六級(jí)已過,英語能力強(qiáng) +3 | 085700資環(huán)調(diào)劑 2026-03-03 | 4/200 |
|
|
[考研] 一志愿東北大學(xué)化學(xué)314分求調(diào)劑 +5 | lr1212.. 2026-03-02 | 5/250 |
|
|
[考研] 化工專碩348,一志愿985求調(diào)劑 +8 | 弗格個(gè) 2026-02-28 | 11/550 |
|
|
[考研] 267求調(diào)劑 +6 | 釣魚佬as 2026-03-02 | 6/300 |
|
|
[考研] 材料類考研調(diào)劑 +6 | gemmgemm 2026-03-01 | 7/350 |
|
|
[考博] 博士自薦 +4 | kkluvs 2026-02-28 | 5/250 |
|
|
[考研] 281求調(diào)劑 +5 | 2026計(jì)算機(jī)_誠心 2026-03-01 | 8/400 |
|
|
[考研] 274求調(diào)劑 +3 | cgyzqwn 2026-03-01 | 7/350 |
|
|
[考研] 311求調(diào)劑 +9 | 南迦720 2026-02-28 | 10/500 |
|
|
[論文投稿]
Optics letters投稿被拒求助
30+3
|
luckyry 2026-02-26 | 4/200 |
|