| 查看: 2629 | 回复: 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 |
回复此楼» 猜你喜欢
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有265人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有19人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有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