文档介绍：H7N9禽流感传播规律模型研究
小组成员: 韩萌朱腾郭轩赫
摘要
本文对H7N9流感病毒的传播规律进行了研究和预测,文章收集了我国各地区H7N9流感病毒感染人数及死亡人数的确切数据,对模型进行了验证,并提出了如何降低流感发病率的可靠方法。
一、问题重述
近期以来h7n9流感病毒正成为人们关注的焦点,通过相关网站获得数据,建立一个模型对流感的走势进行预测。
二、问题分析
h7n9流感的传播是一道传染病问题。在数学建模领域已经有很多关于这方面的研究。本文选取我国各地区的感染者确认病例数量,对h7n9流感的传播规律进行预测。
三、建立模型
1、模型假设
在h7n9流感传播期内,中国境内的总人数不变,既不考虑生死,也不考虑迁移。
人类对h7n9的控制按照正常计划进行,中间不出现任何意外情况,没有其他重大传染病的干扰及政治、经济因素的干扰。
病人之间的接触对感染病的传播没有影响,日治愈率(每天被治愈的病人占总病人数的比例)为常数。为该模型的缺陷是结果常与实际有一定程度差距,这是因为模型中假设有效接触率传染力是不变的。
模型构成
以时间x为自变量,感染人数y1和死亡人数y2为因变量,建立以下多项式模型:a、b为系数,c1、c2为常数。
Y1=a1x^n+a2x^(n-1)+a3x^(n-2)+a……+anx^n+c1
Y2=b1x^n+b2x^(n-1)+b3x^(n-2)+b……+bnx^n+c2
模型求解
用matlab对模型进行求解:
1:多项式求解的符号:
Y1=a1x^n+a2x^(n-1)+a3x^(n-2)+a……+anx^n+c1(令n=10)
Y2=b1x^n+b2x^(n-1)+b3x^(n-2)+b……+bnx^n+c2(令n=10)
X=time()(说明:time()是将时间转换成数字的函数)
function [x]=time(a,b,c)%%一天为单位的,计算与2013/4/5的时间差.%%以2013/4/5为准
if(a<2013 )
disp('你输入的数据不符合要求!!!');
elseif(b<1||b>12||(a==2013&&b<4))
disp('你输入的数据不符合要求!!!');
elseif (c>31||c<1||(a==2013&&b==4&&c<5))
disp('你输入的数据不符合要求!!!');
else
switch b %%1;3;5;7;8;10;12
case 1
z=31-c;
case 2
z=c;
case 3
z=28+c;
case 4
z=28+31+c;
case 5
z=28+31+30+c;
case 6
z=28+31+30+31+c;
case 7
z=28+31+30+31+30+c;
case 8
z=28+31+30+31+30+31+c;
case 9
z=28+31+30+31+30+31+31+c;
case 10
z=28+31+30+31+30+31+31+30+c;
case 11
z=28+31+30+31+30+31+31+30+31+c;
case 12
z=28+