程序:
clear
tspan=0:0.01:100;
y0=[6.60,0.02,4.00]; %给定稳定后的初值
[T,Y]=ode45(@odefun,tspan,y0); %求解常微分方程组(式 9.11)
figure(1);
plot3(Y(:,1),Y(:,2),Y(:,3));
figure(2);
plot(Y(:,1),Y(:,2));
figure(3);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.');
legend('x','y','z');
%子程序:编写混沌方程组
function dy=odefun(t,y)
dy=zeros(3,1);
dy(1)=-y(2)-y(3);
dy(2)=y(1)+0.2*y(2);
dy(3)=0.2+y(3)*(y(1)-3.5);
End