example4_1.m

来自「数字信号处理的比较全面的程序代码」· M 代码 · 共 13 行

clear
Nmax=2048*4;
fft_time=zeros(1,Nmax);
for n=1:1:Nmax
    x=rand(1,n);
    t=clock;
    fft(x);
    fft_time(n)=etime(clock,t);
end
n=[1:1:Nmax];
plot(n,fft_time,'.');
xlabel('N');ylabel('time(ms)');
title('FFT runtime');