需要帮助来运行旧的Matlab代码

Question

我想使用Matlab数组实现FFT（1d）。 谷歌搜索后，我遇到了这段代码 。 但是，这似乎是2007年编写的旧代码，而我正在使用Matlab 2012。

我遇到类似“ 索引超出矩阵尺寸 ”的错误。 那么有人可以更改代码，以便我可以在计算机上运行它吗？

PS：我知道Matlab具有内置的fft例程。 但是我想使用Matlab数组来实现它，因为将来我会将这段代码用于其他目的。 另外，我还需要fft来处理大小为44k的数组，因此基本的DFT实现将无法使用。

Answer 1

代码的问题在第3行：

 p=ceil(l)

该代码假定您的输入信号的长度为2的幂（例如4096个样本或9192）。 如果使用小于2的幂的数字，则此处的p将创建一个比输入大的矩阵。

有2个可能的修复：

1，使用padarray用额外的零填充输入，直到达到下一个2的幂。

2，将ceil改为floor以减少评估数据。 但是，这意味着只对部分输入进行评估。

如果要使用方法1，则需要执行以下操作：

将函数更改为此：

size=length(input);
l=log2(size);
p=ceil(l); 
input(2^p) = 0;
Y=input'; %//NOTE HERE!!!!! the `'` is used if your data is a COLUMN vector, if it is a ROW VECTOR remove this `'`
N = 2^p;
N2=N/2; 
YY = -pi*sqrt(-1)/N2; 
WW = exp(YY);  
JJ = 0 : N2-1;  
 W=WW.^JJ;
for L = 1 : p-1
   u=Y(:,1:N2);
   v=Y(:,N2+1:N);
   t=u+v;
   S=W.*(u-v);
   Y=[t ; S];
   U=W(:,1:2:N2);
   W=[U ;U];
   N=N2;
   N2=N2/2;
end;
u=Y(:,1);
v=Y(:,2);
Y=[u+v;u-v];
Y