% 計算混沌時間序列 Lyapunov 指數 - Lorenz 吸引子% 使用平臺 - Matlab6.5 / Matlab7.0% 作者：陸振波，海軍工程大學% 歡迎同行來信交流與合作，更多文章與程序下載請訪問我的個人主頁% 電子郵件：luzhenbo@yahoo.com.cn% 個人主頁：http://luzhenbo.88uu.com.cnclcclearclose all%-----------------------------------------------------------------% 產生 Lorenz 混沌時間序列sigma = 16;             % Lorenz 方程參數 ab = 4;                  %                 br = 45.92;              %                 c            y = [-1,0,1];           % 起始點 (1 x 3 的行向量)h = 0.01;               % 積分時間步長k1 = 30000;             % 前面的迭代點數k2 = 6000;              % 后面的迭代點數z = LorenzData(y,h,k1+k2,sigma,r,b);x = z(k1+1:end,1);      % 時間序列%-----------------------------------------------------------------% Lorenz 入口參數tau = 10;               % 時延m = 3;                  % 嵌入維taumax = 300;           % 最大離散步進時間P = 100;                 % 序列平均周期fs = 1/h;               % 采樣頻率%-----------------------------------------------------------------Y1 = Lyapunov_rosenstein_2(x,tau,m,taumax,P,fs);figuresubplot(211); plot(Y1); grid; xlabel('i'); ylabel('y(i)'); title('Lorenz'); axis tight;subplot(212); plot(diff(Y1)); grid; xlabel('n'); ylabel('slope'); axis tight;linear_zone = [50:200]';  % 線性區域F1 = polyfit(linear_zone,Y1(linear_zone),1);Lyapunov_2 = F1(1)                  % 以 2 為底的指數Lyapunov_e = Lyapunov_2/1.4427      % 以 e 為底的指數