生物种群动态模拟.doc

《生物种群动态模拟.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、西北农林科技大学实验报告学院名称：理学院专业年级：2011级信计1班姓名：学号：课程：数学模型与数学建模报告日期：2013年12月14日1实验题目:生物种群动态模拟2实验问题陈述:在这个问题中我们采用较复杂的增长模型g（p）=(r/c)(p-c)(1-p/k),其中参数c表示种群水平的最小值，低于这个水平出现的负增长，假设参数a=10∧-7，最小值的种群水平是蓝鲸3000条，长须鲸15000条。(a)两种鲸鱼种群能否共存？建一个动力系统模型描述种群的动态。(b)画这个模型的向量场。确定每个平衡态的位置。(c)确定状态

2、空间每个平衡态是否稳定。d)假设存在5000条蓝鲸和70000条长须鲸。关于这两种鲸鱼的将来这个模型能预测些什么？3实验目的：(1)建立恰当的动力系统模型来描述待求解问题；(2)画出模型向量场并确定平衡态位置；(3)确定状态空间中的平衡态是否稳定；(4)预测这两种鲸鱼的将来情况。4实验内容对于单个种群的增长，有，.当竞争种群存在时，有模型（4.2.1）这就是两个种群的竞争模型。在鲸鱼种群竞争的问题中，已知：，；，；由假设5，，但未知。由此可得模型为（4.2.2）这是一个有两个状态变量的微分方程组。主要依赖于定性分析的

3、方法了解这两个种群的动态特征。这是一个二维的自治系统，相空间应该是二维的平面。针对实际问题的要求，应该在这个平面的第一象限中讨论系统的动态。为叙述方便，将模型（4.2.1）记为：，，其中，，。5实验结果分析与讨论5.1分析Ⅰ——种群的动态特征（一）平衡点：，（二）等倾线：（三）定性分析：不稳定。关于有：，在附近的轨线为同心圆。（四）周期性：系统的解为在相平面的第Ⅰ象限内的封闭轨线。①轨线方程有解②系统的解是相平面的第Ⅰ象限内的封闭曲线。③周期：给定，存在，使得，且对于任何都不成立上述关系，则称为种群波动的周期。5.2

4、分析Ⅱ——人类控制活动的影响①周期平均：为周期为T的函数，则称为的周期平均。由方程可知，于是有②人类控制行为的影响按比率对进行捕杀，则有模型（5.2.1）模型的周期平均（即平衡点）为（5.2.2）6实验程序（Matlab或者其它软件语言陈述）6.1绘制向量场clear;clcsymsBFalpha=0.00001;%1*10^(-7);DB=0.05*B-alpha*B*F;DF=0.08*F-alpha*B*F;[Bsteady,Fsteady]=solve(DB,DF);disp('平衡点为：')disp([Bs

5、teadyFsteady])M=10;Bmin=0;Bmax=5000;Fmin=0;Fmax=70000;[X1,X2]=meshgrid(Bmin:(Bmax-Bmin)/M:Bmax,Fmin:(Fmax-Fmin)/M,Fmax);dX1=0.05*X1-alpha*X1.*X2;dX2=0.08*X2-alpha*X1.*X2;quiver(X1,X2,dX1,dX2,'r');axis([BminBmaxFminFmax]);holdonezplot(DB,[BminBmaxFminFmax]),hold

6、onezplot(DF,[BminBmaxFminFmax])xlabel('蓝鲸');ylabel('长须鲸');运行程序可得：平衡点为：[0,0][8000,5000]使用MATLAB绘制向量场functionzqjh()x=input('x=');y=input('y=');[x,y]=meshgrid(0:1:15,30:1:45);u=x.*(0.05-0.01*y);%方程B'=Bf1(B,F)v=y.*(0.08-0.01*x);%方程F'=Ff2(B,F)quiver(x,y,u,v);grid;ho

7、ld;x=2:0.01:6;plot(x,50-(10/3).*x);holdon%水平集LBx=0:0.01:15;plot(x,40-0.5.*x);%水平集LFend在MATLAB界面中，输入x=3.5y=38.2Currentplotheld其中，为了方便，对模型表达式做了一些处理：。运行程序可得：注：上图即为鲸鱼竞争种群动态的向量场。其中，横坐标为B，纵坐标为F，较陡的直线是水平集的部分，较平缓的直线是水平集的部分。图中的箭线指示该点相轨线前进的方向，箭线的长短表明该点相轨线前进速度的大小。7．实验总结分析

8、不同竞争力下的蓝鲸与长须鲸的种群共存动态状况。本次实习模型涉及到了平衡态、水平集等知识，并得到在竞争力为时的鲸鱼竞争种群动态的向量场图，使得两种群的竞争状态表现的更加直观，这类模型的建立有利于对珍惜生物品种的生存状况及发展前景的分析，可以更好的保护它们。