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

admin

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 三维曲线')
msgbox('MATLAB编程答疑，请加QQ: 1530497909','MATLAB答疑','help')
web http://url.cn/TKcdXk -browser
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