24小时热门版块排行榜    

查看: 1275  |  回复: 3

cigarette

银虫 (小有名气)

[求助] 求高版本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
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★
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

2楼2012-11-30 14:16:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forcal

禁虫 (初入文坛)


csgt0: 金币+1, 应助指数+1, 谢谢 2012-12-03 11:10:54
本帖内容被屏蔽

3楼2012-12-01 07:15:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forcal

禁虫 (初入文坛)

★ ★
csgt0: 金币+2, 谢谢 2012-12-03 11:11:03
本帖内容被屏蔽

4楼2012-12-01 07:21:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cigarette 的主题更新
信息提示
请填处理意见