| 查看: 920 | 回复: 4 | ||
m512990444新虫 (初入文坛)
|
[求助]
求助1 stopt 求出参数，已知编程代码
|
|
Parameters k1[0,],k2[0,],k3[0,],k4[0,],k5[0,],k6[0,]; Constant E0=1.729; Variable t,F,S,M,GOS3,GOS4,GOS5; ODEFunction S'=-(k1*k2*(k3*M+k4*GOS3+k5*GOS4)*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S); M'=-(k1*k2*k3*S*M*E0)/((k1+k2+k6 )*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S); F'=(k1*k2*k8*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S); GOS3'=((k3*M-k4*GOS3)*k1*k2*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S); GOS4'=((k4*GOS3-k5*GOS4)*k1*k2*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S); GOS5'=(k5*GOS4*k1*k2*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S); Data; 0 0 146.64 78.45 0 0 0 60 16.85 129.79 65 10.06 1.53 0 240 53.69 92.95 46.97 22.65 12.39 0 480 83.3 63.34 34.79 25.78 23.75 3.02 660 107.47 39.17 29 25.78 32.36 5.04 Data; 0 0 147.99 148.86 0 0 0 60 22.79 125.2 128.16 16.25 2.88 0 240 68.22 79.77 99.47 41 12.53 0.65 480 99.91 48.08 81.77 51.47 25.02 2.53 660 129.82 18.17 66.9 57.22 35.48 4.85 Data; 0 0.00 153.17 314.62 0.00 0.00 0.00 60 21.54 122.86 292.67 20.46 2.15 0.00 240 74.42 69.98 235.98 58.22 15.92 0.00 480 111.48 32.92 202.04 82.67 25.45 1.02 660 132.25 9.23 177.58 91.32 35.57 2.90 |
回复此楼» 猜你喜欢
光学工程学硕调剂信息
已经有32人回复
-
欢迎加入课题组
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有272人回复
-
散金币，求好运，祝面上顺利！
已经有78人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
m512990444: 金币+10, ★有帮助 2016-03-24 10:13:59
m512990444: 金币+10, ★有帮助 2016-03-24 11:52:05
感谢参与，应助指数 +1
m512990444: 金币+10, ★有帮助 2016-03-24 10:13:59
m512990444: 金币+10, ★有帮助 2016-03-24 11:52:05
|
供参考： 均方差(RMSE): 13.9366245195253 残差平方和(SSE): 13984.5242158729 相关系数(R): 0.981929994660137 相关系数之平方(R^2): 0.964186514413256 确定系数(DC): 0.730266057946131 F统计(F-Statistic): -0.656005330812249 参数 最佳估算 -------------------- ------------- k8 0.136195429337862 k1 0.0202551743597175 k2 0.0302696407861668 k3 0.000850831443242199 k4 0.00106431484597736 k5 0.000284975127225702 k6 0.000333626838395123 |
2楼2016-03-21 16:03:20
m512990444
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 338.6
- 帖子: 27
- 在线: 20.7小时
- 虫号: 1988090
- 注册: 2012-09-09
- 性别: MM
- 专业: 食品科学基础
|
F'=(k1*k2*k8*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S)里面的k8不是参数，是 k3*M+k4*GOS3+k5*GOS4，我之前简写了，没有改过来， 正确的应该是F'=(k1*k2* (k3*M+k4*GOS3+k5*GOS4) *S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S)，能再帮忙拟合下吗？ 发自小木虫Android客户端 |
3楼2016-03-21 21:11:41
m512990444
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 338.6
- 帖子: 27
- 在线: 20.7小时
- 虫号: 1988090
- 注册: 2012-09-09
- 性别: MM
- 专业: 食品科学基础
|
不好意思，F'=(k1*k2*k8*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S)里面的k8不是参数，是 k3*M+k4*GOS3+k5*GOS4，我之前简写了，没有改过来， 正确的应该是F'=(k1*k2* (k3*M+k4*GOS3+k5*GOS4) *S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S)，能再帮忙拟合下吗？ 发自小木虫Android客户端 |
4楼2016-03-22 21:21:26
m512990444
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 338.6
- 帖子: 27
- 在线: 20.7小时
- 虫号: 1988090
- 注册: 2012-09-09
- 性别: MM
- 专业: 食品科学基础
5楼2016-03-24 16:17:44
