数学建模(方红)教学课件18.微分方程.pptVIP

华中农业大学数学建模基地运行结果黑色星号--Logistic模型预测值，绿色圆圈--Malthus模型预测值，蓝色曲线为实际统计值。华中农业大学数学建模基地传染病模型随着卫生设施的改善，医疗水平的提高及人类文明的不断发展，诸如霍乱、天花等曾经肆虐全球的传染性疾病已经得到了有效的控制。但是一些新的、不断变异着的传染病毒却悄悄地向人类袭来，20世纪80年代十分险恶的艾滋病毒开始肆虐全球，至今仍在蔓延；2003年春来历不明的SARS病毒突袭人间，给人们的生命财产带来了极大的危害。长期以来，建立传染病的数学模型来描述传染病的传播过程、分析受感染人数的变化规律、探索制止传染病蔓延的手段等，一直是有关专家关注的一个热点问题。华中农业大学数学建模基地已感染人数(病人)i(t)每个病人每天有效接触(足以使人致病)人数为?模型1假设若有效接触的是病人，则不能使病人数增加必须区分已感染者(病人)和未感染者(健康人)建模?传染病模型华中农业大学数学建模基地模型2区分已感染者(病人)和未感染者(健康人)假设1）总人数N不变，病人和健康人的比例分别为2）每个病人每天有效接触人数为?,且使接触的健康人致病建模?~日接触率SI模型传染病模型华中农业大学数学建模基地模型2Logistic模型当时，用Matlab程序为:y=dsolve(Dy=0.1*y*(1-y),y(0)=0.09,x)%求解此微分方程ezplot(y,[0,60])%画出微分方程的图像ezplot(0.1*y*(1-y),[0,1])%画出y的导数的图像传染病模型华中农业大学数学建模基地传染病模型tm~传染病高潮到来时刻?(日接触率)??tm?t=tm,di/dt最大华中农业大学数学建模基地传染病模型II的函数图像华中农业大学数学建模基地模型3传染病无免疫性——病人治愈成为健康人，健康人可再次被感染增加假设SIS模型3）病人每天治愈的比例为??~日治愈率建模?~日接触率1/?~感染期?~一个感染期内每个病人的有效接触人数，称为接触数。传染病模型华中农业大学数学建模基地编写MATLAB程序如下：y=dsolve(Dy=0.01*y*(1-y)-0.05*y,y(0)=0.7,x);%ezplot(y,[0,120])y2=dsolve(Dy=0.3*y*(1-y)-0.15*y,y(0)=0.7,x);y3=dsolve(Dy=0.3*y*(1-y)-0.15*y,y(0)=0.3,x);figure,ezplot(y2,[0,25]);figure,ezplot(y3,[0,25])传染病模型#2022华中农业大学数学建模基地传染病模型不难看出，接触数?=1~阈值华中农业大学数学建模基地模型4传染病有免疫性——病人治愈后即移出感染系统，称移出者SIR模型假设1）总人数N不变，病人、健康人和移出者的比例分别为2）病人的日接触率?,日治愈率?,接触数?=?/?建模需建立的两个方程传染病模型华中农业大学数学建模基地模型4SIR模型无法求出的解析解求数值解华中农业大学数学建模基地MATLAB程序如下：ts=0:50;x0=[0.02,0.98];[t,x]=ode45(ill,ts,x0)%调用ode45求解ill方程组plot(t,x(:,1),t,x(:,2)),grid,%画出健康者和病人的变化曲线figure,plot(x(:,2),x(:,1)),grid%画出相图functiony=ill(t,x)%函数ill，表示模型IVa=1;b=0.3;y=[a*x(1)*x(2)-b*x(1),-a*x(1)*x(2)];传染病模型华中农业大学数学建模基地画出健康者和病人的变化曲线华中农业大学数学建模基地传染病模型结论：在初始时刻健康者和病人百分比的总和为1；病人的数量先增加然后下降，说明在某时刻传染病得到抑制；而治愈的人群退出此系统，所以最后系统的人群数量为0；这时所有的人群均是免疫者。华中农业大学数学建模基地地中海鲨鱼问题意大利生物