
草原生态平衡的建模分析(共8页).docx
7页精选优质文档-----倾情为你奉上草原生态平衡的建模分析王X(西安电子科技大学 电子工程学院,陕西 西安 )摘要: 本文研究了草原生态系统中黄羊,狼群和草场之间的制约与促进的关系.文章首先把草原生态系统看作是一个离散型动力系统,并将其因素简化为黄羊,狼群和草场.再运用Logistic 模型和Volterra 模型对该动力系统进行建模,分析其三个因素之间的具体关系,进而对草原生态系统的平衡进行评价和预测.关键词: 动力系统;Logistic模型;Volterra 模型Modeling and analysis of grassland ecosystemQING Wang(School of Electronic Engineering, Xidian Univ., Xi’an , China;)Abstract: This paper studies the relationship and advance grassland ecosystems Mongolian gazelle, wolves and pasture between article first prairie ecosystem as a discrete dynamical systems, and simplify their factors Mongolian gazelle, wolves and pasture Logistic model and then use the power system Volterra model to model, analyze the specific relationship between the three factors, and thus balance the prairie ecosystem evaluation and forecasting.Key words: Power systems; Logistic model; Volterra model高原草原生态环境中主要生活着黄羊群和狼群.黄羊的过度繁殖会引起草原退化和狼群数量增加,狼群的增多又会导致黄羊群数量减少.草原,黄羊群,狼群这三者之间有着相互制约又相互促进的关系.试建立草原上“草场—黄羊—狼群”的生态模型,并给出保持生态平衡的建议.1 问题假设1.1草场的基本假设a.草场总面积为1000 平方公里.b.每平方公里的草场在供养50 只以下的黄羊时,草场不退化.当黄羊数量超过每平方公里50 只时,草场衰减,衰减率与黄羊超过50 只的数量成正比,比例系数为0.0001.1.2黄羊种群基本假设a.当前黄羊种群数量为60000 只.b.草场充足,没有狼群的情况下,黄羊净增长率为0.1.c.草场不充足会导致黄羊种群的繁殖率下降,下降率与每平方公里平均黄羊数量减50 只成正比,比例系数为0.001.d.狼群会减少黄羊的数量.1.3狼群的基本假定a.当前狼群总量为100 只.b.黄羊总群数量与狼群数量之比超过300/1 时,狼群的净增长率为0.01.低于300/1 时,狼群的繁殖率下降,下降率与狼群总量与黄羊总量的比值成比例,比例系数为0.5.c.每只狼平均每年吃掉20 只黄羊.2 符号说明n 年份.S(n) 第n 年草场的总面积.X(n) 第n 年黄羊的种群数量.Y(n) 第n 年狼群的种群数量.r 草场面积减少率与每平方公里平均黄羊数量减50 只的数量的比例系数.m 草场充足,没有狼群的情况下,黄羊净增长率.t 黄羊总群数量与狼群数量之比超过300/1 时,狼群的净增长率.a 每只狼平均每年吃掉的黄羊数.3 问题分析与建模第一步:按照对草场的基本假设,易知当 xnsn≤50 时,有sn+1=sn;当xnsn>50时,sn+1=sn*1-rxnsn-50;第二步:按照对黄羊种群的假设,利用Logistic 模型和Volterra 模型,可以得出黄羊数量的状态方程为:当xnsn≤50 时,有xn+1=1+m*xn-a*Yn当xnsn>50 时,有xn+1=xn*1-0.001*xnsn-50其中,Logistic 模型刻画了草场的退化对黄羊数量增长率的影响,0.001*xnsn-50就是Logistic 模型的一项.而a*Yn则刻画了狼群对黄羊数量增长的影响.第三步:按照对狼群的假设,利用Logistic 模型和Volterra模型,可以得出狼群的状态方程为:当Ynxn<1300 时,Yn+1=1+t*Yn ,在食物充足时,狼群将呈一次函数的形式增长.当Ynxn≥1300 时,Yn+1=1-0.5*xnYn*Yn ,其中,0.5*xnYn是食物对狼群增长率的影响.到此,得到了描述“草原—黄羊—狼群”生态系统模型的完整数学表达式.用Matlab 画出各个量之间的关系图形.Matlab 编程如下:clc;clear;s(1)=1000;x(1)=60000;y(1)=100;i=1;k=200;t=1:1:k;while i
