24小时热门版块排行榜    

查看: 1521  |  回复: 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的回帖

forcal

禁虫 (初入文坛)

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

4楼2012-12-01 07:21:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +4 myrtle 2026-08-26 4/200 2026-08-31 08:43 by 逸域隐帆
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 能否申诉? +6 echo8914667 2026-08-30 7/350 2026-08-31 08:41 by kingkocxr
[基金申请] 面上意见出来了 +10 黄鸟于飞Chao 2026-08-29 19/950 2026-08-31 08:01 by blueearth171
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见