《鱼类物种的数量模型》.docVIP

鱼类物种的数量模型 摘要 本文针对种鱼以及其捕食者种鱼各自随时间的数量变化，估计出种鱼和种鱼之间的相互影响，在此基础上，依据问题中给出的微分方程组，进行数据拟合以及修改，从而最终得出两物种的数量模型。 针对问题一，根据已给数据，运用绘制图形，发现种鱼的数量符合阻滞增长模型（模型），由此建立与的数学关系为，于是运用拟合，计算出种鱼的自然增长率以及环境资源容许的种群最大数量，总体拟合程度较高。 针对问题二，将数值以及初始值带入，运用计算出微分方程的数值解和并作图。 针对问题三，根据已知数据，首先分别绘制及的图像，再将问题二中数值带入，按照问题二的方法计算出两组数据, ,并与原数据比较，发现图形走势基本一致，再根据数据初步拟合的结果显示，拟合后的图像与原图像横坐标相差倍，故通过遍历()邻域内的值优化求解。 然后通过控制变量，采用最小二乘拟合法循环逐个拟合四个参数，用平方和残差值的最小值来评价参数的优劣程度，从而筛选出最优解，并代入微分方程组，将再次计算出的数据与原数据比较，发现拟合程度比较理想。拟合过程中还采用微元的思想，将整个时间区间分成无限个小区间进行积分求和运算，得到的数值解为 ，最后根据求出两个初值分别为。 针对问题四，通过观察问题三绘制出的拟合图像，发现计算出的数据比原数据略大，说明原微分方程组是理想模型，没有考虑起物种自身阻滞作用的项，于是在原有的微分方程基础上增加了项，提高了数据的拟合程度。 针对问题五，分析了微分方程组中的各参数代表的含义，得出不依赖于数据的结论，即无论y取何种初值，图像的形状不变，只是在原有基础上进行平移，基于这个结论， 的数值可解，但对于y的初值无法给定。 本文亮点首先在于问题三中的数据拟合过程，通过对拟定的未优化数据与原数据各自形成的图像特征的比较，确定准确值存在的大致范围，运用最小二乘拟合法逐个拟合优化，最终取得最优解。然后在问题五中，清楚地分析了与、的关系中最重要的一点就是和各时刻对应的数值无关，而是和图像的形状、周期有关。 关键词：阻滞作用、最小二乘拟合、 问题重述 在某水域生长的一种鱼类（A种鱼），该水域出现了一种以A种鱼为食的另一种鱼类（B种鱼），记时刻A种鱼的数量为，B种鱼的数量为。A种鱼与B种鱼的数量相互影响，本题要求通过建立数学模型解决以下五个问题： 根据表格中记录的一段时间内所观测到的A种鱼在若干观测观测的数量，描述A种鱼的数量随时间变化的情况。 如果这两种鱼的数量变化符合以下的微分方程组： 根据题中给出的及初始值，计算t=10,20,…,100时两种鱼的数量。 根据在一段时间内观测到的两种鱼的数量，如果两种鱼的数量变化符合问题二中的微分方程组，求其中的常数以及初始时刻两种鱼的数量。 问题二中的微分方程组能否很好的反映问题三数据？请你修改微分方程组的形式，使得改进后的微分方程组更好的反映问题三中的数据。 如果在问题三中，无法观测到B种鱼的数据（），只能观测到A种鱼的数据（），是否可以求解与问题三与问题四相应的问题。 模型假设 A种鱼的数量变化符合阻滞增长模型； 除了食饵、捕食者关系以及自身阻滞作用，其他因素对两种鱼的数量无影响； 符号说明 N：环境资源所能容纳的最大种群数量； r：种群数量的自然增长率； ：A种鱼在初始时刻的种群数量； ：B种鱼在初始时刻的种群数量； ：平方和残差值的最小值； ：表示的是单位数量A种鱼提供的供养B种鱼的食物量为单位数量的B种鱼消耗的供养B的食物量的倍； ：表示单位数量B种鱼消耗A种鱼作为食物量为单位数量的A种鱼自身减少数量的倍。 模型的建立与求解 4.1问题一 4.1.1问题一的分析 当水域中只有A种鱼时便构成了单种群模型，为了清楚的认识问题一中A种鱼的数量变化规律，根据A种鱼在从0到19之间不同时刻的数量，首先运用绘制出A种鱼数量随时间变化的图像如下（如图4.1.1-1），它大体上描述了该种群数量随时间变化的规律，从而判断出A种鱼符合的模型，并根据数据拟合出模型中的未知参数。 图4.1.1-1 4.1.2模型Ⅰ的建立 根据图4.1.1-1，不难看出A种鱼的增长规律符合阻滞增长模型（logistic模型）。将t时刻的种群数量记为，并将视为连续、可微函数，记初始时刻的数量为,假设A种鱼的自然增长率为常数r，可以注意到r随着数量的增加而下降，即表现为阻滞越来越大。若将r表示为x的函数，则它是减函数，有 ， 对作一个简单的假定，为此，引进环境资源所能容纳的最大种群数量为N，设 带入上述公式，得 上述式右端的因子体现了种群数量自身的增长趋势，因子则体现了种群数量增长的阻滞作用，式构成阻滞增长模