| 查看: 2199 | 回復(fù): 11 | |||
蘇拉鐵蟲(chóng) (小有名氣)
|
[求助]
matlab擬合反應(yīng)動(dòng)力學(xué)參數(shù) 已有2人參與
|
||
|
各位大佬,可以幫我擬合一下動(dòng)力學(xué)參數(shù)嗎,可以有償。先謝過(guò)各位大佬。 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 |

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

鐵蟲(chóng) (小有名氣)
至尊木蟲(chóng) (知名作家)
|
https://ww2.mathworks.cn/help/matlab/ref/polyfit.html https://ww2.mathworks.cn/products/curvefitting/features.html https://zhuanlan.zhihu.com/p/21754490 |
鐵蟲(chóng) (小有名氣)
鐵蟲(chóng) (小有名氣)
鐵蟲(chóng) (小有名氣)
|
大哥,數(shù)據(jù)和模型發(fā)了,你抽空幫我用程序跑一下吧,非常感激。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 283求調(diào)劑 +10 | 鹿沫笙 2026-03-02 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 材料與化工304求調(diào)劑 +3 | 邱gl 2026-03-05 | 6/300 |
|
|
[考研] 材料328求調(diào)劑 +10 | 一個(gè)蘿卜02 2026-03-03 | 10/500 |
|
|
[考研] 376求調(diào)劑 +3 | 王浩然sam 2026-03-04 | 3/150 |
|
|
[考研] 調(diào)劑材料學(xué)碩 +6 | 詞凝Y 2026-03-02 | 6/300 |
|
|
[考研] 材料復(fù)試調(diào)劑 +7 | 學(xué)材料的點(diǎn) 2026-03-01 | 8/400 |
|
|
[考研] 070300調(diào)劑一志愿南京大學(xué)化學(xué)專(zhuān)業(yè),初試273分,有論文有專(zhuān)利有競(jìng)賽,求調(diào)劑B區(qū)211 +3 | 82206202 2026-02-28 | 8/400 |
|
|
[考研] 材料專(zhuān)碩346求調(diào)劑 +3 | 旺一下 2026-03-04 | 3/150 |
|
|
[考研] 281求調(diào)劑 +3 | 我是小小蔥蔥 2026-03-03 | 5/250 |
|
|
[考研] 一志愿西工大材料學(xué),英一數(shù)二總分321分,求調(diào)劑。 +3 | zz05zz 2026-03-02 | 5/250 |
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
[考研] 266材料化工求調(diào)劑 +3 | 哇塞王帥 2026-03-03 | 3/150 |
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
[考研] 275求調(diào)劑 +7 | 明遠(yuǎn)求學(xué) 2026-03-01 | 7/350 |
|
|
[考研] 265分求調(diào)劑不調(diào)專(zhuān)業(yè)和學(xué)校有行學(xué)上就 +6 | 禮堂丁真258 2026-02-28 | 9/450 |
|
|
[考研] 295求調(diào)劑 +8 | 19171856320 2026-02-28 | 8/400 |
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|
|
[基金申請(qǐng)] 成果系統(tǒng)訪問(wèn)量大,請(qǐng)一小時(shí)后再?lài)L試。---NSFC啥時(shí)候好哦,已經(jīng)兩天這樣了 +4 | NSFC2026我來(lái)了 2026-02-28 | 4/200 |
|