matlab Lorenz混沌系统求解

Posted 2023-03-04 studyer_domi

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

4、参考论文

略