24小时热门版块排行榜    

查看: 1624  |  回复: 2

xwndf250

银虫 (小有名气)

[求助] 应用matlab已知数据求参数问题 已有2人参与

流行病模型:
dI/dt=a*S*I-b*I
dS/dt=-a*S*I
dR/dt=b*I
已知
第1天:S=100 I=5 R=0
第2天:S=98  I=7 R=0
第3天:S=90 I=15 R=0
第4天:S=85 I=19 R=1
第5天:S=80 I=20 R=5
第6天:S=80 I=15 R=10
第7天:S=75 I=15 R=15
第8天:S=75 I=10 R=20
第9天:S=70 I=14 R=21
第10天:S=68 I=16 R=21
想求一下每天的参数 a b 请问怎么求解呢?(matlab)
另外,如何应用matlab 将实际点图以及公式图像 放到一个图上。观察他的模拟结果情况。
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
xwndf250: 金币+15, 有帮助, 我只想要matlab的的程序,不过还是要谢谢你。另外对于求参数问题,这个1stOpt有什么要求吗?比如像求解2个参数,至少知道几个已知参数为好? 2014-03-24 10:06:34
1stOpt求解:
CODE:
Variable t,S,I,R;
ODEFunction
I'=a*S*I-b*I;
S'=-a*S*I;
R'=b*I;
Data;
//t,S,I,R
1,100,5,0
2,98,7,0
3,90,15,0
4,85,19,1
5,80,20,5
6,80,15,10
7,75,15,15
8,75,10,20
9,70,14,21
10,68,16,21

均方差(RMSE):4.34441684362604
残差平方和(SSE):509.596858201905
相关系数(R): 0.672258204473979
相关系数之平方(R^2): 0.451931093482578
确定系数(DC): 0.147995676561136
F统计(F-Statistic): 34.0062120157514

参数                  最佳估算
--------------------        -------------
a        0.00461118762442189
b        0.252793800382281

====== 结果输出 ======

文件: 数据文件-1
No        t        目标 s        计算 s        目标 i        计算 i        目标 r        计算 r
1        2        98        97.4741145006773        7        6.12335991909717        0        1.40252558022548
2        3        90        94.4863088186581        15        7.40445719340253        0        3.10923398793936
3        4        85        91.0210096991185        19        8.82136082391591        1        5.15762947696557
4        5        80        87.0912508279901        20        10.3316165847728        5        7.57713258723717
5        6        80        82.7448628642643        15        11.8714279459937        10        10.383709189742
6        7        75        78.0658734111211        15        13.3592965817125        15        13.5748300071664
7        8        75        73.1693543000699        10        14.7046582202379        20        17.1259874796922
8        9        70        68.1898075656599        14        15.8202753349948        21        20.9899170993453
9        10        68        63.2654522619756        16        16.6354209350298        21        25.0991268029947
2楼2014-03-23 22:34:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
xwndf250: 金币+15, ★★★很有帮助, 还是1stopt...不过最起码给了个matlab 参考了 谢谢 2014-03-26 22:51:04
引用回帖:
2楼: Originally posted by dingd at 2014-03-23 22:34:57
1stOpt求解:

Variable t,S,I,R;
ODEFunction
I'=a*S*I-b*I;
S'=-a*S*I;
R'=b*I;
Data;
//t,S,I,R
1,100,5,0
2,98,7,0
3,90,15,0
4,85,19,1
5,80,20,5
6,80,15,10
7,75,15,15
8,75,10,20
9,70,14 ...

常微分方程组参数的拟合,MATLAB参见:
http://muchong.com/bbs/viewthread.php?tid=7151341&authorid=1122189
在拟合方面,1stopt更强大哈。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
3楼2014-03-24 14:28:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xwndf250 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 08:59 by 6vVgjDL4CnGu
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 08:27 by 6vVgjDL4CnGu
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 KXLV3nuBVBY7 2026-08-07 4/200 2026-08-08 08:21 by 6vVgjDL4CnGu
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见