matlab 混沌函数工具箱

这是常用的混沌函数工具箱，有现成的m文件使用，欢迎大家下载

function s=AAFT（xc）
%Syntax: s=AAFT（xc）
%___________________
%
% Makes c Amplitude Adjusted Fourier Transformed （AAFT） surrogates of a time
% series x.
%
% s is the AAFT time series.
% x is the original time series.
% c is the number of surrogates.
%
%
% References:
%
% Theiler J Galdrikian B Longtin A Eubank S Farmer D J （1992）: Using 
% Surrogate Data to Detect Nonlinearity in Time Series. In Nonlinear Modeling
% and Forecasting eds. Casdagli M & Eubank S. 163-188. Addison-Wesley
%
% Theiler J Eubank SGaldrikian B Longtin A  Farmer D J （1992）: Testing
% for nonlinearity in time series: the method of surrogate data. Physica D
% 58: 77-94
%
% 
% Alexandros Leontitsis
% Department of Education
% University of Ioannina
% 45110 - Dourouti
% Ioannina
% Greece
%
% University e-mail: me00743@cc.uoi.gr
% Lifetime e-mail: leoaleq@yahoo.com
% Homepage: http://www.geocities.com/CapeCanaveral/Lab/1421
%
% 12 Apr 2001.

if nargin<1 | isempty（x）==1
   error（‘You should provide a time series.‘）;
else
   % x must be a vector
   if min（size（x））>1
      error（‘Invalid time series.‘）;
   end
   x=x（:）;
end

if nargin<2 | isempty（c）==1
   c=1;
else
   % c must be scalar
   if sum（size（c））>2
      error（‘c must be scalar.‘）;
   end
   % c must be greater or equal than 1
   if c<1
      error（‘c must be greater or equal than 1.‘）;
   end
end

for i=1:c
    % Initialize
    y=x;
    % Make n normal random devaiates
    normal=sort（randn（size（y）））;
    % Sort y and extract the ranks
    [yT]=sort（y）;
    [Tr]=sort（T）;
    % Assign the ranks of y to the normal deviates and apply the phase
    %  randomization
    normal=phaseran（normal（r））;
    % Extract the ranks of the phase randomized normal deviates
    [normalT]=sort（normal）;
    [Tr]=sort（T）;
    % Assign the ranks of the phase randomized normal deviates to y and
    %  obtain the AAFT surrogates
    s（:i）=y（r）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1922  2003-10-29 21:29  logistic.m
     文件        1974  2003-10-29 21:28  AAFT.m
     文件        2216  2003-10-27 00:38  Contents.m
     文件        1687  2003-10-29 21:29  derivs.m
     文件        5487  2003-10-29 21:23  gencorint.m
     文件        2409  2003-10-29 21:29  henon.m
     文件        2706  2003-10-29 21:24  IAAFT.m
     文件        2434  2003-10-29 21:29  ikeda.m
     文件         787  2004-06-27 22:03  Knearest.m
     文件        1366  2002-11-20 02:31  vr.m
     文件        3299  2003-10-29 21:30  lorentz.m
     文件        2663  2003-10-29 21:30  mackeyglass.m
     文件        2869  2003-10-29 21:38  noiseest.m
     文件        6307  2004-07-06 15:48  noisergeo.m
     文件        4360  2004-06-30 22:21  noiserSchreiber.m
     文件        2029  2003-10-29 21:31  phaseran.m
     文件        2013  2003-10-29 21:27  phasespace.m
     文件        1941  2003-10-29 21:31  quadratic.m
     文件         961  2004-06-30 23:50  radnearest.m
     文件        3246  2003-10-29 21:31  rossler.m
     文件        1606  2003-10-29 21:31  shuffle.m
     文件        1111  2003-10-29 21:32  SSA.m
     文件        2523  2003-11-08 12:24  SSAeye.m
     文件        2297  2003-11-03 19:15  SSAforeiter.m
     文件        2012  2003-10-29 21:32  SSAinv.m