[求助] 求高版本1stopt，解决一个问题，微分方程组参数拟合

dy0/dt=-(k1+k2)*y0+0.13*k2*y1;
dy1/dt=-k2*y0-0.13*k2*y1;

求k1，k2；

数据有：

t y0+y1 y0
1 0.962625393 0.108646918
2 0.844060774 0.101976843
4 0.733472444 0.074286987
18 0.377241439 0.050813504
24 0.316338229 0.052819116

回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:

麻烦大神帮忙用1stopt拟合得到未知参数已经有11人回复
求助大神，关于1stopt带积分的函数拟合已经有7人回复
这样的动力学方程能够有1stopt拟合出参数来么已经有4人回复
求助 1stopt低版本方程的拟合（含若干未知参数）已经有10人回复
动力学方程组参数求解问题已经有9人回复
求大哥大姐帮跑一下1stopt~新手啊，版本低啊已经有6人回复
1stopt自变量超多时是否有简单语法？AND 自变量多数据量大的话1stopt都会算的特慢么已经有7人回复
求助1stopt拟合动力学参数已经有4人回复
求大神帮忙拟合一个非线性方程，求出模型参数已经有15人回复
1stopt 计算问题已经有7人回复
求助，用1stopt解方程组，两次得到的解不唯一？什么原因，大神帮忙已经有6人回复
微分方程组求参数问题，求高人指教，感谢已经有11人回复
微分方程组参数拟合的问题（多参数动力系统参数估计）已经有11人回复
求大哥大姐帮跑一下1stopt~新手啊，版本低啊已经有9人回复
MATLAB微分方程参数拟合问题，求大神已经有7人回复
大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？已经有8人回复

1楼 2012-11-30 12:00:33

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

forcal

禁虫 (初入文坛)

★ ★
csgt0: 金币+2, 谢谢 2012-12-03 11:11:03

本帖内容被屏蔽

4楼2012-12-01 07:21:40

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

查看全部 4 个回答

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★
csgt0: 金币+3, 应助指数+1, 多谢应助 2012-11-30 14:34:09

先把数据“y0+y1”中的y1单独分离出来：

CODE:

Variable t,y1,y0;

ODEFunction y0'=-(k1+k2)*y0+0.13*k2*y1;

            y1'=-k2*y0-0.13*k2*y1;

Data;

//t        y1                y0

1        0.853978475        0.108646918

2        0.742083931        0.101976843

4        0.659185457        0.074286987

18        0.326427935        0.050813504

24        0.263519113        0.052819116

均方差(RMSE): 0.0328672483772618
残差平方和(SSE): 0.00864204812714095
相关系数(R): 0.998154499853223
相关系数之平方(R^2): 0.996312405577239
决定系数(DC): 0.99269172179015
F统计(F-Statistic): 49.4425008230623

参数最佳估算
-------------------- -------------
k1 0.195010932089137
k2 0.255123230578225

jv2.JPG

jv1.jpg