MATLAB求解混沌系统微分方程组
function main()clc
close all
global c u%定义全局变量
c=3;
u=7;
=ode45(@fun,,); %根据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
页:
[1]