Цикл преобразования Фурье.

[Dusty]

Добрый день. Стоит задача написать вручную быстрое преобразование фурье. Зашел в скрипт функции fft и попытался сделать в цикле эту формулу дабы преобразовывала изображение. К сожалению либо нехватает цикла либо я запутался в индексации. Прошу выручайте)
Мой код

Код: Выделить всё

close all
f = imread('lena.bmp');
imshow(f);
size(f);
[M,N] = size(f);
 
for i = 1:length(M)
 
        for j = 1:length(N)
       
        alpha(i,j) = f(i,j)*exp((-1i*2*pi*(i-1)*(j-1))/length(M));
   
        end
   
FFT(i,j) = sum(alpha);
   
end
 
 
y = fft(f);
figure, imshow(FFT);
figure, imshow(y);

[Sniper53]

Мой код

Плагиат