 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 745 | 回復(fù): 1 | ||
追星的吐司金蟲 (小有名氣)
|
[求助]
1stopt微分方程動(dòng)力學(xué)參數(shù)擬合問題，急求高手指點(diǎn)�。�！
|
|
1stopt微分方程參數(shù)擬合問題，以下是我的運(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í)，想知道問題如何解決 運(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)有209人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有28人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有26人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
UJN物理學(xué)專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(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ā)表 | |
|---|---|---|---|---|
|
[考研] 385分 生物學(xué)（071000）求調(diào)劑 +3 | qf626 2026-04-01 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085600 295分求調(diào)劑 +15 | W55j 2026-03-30 | 18/900 |
|
|
|
[考研] 一志愿上海海洋大學(xué)083200食品學(xué)碩，求調(diào)劑，接受其他專業(yè) +5 | what張 2026-04-01 | 5/250 |
|
|
|
[考研] 26考研調(diào)劑 +3 | Wnz.20030617 2026-04-01 | 3/150 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +11 | 傾____蕭 2026-03-27 | 12/600 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +5 | 大饞小子 2026-03-27 | 7/350 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-30 | 7/350 |
|
|
|
[考研] 282求調(diào)劑 +6 | 呼吸都是減肥 2026-04-01 | 6/300 |
|
|
|
[考研] 309分085801求調(diào)劑 +7 | 學(xué)員Gtwj7W 2026-03-31 | 7/350 |
|
|
|
[考研] 311（085601）求調(diào)劑 +12 | liziyeyeye 2026-03-28 | 13/650 |
|
|
|
[考研] 一志愿西交大080500材料學(xué)碩349 +6 | jqx1258 2026-03-31 | 7/350 |
|
|
|
[考研] 求調(diào)劑 生物學(xué) 377分 +6 | zzll03 2026-03-31 | 6/300 |
|
|
|
[考研] 315求調(diào)劑 +6 | akie... 2026-03-28 | 7/350 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +16 | kikiki7 2026-03-30 | 16/800 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +25 | 稻妻小編 2026-03-26 | 25/1250 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +16 | 李李不服輸 2026-03-26 | 16/800 |
|
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|