实验六机理模型与平衡原理解析.docVIP

西北农林科技大学实验报告 学院名称：理学院 专业年级：2013级信计1班 课 程：数学模型与数学建模 报告日期：2015年12月15日 实验六 机理模型与平衡原理 实验目的 如果对所研究的问题了解的比较深入，知道产生现象的内在的机理，那么依据机理建模，则模型具有更好的可靠性和广泛性。不考虑随机因素，假设每一时刻是确定的如果对系统状态的观测和描述只在离散的时间点上，则构成差分方程模型；如果考虑系统随时间连续变化，则是微分方程模型。本节主要以这两类方程为例，介绍用MATLAB软件求解机理模型的基本方法。 一、差分方程模型 实验题目 由一对兔子开始，一年可以繁殖出多少只兔子？如果一对兔子每个月可以生一对小兔子，兔子在出生两个月后就具有繁殖能力，由一对刚出生一个月的兔子开始，一年内兔子种群数量如何变化。求这个种群的稳定分布和固有增长率。 实验内容 解 假设 兔子每经过一个月底就增加一个月龄； 月龄大于等于2的兔子都具有繁殖能力； 具有繁殖能力的兔子每一个月一定生一对兔子； 兔子不离开群体（不考虑死亡） 记第n个月初的幼兔（一月龄兔）数量为a0(n）,成兔（月龄大于等于2）数量为a1(n），则兔子总数为a(n)= a0(n）+a1(n）,平衡关系为： 建立模型： 这个一阶差分方程的矩阵表达式为 其中 ， 利用迭代方法求数值解，也就是按时间步长法仿真种群增长的动态过程，模拟幼兔和成兔占整体比例随时间的变化。 ?a=[0?1;1?1];?x=[1?0];for?k=2:12y=a*x(:,k-1);x=[x?y];end?zz=repmat(sum(x),[2?1]);z=x./zz;t=1:12;?plot(t,x(1,:),r*,t,x(2,:),b*),grid;?plot(t,z(1,:),r*,t,z(2,:),b*),grid; 由数值模拟结果可见，兔子数量递增，但是幼兔和成兔在种群中所占比例很快会趋于一个极限。由线性差分方程的定性理论可知，这个极限就是对应于差分方程系数矩阵A主特征值得归一化特征向量。因为A是非负矩阵，由矩阵理论知，A主特征值是正实数，是最大的特征值。 ?[v,d]=eig(a)v?=???-0.8507????0.5257????0.5257????0.8507d?=???-0.6180?????????0?????????0????1.6180 ?max(max(d))ans?=????1.6180 ?v(:,2)/sum(v(:,2))ans?=????0.3820????0.6180 由数值计算可知，系数矩阵模A最大的特征值r=1.618，生物上称之为种群的内禀增长率，是个大于一的实数。因此种群数量随时间递增。相应的归一化的特征向量的两个分量0.382和0.618正是幼兔和成兔在种群中所占比例趋近的稳定值，生物上称之为种群的稳定分布。 从这个例子的讨论可见，数值模拟能帮我们对系统的变化有更直观的认识。但是只有通过计算方程组系数矩阵的特征值和特征向量，运用差分方程定性分析才能对解得渐进性质给出确定的结论。 一、微分方程模型 实验题目 蓝鲸的內禀增长率每年估计为5%，估计蓝鲸的最大环境载量为150000条，磷虾是蓝鲸喜欢的一种食物。磷虾的最大饱和种群为500吨/英亩*。当缺少捕食者，环境不拥挤时，磷虾种群以每年25%的速率增长。磷虾500吨/英亩可以提高蓝鲸2%的年增长率，同时150000条蓝鲸将减少磷虾10%的年增长率。 组建一个蓝鲸和磷虾的动态模型，模拟两个种群随时间的变化。假设初始状态为蓝鲸5000条，磷虾750吨/英亩； 确定蓝鲸与磷虾是否可以长期共存； 假设捕捞使得蓝鲸只剩下它的平衡态的5%，而磷虾保持平衡态的数量。描述一旦停止捕捞将发生什么情况。蓝鲸恢复需要多长时间？磷虾群体将发生什么变化？进一步，给出蓝鲸种群恢复时间对它所受伤害程度的依赖关系。 实验内容 解 （1）假设 蓝鲸和磷虾单独存在时，两种群都依靠有限的自然资源增长，即遵循Logistic模型； 蓝鲸捕食磷虾，蓝鲸平均增长率的增加量正比于单位区域内磷虾数量，磷虾平均增长率的减少量正比于蓝鲸数量 记N1(t)为蓝鲸在t时刻的种群数量（条），N2(t)为磷虾在t时刻的种群质量（吨/英亩），于是，依据假设，建立蓝鲸和磷虾两个种群的动态模型 其中r1=0.05和r2=0.25分别表示蓝鲸和磷虾种群的内禀增长率，k1=150000(条）和k2=500（吨/英亩）分别表示蓝鲸和磷虾种群环境载量，a12=0.02/500表示每英亩每吨磷虾对蓝鲸平均增长率的改变，a21=0.10/150000表示每条蓝鲸对磷虾平均增长率的