 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 9 | 1/1 | 返回列表 |
| 查看: 1943 | 回復(fù): 8 | ||||
winterhao鐵桿木蟲 (正式寫手)
|
[求助]
matlab常微分方程多參數(shù)擬合的參數(shù)之間cross correlations 求助@月只藍(lán)以及matlab高手 已有1人參與
|
|
用lsqnonlin進(jìn)行參數(shù)擬合后 怎樣求參數(shù) k1 k2 k3 之間的相關(guān)性 cross correlations lsqnonlin有這個(gè)功能嗎 還是需要我用別的函數(shù)去評(píng)價(jià) 我看了您以前列出的例子 如果能結(jié)合任何一個(gè)簡(jiǎn)單的例子的code給我講講 那真是太感謝了 |
回復(fù)此樓» 猜你喜歡
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有97人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有5人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有29人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有53人回復(fù)
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
2樓2015-03-10 08:40:23
winterhao
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 4492.7
- 散金: 20
- 紅花: 2
- 帖子: 504
- 在線: 187.2小時(shí)
- 蟲號(hào): 666604
- 注冊(cè): 2008-12-02
- 性別: GG
- 專業(yè): 流體力學(xué)
3樓2015-03-10 18:33:20
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
4樓2015-03-10 19:01:13
winterhao
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 4492.7
- 散金: 20
- 紅花: 2
- 帖子: 504
- 在線: 187.2小時(shí)
- 蟲號(hào): 666604
- 注冊(cè): 2008-12-02
- 性別: GG
- 專業(yè): 流體力學(xué)
5樓2015-03-10 19:27:31
winterhao
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 4492.7
- 散金: 20
- 紅花: 2
- 帖子: 504
- 在線: 187.2小時(shí)
- 蟲號(hào): 666604
- 注冊(cè): 2008-12-02
- 性別: GG
- 專業(yè): 流體力學(xué)
|
這是一個(gè)關(guān)于Mg(OH)2顆粒與水和CO2反應(yīng)生成 MgCO3晶體的 多相流反應(yīng) 包含了Mg(OH)2顆粒溶解模型 CO2吸收模型 以及MgCO3結(jié)晶模型 前兩個(gè)模型參數(shù)已經(jīng)確定 現(xiàn)在只需要用實(shí)驗(yàn)數(shù)據(jù) 擬合 MgCO3結(jié)晶模型的參數(shù) 我已經(jīng)寫好了 lsqnonlin的code 請(qǐng)見matlab文件 可以直接運(yùn)行 運(yùn)行發(fā)現(xiàn)： B = k1*(S^k3); % 成核速率 G = k2*(S^k4)；% 生長(zhǎng)速率 第一個(gè)參數(shù) k1 由于數(shù)量級(jí)很大（物理模型原因） 在擬合時(shí)變化很小 單獨(dú)擬合這個(gè)參數(shù)變動(dòng)也不大 擬合結(jié)果R2=0.90 我現(xiàn)在不確定是模型本身原因 還是我的初值選的不好的原因 我正在看您寫的遺傳算法 請(qǐng)您運(yùn)行一下我的文件 幫我找找原因 非常感謝 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : fitting_lsqnonlin_Crystallization_kinetics.m
2015-03-10 20:11:18, 13.29 K
6樓2015-03-10 20:12:09
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
winterhao: 金幣+50, ★★★很有幫助, 遺傳算法我也試過(guò)了 R2能到0.93 我再試試1stopt 非常感謝 2015-03-10 22:30:32
winterhao: 金幣+50, ★★★很有幫助, 遺傳算法我也試過(guò)了 R2能到0.93 我再試試1stopt 非常感謝 2015-03-10 22:30:32
|
您好，m文件我看了，程序是可以運(yùn)行的，計(jì)算結(jié)果不太好。 我覺得還是初值的問(wèn)題，4個(gè)k初值之間的數(shù)量級(jí)相差太大，MATLAB在比較復(fù)雜的常微分方程擬合方面，初值的選取一直是一個(gè)問(wèn)題。 如果遺傳算法效果還是不理想，試試用1stopt軟件，由于只有4個(gè)參數(shù)，你試試高版本的1stopt軟件吧：http://www.gaoyang168.com/bbs/viewthread.php?tid=8050245 |
7樓2015-03-10 22:01:33
winterhao
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 4492.7
- 散金: 20
- 紅花: 2
- 帖子: 504
- 在線: 187.2小時(shí)
- 蟲號(hào): 666604
- 注冊(cè): 2008-12-02
- 性別: GG
- 專業(yè): 流體力學(xué)
|
我有一個(gè)想法不知道可不可行 因?yàn)楣ぷ髁勘容^大 所以先和您討論一下 如果對(duì)參數(shù)和數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化或者歸一化 使他們都限制在{0,1} 或者{-1,1}會(huì)不會(huì)對(duì)參數(shù)的擬合有幫助？ 因?yàn)榭吹揭粋�(gè)帖子說(shuō) 如果兩個(gè)參數(shù)差三個(gè)數(shù)量級(jí)以上 不是大的不動(dòng) 就是小的不動(dòng) 不可能找到最優(yōu) 這種級(jí)數(shù)的差距即使1stopt可能也解決不了 我現(xiàn)在也不太明白 所以問(wèn)問(wèn)您 |
8樓2015-03-11 04:15:46
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
9樓2015-03-11 08:32:18
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ICPP2025 的主題更新
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專碩 +5 | 陌の森林 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 321求調(diào)劑 +3 | 何潤(rùn)采123 2026-03-18 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +8 | llllkkkhh 2026-03-18 | 8/400 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +6 | d如愿上岸 2026-03-12 | 9/450 |
|
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +3 | FY_yy 2026-03-14 | 3/150 |
|
|
|
[考研] 314求調(diào)劑 +8 | 無(wú)懈可擊的巨人 2026-03-12 | 8/400 |
|
|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂(lè) 2026-03-14 | 5/250 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +5 | 超級(jí)伊昂大王 2026-03-12 | 5/250 |
|
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。，。，。，。i 2026-03-12 | 4/200 |
|
|
|
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|