matlab第2次课.ppt

实验四 函数的迭代、混沌与分形 实验目的 理解迭代的基本含义 掌握迭代数列的系列图形表示方法 以一类特殊二次函数（Logistic函数）为例，掌握二次函数迭代数列的收敛性分析方法 熟悉编写函数迭代的Matlab程序 了解二元函数迭代的方法及其图形特征 作业说明 1、分段函数的定义 2x , 0=x=1/2, f(x)= 2(1-x), 1/2x=1; f=inline(‘2*x*(x=0x=1/2)+2*(1-x)*(x1/2x=1)’); 但画蜘蛛网图时不能直接使用。 2、迭代时，xn+1与xn的前8位有效数字一致时终止计算 用WHILE(vpa(xn,8)~=vpa(f(xn),8)) 输出前8位有效数字，fprintf(‘%.8g’,xn) syms x; f=inline(3.5*x*(1-x)); x0=0.12; for i=1:1:100 if i50 plot(i,f(x0),.); end; x0=f(x0); hold on; end; hold off syms x; f=inline(3.5*x*(1-x)); x0=0.12; for i=1:1:100 if i50 plot(i,f(x0),.); end; x0=f(x0); hold on; end; hold off * * 实验四 函数的迭代、混沌与分形 1、 定义 给定某个初值，反复作用以同一个函数的过程称为迭代 ，一般形式为 它生成了一个序列{ }，称为迭代序列． 2、迭代序列的收敛性 设函数 满足： (1)对任意 ； 在( )内可导，且存在常数 使得 则当初值 时，由 生成的迭代序列收敛． 问题1：如果迭代序列收敛，收敛点会满足怎样的条件？ 3.分式线性函数的迭代 例： 先取初值x0=5.5 f=inline((25*x-85)/(x+3));%先定义函数 x0=5.5; for i=1:1:20 x0=f(x0); fprintf(%g,%g\n,i,x0); end 迭代次数 迭代序列 迭代次数 迭代序列 17. 20 16.9711 10 17. 19 16.9281 9 17. 18 16.8218 8 17. 17 16.5642 7 16.9999 16 15.9668 6 16.9997 15 14.7125 5 16.9993 14 12.5529 4 16.9981 13 9.85437 3 16.9954 12 7.5641 2 16.9884 11 6.17647 1 取其它的初值做试验 14 收敛于17 1000 14 收敛于17 100 12 收敛于17 20 17 收敛于17 6 19 收敛于17 5.1 0 收敛于5 5 19 收敛于17 4.9 17 收敛于17 4 17 收敛于17 0 16 收敛于17 -20 16 收敛于17 -500 16 收敛于17 -40000 得到收敛点的迭代次数 收敛性 初值 结论：只要初值不取为5，迭代序列总收敛于17。 易知，f(x)的不动点恰好是17与5。5称为排斥点，17称为吸引点。 问题2 为何17是吸引点，5是排斥点？ 例1 用分式函数的迭代法近似计算 4.迭代的可视化（蜘蛛网图） (x1 y1) (x2 y2) (x3y3) (x4 y4) f=inline((25*x-85)/(x+3)); x=[];y=[]; x(1)=5.5; y(1)=0;x(2)=x(1);y(2)=f(x(1)); for i=1:100 x(1+2*i)=y(2*i); x(2+2*i)=x(1+2*i); y(1+2*i)=x(1+2*i); y(2+2*i)=f(x(2+2*i)); end plot(x,y,r); hold on; syms x; ezplot(x,[0,20]); ezplot(f(x),[0,20]); axis([0,20,0,20]); hold off 5. 认识混沌 迭代序列如果不收敛，会出现什么情况？ 1. 迭代次数充分大时，迭代序列出现周期性重复 k称为该序列的周期 2. 序列没有规律、杂乱无