差分方程模型.pptVIP

鲑鱼数量的变化问题 * * 海洋中鱼的数量通常是按繁殖期的长短呈周期变化的。以太平洋里的鲑鱼为例，其生长繁殖过程大致是，成年的鲑鱼产下大量的卵，在卵成长为幼鱼和幼鱼长大的过程中，相当大的部分被成年的鱼吃掉，剩下来的还要被恶劣的环境淘汰一些，而成年的鱼在产卵后则活不了多久就会死掉。 一 问题的提出 试建立鲑鱼产卵期到来之前，其数量变化规律的数学模型。 二 生长特点 1 呈周期性变化； 2 在每个周期里，经过了从卵、幼鱼到成年鱼 的变化过程。 一般地，长期观察是呈离散变化，而在每个离散时间段里呈连续变化。 如：树木的生长、冰箱温度的变化等， 嵌入式模型 嵌入式模型 它把一个个短期内描述连续变化过程的微分方程，嵌入到一个长期的描述离散变化规律的差分方程中，而那些描述短期演变过程的微分方程在定性上应该是相同的，只是在定量上参数与初始条件有所改变。 三 符号的说明 ：第n个繁殖期（周期）开始时成年鲑鱼（鲑鱼）的数量，以条数计，n=1,2,…； ：在每个周期内，时刻t 幼鱼的数量； 为了反映每个周期末和下个周期开始时的突变性，引入下列记号： 可以很小。在区间 内允许数量上的突变 四 模型的假设 1 成正比，比例系数为 ，表示每 条鱼的产卵量； 2 单位时间内 减少的比例与 成正比， 比例系数为 ，表示鲑鱼吞食幼鱼的能力； 3 成正比，比例系数为 ，表示在 繁殖期末幼鱼存活长成鲑鱼的比例。 五 模型建立 根据假设条件容易写出 （1） （2） （3） 方程（2）的解为 （4） 将（1），（4）代入（3）式得 （5） 若记 （6） 则方程（5）可写作 （7） 差分方程（7）是将每周期内的微分方程（2）嵌入 （1）、（3）得到的。这种嵌入式模型的一般形式 可以表为 （8） 差分方程（7）无法求出 的显式表达式，只能递推 求数值解。例如设 （表示1个数量单位， 比如 条），第1代（ ） 鲑鱼吞食掉90%的幼鱼 即 ，代入（4），（6）是可以算出 若 分别取 ，则由（6）式 将 代入（7）式递推计算 ，考察鲑鱼数量的 周期变化规律，结果见表。 2.396 1.725 21 0.1983 0.3744 0.7007 10 0.4308 0.3562 0.6990 20 2.226 1.698 0.6961 9 1.796 1.727 0.6989 19 0.2882 0.5584 0.7037 8 0.1821 0.3581 0.6990 18 2.022 1.458 0.6912 7 2.272 1.724 0.6989 17 0.2278 0.7562 0.7117 6 0.5919 0.3545 0.6991 16 2.151 1.265 0.6778 5 1.611 1.730 0.6988 15 0.2625 0.8876 0.7329 4 0.1526 0.3614 0.6992 14 2.074 1.156 0.6402 3 2.367 1.719 0.6986 13 0.7115 0.9611 0.7906 2 0.3691 0.3488 0.6996 12 1.500 1.100 0.500 1 1.884 1.739 0.6979 11 1.000 1.000 1.000 0 1.676 1.726 31 0.3225 0.3568 40 0.1771 0.3569 30 1.960 1.726 39 2.287 1.726 29 0.2137 0.3568 38 0.3147 0.3568 28 2.186 1.726 37 1.973 1.726 27 0.6481 0.3568 36 0.2167 0.3569 26 1.557 1.726 35 2.178 1.726 25 0.1449 0.3568 34 0.6526 0.3567 24 2.394 1.726 33 1.552 1.726 23 0.4531 0.3568 32 0.1443 0.3572 22 按（7）式（b=2.3和不同的a）计算的 由表可知，对于 趋向稳态值0.699，即初值 的70%；对于 交替地趋向两个稳态值0.3568 和1.726，对于 则难以看出什么规律。 六 平衡点及稳定性分析 为了研究对于不同的 ，鲑鱼数量 的变化 规律，我们利用差分方程求解的有关结果讨论（7） 的平衡点及稳定性。 方程（7）的平衡点 满足 （9） 注： 也是方程（7）的平衡点，但容易验证它不 是稳定的 ，不再讨论，以后平衡点均指非零的。 （9）的非零解为 （10） 平衡点稳定的条件是 这里 因为 所以当 ，即