=1:1:100;Р t=clock;Р a1=sin(n)+cos(n);Р fft(a1);Р fft_time(n)=etime(clock,t);Рend;Рn=1:100;Рsubplot(1,2,2);Рplot(n,fft_time,'.');Рxlabel('N');ylabel('time/s');Рtitle('FFT');Р*由图可以看出FFT计算比DFT计算更省时间。Р1.2 对离散确定信号作如下谱分析:Р(1)截取使成为有限长序列N(),(长度N自己选)写程序计算出的N点DFT ,画出时域序列图xn~n和相应的幅频图Р。Р假设N取32, Р即 0 ≤n ≤ 31 时, 编写程序,计算出X(n)的32点DFT X(k)。Р源程序如下:РN=32;Рn=0:1:N-1;Рxn=cos(0.48*pi*n)+cos(0.52*pi*n);Рsubplot(3,1,1)Рstem(n,xn,'.k');Рtitle('时域序列图xn');Рxlabel('n');Рaxis([0,10,-2.5,2.5]);Рw=2*pi*(0:1:2047)/2048;РXw=xn*exp(-j*n'*w);Рsubplot(3,1,2);Рplot(w/pi,abs(Xw));Рtitle('幅频特性曲线X(ejw)');Рxlabel('w');Рaxis([0,1,0,20]);Рsubplot(3,1,3)Рk=0:1:N-1;РWN=exp(-j*2*pi/N);Рnk=n'*k;РWNnk=WN.^nk;РXK=xn*WNnk;Рk1=0:1:31;Рw1=2*pi/32*k1;Рstem(w1/pi,abs(XK(1:1:32)),'.k');Рtitle('频域序列图XK');Рxlabel('频率(单位:pi)');Рaxis([0,1,0,20]);
全文预览
