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

登录注册

应《网络安全法》要求，自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

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有227人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复

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

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

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 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 的主题更新

返回列表