% Cao方法計算混沌時間吸引子的嵌入維 - 主函數% 使用平臺 - 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,0];            % 起始點 (1 x 3 的行向量)h = 0.01;               % 積分時間步長k1 = 10000;             % 前面的迭代點數k2 = 5000;              % 后面的迭代點數z = LorenzData(y,h,k1+k2,sigma,r,b);    % 用四階 Runge-Kutta 法產生 k1+k2 個點x = z(k1+1:end,1);                      % 去掉前面 k1 個過渡點, 時間序列 - 列向量%-----------------------------------------------------------------% 入口參數tau = 11;               % 時延DimensionMax = 8;       % 最大嵌入維k = 3;                  % 需要計算的近鄰點個數%-----------------------------------------------------------------% 調用cao函數[E1,E2] = cao_luzhenbo(x,tau,DimensionMax,k);%-----------------------------------------------------------------% 結果作圖figure(1); plot(1:DimensionMax,E1,'b+-',1:DimensionMax,E2,'kx:');xlabel('Dimension(d)');ylabel('E1(d) & E2(d)');title('Cao方法求嵌入維');legend('E1','E2',0);