2009年8月吉林师范大学学报(自然科学版)
.3第3期Journal of Jilin Normal University (Natural Science Edition)Aug.2009
收稿日期:2009 06 23 基金项目:四平市科技发展计划项目(2009 2010)作者简介:张双红(1978 ),女,吉林省四平市人,现为吉林师范大学数学学院讲师,硕士.研究方向:数学研究与教学.
Logistic 模型的Matlab 计算与可视化
张双红
(吉林师范大学数学学院,吉林四平136000)
摘 要:探讨了Logistic 模型的Matlab 计算及其可视化问题,尤其是迭代收敛过程蛛网图形的可视化,将Logis tic
模型差分形式从收敛、分叉到混沌的过程直观的展现出来.
关键词:Logistic 模型;差分方程;Matlab;可视化
中图分类号:O242 文献标识码:A 文章编号:1674 3873 (2009)03 0097 03
0 引言
Logistic 模型(1)在种群系统中有着广泛的应用.实际问题若涉及时间段时则其差分形式(2)比连续时间更符合实际.
d x d t =rx (1-x N
)(1)x k +
1-x k =rx k (1-x k
N ),k =0,1,
(2)
取N =1,x 0=0.1,r =0.3,1.8,2.5,按照迭代方程(2)利用Matlab 计算x k 并画图(如图1),可以看到序列{x k }并不是总收敛的,并且对初值具有敏感性[1].
图1 Logistic 模型差分形式的Matlab 模拟
1 Logistic 模型的简化差分形式与迭代可视化
若令b =r +1,y k =
r
(r +1)N x k
,则(2)式简化为如下非线性差分方程
y k +1=by k (1-y k ),k =0,1,
(3)
取初始值为0.2,并选取不同的b 值利用Matlab 进行数值计算,结果我们可以看到当b =3.3,b =3 45,b =3.55时迭代序列分别收敛到2个,4个,8个极限点(如表1),呈现2n
倍周期收敛的情形.
对于Logistic 模型简化差分形式(3)的收敛性,可以利用计算机直接进行模拟[2],也可以通过蛛网图形来直观展现这一过程.编写Matlab 程序对于相同初始值,不同的b 和迭代次数进行蛛网迭代(如图2).显然图2中(a)和(b)迭代序列收敛到一个不动点,(c)收敛到两个不动点,(d)收敛到4个不动点,(e)收敛到8个不动点,与数值计算结果是一致的.
!
97!