 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2214 | 回復(fù): 11 | |||
蘇拉鐵蟲 (小有名氣)
|
[求助]
matlab擬合反應(yīng)動力學(xué)參數(shù) 已有2人參與
|
||
|
各位大佬，可以幫我擬合一下動力學(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 |
回復(fù)此樓» 猜你喜歡
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語口語老師找學(xué)生
已經(jīng)有4人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有87人回復(fù)
-
計算三階力常數(shù)用thirdorder擴(kuò)胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(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)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(fù)
獨(dú)孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機(jī)械動力學(xué)
- 管轄: 計算模擬
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
均方差(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)計(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 計算 b 目標(biāo) c 計算 c 目標(biāo) a 計算 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 |
2樓2019-03-23 20:10:57
蘇拉
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 247.9
- 帖子: 77
- 在線: 14.5小時
- 蟲號: 6170898
- 注冊: 2017-03-28
- 性別: GG
- 專業(yè): 能源化工
|
你好，這里求出來的是全局最優(yōu)解還是局部最優(yōu)解，因為我用不同的初值擬合也能收斂，只是每次求出來的參數(shù)可能都不一樣如(k1 = 17.932955894914， k2 = 801.597983738672, kA = 8.49463586176274, KB = 0.0324248271839244, kc = 2.11449075468859, kH2 = 0.537868687305816）。另外，可以加上帶5%誤差分布的擬合參數(shù)嗎，是用MATLAB計算的嗎，可以的話還想求一份計算的代碼，麻煩你了，我可以向你付費(fèi)咨詢的。 |
3樓2019-03-24 10:18:32
獨(dú)孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機(jī)械動力學(xué)
- 管轄: 計算模擬
【答案】應(yīng)助回帖
|
這是用1stopt擬合的，采用通用全局優(yōu)化算法，至于每次結(jié)果有些許差別那是正常的，只要總體趨勢滿足，擬合的統(tǒng)計參數(shù)符合要求就可以了。擬合只是尋找理想的最優(yōu)解，結(jié)果往往只有更優(yōu)，沒有最優(yōu)。 發(fā)自小木蟲Android客戶端 |
» 本帖已獲得的紅花（最新10朵）
蘇拉4樓2019-03-24 13:03:01
5樓2019-03-24 13:19:11
蘇拉
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 247.9
- 帖子: 77
- 在線: 14.5小時
- 蟲號: 6170898
- 注冊: 2017-03-28
- 性別: GG
- 專業(yè): 能源化工
6樓2019-03-24 16:17:01
FMStation
至尊木蟲 (知名作家)
- 應(yīng)助: 591 (博士)
- 貴賓: 0.03
- 金幣: 18311.2
- 紅花: 97
- 帖子: 8891
- 在線: 1489.3小時
- 蟲號: 2400059
- 注冊: 2013-04-04
- 專業(yè): 計算機(jī)應(yīng)用技術(shù)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
https://ww2.mathworks.cn/help/matlab/ref/polyfit.html https://ww2.mathworks.cn/products/curvefitting/features.html https://zhuanlan.zhihu.com/p/21754490 |
7樓2019-03-24 18:22:24
蘇拉
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 247.9
- 帖子: 77
- 在線: 14.5小時
- 蟲號: 6170898
- 注冊: 2017-03-28
- 性別: GG
- 專業(yè): 能源化工
8樓2019-03-25 09:51:46
蘇拉
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 247.9
- 帖子: 77
- 在線: 14.5小時
- 蟲號: 6170898
- 注冊: 2017-03-28
- 性別: GG
- 專業(yè): 能源化工
9樓2019-03-25 11:04:13
蘇拉
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 247.9
- 帖子: 77
- 在線: 14.5小時
- 蟲號: 6170898
- 注冊: 2017-03-28
- 性別: GG
- 專業(yè): 能源化工
|
大哥，數(shù)據(jù)和模型發(fā)了，你抽空幫我用程序跑一下吧，非常感激。 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : Desktop.rar
2019-03-25 15:59:38, 24.88 K
10樓2019-03-25 16:03:01
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 蘇拉 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 315食品工程求調(diào)劑 +3 | Oreov0 2026-03-06 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 家人們 調(diào)劑不迷路 看這里 +8 | likeihood 2026-03-09 | 10/500 |
|
|
|
[考研]
 295復(fù)試調(diào)劑
+4
|
簡木ChuFront 2026-03-09 | 4/200 |
|
|
|
[考研] 一志愿鄭州大學(xué)，學(xué)碩，物理化學(xué)， 333求調(diào)劑 +7 | 李魔女斗篷 2026-03-04 | 8/400 |
|
|
|
[考研] 320材料與化工，求調(diào)劑 +11 | 鶴遨予卿 2026-03-04 | 13/650 |
|
|
|
[考研] 327分求調(diào)劑086 +3 | 西紅柿?小帥 2026-03-09 | 6/300 |
|
|
|
[考研] 269求調(diào)劑 +3 | 朔朔話 2026-03-08 | 4/200 |
|
|
|
[考研] 347求調(diào)劑 +4 | 浮云滿足 2026-03-07 | 4/200 |
|
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
|
[考研] 材料工程專碩調(diào)劑自薦信，初試323 +4 | 李白26 2026-03-07 | 4/200 |
|
|
|
[考研] 求調(diào)劑，不管什么專業(yè)，我是可塑造的人才一枚，希望遇到知己老師撈撈我 +4 | 13102137290 2026-03-06 | 5/250 |
|
|
|
[考研] 化學(xué)290求調(diào)劑 +3 | 兩顆 西柚 2026-03-07 | 4/200 |
|
|
|
[考研] 0856材料與化工求調(diào)劑！ +5 | 化工考生111 2026-03-04 | 11/550 |
|
|
|
[考研] 282求調(diào)劑 +7 | 夕～日 2026-03-05 | 8/400 |
|
|
|
[考研] 274求調(diào)劑 +9 | 一個學(xué)習(xí)者 2026-03-04 | 9/450 |
|
|
|
[考研] 紡織、生物、化學(xué)、材料等專業(yè) +3 | Eember. 2026-03-05 | 7/350 |
|
|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級已過 求調(diào)劑 +5 | ccc! 2026-03-03 | 6/300 |
|
|
|
[考研] 材料085600 303求調(diào)劑 +7 | 1bygone 2026-03-04 | 7/350 |
|
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
|
[考研] 281求調(diào)劑 +3 | 我是小小蔥蔥 2026-03-03 | 5/250 |
|