差分形式的阻滞增长模型matlab,（精）差分形式的阻滞增长模型.ppt

差分方程模型 -----差分形式的阻滞增长模型 平衡点及其稳定性 倍周期收敛 * 前面我们介绍了人口数量增长的阻滞型模型: 它还可以用来近似描述其他受环境约束的事物的增长规律,如:种群数量的增长,传染病的传播,耐用消费品在有限市场上的销售等. 此模型的意义是:在t时刻单位时间内的人口数量的变化量仅仅与此时的人口数量x有关(等于右边的值),其中的r表示人口的固有增长率,N表示能容纳的最大人口数. 有时我们将时间离散化来研究可能方便些.例如:有些生物(比如鱼)每年在固定的时间繁殖,我们用繁殖周期作为时段来研究其增长规律比我们简单地以连续时间处理应该更好些.我们类似认为: 经过 单位时间，即一个繁殖周期的种群数量的增长量仅仅与前一个时期的种群数量有关,且有类似于上面的表达式. 于是模型为 即 这是一个一阶非线性差分方程.对于给定的初值,我们可以从这个递推公式运用计算机很容易地计算出一些xk.这是在计算机出现以后的一个新的特点. 但是我们更关心的是当时间趋于无穷时,即k趋于无穷时, xk的极限如何. 数值计算 .7357 .6418 .7223 .6558 .7123 .6655 .705 .672 0.3 3.2 .6555 .6260 .6626 .6155 .6737 .597 .692 .448 0.2 2.8 .5834 .5831 .584 .582 .586 .576 .600 .504 0.3 2.4 .4444 .4444 .444 .443 .438 .419 .369 .288 0.2 1.8 .171 .172 .174 .176 .179 .182 .186 .192 0.2 1.2 8 7 6 5 4 3 2 1 0 b\k .5210 .7953 .5384 .7859 .5664 .7702 .5965 .7522 .6222 3.2 .6411 .6450 .6402 .6462 .6386 .6481 .6362 .6510 .6323 2.8 .58333 .58333 .58333 .58333 .58333 .58333 .58333 .58335 .5833 2.4 .4444 .4444 .4444 .4444 .4444 .4444 .4444 .4444 .4444 1.8 .1671 .1672 .1673 .1675 .1677 .168 .1685 .169 .170 1.2 17 16 15 14 13 12 11 10 9 b\k 归纳看到的现象 我们很容易求得差分方程(2)的平衡点为0和(b-1)/b.它们分别对应于差分方程(1)的平衡点0和N. 我们将这个差分方程(2)在平衡点附近展开,有 注意到b=1+r>1,平衡点0是不稳定的. 略去高阶项得 因此当|2-b|<1时,平衡点x*＝ (b-1)/b是稳定的. |2-b|>1即 b>3时,平衡点x*是不稳定的. 注:由于b=r+1,结论表明只有r<2时y*=N才是差分方程(1)的平衡点;这与微分方程不同,微分方程中y=N是稳定的平衡点(没有条件). 事情至此好像结束了.当然,我们还可以进一步判断稳定的平衡点是否为全局稳定的.但是,数值计算表明,对于有些b值,平衡点不稳定,但是xk好象在某几个值附近循环摆动. 我们只需要用计算器多迭代计算几次即可. 当b<3时,平衡点x*是稳定的. 当b>3时,平衡点x*是不稳定的.如果序列xk存在两个收敛的子列我们就称之为2倍周期收敛. 我们称之为单周期收敛. 一般地,我们记(2)式为 所谓2倍周期收敛的点就是(4)式的平衡点:即满足 本例 (2)的平衡点为0和(b-1)/b.仍然满足上式,现在我们求的是另外两个根: 或 故 故这两个平衡点具有相同的稳定性.且 结论: *