版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

m512990444

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 338.6
帖子: 27
在线: 20.7小时
虫号: 1988090
注册: 2012-09-09
性别: MM
专业: 食品科学基础

[求助] 求助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

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有94人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有0人回复

1楼 2016-03-21 14:42:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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

供参考：

均方差(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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

m512990444

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 338.6
帖子: 27
在线: 20.7小时
虫号: 1988090
注册: 2012-09-09
性别: MM
专业: 食品科学基础

引用回帖:

2楼: Originally posted by dingd at 2016-03-21 16:03:20
供参考：

均方差(RMSE): 13.9366245195253
残差平方和(SSE): 13984.5242158729
相关系数(R): 0.981929994660137
相关系数之平方(R^2): 0.964186514413256
确定系数(DC): 0.730266057946131
F统计(F-Statist ...

不好意思，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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

m512990444

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 338.6
帖子: 27
在线: 20.7小时
虫号: 1988090
注册: 2012-09-09
性别: MM
专业: 食品科学基础

引用回帖:

大神，能不能帮忙再模拟一下，谢谢

赞一下

回复此楼

5楼2016-03-24 16:17:44

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 m512990444 的主题更新

返回列表