差分方程模型说课讲解.ppt

《差分方程模型说课讲解.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、差分方程模型二生长特点1呈周期性变化；2在每个周期里，经过了从卵、幼鱼到成年鱼的变化过程。一般地，长期观察是呈离散变化，而在每个离散时间段里呈连续变化。如：树木的生长、冰箱温度的变化等，嵌入式模型嵌入式模型它把一个个短期内描述连续变化过程的微分方程，嵌入到一个长期的描述离散变化规律的差分方程中，而那些描述短期演变过程的微分方程在定性上应该是相同的，只是在定量上参数与初始条件有所改变。三符号的说明：第n个繁殖期（周期）开始时成年鲑鱼（鲑鱼）的数量，以条数计，n=1,2,…；：在每个周期内，时刻t幼鱼的数量；为了反映每个周期末和下个周期开始时的突变性，引入下列记号：可以很小。在区间内

2、允许数量上的突变四模型的假设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

3、）是可以算出若分别取，则由（6）式将代入（7）式递推计算，考察鲑鱼数量的周期变化规律，结果见表。01.0001.0001.000110.69791.7391.88410.5001.1001.500120.69960.34880.369120.79060.96110.7115130.69861.7192.36730.64021.1562.074140.69920.36140.152640.73290.88760.2625150.69881.7301.61150.67781.2652.151160.69910.35450.591960.71170.75620.2278170.6989

4、1.7242.27270.69121.4582.022180.69900.35810.182180.70370.55840.2882190.69891.7271.79690.69611.6982.226200.69900.35620.4308100.70070.37440.1983211.7252.396220.35720.1443320.35680.4531231.7261.552331.7262.394240.35670.6526340.35680.1449251.7262.178351.7261.557260.35690.2167360.35680.6481271.7261

5、.973371.7262.186280.35680.3147380.35680.2137291.7262.287391.7261.960300.35690.1771400.35680.3225311.7261.676按（7）式（b=2.3和不同的a）计算的由表可知，对于趋向稳态值0.699，即初值的70%；对于交替地趋向两个稳态值0.3568和1.726，对于则难以看出什么规律。六平衡点及稳定性分析为了研究对于不同的，鲑鱼数量的变化规律，我们利用差分方程求解的有关结果讨论（7）的平衡点及稳定性。方程（7）的平衡点满足（9）注：也是方程（7）的平衡点，但容易验证它不是稳定的，不再讨

6、论，以后平衡点均指非零的。（9）的非零解为（10）平衡点稳定的条件是这里因为所以当，即时是稳定平衡点，而当时不稳定。这个结果表明，是否稳定只取决于与无关。而，注意到和的含义可知表示的是鲑鱼从一个周期到下一周期增长关系的一个因素（增长率还与有关），正是这因素决定了的稳定状态情况。根据上述分析，当时，稳定，且若由（10）可得，而当时不稳定这与前面的数值结果（见表）是一致的。为了进一步研究（如）时的变化情况，应该考察方程（参考倍周期收敛的相关内容）（11）其中的具体形式由方程（7）给出。首先用无量纲化方法简化方程，令（12）则方程化为（7）化为（13）（14）显然，当时是方程（13）的

7、稳定平衡点，而时不稳定。下面讨论的情况。考察方程（15）其中由（13）给定。方程（15）的平衡点除以外还有和，满足和（16）由（16）不难得到（17）于是是方程（18）的两个根。若记函数（19）则曲线和直线有3个交点，其横坐标是，1和（见图）。当时用数值方法可以算出（20）是方程（15）稳定平衡点的条件是O12经过比较精密的计算得到，当（21）时上述条件成立。这个结果表明，在条件（21）下方程（13）给出的序列是2倍周期稳定的，即子序列和当时分别趋向于和代回到变量，由（14）式可