php中文网 | cnphp.com

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 461|回复: 0

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

[复制链接]

2571

主题

2578

帖子

9095

积分

管理员

Rank: 9Rank: 9Rank: 9

UID
1
威望
0
积分
6402
贡献
0
注册时间
2021-4-14
最后登录
2024-3-24
在线时间
652 小时
QQ
发表于 2022-6-16 17:56:59 | 显示全部楼层 |阅读模式
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





上一篇:python 绘制时钟
下一篇:python爬取京东商品图片
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|php中文网 | cnphp.com ( 赣ICP备2021002321号-2 )51LA统计

GMT+8, 2024-3-28 18:36 , Processed in 0.174222 second(s), 36 queries , Gzip On.

Powered by Discuz! X3.4 Licensed

Copyright © 2001-2020, Tencent Cloud.

申明:本站所有资源皆搜集自网络,相关版权归版权持有人所有,如有侵权,请电邮(fiorkn@foxmail.com)告之,本站会尽快删除。

快速回复 返回顶部 返回列表