孤岛疾病问题的探讨.docVIP

孤岛疾病问题的探讨 摘要 建立传染病的数学模型来描述传染病的传播过程，分析受感染人数的变化规律，探索制止传染病蔓延的手段等，一直是各国有关专家和官员关注的课题。 本文是一个孤岛疾病传播问题。首先根据孤岛的自身特征方面的条件解决问题a:列出假设1、2，问题b:根据模型建立的被感染人数X与时刻t的关系，利用matlab软件画出dX/dt关于X的图像。问题c:又假设感染的人数X1N/2,画出X关于t的图形。在此问题的基础上取X1=N/4和X1=3N/4；问题d:求解出微分方程，最后把X作为t的函数，解出前面给的模型。问题e:此模型方程是logistic模型，求出X的表达式,得出当t趋于无穷大时X的极限为N。问题f:有代表的数据列出表格求解出这些数据，得出的这些数据支持该模型，估计模型中的常数并预测t=12天被感染的人数，列出关系模型。问题h:分析上述模型的优缺点，试给出改进方案。 【关键词】孤岛疾病 微分方程 关系模型 改进方案 一、问题重述 如何确定孤岛的疾病传播速度问题通过正确的措施进行控制并对人们的身心健康有很大的帮助，也是必不可少的。 考虑在一个人口数量为N的孤岛上，一部分到岛外旅游的居民回来使该岛感染了一种高传染性的疾病。请预测在某时刻t将会被感染的人数X。考虑一下模型，其中k0为常数： ??? 列出这个模型所隐含的两条主要假设，说明这些假设有什么依据？ (b)?? 画出关于X的图形 (c)??? 若初始被感染的人数，画出X关于t的图形；若初始被感染人数为，画出X关于t的图形。 (d)?? 把X作为t的函数，解出前面给出的模型。 (e)???由(d)，当t趋于无穷时求X的极限。 (f)????设岛上的人口有5000人，在传染期的不同时刻被感染人数如下表 天数t 2 6 10 被感染人数X 1887 4087 4853 ln(X/(N-X)) -0.5 1.5 3.5 问这些数据能否支持所给的模型？ (g) 利用(f)的结果估计模型中的常数，并预测t=12天时被感染的人数。 (h) 分析上述模型的优缺点，试给出改进方案。 二、模型假设 1、人口数量N不变，因为是孤岛。 2、人口分为健康人和被感染的病人，数量分别为X ,N-X。 3、在规定的时刻内人口变化X取整数，因为人口为整数 三、符号定义说明 X（t）：t时刻被感染的人数 S（t）：t时刻未被感染的人数 k ：常数 C1、C2：皆为参数 N ：孤岛上的人口总数,即N=X（t）+S（t） 四、模型的建立与求解 4.1模型的建立依据 列出这个模型所隐含的两条主要假设，说明这些假设有什么依据？ 假设1人口数量N不变，因为是孤岛。 假设2人口分为健康人和被感染的病人数量分别为X ,N-X。 4.2模型的建立与求解 (b) 画出关于X的图形 设k=0.1 X=0:0.1:1; ezplot(0.1*X*(1-X),[0,1]); (c)若初始被感染的人数，画出X关于t的图形；若初始被感染人数为，画出X关于t的图形。 1、取，则 y=dsolve(Dy=0.1*y*(1-y),y(0)=0.25,x); ezplot(y,[0,50]); 2、，则 y=dsolve(Dy=0.1*y*(1-y),’y(0)=0.75,x); ezplot(y,[0,50]); (d)把X作为t的函数，解出前面给出的模型。 syms N k:X=dsolve(DX=k*X*（N-X）) 即 (e) ?由(d)，当t趋于无穷时求X的极限。 (f)????设岛上的人口有5000人，在传染期的不同时刻被感染人数如下表 天数t 2 6 10 被感染人数X 1887 4087 4853 ln(X/(N-X)) -0.5 1.5 3.5 问这些数据能否支持所给的模型？ 由 可以看出为线性变化的，所以可以认为这些数据支持该模型。 (g)利用(f)的结果估计模型中的常数，并预测t=12天时被感染的人数。 得到k=0.5 再代入k值，得C2=1.5，而N=5000， 当t=4+10=14时，可以解出X的值： X= 4979 确定参数C1=N/X0，则有 五、模型评估及改进 5.1 模型的优点 （1）用假设分析、微分求导导数、线性相关等方法，解决问题，方法简单易懂，过程清晰且准确。 （2）又利