王智LogisticEquation的feigenhaum分岔图研究实验报告.docVIP

Logistic Equation的feigenhaum分岔图研究 脑所 201121430048 王智 实验目的： 由Logistic Equation做出Feigenhaum分岔图，并通过分岔图研究该方程的特性。 实验内容： 1、 用计算机实现分叉图； 2、 计算Feigenham’s number； 3、 进一步找到周期3的R值； 4、 研究自相似性。 实验原理： 非线性有限差分方程（Logistic Equation）： （1） 其中R为系统的参数，方程初始条件为; 固定点（fixed point）：满足的点，可解得： （2） 分岔图的基本特性 (1)系统参数(R)的不同给系统带来的影响 初始条件取值满足01,随着R值的增加，系统分别表现出稳定，周期，混沌等不同的行为。 1) 当0R1时，系统只有一个稳定的点(零点)； 2) 当1R3时，Logistic方程迭代的结果趋向于一个稳定的点，此时为周期1的解，系统处于稳定状态； 3）当3R3.449时，方程迭代的结果逐渐趋向于两个稳定的点，此时系统为周期2状态，满足； 4)当3.449R3.556时，系统呈现周期4,8,16…周期倍增现象（period-doubling）； 5）当3.,556 R4时，系统在大部分R值呈现出混沌（chaos）的状态，在部分R值呈现出周期3的状态； 6）当R4时，轨线最终逃逸（escape）到无穷。 (2)初始状态( x0)的不同对系统的影响 1) 系统处于周期等稳定状态时对初始值不敏感。 2) 系统处于混沌状态对初始值非常敏感。 (3) feigenbaum’s number 利用分叉点的R值，计算相邻分岔点的间距之比，所得的极限值即为feigenbaum’s number，计算公式为： （3） 方法实现： 实验平台：MATLAB 1、 计算机实现分叉图 程序运行首先生成Logistic方程的Feigenhaum分岔图。横轴是参数R，取值范围是2.6—4，取样间隔为1e-4；纵轴是变量x，取值范围是0—1，方程迭代次数选择250。 当R从2.6增加到4的时候，迭代Logistic方程 ，将迭代150次至迭代250次的x的取值在屏幕上输出最终得到Logistic方程的Feigenhaum分岔图。 2、计算feigenbaum’s number 将每次迭代计算出的x值存入X矩阵，X矩阵为250行，代表程序迭代250次，14001列，代表14001个R值，因此该矩阵的每一列即代表了选取R为定值后方程的迭代过程。由于所选取的R值的精度有限，因此选择8周期分叉点与16周期分叉点近似计算feigenbaum’s number。满足的R值即认为是周期8对应的R值，考察迭代160次~180次的x值（共21个x值），超过14个x值满足上述关系，程序则认定该R值满足周期8关系，根据该方法同样可以找出周期16对应的R值，进而计算得到feigenbaum’s number。 3、找到周期3对应的R值 周期3 满足关系，因此当满足的R值即认为是周期3对应的R值，同样考察迭代160次~180次的x值（共21个x值），超过14个x值满足上述关系，程序则认定该R值满足周期3关系。 4、自相似性研究 在得到Logistic方程的Feigenhaum分岔图后，通过对X,Y坐标轴取值范围的逐渐缩小，可以逐渐放大分岔图的局部，进而研究自相似性。 实验结果呈现： 1、系统参数(R)的不同给系统带来的影响 初始值选择x=0.6，首先分别作出R=0.5, R=1.,4, R=2.8, R=3.2, R=3.5, R=3.9时x随迭代次数增加而表现出的状态，横轴代表迭代次数t，纵轴代表x： 1）R=0.5时，x逐渐趋于稳定值x=0； 图1. 2）R=1.4时，x逐渐趋于稳定值x=0.28； 图2. 3）R=2.8，x交替趋近于稳定值x=0.64 图3. 4）R=3.2时，x呈现周期2，两个稳定值分别为x=0.8和x=0.52 图4. 5）R=3.5时，x呈现周期4，四个稳定值大致分别为x=0.88,x=0.82,x=0.51,x=0.39 图5. 6）R=3.9时，x呈现混沌状态 图6. 2、计算机实现分叉图 R从2.6增加到4，，迭代Logistic方程 ，把对应的x取值在屏幕上输出作图，横轴代表R，纵轴代表x 图7. Logistic方程的Feigenhaum分岔图 3、计算Feigenhaum’s number 计算Feigenhaum’s number所得结果的准确