尝试对长度幂为2的向量实现快速傅立叶变换

如果操作正确,则应具有运行时间O(n log n)。

function d = ffTU(f)

n = [size(f)](1);

if n==1
  d = f;
  return;

else

  even=ffTU(f(2:2:end,:));
  odd =ffTU(f(1:2:end,:));

  for k=0:(n/2)-1
    T(k+1)= exp(-2i*pi*k/n)*odd(k+1);

  end

  for k=0:(n/2)-1
    d(k+1) = even(k+1) + T(k+1) + even(k+1) - T(k+1);
  end

end

我不断收到以下错误:

error: ffTU: A(I): index out of bounds; value 2 out of bound 1
error: called from
    ffTU at line 16 column 11

我知道我必须调整函数的索引,但我认为我已经解决了这个问题。