| 查看: 739 | 回复: 1 | ||
追星的吐司金虫 (小有名气)
|
[求助]
1stopt微分方程动力学参数拟合问题，急求高手指点！！！
|
|
1stopt微分方程参数拟合问题，以下是我的运行代码： 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 拟合结果k4特别大，不符合事实，想知道问题如何解决 运行结果： 模型公式: 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 迭代数: 114 常微分方程算法: 龙格-库塔-费尔博格法(Runge-Kutta-Fehlberg Method) 优化算法: 通用全局优化算法(UGO1) 计算结束原因: 达到收敛判断标准 计算用时(时:分:秒:微秒): 00:01:54:998 均方差(RMSE): 0.0810597491269555 残差平方和(SSE): 0.275968682998049 参数 最佳估算 -------------------- ------------- k1 1788.61544280086 k2 1789.29764640257 k3 167.08336428024 k4 9999.96896608683 k5 1317.49753375156 k6 4724.79409921573 |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有118人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有26人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有34人回复
追星的吐司
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1563.4
- 散金: 15
- 红花: 2
- 帖子: 214
- 在线: 146.1小时
- 虫号: 3777120
- 注册: 2015-03-31
- 专业: 法学其他学科
2楼2018-01-16 10:50:33
