compare my-fft.m

来自「自制的快速傅立叶变换函数」· M 代码 · 共 20 行

L=[2^1,2^2,2^3,2^4,2^5,2^6,2^7,2^8]
for i=1:8
    tic;
    dft(ones(1,L(i)));
    T1(i)=toc;
end
T1

for j=1:8
    tic;
    for (k=1:1000)
    my_fft(ones(1,L(j)));
    end;
    T3(j)=toc/1000;
end
T3

loglog(L,T1)
hold on,loglog(L,T3)