How to do a Amplitude Adjusted Fourier Transform on time series?
10 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
I need to use MATLAB (R2022b) to surrogate my data (1x500 excel) using Amplitude Adjusted Fourier Transform.
Any advice how to do that?
0 Kommentare
Antworten (1)
Jaimin
am 3 Jan. 2025
Hi @Mir Sahand
To surrogate your data using the Amplitude Adjusted Fourier Transform (AAFT) in MATLAB, you will need to perform the AAFT algorithm, which involves randomizing the phases of the Fourier transform of your data and then adjusting the amplitude.
For understanding kindly refer following code.
% Step 1: Import Data
filename = 'your_data.xlsx'; % Replace with your file name
data = readtable(filename);
% Ensure data is a row vector
if size(data, 1) > 1
data = data';
end
% Step 2: AAFT Algorithm
% Sort the original data
sortedData = sort(data);
% Randomize the original data
randomizedData = data(randperm(length(data)));
randomizedFFT = fft(randomizedData);
% Randomize the phases
randomizedPhases = angle(randomizedFFT);
amplitudes = abs(randomizedFFT);
randomizedFFT = amplitudes .* exp(1i * randomizedPhases);
% Inverse Fourier Transform
surrogateData = ifft(randomizedFFT, 'symmetric');
% 2.6: Adjust the amplitude by sorting
[~, sortIndex] = sort(surrogateData);
surrogateData(sortIndex) = sortedData;
For more information kindly refer following MathWorks documentation.
1 Kommentar
Paul
am 4 Jan. 2025
Do you have source for this code/algorithm? I'm curious about a few things.
rng(100);
% Step 1: Import Data
data = 1:10;
data = data(randperm(10));
% Step 2: AAFT Algorithm
% Sort the original data
sortedData = sort(data);
% Randomize the original data
randomizedData = data(randperm(length(data)));
randomizedFFT = fft(randomizedData);
temp = randomizedFFT;
Nothing is being randomized here. All this code does is recreate the randomizedFFT (to withing numerical precision).
% Randomize the phases
randomizedPhases = angle(randomizedFFT);
amplitudes = abs(randomizedFFT);
randomizedFFT = amplitudes .* exp(1i * randomizedPhases);
norm(temp-randomizedFFT)
% Inverse Fourier Transform
surrogateData = ifft(randomizedFFT, 'symmetric');
% 2.6: Adjust the amplitude by sorting
[~, sortIndex] = sort(surrogateData);
surrogateData(sortIndex) = sortedData;
Consequently, surrogateData is identical to randomizedData.
randomizedData
surrogateData
It's unclear as to why any more code is needed after randomizedData is computed.
Siehe auch
Kategorien
Mehr zu Discrete Fourier and Cosine Transforms finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!