混沌的数值表征与普适结构.docVIP

实验报告 第七次数学试验 实验人 实验题目 混沌的数值表征与普适结构 实验类别 A 理论满分值 100分 实验报告正文 一. 实验题目(题目使用四号, 黑体)clear all; clc; x(1)=1; af=7*rand(1) %af=8+2*rand(1) %af=13+1*rand(1) %af=14.5+0.3*rand(1) %af=15+1*rand(1) for i=1:100 x(i+1)=af*x(i)*exp(-x(i)); end plot(x) 对其中的af分别运行，得到在各区域的参数的值分别是6.7025,9.9412,13.4218， 14.5170,15.9340，且其迭代示意图分别如下： 2、判断进入混沌区域的相关程序如下： clear all; clc; x(1)=1 for af=14.7:0.001:14.8 for i=1:3000 x(i+1)=af*x(i)*exp(-x(i)); end y=round(x(1500:1:end).*1e+6)./1e+6; A=unique(round(x(1500:1:end).*1e+6)./1e+6); a=length(A) if a1000 af; a; break end end fprintf(当af=%2.4f时进入混沌区域\n,af); 运行得结果为：当af=14.7410时进入混沌区域。 根据第1题，分别取其中的af为14.6、14.7、14.7410、14.8进行比较图形可知当af为14.7410时，确实是开始进入混沌区域。 3、当分别取参数的值是6.7025,9.9412,13.4218，14.5170,15.9340时其蛛网迭代图的matlab程序如下： clear all; clc; af=6.7025 %af=9.9412 %af=13.4218 %af=14.5170 %af=15.9340 x2=0:0.01:7 y2=x2 y3=af.*x2.*exp(-x2); x=0.3 for i=400:500 x=af*x*exp(-x) y(i)=x end x1(1)=0.3 y1(1)=0 for i=400:500 x1(2*i+1)=y(i); x1(2*i)=x1(2*i-1); y1(2*i)=y(i); y1(2*i+1)=y1(2*i); plot(x2,y2,r-,x2,y3,k-,x1,y1,b-) pause(0.01) title(af=6.7025的蛛网迭代示意图) end 当af取以上不同值时其蛛网迭代示意图分别如下： 四. 实验结果 1. 参数表 Regions region 1 region 2 region 3 region 4 region 5 Values 6.7025 9.9412 13.4218 14.5170 15.9340 2. 当取何值时开始进入混沌区域, 为什么? 答：当取14.7410时进入混沌区域，根据第1题，分别取其中的af为14.6、14.7、14.7410、14.8进行比较图形可知当af为14.7410时，确实是开始进入混沌区域。 3. 蛛网迭代示意图的作图思想是什么? 答：由初始点A0出发作横轴的垂线交抛物线L于一点B0，再由这一点出发作横轴的平行线，与直线l的交点将为A1，再有A1出发作横轴的垂线交L于B1，再有B1出发作横轴的平行线与l的交点为A2，用此方法，可依次得到点B2，A3，B3，…，直至所有的迭代点。 实验日期: 2012 年6 月14日 实验报告 续表 实验日期: 2012 年6月14日 数学系 10 信息与计算科学 实验报告