关于澳洲瓢虫和介壳虫在使用DDT前后数量变化的分析.docxVIP

关于澳洲瓢虫与介壳虫在使用DDT前后数量变化的分析摘要 鉴于澳洲柠檬园中因介壳虫虫害，引进瓢虫，又由于后来农药DDT的普遍应用，当柠檬园主用DDT除害虫时，由于农药对两者均有作用，而使瓢虫大量死亡，介壳虫数量反而数量增加，本论文将从简单模型入手，进而通过建立Lotka-Volterra模型，最后优化考虑DDT浓度随时间变化的影响，完善模型，最终解释此现象。 关键词：Lotka-Volterra模型 自然阻滞 农药浓度 问题的提出 1968年，澳洲柠檬园中因介壳虫虫害，引进瓢虫，又由于后来农药DDT的普遍应用，柠檬园主用DDT进一步除害虫时，却发现由于农药对两者均有作用，而使瓢虫大量死亡，介壳虫数量反而数量增加，使虫害没有降低反而增加。为了解释此现象，必须建立相应的数学模型，找出害虫增长受DDT、瓢虫影响和瓢虫增长受害虫、DDT影响及其相互之间的关系。问题的分析此题在寻求DDT对介壳虫和瓢虫的影响关系，我们知道介壳虫和瓢虫之间构食饲—捕食者模型，我们可以先用Volterra模型描述介壳虫和瓢虫间的关系，再分析加入农药后，两虫数量的变化情况。在能解释题目所给的关系后，我们建立模型2增加logistics项，从而对Volterra的局限性进行了修改。在模型3中，我们加入考虑了时间变化对农药浓度的影响。建立数学模型过程（一）基本模型1、模型假设1）介壳虫与瓢虫构成一个封闭系统，不考虑其他生物的影响；2）DDT能在很快的时间发挥效用（不记发效时间）；3）农药对两虫的效用是相同的；2、定义符号说明 1）函数变量 ：时刻介壳虫的数量：时刻瓢虫的数量：介壳虫的增长率：瓢虫的增长率：单位时间内瓢虫对介壳虫的捕食率：单位时间内介壳虫对瓢虫的增长贡献率 农药对介壳虫和瓢虫的分别效用3、模型建立模型1：对于介壳虫来说，其增长来源于自身的生长繁殖，假设这时对于介壳虫生长环境是非常优良的，不用考虑环境最大承载量，造成介壳虫死亡的原因只有瓢虫的捕食，其减少的量和瓢虫的相互接触程度成正比。对于瓢虫来说，一旦没有了介壳虫，他就没有了食物，我们假设瓢虫的死亡和他的数量成正比，介壳虫的存在使瓢虫的数量增长，增长的速率与瓢虫与介壳虫的接触成正比。所以可得到介壳虫和瓢虫的微分方程。模型求解：程序如下：function dydt=ican(t,y)a=1;b=0.1;c=0.5;d=0.02;dydt=[(a-b*y(2))*y(1) (-c+d*y(1)).*y(2)];y(1):x y(2):y(t)4、模型求解根据上面提出的微分方程，能得到：这是个可分离变量的一阶微分方程，可改写为： 对方程求积分得： ， 为常数。运用自然对数和指数的性质，可将方程改写为：其中为常数。方程可确定相平面上的解轨：图形显示其相轨线是闭合的，它表示周期运动。模型的解释：这两个物种的曲线可知，在一个完整的循环内，介壳虫的数量在其最大和最小值之间波动。注意到，当介壳虫充沛时，瓢虫有最大的增长率，但是介壳虫的数量减少之后，瓢虫的数量才达到最大值。瓢虫的变化落后于介壳虫。加入农药后 给定初始量级介壳虫和瓢虫会随时间沿一条闭轨线运动，设表示经过一个完整时间周期返回初始点的时间，介壳虫和瓢虫一个时间周期的平均量分别由以下积分给出： 对方程两边同时积分，结合得：，由相轨的周期性，求出平均值，我们假设，农药使介壳虫和瓢虫的数量以减少，常数表示农药的效用。改进的模型可变为：改进方程组与原来方程组有同样的形式，所以：当使用农药时事实上增加了介壳虫的数量，减少了瓢虫的数量。如图：模型2： 自然阻滞增长模型自然环境中，介壳虫和瓢虫都不会无限制的增长，所以我们在下面的模型中考虑介壳虫和瓢虫的环境最大承载量。设分别为介壳虫和瓢虫的环境最大承载量。 建立方程有： 模型的求解：程序如下：function dydt=ican(t,y)a=1;b=0.1;c=0.5;d=0.02;m=200;n=100;dydt=[a*y(1)*(1-y(1)/m-b*y(2)) c*y(2)*(-1-y(2)/n+d*y(1))];两虫曲线如下：相轨线如下：加入农药模型程序如下：function dydt=ican(t,y)a=1;b=0.1;c=0.5;d=0.02;m=200;n=100;r=0.5;dydt=[a*y(1)*(1-y(1)/m-b*y(2))-r*y(1) c*y(2)*(-1-y(2)/n+d*y(1))-r*y(2)];曲线关系如下：模型评价：自然界长期存在的呈周期变化的生态平衡系统应该是稳定的，即系统受到不可避免的干扰而偏离原来的周期轨道后，其内部制约作用会使系统自动恢复原状，而Volterra模型描述的周期变化却不是稳定的，一旦离开某一条闭轨线，就进入另一条闭轨线，加入考虑自身阻滞作用的logi