example4_1.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

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');