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