 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 746 | 回復(fù): 1 | ||
追星的吐司金蟲 (小有名氣)
|
[求助]
1stopt微分方程動(dòng)力學(xué)參數(shù)擬合問(wèn)題，急求高手指點(diǎn)�。�！
|
|
1stopt微分方程參數(shù)擬合問(wèn)題，以下是我的運(yùn)行代碼： ODEAlgorithm = RKF45; Parameters k(6)[0,10000]; Variable x,y(7); ODEStep = 1; ODEFunction y1'=-0.01*k1*y1; y2'=0.01014*k1*y1-0.01*k2*y2; y3'=0.006389*k2*y2+0.00697*k5*y6-0.01*k3*y3; y4'=0.010217*k3*y3-0.01*k4*y4; y5'=0.004681*k4*y4+0.004783*k2*y2+0.003333*k5*y6+0.003099*k6*y7; y6'=-0.01*k5*y6; y7'=-0.01*k6*y7; Data; 0,0.415,0,0,0,0,0.301,0.284 0.4174,NAN,0.0000,0.3415,0.0021,0.6131,NAN,NAN 0.4826,NAN,0.0000,0.3733,0.0023,0.5855,NAN,NAN 0.5720,NAN,0.0000,0.3329,0.0012,0.6356,NAN,NAN 0.6176,NAN,0.0000,0.3268,0.0016,0.6527,NAN,NAN 0.7765,NAN,0.0000,0.2512,0.0012,0.7368,NAN,NAN 1.0455,NAN,0.0000,0.0672,0.0019,0.9201,NAN,NAN 擬合結(jié)果k4特別大，不符合事實(shí)，想知道問(wèn)題如何解決 運(yùn)行結(jié)果： 模型公式: y1' = -0.01*k1*y1 y1' = -0.01*1788.61544280086*y1 模型公式: y2' = 0.01014*k1*y1-0.01*k2*y2 y2' = 0.01014*1788.61544280086*y1-0.01*1789.29764640257*y2 模型公式: y3' = 0.006389*k2*y2+0.00697*k5*y6-0.01*k3*y3 y3' = 0.006389*1789.29764640257*y2+0.00697*1317.49753375156*y6-0.01*167.08336428024*y3 模型公式: y4' = 0.010217*k3*y3-0.01*k4*y4 y4' = 0.010217*167.08336428024*y3-0.01*9999.96896608683*y4 模型公式: y5' = 0.004681*k4*y4+0.004783*k2*y2+0.003333*k5*y6+0.003099*k6*y7 y5' = 0.004681*9999.96896608683*y4+0.004783*1789.29764640257*y2+0.003333*1317.49753375156*y6+0.003099*4724.79409921573*y7 模型公式: y6' = -0.01*k5*y6 y6' = -0.01*1317.49753375156*y6 模型公式: y7' = -0.01*k6*y7 y7' = -0.01*4724.79409921573*y7 迭代數(shù): 114 常微分方程算法: 龍格-庫(kù)塔-費(fèi)爾博格法(Runge-Kutta-Fehlberg Method) 優(yōu)化算法: 通用全局優(yōu)化算法(UGO1) 計(jì)算結(jié)束原因: 達(dá)到收斂判斷標(biāo)準(zhǔn) 計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:01:54:998 均方差(RMSE): 0.0810597491269555 殘差平方和(SSE): 0.275968682998049 參數(shù) 最佳估算 -------------------- ------------- k1 1788.61544280086 k2 1789.29764640257 k3 167.08336428024 k4 9999.96896608683 k5 1317.49753375156 k6 4724.79409921573 |
回復(fù)此樓» 猜你喜歡
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有279人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院
已經(jīng)有29人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院
已經(jīng)有26人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
UJN物理學(xué)專業(yè)調(diào)劑
已經(jīng)有5人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復(fù)
-
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名
已經(jīng)有0人回復(fù)
追星的吐司
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1563.4
- 散金: 15
- 紅花: 2
- 帖子: 214
- 在線: 146.1小時(shí)
- 蟲號(hào): 3777120
- 注冊(cè): 2015-03-31
- 專業(yè): 法學(xué)其他學(xué)科
2樓2018-01-16 10:50:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 追星的吐司 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 366求調(diào)劑一志愿東北大學(xué) +5 | 運(yùn)氣來(lái)得若有似?/a> 2026-04-02 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑！生物與醫(yī)藥專碩 +3 | 逆轉(zhuǎn)陸先生 2026-04-01 | 3/150 |
|
|
[考研] 一志愿346上海大學(xué)生物學(xué) +3 | 上海大學(xué)346調(diào)劑 2026-04-01 | 3/150 |
|
|
|
[考研] 材料工程322分 +7 | 哈哈哈吼吼吼哈 2026-04-01 | 7/350 |
|
|
|
[考研] 安全工程 285 求調(diào)劑 +3 | Xinyu56 2026-04-01 | 4/200 |
|
|
|
[考研] 070300化學(xué)279求調(diào)劑 +15 | 哈哈哈^_^ 2026-03-31 | 17/850 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +21 | 荷包蛋hyj 2026-03-26 | 21/1050 |
|
|
|
[考研] 303分 0807學(xué)碩求調(diào)劑 +3 | TYC3632 2026-04-01 | 3/150 |
|
|
|
[考研] 085600 一志愿9 總分351 求調(diào)劑學(xué)校 +7 | czhcz 2026-03-31 | 9/450 |
|
|
|
[考研] 070300求調(diào)劑306分 +5 | 26要上岸 2026-03-27 | 5/250 |
|
|
|
[考研] 一志愿華南師范361分，化學(xué)求調(diào)劑 +4 | Nicole88888 2026-04-01 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +4 | 阿科逸 2026-03-31 | 4/200 |
|
|
|
[考研] 調(diào)劑申請(qǐng) +8 | 張張張張zy 2026-03-31 | 9/450 |
|
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +13 | 深V宿舍吧 2026-03-29 | 13/650 |
|
|
|
[考研] 0856求調(diào)劑 +9 | 楒桉 2026-03-28 | 9/450 |
|
|
|
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +5 | 小卡不卡. 2026-03-29 | 5/250 |
|
|
|
[考研] 071010 323 分求調(diào)劑 +3 | Baekzhy 2026-03-27 | 3/150 |
|
|
|
[考研] 312，生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|