| 查看: 2705 | 回复: 3 | |||
[求助]
1st opt拟合方程——————求助！！ 已有1人参与
|
|
有高版本1st opt的大神帮忙跑下，谢谢！！！ Test1; Parameters a,b,c,d,e,f,g,h; SharedModel; Variable t,X,P,S; ODEFunction X'=a*X*(1-X/b)*(1+S/c)^(-1); P'=d*X'+e*X; S'=f*X'+g*P'+h*X; Data; t,X,P,S; 0 0.35 0 99.8 12 2.75 0.55 88.96 24 4.6 3.72 73.49 36 6.27 8.19 61.57 48 8.07 14.12 47.38 60 9.73 19.37 33.63 72 10.41 22.54 24.65 84 10.7 27.61 11.19 96 10.53 32.49 3.2 108 10.56 35.62 0 120 10.59 37.58 0 132 10.4 38.44 0 144 10.77 39.88 0 |
回复此楼» 猜你喜欢
光学工程学硕调剂信息
已经有32人回复
-
欢迎加入课题组
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有270人回复
-
散金币，求好运，祝面上顺利！
已经有79人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
皓小天
木虫之王 (职业作家)
- 应助: 19 (小学生)
- 贵宾: 0.1
- 金币: 76016.6
- 散金: 502
- 红花: 137
- 帖子: 3801
- 在线: 869小时
- 虫号: 1864745
- 注册: 2012-06-19
- 性别: GG
- 专业: 力学中的基本问题和方法
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zhl717: 金币+10, ★★★很有帮助 2015-09-06 09:53:54
zhl717: 金币+10, ★★★很有帮助 2015-09-06 09:53:54
|
首先这个不是sharemodel，而且你代码写的也有错。你的ode方程可以简化形式的，把第一个X’代入第二个ode方程，把第二个新的方程和第一个方程代入第三个方程，那么整个问题就简化成一个微分方程了。方程如下： S'=f*(a*X*(1-X/b)*(1+S/c)^(-1))+g*(d*(a*X*(1-X/b)*(1+S/c)^(-1))+e*X)+h*X 代码改成如下： Parameters a,b,c,d,e,f,g,h; Variable t,X,P,S; ODEFunction S'=f*(a*X*(1-X/b)*(1+S/c)^(-1))+g*(d*(a*X*(1-X/b)*(1+S/c)^(-1))+e*X)+h*X; Data; 0 0.35 0 99.8 12 2.75 0.55 88.96 24 4.6 3.72 73.49 36 6.27 8.19 61.57 48 8.07 14.12 47.38 60 9.73 19.37 33.63 72 10.41 22.54 24.65 84 10.7 27.61 11.19 96 10.53 32.49 3.2 108 10.56 35.62 0 120 10.59 37.58 0 132 10.4 38.44 0 144 10.77 39.88 0 结果如下： 均方差(RMSE):1.70125017832388 残差平方和(SSE):34.7310260309643 相关系数(R): 0.998920481960877 相关系数之平方(R^2): 0.997842129280951 确定系数(DC): 0.996942448980196 F统计(F-Statistic): 186.89093582654 参数 最佳估算 -------------------- ------------- a -0.307553362555249 b 11.7009599853419 c -0.0108960231582921 d 17354.8847310689 e -0.526630322610445 f 0.0237395701129598 g -0.824893382148287 h -0.326710215811974 ====== 结果输出 ====== 文件: 数据文件-1 No t 目标 s 计算 s 1 12 88.96 90.6377657133101 2 24 73.49 77.4840744651627 3 36 61.57 61.5302531902025 4 48 47.38 45.0285350315543 5 60 33.63 33.8481939424802 6 72 24.65 25.1096267702726 7 84 11.19 11.1348638405025 8 96 3.2 0.0109328172964804 9 108 0 0.0109000671739397 10 120 0 0.0109000671739397 11 132 0 0.0109000671739397 12 144 0 0.0109000671739397 |
2楼2015-08-30 22:34:38
3楼2015-09-06 10:24:28
4楼2015-09-06 10:49:50
