www.pudn.com > Chaos_Prediction.rar > Main_Lorenz.m

 
% Lorenz 吸引子三维相空间图，这里用四阶 Runge-Kutta 法得到微方程的离散序列 
% 使用平台 - Matlab7.0 
% 作者：陆振波，海军工程大学 
% 欢迎同行来信交流与合作，更多文章与程序下载请访问我的个人主页 
% 电子邮件：luzhenbo@yahoo.com.cn 
% 个人主页：http://luzhenbo.88uu.com.cn 
 
clc 
clear all 
close all 
 
%-------------------------------------------------------------------------- 
% 方程表达式 
% dx/dt = sigma*(y-x) 
% dy/dt = r*x - y - x*z 
% dz/dt = -b*z + x*y 
 
sigma = 16;             % Lorenz 方程参数 a = 16 | 10 
b = 4;                  %                 b = 4 | 8/3 
r = 45.92;              %                 c = 45.92 | 28     
 
y = [-1,0,1];        % 起始点 (1 x 3 的行向量) 
h = 0.01;            % 积分时间步长 
 
k1 = 8000;           % 前面的迭代点数 
k2 = 3000;           % 后面的迭代点数 
 
data = LorenzData(y,h,k1+k2,sigma,r,b); 
data = data(k1+1:end,:); 
 
%-------------------------------------------------------------------------- 
 
X = data(:,1); 
Y = data(:,2); 
Z = data(:,3); 
 
figure 
plot3(Z,Y,X); 
xlabel('Z');ylabel('Y');zlabel('X'); 
title('Lorenz attractor');