 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 9 | 1/1 | 返回列表 |
| 查看: 1945 | 回復(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ì)功能材料國際會(huì)議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有214人回復(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)問還有沒有用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)問題，其責(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, ★★★很有幫助, 遺傳算法我也試過了 R2能到0.93 我再試試1stopt 非常感謝 2015-03-10 22:30:32
winterhao: 金幣+50, ★★★很有幫助, 遺傳算法我也試過了 R2能到0.93 我再試試1stopt 非常感謝 2015-03-10 22:30:32
|
您好，m文件我看了，程序是可以運(yùn)行的，計(jì)算結(jié)果不太好。 我覺得還是初值的問題，4個(gè)k初值之間的數(shù)量級(jí)相差太大，MATLAB在比較復(fù)雜的常微分方程擬合方面，初值的選取一直是一個(gè)問題。 如果遺傳算法效果還是不理想，試試用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è)帖子說 如果兩個(gè)參數(shù)差三個(gè)數(shù)量級(jí)以上 不是大的不動(dòng) 就是小的不動(dòng) 不可能找到最優(yōu) 這種級(jí)數(shù)的差距即使1stopt可能也解決不了 我現(xià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ā)表 | |
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 +3 | 何潤(rùn)采123 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +8 | SUICHILD 2026-03-12 | 8/400 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +4 | 芬達(dá)46 2026-03-12 | 4/200 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
|
[論文投稿]
投稿問題
 5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|