| 查看: 1509 | 回复: 4 | |||
| 【悬赏金币】回答本帖问题，作者happy_rolly将赠送您 5 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
happy_rolly新虫 (初入文坛)
|
[求助]
用matlab求动力学参数代码 求指教已有1人参与
|
||
|
附件中的文章 跟我的研究方向比较相关 所以想重现一下作者的结果, 首先建立传染病动力学SEIQR模型（分别表示易感者、潜伏期、感染未住院、感染并住院、康复），然后写出了模型的常微分方程组，共有8个参数，能给出数值的有μ和δ，另外六个参数给出了初始值范围，并且需要估计 S E I Q R五个未知数中，S Q R的初始值可以给出，E I是给定范围 ，要估计 最优参数的估计采用的是fminsearch，目标函数是MSS=(log2 (1+data per month) 2-log2 (1+ infectious and hospitalized on day 1+...+on day 28(or 29,30,31)))2 目标函数表达式中data per month就是表1中 每月发病人数/10, 另外 infectious and hospitalized on day 1+...+on day 28(or 29,30,31)应该是来源于常微分方程组的解 目前为止我只写出了常微分方程组的代码，见li.m 至于如何构建MSS这样的目标函数，以及如何利用常微分方程组的解和每月发病数据来求MSS最小值 并且要得出另外六个参数值和E I的初始值 我一点头绪也没有 求大神指教 必有重谢 最好能给出matlab代码 医学生学matlab好困难，整个课题组就我一个人在做这个 *li.m function f=li(t,y) u=3.9139*10^-5; f=[u*(y(1)+y(2)+y(3)+y(4)+y(5))-b*(y(3)+y(4))*y(1)/(y(1)+y(2)+y(3)+y(4)+y(5))+r*y(5)-u*y(1); b*(y(3)+y(4))*y(1)/(y(1)+y(2)+y(3)+y(4)+y(5))-a*y(2)-u*y(2); a*(1-p)*y(2)-m*y(3)-u*y(3); a*p*y(2)-n*y(4)-u*y(4); m*y(3)+n*y(4)-r*y(5)-u*y(5)];     常微分方程组.png  表1 每月发病数据.png |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com - 附件 1 : Li-2014-Modeling_and_Prevent.pdf
2018-04-18 23:20:19, 315.08 K
» 收录本帖的淘贴专辑推荐
 常微分方程 |
» 猜你喜欢
Error in routine lambda (100): wrong # or too many modes
已经有0人回复
-
祈福国自然面上
已经有61人回复
-
物理学I论文润色/翻译怎么收费?
已经有143人回复
-
请问castep计算出的能带图横轴数据缩放因子怎么计算
已经有1人回复
-
Result in Physics 从见刊开始 一般多久可以被SCI收录
已经有26人回复
-
面青地会评时间？？？
已经有8人回复
-
科研求助
已经有1人回复
-
深圳大学高等研究院暑期夏令营
已经有10人回复
-
深圳大学高等研究院暑期夏令营
已经有8人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15356.8
- 散金: 101
- 红花: 234
- 帖子: 3409
- 在线: 1222.7小时
- 虫号: 291104
- 注册: 2006-10-28
3楼2018-04-20 17:09:41
20