diesenbacher.net

Completely needless and trivial – my own implementation of fft (fast fourier transformation) in matlab. Yes, I know, there ist already fft, fft2, and fftn – but as a proof of concept:

function myfft(vec)
%MYFFT My own simple version of FFT 
  vec * fftmatrix(length(vec))
end

function thematrix = fftmatrix(n)
% FFTMATRIX Generates a fft-matrix of the size n x n
  thematrix = zeros(n,n);
  for izeile = 1:n
    for ispalte = 1:n
      thematrix(izeile,ispalte) = exp(i*(2*pi/n))^((izeile-1)*(ispalte-1));
    end
  end
end