| 5 | 1/1 | 返回列表 |
| 查看: 2209 | 回復(fù): 11 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
蘇拉鐵蟲 (小有名氣)
|
[求助]
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é)果有些許差別那是正常的,只要總體趨勢(shì)滿足,擬合的統(tǒng)計(jì)參數(shù)符合要求就可以了。擬合只是尋找理想的最優(yōu)解,結(jié)果往往只有更優(yōu),沒有最優(yōu)。 發(fā)自小木蟲Android客戶端 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)求調(diào)劑學(xué)校 +3 | 一志愿京區(qū)211 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 0703求調(diào)劑 +5 | c297914 2026-03-06 | 6/300 |
|
|
[考研] 0307化學(xué)求調(diào)劑 +6 | 0ok0k 2026-03-07 | 6/300 |
|
|
[考研] 288求調(diào)劑085600材料與化工 +9 | Daunrin 2026-03-07 | 10/500 |
|
|
[考研] 304求調(diào)劑 +4 | 52hz~~ 2026-03-05 | 5/250 |
|
|
[考研] 306求調(diào)劑 +7 | Bahati 2026-03-05 | 7/350 |
|
|
[考研] 一志愿211,0860總分286食品工程專業(yè)求調(diào)劑 +6 | 林林Winnie 2026-03-05 | 6/300 |
|
|
[考研] 324求調(diào)劑 +5 | wxz2 2026-03-03 | 8/400 |
|
|
[考研] 301求調(diào)劑 +7 | 朝天椒的雙馬尾 2026-03-05 | 7/350 |
|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來 +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|
|
[考研] 材料與化工354調(diào)劑 +4 | Lucy-xiao 2026-03-06 | 7/350 |
|
|
[考研] 332材料求調(diào)劑 +6 | zjy101327 2026-03-05 | 7/350 |
|
|
[考研] 282求調(diào)劑 +7 | 夕~日 2026-03-05 | 8/400 |
|
|
[考研] 化工282求調(diào)劑一志愿211 +5 | NA0912 2026-03-05 | 6/300 |
|
|
[考研] 347求調(diào)劑 +6 | 啊歐歐歐 2026-03-03 | 8/400 |
|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +7 | 1一11 2026-03-03 | 7/350 |
|
|
[考研] 一志愿314求調(diào)劑 +7 | 202111120625 2026-03-03 | 7/350 |
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|
|
[考研] 322,求調(diào)劑 +3 | 菜菜愛玩 2026-03-04 | 3/150 |
|