| 查看: 1656 | 回复: 2 | ||||
[求助]
应用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 将实际点图以及公式图像 放到一个图上。观察他的模拟结果情况。 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有82人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求回归方程参数的最优化算法matlab代码
已经有17人回复
- 求大神用matlab给我处理一组数据，有人没有啊？
已经有15人回复
- 怎样拟合数据求出x=f(A,B,C),y=g(A,B,C,D)中的A,B,C,D参数
已经有18人回复
- 通过实验数据拟合，求解公式中的参数
已经有14人回复
- 紧急求助，利用Matlab对实验数据进行拟合求解参数。
已经有27人回复
- 关于matlab的参数估计
已经有15人回复
- 重金请教-如何用MAtlab对力学模型进行参数识别
已经有12人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- matlab数据拟合问题求助
已经有7人回复
- 用MATLAB或者ORIGIN怎么拟合复数数据
已经有10人回复
- matlab求解方程中的参数
已经有21人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- matlab反求未知输入参数
已经有4人回复
- 已知一随机信号的功率谱密度函数，如何利用Matlab产生其单次时域样本
已经有6人回复
- matlab非线性参数拟合问题
已经有7人回复
- 【求助】MATLAB导入数据精度问题。。。
已经有4人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】问一个关于matlab中legend的问题
已经有11人回复
- 【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】
已经有15人回复
- 【求助】matlab中保存计算数据的问题
已经有8人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
xwndf250: 金币+15, ★有帮助, 我只想要matlab的的程序，不过还是要谢谢你。另外对于求参数问题，这个1stOpt有什么要求吗？比如像求解2个参数，至少知道几个已知参数为好？ 2014-03-24 10:06:34
感谢参与，应助指数 +1
xwndf250: 金币+15, ★有帮助, 我只想要matlab的的程序，不过还是要谢谢你。另外对于求参数问题，这个1stOpt有什么要求吗？比如像求解2个参数，至少知道几个已知参数为好？ 2014-03-24 10:06:34
|
1stOpt求解： 均方差(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
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
xwndf250: 金币+15, ★★★很有帮助, 还是1stopt...不过最起码给了个matlab 参考了 谢谢 2014-03-26 22:51:04
感谢参与，应助指数 +1
xwndf250: 金币+15, ★★★很有帮助, 还是1stopt...不过最起码给了个matlab 参考了 谢谢 2014-03-26 22:51:04
|
常微分方程组参数的拟合，MATLAB参见： http://muchong.com/bbs/viewthread.php?tid=7151341&authorid=1122189 在拟合方面，1stopt更强大哈。 |
3楼2014-03-24 14:28:45
30