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

% Duffing2 吸引子三维相空间图，这里用四阶 Runge-Kutta 法得到微方程的离散序列 
% 使用平台 - Matlab6.5 / Matlab7.0 
% 作者：陆振波，海军工程大学 
% 欢迎同行来信交流与合作，更多文章与程序下载请访问我的个人主页 
% 电子邮件：luzhenbo@yahoo.com.cn 
% 个人主页：http://luzhenbo.88uu.com.cn 
% 
% 方程表达式 
% dx/dt = y 
% dy/dt = -delta*y + a*x*(1-x^2) + f*cos(z) 
% dz/dt = omega 
 
clc 
clear 
close all 
 
delta = 0.05; 
a = 0.5; 
f = 7.5; 
omega = 1; 
 
y = [-1,0,1];          % 起始点 (1 x 3 的行向量) 
h = 0.05;             % 积分时间步长 
 
k1 = 8000;            % 前面的迭代点数 
k2 = 6000;             % 前面的迭代点数 
 
z = DuffingData2(y,h,k1+k2,delta,a,f,omega); 
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('Duffing attractor');