www.pudn.com > ChaosToolbox1p0_trial.rar > Main_Rossler.m

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