大学生数学建模有关传染病论文.docVIP

PAGE PAGE 1 2013上学期数学论 2013上学期数学论 姓 名：杨丽香、涂蓉 学 号：（02） 、（04） 学 院：湖南信息职业技术学院 专 业：计算机网络 指 导 教 师:祝文达 2013 年 06月 06日 邮箱：ZWD233@SINA.COM 传染病 一、摘要: 描述传染病的传播过程，分析受感染人数的变化规律，预报传染病高潮到来的时刻，预防传染病蔓延的手段，按照传播过程和一般规律，建立模型。利用了数学、力学、物理等学科中的定理来建立微分方程模型。利用已知的定理与规律寻找微元之间的关系式，与第一种方法不同的是对微元而不是直接对函数及其导数应用规律。在我们的实际问题中，许多现象的规律性并不是很清楚，如果有所了解也是极其复杂的，建模时在不同的假设中去模拟实际的现象，建立能近似反映问题的微积分方程，然后从数学上去求解或分析所建的方程及其解的性质，再去与实际情况相对比，检验此模型能否刻画模拟了某些实际现象。 二、问题重述 问题: 有一种传染病（如SARS、甲型H1N1）正在流行。现在希望建立适当的数学模型，利用已经掌握的一些数据资料对该传染病进行有效地研究，以期对其传播蔓延进行必要的控制，减少人民生命财产的损失。考虑如下的几个问题，建立适当的数学模型，并进行一定的比较分析和评价展望。 1、不考虑环境的限制，设单位时间内感染人数的增长率是常数，建立模型求t时刻的感染人数。 2、假设环境条件下所允许的最大可感染人数为 。单位时间内感染人数的增长率是感染人数的线性函数，最大感染时的增长率为零。建立模型求t时刻的感染人数。 3、现有卫生防疫部门采集到的某地区一定时间内一定间隔区间的感染人数数据（见下表），利用该数据确定上述两个模型中的相关参数，并将它们的预测值与实际数据进行比较分析（计算仿真偏差）并对两个模型进行适当的评价。（注：该问题中，设最大可感染人数为2000人） 4、假设总人口可分为传染病患者和易感染者，易感染者因与患病者接触而得病，而患病者会因治愈而减少且对该传染病具有很强的免疫功能，建立模型分析t时刻患病者与易感染者的关系，并对传染情况（如流行趋势，是否最终消灭）进行预测。 三、模型假设 模型一： 1)、感染人数是时间的连续可微函数； 2)、单位时间内感染人数的增长是常数，或单位时间内感染人数的增长量与当时的感染人数成正比。 模型二： 1)、感染人数是时间的连续可微函数； 感染人数受环境条件的限制，有一个最大的可感染人数 。 3)、单位时间内感染人数的增长率和感染人数有关，是其线性函数，最大感染时对应增长率为零。 四、模型建立与求解 模型建立 模型一：设t时刻的感染人数为，初始时刻(t=0 )的感染者人数为，感染者的增长率为r，根据单位时间内感染人数的增长率是常数的假设，t到 时间内感染人数的增量为: 因此，满足如下的微分方程： MATLAB dsolve(‘Dx=r0*(1-x/xm)*x’,’x(0)=x0’) ans=xm/(1+exp(-r0*t)*(xm-x0)/x0) X=[39,53,72,96,129,171,232,314,386,502,629,760,920,1065,1232] x1=dsolve(Dx=r1*x,x(0)=39) r1=log(53/39) r1=0.3067 模型二：仍然设t时刻的感染人数为 ，初始时刻 的感染者人数为0 ，感染者人数为0时，感染人数的增长率为 。根据单位时间内感染人数的增长率和感染人数有关，是其线性函数的假设，可得增长率关于感染者人数的线性函数关系式： 进一步，由最大感染时对应的增长率为零可确定参数k的值为： 因此，在该模型的假设下，感染人数 应满足如下的微分方程： MATLAB dsolve(‘Dx=r0*(1-x/xm)*x’,’x(0)=x0’ ans=xm/(1+exp(-r0*t)*(xm-x0)/x0) 模型求解 这是一个非线性微分方程，利用微分方程中的分离变数法，求得其解为： = X=[39,53,72,96,129,171,232,314,386,502,629,760,920,1065,1232] ; r1=0.3067;t=0:14;X