数学建模细菌繁殖问题.pdf

《数学建模细菌繁殖问题.pdf》由会员分享，可在线阅读，更多相关《数学建模细菌繁殖问题.pdf（13页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、细菌繁殖细菌繁殖摘要摘要本文针对酵母菌种群繁殖的基本特点，为达到解决所列出的三个问题的目的，建立了符合实际情况的预测模型。预测模型： 根据题目给出的已知条件， 最终建立了符合本题的 Logistic 模型。综合考虑了各种因素，利用计算机 MATLAB编程分别对问题进行求解，并分别绘制出本题的Logistic数学模型和问题三中所列的二次多项式的曲线， 以供对比。对于问题一得出，本文建立了种群预测的Malthus 模型以及符合本题的Logistic 模型，模型中参数 K 的值为：0.00081411，参数 M 的值为：663.97。对于问题二得出，自初始时刻起，20 小时时酵母菌的数量为：663.

2、06。该种群的增长呈现出 S 型，前期呈指数型增长，中后期增长缓慢，种群数量最终达到最大值：663.97。对于问题三得出，根据计算机 MATLAB程序绘制出的本题 Logistic 数学模型以及问题三中所列的二次多项式的曲线。 对两条曲线进行对比，易知符合本题的 Logistic 模型具有更好的预测能力。关键词：关键词：Malthus模型；Logistic 模型；MATLAB；预测1 1 问题重述问题重述已知酵母菌种群在培养物中的增长情况，见附录中表 a 所示。现根据已有的数据来预测酵母菌的数量，要求尽量与实际相符。根据以上题目所给的条件及数据，回答以下问题：问题一：建立酵母菌数量的数学模型，

3、确定模型中的未知参数；问题二：利用问题一中的模型，预测 20 小时时酵母菌的数量；问题三：若用二次多项式N(t) k0 k1t k2t2（其中ki(i 0,1,2)为常数）作为新模型，试从误差角度说明新模型与问题一中的模型哪个具有更好的预测能力，并画出对比曲线。2 2 问题的基本假设与说明问题的基本假设与说明1）假设题目所给的数据全部真实可靠，可以作为检验所建立的数学模型的准确性的事实依据。2）在自然环境中，细菌繁殖增长会受到各方面复杂因素的影响，为简化模型， 本文以题目中给出的实测数据，作为衡量所建立的数学模型准确度的主要因素。3）本文中该酵母菌种群的繁殖方式不随时间变化。3 3 符号说明符

4、号说明符号tN(t)MK、kN0表示意义时刻（单位为小时）t 时刻时酵母菌的数量酵母菌数量的最大值常系数酵母菌的初始数量：9.64 4 问题的分析问题的分析自然界中某酵母菌种群数量的变化和随着时间的发展过程， 是由很多因素决定的，自然环境、资源制约、种群的繁衍能力、种群的存活能力等，都能严重的影响种群的繁衍过程。然而，自然环境、资源制约却是决定该种群数量变化的直接原因。综合考虑这些因素成为构建符合本题中酵母菌种群繁殖预测模型的关键。建立模型对该酵母菌种群发展过程进行定量预测， 就是根据现有的统计资料和初始数据，从当前实际出发，并对未来的种群发展过程，提出合理的控制要求和假设说明，应用科学的方法

5、，预测该种群数量的发展趋势。为此，本文建立了具有预测性的 Malthus 模型，在综合考虑各影响因素后，建立了符合本题的Logistic 模型。Logistic 模型相比 Malthus 模型以及题中所述的二次多项式模型，更符合题目要求，用题中所给的实测数据检验后发现在误差允许范围内，是十分准确的；从误差的角度分析，Logistic 模型具有更好的预测能力。5 5 模型的建立与求解模型的建立与求解5.15.1 数据预处理数据预处理由于题中所给数据的不完备性， 并不能由它来预测未来种群的发展情况，但是基于抽样调查的等概率性，可以认为它反应的种群增长情况是符合实际情况的，因此认为，根据题中数据，结

6、合所建立的合理的数学模型，准确地对该酵母菌种群的繁殖增长数量作出合理预测。题中所给数据见附录中表 a。建模初始，本文将题中所给数据分为两部分考虑，其中前八组数据为第一部分，其余数据为第二部分。5.25.2 模型一：种群预测的模型一：种群预测的 MalthusMalthus 模型模型5.2.1 模型的建立在任意时刻 t，细菌的繁殖速度显然可以用表达式v(t ) limtt0N(t) N(t0)t t0来表示，设t0时刻细菌数量为N(to) N0。我们将时间间隔0, t 分成 n 等份。由于细菌的繁殖是连续变化的，在很短的一段时间内细菌数量的变化是很小的，繁t殖速度可近似看成是不变的。因此，在第一

7、段时间0,内，细菌数量满足关系n式N(t n) N0 kN0t n0,t n时段内细菌的增量为t t N N0 kN0nnt故时刻细菌数量为nt t N N01knnt 2t同理，第二时段,末细菌的数量为n nt 2t N N01knnn1nt依次类推，可以得到，最后一时段t,末细菌的数量为nn2t Nt N01k（1）n由 于 这 是 一 个 近 似 值 。 因 为 我 们 假 设 了 在 每 一 小 段 时 间i1it且等于该时段初始时刻的变nt,n(i 1,2,n)内细菌的繁殖速度是不变的，n化速度。 但这种近似程度将随着小区间的长度的缩小精度越高。若对时间间隔无限细分，就可以得到精确值

8、。所以，经过时间t后细菌总数为tnktktktN(t) limN0(1k) N0lim(1) N0ekt（2）nnnnn即种群预测的 Malthus 模型为：N(t) N0ekt。5.2.2 模型的求解本文结合题中所给实测数据的第一部分，运用计算机 MATLAB程序对（2）式进行求解，程序见附录中程序一，得到：k=0.4580即：N(t) 9.6e0.4580t，将其拟合的数据与题中酵母菌实测数据进行比对，见图一：400350300250已 知 酵 母 菌 的 观 测 数 据Malthus拟 合 的 数 据数量n20015010050001234时 间 t/ 小 时5678图一 酵母菌部分实测

9、数据与 Malthus 数学模型曲线为验证（ 2）式的准确性，结合题中所给的全部实测数据，运用计算机MATLAB程序对（2）式进行求解，程序见附录中程序二，得到：k=0.25781即：N(t) 9.6e0.25781t，将其拟合的数据与题中酵母菌实测数据进行比对，见图二：1000900800700600已 知 酵 母 菌 的 观 测 数 据Malthus拟 合 的 数 据数量n50040030020010000246810时 间 t/ 小 时12141618图二 酵母菌实测数据与 Malthus 数学模型曲线根据图二可以看出，在种群繁殖增长前期，Malthus 数学模型可以较为准确的表示出其增

10、长规律；在种群增长中后期却有很大偏差。本文结合 Malthus 数学模型拟合的数据与题中酵母菌实测数据对比，得各个时刻种群数量的误差见表 1表 1 各个时刻种群数量的误差表时间 t/小时误差时间 t/小时误差时间 t/小时误差00714181529163101741118512613-5.8767-12.923-26.395-44.177-84.759-129.51-198.95-275.19-343.29-386.85-396.06-383.04-355.36-286.16-192.17-62.002108.96332.79其误差偏大，不利于中长期预测。算得其误差平方和为： 1.084810

11、6，算得第 20 小时时，酵母菌的数量为：1665.6，是不符合实际情况的。又由图二可知，Malthus 数学模型预测该种群呈现出无限增长的趋势，显然不符合实际。因此可以说明，Malthus 数学模型不能准确的表示出该种群的发展趋势，不具备预测该种群增长数量的能力。为此，本文对该模型进行了改进与优化，充分考虑影响该种群繁殖增长的各种符合实际情况的因素后， 建立了能准确表示出该种群的发展趋势，具备更好的预测能力的 Logistic 模型。5.35.3 模型二：符合本题的模型二：符合本题的 LogisticLogistic 模型模型5.3.1 模型的建立结合 Malthus 数学模型的推导，本文建

12、立了符合本题的Logistic 模型，其数学表达式为dN(t) KN(t)(M N(t)（3）dt该式中 K 为常系数，M 为酵母菌数量的最大值，N(t)为任意时刻 t 时酵母菌的数量。整理该式得到 N(t)的表达式为N(t)=当t 时，N(t)的表达式为N(t)=N0MekM (tt0)M N0 N0ekM (tt0)（4）N0MN0(M N0)ekM (tt0)即符合本题的 Logistic 模型为：N(t)=5.3.2 模型的求解N0MN0(M N0)ekM (tt0)结合题中所给实测数据，运用计算机 MATLAB 程序对（4）式进行求解，程序见附录中程序三，得到K=0.00081411

13、，M=663.97即：N(t)=6374.1129.6 654.37ekM (tt0)，将其拟合的数据与题中酵母菌实测数据进行比对，见图三：700已 知 酵 母 菌 的 观 测 数 据Logistic模 型 拟 合 的 数 据600500400数量n30020010000246810时 间 t/ 小 时12141618图三 酵母菌实测增长数据与 Logistic 数学模型曲线由图三可以看出，在种群繁殖增长的整个过程中，Logistic 数学模型可以准确的表示出其增长规律，题中所给的实测数据与Logistic 数学模型拟合数据，在误差允许的范围内，几乎一致。本文结合Logistic 数学模型拟合

14、的数据与题中酵母菌实测数据对比，得各个时刻种群数量的误差见表 2表 2 各个时刻种群数量的误差表时间 t/小时误差时间 t/小时误差时间 t/小时误差0073.107914181529163-1.305710-4.971174113.7856185126.698666.772613-3.3014-1.9864-1.4754-1.7208-5.81023.9774-0.362410.57621-0.477080.23238-0.21553-0.50652在误差允许范围内，预测的数据是合理的。算得其误差平方和为：211.75，算得第 20 小时时，酵母菌的数量为：663.06，又由图像知，Logi

15、stic 数学模型预测该种群呈现出前期、中期增长较快，呈现出J 型曲线，在后期随着种群数量不断接近其最大值 M=663.97 种群增长缓慢， 种群结构趋于稳定， 其数量变化很小，显然是符合实际的。因此可以说明，Logistic 数学模型能准确的表示出该种群的发展趋势，具备预测该种群增长数量的能力。5.45.4 二次多项式二次多项式N(t) k0 k1t k2t2模型模型根据题目所述， 结合题中所给的实测数据， 本文利用计算机 MATLAB程序，运行程序见附录中程序四，计算得出该模型的参数值为k0=93.785，k1=65.706，k2=1.1232则该二次多项式N(t) k0 k1t k2t2

16、可化为N(t) 93.78565.706t 1.1232t2整理得N(t) 1.1232（5）t265.706t 93.785为达到在误差角度下，比较 Logistic 数学模型与该模型的预测能力的目的。本文利用计算机 MATLAB 程序，将题中的实测数据、Logistic 数学模型以及二次多项式模型，绘制成如下图四的曲线：8007006005004003002001000-1000已 知 酵 母 菌 的 观 测 数 据Logistic模 型 拟 合 的 数 据二 次 拟 合 的 数 据数量n246810时 间 t/ 小 时12141618图四 酵母菌实测数据、Logistic 模型以及二次多

17、项式的对比曲线由图四可以看出，题中所给的实测数据均匀的分布在 Logistic 模型的曲线上，相比之下，该二次多项式曲线误差较大，与实际情况不符。本文结合该二次多项式模型拟合的数据与题中酵母菌实测数据对比， 得各个时刻种群数量的误差见表 3表 3 各个时刻种群数量的误差表时间 t/小时误差时间 t/小时误差时间 t/小时误差0753.81514189.27141524.133491614.052346.023101738.99479.965111854.43587.06212685.41113-103.38-47.502-34.418-62.355-66.637-61.867-58.842-3

18、4.864-12.033其误差偏大，预测的数据准确度偏低。算得其误差平方和为： 59212。显然，Logistic 数学模型的误差更小，具有更好的预测能力。6 6 模型的优缺点分析模型的优缺点分析Malthus 数学模型在短期预测中具有准确度高， 操作简便， 容易实现等优势，但是不能准确的预测出该种群长期的发展趋势， 不具备预测该种群长期增长数量的能力。为此，本文对该模型进行了改进与优化，充分考虑影响该种群繁殖增长的各种符合实际情况的因素后，建立了能准确表示出该种群的发展趋势，具备更好的预测能力的 Logistic 模型。在种群繁殖增长的整个过程中，Logistic 数学模型可以准确的表示出其

19、增长规律， 题中所给的实测数据与 Logistic 数学模型拟合数据， 在误差允许的范围内，几乎一致。相比较二次多项式，Logistic 数学模型误差小，预测更准确。但 Logistic 数学模型操作繁杂，更适用于中长期预测；在短期种群增长预测中通常选用Malthus 数学模型。7 7 模型的改进及其推广模型的改进及其推广Logistic 数学模型在中长期种群繁殖增长预测中，具有误差小、准确度高的优势，可将其用于卫生防疫部门监测细菌的繁殖状况，提出预警机制。在疾病防治等公共医疗卫生中，Logistic 数学模型可以很好的观察流行疾病的发展态势，为更好的预防疾病防治，提供可靠保证和理论事实依据。

20、8 8 参考文献参考文献1 陈恩水，王峰，朱道元. .数学建模与实验. .北京：科学出版社，20082 熊启才，张东升. .数学模型方法及应用. .重庆：重庆大学出版社，20053 秦新强，赵凤群. .线性代数学习指导. .北京：机械工业出版社，20064 刘卫国. .MATLAB程序设计与应用. .北京：高等教育出版社，20065 邬学军，周凯. .数学建模竞赛铺导教程. .杭州：浙江大学出版社，2009附录：附录：表 a 酵母菌在培养物中的增长时间 t（小时）数量 N（t）时间 t（小时）数量 N（t）时间 t（小时）数量 N（t）时间 t（小时）数量 N（t）09.65119.61051

21、3.315651.1118.36174.611559.716655.92297257.312594.817659.6347.28350.713629.418661.8471.1944114640.8程序一：程序一：种群预测的 Malthus 模型（部分数据）function f=curvefun3(c,t)f=9.6*exp(c*t);clct=012345678;n=9.618.32947.271.1119.6174.6257.3350.7;c0=1;c, options,funval, Jacob =lsqcurvefit(curvefun3,c0,t,n)m= curvefun3(c,t

22、)plot(t,n,o,t,m,*)xlabel(时间 t/小时);ylabel(数量 n);legend(已知酵母菌的观测数据,Malthus 拟合的数据)程序二：程序二：种群预测的 Malthus 模型（全部数据）function f=curvefun3(c,t)f=9.6*exp(c*t);k= 0.25781function f=curvefun3(c,t)f=9.6*exp(c*t);clct=0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 1516 17 18;n=9.6 18.3 29 47.2 71.1 119.6 174.6 257.3 350.7441

23、 513.3 559.7 594.8 629.4 640.8 651.1 655.9 659.6661.8;c0=1;x, options,funval, Jacob =lsqcurvefit(curvefun3,c0,t,noptions = 1.0848e+006Jacob = 3clct=012345678910 11 12 13 14 15 16 17 18;n=9.6 18.329 47.271.1119.6174.6257.3350.7441 513.3559.7594.8629.4640.8651.1655.9659.6661.8;m=9.6*exp(0.25781.*t);pl

24、ot(t,n,o,t,m,*)xlabel(时间 t/小时);ylabel(数量 n);legend(已知酵母菌的观测数据,Malthus 拟合的数据);程序三：符合本题的 Logistic 模型function f=curvefun4(c,t)f=-9.6.*c(1)./(-9.6-exp(-c(1).*c(2).*t)*c(1)+exp(-c(1).*c(2).*t).*9.;clct=0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 1516 17 18;n=9.6 18.3 29 47.2 71.1 119.6 174.6 257.3 350.7441 513.3

25、559.7 594.8 629.4 640.8 651.1 655.9 659.6661.8;c0=700 0.001; x, options,funval, Jacob =lsqcurvefit (curvefun4,c0,t,n) f= curvefun4(c,t)plot(t,n,o,t,f,*)xlabel(时间 t/小时);ylabel(数量 n);legend(已知酵母菌的观测数据,Logistic 模型拟合的数据);运行结果x = 663.97 0.00081411options = 211.75Jacob = 3程序四：程序四：酵母菌实测数据及 Logistic 模型和题目中二

26、次多项式对比曲线t=0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 1516 17 18;n=9.6 18.3 29 47.2 71.1 119.6 174.6 257.3 350.7441 513.3 559.7 594.8 629.4 640.8 651.1 655.9 659.6661.8;a=polyfit(t,n,2)z=polyval(a,t);plot(t,n,o,t,z,*)xlabel(时间 t/小时);ylabel(数量 n);legend(已知酵母菌的观测数据,二次拟合的数据);for i=1:18 y(i)=z(1,i)-n(1,i); k(i)=y(i)*y(i);endyk=sum(k)a = -1.1232 65.706 -93.785k = 59212