matlab Lorenz混沌系统求解

最新推荐文章于 2023-09-10 01:14:28 发布

studyer_domi

最新推荐文章于 2023-09-10 01:14:28 发布

阅读量931

点赞数 3

分类专栏： Matlab系列案例文章标签： matlab Lorenz 混沌系统求解 Powered by 金山文档

本文链接：https://blog.csdn.net/qingfengxd1/article/details/129091835

版权

Matlab系列案例专栏收录该内容

494 篇文章 266 订阅

订阅专栏

该文使用MATLAB的ode45函数求解Lorenz混沌系统的动态行为，通过绘制X、Y、Z三个维度的轨迹以及3D图展示混沌特性。此外，还进行了Logistic映射的迭代过程，探讨了参数r对系统稳定性的影响。

摘要由CSDN通过智能技术生成

1、内容简介

matlab Lorenz混沌系统求解

649-可以交流、咨询、答疑

2、内容说明

略

3、仿真分析

clc
close all
clear
sigma = 10; alpha = 30; beta = 8/3;
x0 = 1; y0 = 1; z0 = 1;
[T1,Y1] = ode45(@(t,x) Lorenzfun(t,x,sigma,alpha,beta),[0,100],[x0;y0;z0]);
%%
figure
plot(T1,Y1(:,1))
xlabel time
ylabel X

figure
plot(T1,Y1(:,2))
xlabel time
ylabel Y

figure
plot(T1,Y1(:,3))
xlabel time
ylabel Z

figure
plot3(Y1(:,1),Y1(:,2),Y1(:,3))
xlabel x1
ylabel x2
zlabel x3

clear;
r=2:0.01:4;%r为横坐标，定义域为-2到4
x=0.5;%对任意的r给定0-1初值？
for i=1:5000 %迭代很多次达到稳定
    x=r.*x.*(1-x);%迭代的式子
end
figure;
hold on;
for i=1:2000%这些循环是要画在图上的
    x=r.*x.*(1-x);
    plot(r,x,'k.','markersize',1);%描点画图
end