% 假近鄰算法求嵌入維 (False Nearest Neighbors, FNN)% 使用平臺 - Matlab6.5 / Matlab7.0% 作者：陸振波，海軍工程大學% 歡迎同行來信交流與合作，更多文章與程序下載請訪問我的個人主頁% 電子郵件：luzhenbo@sina.com% 個人主頁：luzhenbo.88uu.com.cn% % 參考文獻: M.B.Kennel, R.Brown, H.D.I.Abarbanel. Determining embedding dimension for phase-space reconstruction using a geometrical construction[J]. Phys. Rev. A 1992,45:3403.clcclearclose all%---------------------------------------------------% 產生 Lorenz 時間序列sigma = 16;             % Lorenz 方程參數 a = 16 | 10b = 4;                  %                 b = 4 | 8/3r = 45.92;              %                 c = 45.92 | 28     y = [-1;0;1];           % 起始點 (3x1 的列向量)h = 0.01;               % 積分時間步長k1 = 10000;             % 前面的迭代點數k2 = 5000;              % 后面的迭代點數z = LorenzData(y,h,k1+k2,sigma,r,b);x = z(k1+1:end,1);      % 時間序列(列向量)%-----------------------------------------------------------------% 函數調用tau = 14;           % 時延d_max = 8;          % 最大嵌入維R_tol = 15;         % 判劇1門限A_tol = 2;          % 判劇2門限[Percent] = FNN(x,tau,d_max,R_tol,A_tol);%-----------------------------------------------------------------% 結果作圖plot(1:d_max,Percent,'-o');axis([1,d_max,0,100]);xlabel('Dimension(d)');ylabel('False NN Percentage');