MATLAB求解混沌系统微分方程组

function main() clc close all global c u %定义全局变量 c=3; u=7; [T,Y]=ode45(@fun,[0 100],[3 2 12]); %根据fun定义的微分方程组，使用ode45函数求解微分方程组 %绘制x随时间t变化曲线 figure plot(T,Y(:,1)) %绘制曲线 xlabel('t') %添加横轴标注t ylabel('x') %添加纵轴标注x title('x的值') %添加标题 %绘制y随时间t变化曲线 figure plot(T,Y(:,2)) %绘制曲线 xlabel('t') %添加横轴标注t ylabel('y') %添加纵轴标注y title('y的值') %添加标题 %绘制z随时间t变化曲线 figure plot(T,Y(:,3)) %绘制曲线 xlabel('t') %添加横轴标注t ylabel('z') %添加纵轴标注z title('z的值') %添加标题 figure plot3(Y(:,1),Y(:,2),Y(:,3)) %plot3绘制三维曲线 xlabel('x') ylabel('y') zlabel('z') title('x y z 三维曲线') end function dy=fun(t,y) %fun定义微分方程组 global c u %全局变量，使用之前需要先声明 dy=zeros(3,1); dy(1)=y(2)-y(1); dy(2)=y(1)*y(3)+c; dy(3)=u-exp(y(1)*y(2)); end