收敛与混沌(迭代).docVIP

第三章 收敛与混沌----迭代 迭代产生的数列可能是：(1)收敛；(2)周期性变化；(3)分岔；(4)混沌。 如：用迭代产生的数列是否收敛？ 答： 0.4000 0.7200 0.6048 0.7171 0.6087 0.7146 0.6119 0.7125 0.6146 0.7106 0.6169 0.7090 0.6190 0.7075 0.6208 0.7062 0.6224 0.7050 0.6239 0.7040 0.6252 不收敛。 §3—1 不动点与迭代 1．什么是迭代 定义：任意给定一个输入，由一个函数表达式得到一个输出；再将作为新的输入，由同一个得到下一个输出；重复。这中对某个函数规则反复将输出作为新输入的重复执行过程就称为迭代。一个迭代过程的数学表示为 其中，称为迭代函数，产生的数列称为迭代数列，称为迭代初值。 迭代函数是关键，迭代数列的变化趋势主要由它决定。 例1：，则迭代式为， 分别取=0, 0.1, 0.8, 1, 2, 99计算，得下面数列： 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, …… 0.1, 0.31, 0.56, 0.74, 0.86, 0.93, 0.96, 0.98 0.99 0.99 …... 0.8 0.89 0.94 0.97 0.98 0.99 0.99 …… 1 1 1 1 1 1 1 ……. 2, 1.41, 1.18, 1.09, 1.04, 1.02, 1.01, 1.00, 1.00, ….. 99, 9.94, 3.15 1.77 1.33 1.15 1.07 1.03 1.01 1.00 1.00 … 结论：迭代函数对于初值保持不动（称为不动点）；对于其余初值（非1的任何正数）都收敛于1 . 问：时，迭代规律？ 2．不动点 定义：若存在点，满足，则称点为迭代函数的一个不动点。 对于同一个，其不动点不一定存在、不一定唯一。 在某个不动点附近取初值，迭代可能收敛于这个不动点，称它为吸引的（如例1中的不动点1）；也可能远离此不动点，称它为排斥的（如例1中的不动点0）。 §3—2 用图示法体现迭代数列的规律 线性迭代： 二次函数迭代： 一种特殊的二次函数迭代：著名的Logistic函数 ，其中，参数在0到4之间取值 对应的Logistic迭代为 本节就以Logistic迭代为例，讨论揭示迭代数列规律的图形方法。你将会看到某写全新的现象和问题，如分岔、混沌等古怪现象。 先求不动点：令，解得两个不动点 与 再分析迭代规律： 若，则 因为“总有”，从而，所以，对于任意初值，都收敛于不动点0； （2）若，则对于任意初值，都收敛于不动点； （3）若，则迭代规律很乱，对于不同的分别呈现诸如收敛、周期性震荡、分岔、混沌之类的从有规律到无规律、又从无规律到有规律等等的非常复杂的、有趣的、怪异的现象。 定义：若迭代数列中，当下标n充分大时，每隔k个数就出现周期性重复，则称此迭代为k—周期。 下面介绍刻画迭代数列规律的三种图形方法 1．线性联结图 横坐标： 纵坐标： 把向邻的点用折线联结。 看书 P26 图 画 图3.1 的程序为： a=3.8;x(1)=0.2; for i=1:20 x(i+1)=a*x(i)*(1-x(i)); end plot(0:20,x) 画 图3.2 的程序为： a=[3.2,3.5,3.5644,3.8284];y =0.2*[1,1,1,1]; for i=1:10000 y=a.*y.*(1-y); end x(1,:)=y; for i=1:20 x(i+1,:)=a.*x(i,:).*(1-x(i,:)); end subplot(2,2,1),plot(0:20,(x(:,1))) subplot(2,2,2),plot(0:20,(x(:,2))) subplot(2,2,3),plot(0:20,(x(:,3))) subplot(2,2,4),plot(0:20,(x(:,4))) 从线性联结图上，看不出当a=3.5644, a=3.8284时，Logistic迭代到底是几---周期的？ 用下面蛛网图就可以看清楚。 2．蛛网图 先计算出，再做下列工作： （1）画曲线和直线； 出发点， 过A做竖直线交曲线于点，过B做水平线交直线于新的点；再过新A做竖直线交曲线于新的点，过新B做水平线交直线于新的点；重复多次。 画 图3.4中第3个图 的程序为： hold on x=