实验三用FFT做谱分析.doc

（1） x(n)=1 0=n=5 构造dft 函数计算的10点和20点的 0 其他 dft并画出图形。 functionxk=dft(xn,N);n=0:N-1;k=0:N-1;w=exp(-j.*2.*pi./N);nk=n.*k;wnk=w.nk;xk=xn.*wnk;N1=10;N2=20;n1=0:N1-1;n2=0:N2-1;xn=ones(1,6);xn1=xn zeros(1,4);xn2=xn zeros(1,14);xk1=dft(xn1,N1);xk2=dft(xn2,N2);subplot(2,1,1);stem(n1,xk1);subplot(2,1,2);stem(n2,xk2);（2）a：n=0:7;x=ones(1,4),zeros(1,4);subplot(2,2,1);stem(n,x);xlabel(N=8);ylabel(x(n);y=fft(x);q=n/8;subplot(2,2,2);stem(q,abs(y);ylabel(fft N=8);n=0:15;x=ones(1,4),zeros(1,12);subplot(2,2,3);stem(n,x);xlabel(N=16);ylabel(x(n);q=n/16;y=fft(x)subplot(2,2,4);stem(q,abs(y);ylabel(fft N=16);b，c在一起：N=8;n=0:N-1;q=n/N;x=cos(pi.*n./4);y=fft(x,N);subplot(2,2,1);stem(q,abs(y);grid on;ylabel(fft N=8); N=16;n=0:N-1;q=n/N;x=cos(pi.*n./4);y=fft(x,N);subplot(2,2,2);stem(q,abs(y);grid on;ylabel(fft N=16);N=8;n=0:N-1;q=n/N;x=sin(pi.*n./8);y=fft(x,N);subplot(2,2,3);stem(q,abs(y);grid on;ylabel(fft N=8);N=16;n=0:N-1;q=n/N;x=sin(pi.*n./8);y=fft(x,N);subplot(2,2,4);stem(q,abs(y);grid on;ylabel(fft N=16);（3）：f1=2,f2=2.05;fs=10;N=128;n=0:N-1;q=n/N;x=sin(2.*pi.*n.*f1/fs)+sin(2.*pi.*f2.*n/fs);y=fft(x,N);subplot(3,1,1);stem(q,abs(y);grid on;ylabel(fft N=128);f1=2,f2=2.05;fs=10;N=512;n=0:N-1;q=n/N;x=sin(2.*pi.*n.*f1/fs)+sin(2.*pi.*f2.*n/fs);y=fft(x,N);subplot(3,1,2);stem(q,abs(y);grid on;ylabel(fft N=512);f1=2,f2=2.05;fs=10;N=512;n=0:N-1;q=n/N;n1=0:128;x=sin(2.*pi.*n1.*f1/fs)+sin(2.*pi.*f2.*n1/fs);y=fft(x,N);subplot(3,1,3);stem(q,abs(y);grid on;ylabel(fft N=128-128);