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

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