43110100122_邵建智_种群模拟模型及模型之混沌现象.doc

实 验 报 告 课程名称 生物系统模拟 姓 名 邵建智 学 号 3110100122 专 业 生物系统工程 实验名称 种群模拟模型及模型之混沌现象 实验一：常规模拟软件练习 实验类型： 上机操作 实验地点： 农生环D-414 指导老师： 方慧 实验日期： 2013 年 10 月 31 日 实验目的 学习种群模拟模型中的Lotak-Volterra竞争模型，连续型捕食者-猎物模型及模型中的混沌现象。 用Matlab软件编写程序，对Lotak-Volterra竞争模型，连续型捕食者-猎物模型进行模拟。并编写May的Logistic模型程序，了解模型混沌现象的产生。 实验仪器设备 计算机 Matlab软件 原理和方法 种群增长模型有两类：一类是与密度无关的种群增长模型；另一类是余密度有关的种群增长模型。一个以内禀增长率增长的种群，其种群数目将以指数方式增加。只有在种群不受资源限制的情况下，这种现象才会发生。尽管种群数量增长很快，但种群增长率不变，不受种群自身密度变化的影响。这类指数增长成为与密度无关的种群增长（density-independent growth）或种群的无限增长。与密度无关的种群增长又可分为两类。如果种群各个世代不相重叠，如许多一年生植物和昆虫，其种群增长是不连续的，称为离散增长，一般用差分方程描述；如果种群的各个世代彼此重叠，如人和多数兽类，其种群增长是连续的，可用微分方程描述。 受自身密度影响的种群增长称为与密度有关的种群增长（density-dependent growth）或种群的有限增长。种群的有限增长同样分为离散的和连续的两类。如长江的连续增长模型-逻辑斯蒂增长模型。 而Lotak-Volterra种群竞争模型中，考虑参与竞争的两个种群在没有竞争时均按逻辑斯蒂增长模型增长，另一竞争种群的出现会对原种群有负面作用。连续型捕食者-猎物模型中，捕食者与被捕食者在同一环境下生存，它们的种群变化速度互相影响，捕食导致被捕食者数量减少，捕食导致捕食者数量增加。进一步分析May’s Logistic模拟模型，则可在模型中观察到混沌系统中常见的分岔图。 实验步骤 Q1： 离散型Lotak-Volterra竞争模型模拟程序： 在命令行对以下4种情况进行模拟 （a）Lotka_Volterra(2,1,0.2,0.1,0.8,0.7,100,150) 两者竞争的结果是：Y物种存活，X物种灭亡 (b) Lotka_Volterra(1,2,0.1,0.2,0.7,0.8,150,100) 两者竞争的结果是：X物种存活，Y物种灭亡 (c) Lotka_Volterra(1,2,0.2,0.1,1.2,1.2,130,150) 两者竞争的结果是：X物种存活，Y物种灭亡 (d) Lotka_Volterra(1,2,0.2,0.1,0.6,0.8,100,100) 两者竞争的结果是：X,Y物种均存活 Q2: 连续型捕食者-猎物模型竞争模型模拟程序： 程序1：predator1: function predator1(X1,Y1) ts=0:0.1:30; x0=[X1,Y1]; [t,x]=ode45(predator,ts,x0); plot(t,x),grid,gtext(x(t)),gtext(y(t)); figure(2); plot(x(:,1),x(:,2)),grid 程序2：predator: function xdot=predator(t,x) r=1; d=0.5; a=0.1; b=0.02; xdot=[(r-a*x(2)).*x(1);(-d+b*x(1)).*x(2)]; 从左图可以看出，随着时间的变化，X与Y物种的数量呈现周期性的变化，从右图可以看出，X物种数量的多少影响Y物种数量的多少，同理，Y物种数量的多少影响X物种数量的多少，两者数量连线形成一个闭环。 Q3: May’s Logistic模型模拟程序： function Logistic() X(1)=0.1; for a=0:0.02:4 for t=1:300 X(t+1)=a*X(t)*(1-X(t)); end n=a*(ones(2