compare.m

自制的快速傅立叶变换函数

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


for j=1:10
    tic;
    for (k=1:1000)
        fft(ones(1,L(j)));
    end;
    T2(j)=toc/1000;
end
T2

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