第五章医药数学模型.ppt

考虑到条件3，模型二的(2)式应改为di(t)/dt=ks(t)i(t)-?·i(t)i(0)=i0(5)将s(t)=n-i(t)代入上式,得di(t)/dt=ki(t)[n-i(t)]-?·i(t)i(0)=i0(6)其解为i(t)=(kt+n/i0)-1(k=?)和i(t)=[k/(k-?)+(n/i0-k/(k-?))exp(-(k-?)t)]-1(k≠?)第53页,共60页，5月，星期六，2024年，5月i0i0接触数?=1~阈值感染期内有效接触感染的健康者人数不超过病人数i00ti?11-1/?i0t??1di/dt0第54页,共60页，5月，星期六，2024年，5月传染病模型四大多数传染病如天花、流感、肝炎、麻疹等治愈后均有很强的免疫力，所以病愈的人既非健康者，也非病人(已感染者)，他们已经退出传染系统，称为移出者我们将人员分成3类:第一类是传染者(i)类;第二类是健康者(易受传染者)(s)类;第三类是移出者(r)类，(包括患病死去的人、痊愈后具有长期免疫力的人以及在痊愈以前被隔离起来的人)令i(t)，s(t)和r(t)分别表示在时刻t的第一、二、三类的人数。第55页,共60页，5月，星期六，2024年，5月模型四的假设1)总人数为常数n，则i(t)十s(t)十r(t)=n(5)2)单位时间内一个病人能传染的人数与当时健康者人数成正比，比例系数为k(传染强度)。3)单位时间内病愈免疫的人数与当时的病人人数成正比，比例系数为a(恢复系数)。第56页,共60页，5月，星期六，2024年，5月模型四的建立由假设3)，有dr/dt=a·i(t)(6)由于引入了r(t)，则模型2的方程(2)式应改为di(t)/dt=k·s(t)i(t)-dr/dt(7)上式表示单位时间内病人人数的增加应等于被传染的人数减去病愈的人数。从(5)~(7)式中消去di/dt，并设r(0)=0，s(0)=s0。得ds/dt=-ks(t)i(t)dr/dt=a·i(t)r(0)=0,s(0)=s0i(t)十s(t)十r(t)=n(8)第57页,共60页，5月，星期六，2024年，5月模型四的求解记p=a/k，称为特征指标，对同一地区，同一种传染病,p是常数。由方程(8)可得ds/dt=-(1/p)s(t)(dr/dt)其解为s(t)=s0e-(1/p)r(t)≈s0[1-(1/p)r(t)+(1/2p2)r2(t)]于是dr/dt=a·i(t)=a[n-r(t)-s(t)]≈a[(n-s0)+(s0/p-1)r(t)-(s0/2p2)r2(t)]再解上式,经过整理简化得limt→∞r(t)≈(2p2/s0)(s0/p－1)上式表明:当s0很小(p)时，根本传染不开；只有s0很大(p)时，才存在传染问题。第58页,共60页，5月，星期六，2024年，5月今设s0=p十h(h0,只有这时才有传染问题)，且hp,则limt→∞r(t)≈[2p2/(p+h)](h/p)≈2h因limr(t)表示最终免疫的人数，也即最终痊愈的人数。显然，这也就是传染病传染的人数，于是由上式可得下列结论第59页,共60页，5月，星期六，2024年，5月感谢大家观看第60页,共60页，5月，星期六，2024年，5月那么(1)的解为如下形式:C1(t)=Ae-at＋Be-bt+E·exp(-k01t)C2(t)=Fe-at＋Ge-bt+H·exp(-k01t)借助待定系数法和初始条件C1(0)=C2(0)=0，可得系数A、B、E、F、G、H的