| 查看: 720 | 回复: 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 |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有141人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
追星的吐司
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1563.4
- 散金: 15
- 红花: 2
- 帖子: 214
- 在线: 146.1小时
- 虫号: 3777120
- 注册: 2015-03-31
- 专业: 法学其他学科
2楼2018-01-16 10:50:33