function chaoticDemoThree %混沌曲线演示 %混沌系统具有初始值敏感性,初始值的微小变化最终会导致混沌系统严重偏离。 %此程序用于演示混沌系统的3个函数形式 % x(i+1)=4*x(i)*(1-x(i)) 值的区间在[0,1] % x(i+1)=1-2*x(i)*x(i) 值的区间在[-1,1] %x(i+1)=cos(k*arccos(x(i)) k阶chebyshev,k为阶数 值的区间在[-1,1] %程序设计:李立宗 %2012年3月3日 % http://blog.csdn.net/superdont % lilizong【at】Gmail clc; clear all; % 设置长度 len=100; %设置初始值 format long; key1=0.7; key2=0.7; key3=0.7; key4=0.7; %第1条曲线 s1=zeros(1,len); s1(1)=key1; for i=2:len s1(i)=4*s1(i-1)*(1-s1(i-1)); end %第2条曲线 s2=zeros(1,len); s2(1)=key2; for i=2:len s2(i)=1-2*s2(i-1)*s2(i-1); end %第3条曲线 3阶 s3=zeros(1,len); s3(1)=key3; for i=2:len s3(i)=cos(3*coth(s3(i-1))); end %第3条曲线 8阶 s4=zeros(1,len); s4(1)=key4; for i=2:len s4(i)=cos(8*coth(s4(i-1))); end subplot(2,2,1),plot(s1,'b-'); title('x(i+1)=4*x(i)*(1-x(i))混沌系统'); subplot(2,2,2),plot(s2,'r-'); title('x(i+1)=1-2*x(i)*x(i)混沌系统'); subplot(2,2,3),plot(s3,'g-'); title('3阶chebyshev函数混沌系统'); subplot(2,2,4),plot(s4,'g-'); title('8阶chebyshev函数混沌系统');