4.2-常微分方程组模型

《4.2-常微分方程组模型》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数学模型安徽大学数学科学学院4.2常微分方程组模型由几个微分方程联立而成的方程组称为微分方程组。本节建立传染病模型、男生追女生模型和种群增长模型，并判断方程组的稳定性。4.2.1传染病模型我国目前的法定传染病有甲、乙、丙三类，共39种。传染病的特点是有病原体，有传染性和流行性，感染后常有免疫性。有些传染病还有季节性或地方性。传染病的分类尚未统一，有人按病原体分类，有人按传播途径分类。传染病的预防应采取以切断主要传播环节为主导的综合措施。传染病的传播和流行必须具备3个环节，即传染源（能排出病原体的人或动物）、传播途径（病原体传染他人的途径）及易感者（对该

2、种传染病无免疫力者）。若能完全切断其中的一个环节，即可防止该种传染病的发生和流行。各种传染病的薄弱环节各不相同。在预防中应充分利用。除主导环节外对其他环节也应采取措施，只有这样才能更好地预防各种传染病。不同类型的传染病，其传播过程有着各自不同的特点，了解这些具体的传染病的传播过程需要了解其病理知识，这里不可能从医学角度一一进行分析，而主要按照一般的传播机理建立几类一般的传染病模型分析受感染人数的变化规律，讨论终止传染病蔓延的方法和手段。模型1.用表示时刻的病人数量，假设病人一旦与健康人群接触就会使健康人群患病，且单位时间内每个病人能够使健康人患病的人数

3、为。初始时刻的病人数为。利用微元法，考虑内病人数的变化，则有，上式两边同时除以，并令，得到（4.2.1）式（4.2.1）的求解程序及结果如下：dsolve('Dx-lambda*x=0','x(0)=x0','t')ans=x0*exp(lambda*t)即式（4.2.1）的解为（4.2.2）式（4.2.2）为指数形式，故称模型（4.2.1）为指数增长模型。根据式（4.2.2），当时，，即所有人都会患病。很显然，这与事实不符。模型2.（SI模型）考虑以下假设：（1）不考虑人口的出生、死亡和迁移等种群动力因素，在疾病传染期所考察地区内总人数始终保持常数不

4、变；（2）人群分健康人群（易感染者，Susceptible）和病人（已感染者，Infective），时刻这两类人在总人数中所占比例分别为和，即有，并设初始时刻患病人数比例为；（3）病人一旦与健康人群接触，就必然具有一定的传染性。假设时刻单位时间内每个病人能传染的易感染者数量与该地区易感染者比例成正比，比例系数为，即时刻单位时间内每个病人可使个易感染者患病，而病人总数为，因而时刻单位时间内共有个易感染者患病。利用微元法，内患病人数的变化量为：上式两段同时除以，并令，有。结合假设（2），有（4.2.3）先观察的图像，建M文件myfun1.m。命令如下：fu

5、nctiony=myfun1(i);y=0.01*i*(1-i);再输入命令：fplot('myfun1',[0,1]);fplot函数的命令格式为fplot(‘fun’,lims)，表示绘制字符串fun制定的函数在lims=[xmin,xmax]或lims=[xmin,xmax,ymin,ymax]上的图形。但是fun必须是M文件的函数名或是独立变量为x的字符串。图形如下：图4.2.1图像由图4.2.1可以看出，当时，达到最大值，此时病人数增长得最快，意味着传染病高潮的到来。模型（4.2.3）的求解程序如下：symsialphat;dsolve('D

6、i-alpha*i*(1-i)','i(0)=i0','t')结果为：ans=1/(1-exp(-alpha*t)*(-1+i0)/i0)这里syms是Matlab中定义多个变量的。根据运行结果，有式（4.2.3）的解：（4.2.4）由式（4.2.4），当时，，这是传染病最高峰时刻，此时与成反比。因为反映了单位时间每个病人的传染数量，所以也被称为单位时间接触率，它直接反映了当地的卫生水平可以推迟传染病高潮的到来。函数的作图命令及图形如下（这里取，）：fplot('1/(1-exp(-0.5*t)*(-1+0.01)/0.01)',[0,30]);xla

7、bel('t');ylabel('i');图4.2.2患病人数比例的图像由图4.2.1可以看出，单调递增，且在时下凸，而在时上凸，且有渐近线，即病人数量一直递增，在病人比例不到时增长速度非常快，而在病人比例超过时增长速度放缓，最终病人的比例将接近于100%。显然这与事实不符，究其原因，是因为模型（4.2.3）中没有考虑病人被治愈的情况。模型3.（SIS模型）增加假设：假设被治愈的人不具有免疫能力，将成为易感染者，且时刻单位时间内治愈者变化率与病人数量成正比，比例系数为，利用微元法，内患病人数的变化量为，故有(4.2.5)显然，当，即没有病人治愈时，模型

8、(4.2.5)变为模型(4.2.3)。令，则有（4.2.6）因为反映了单位时间每个病人的传染数