实验三用fft做谱分析

《实验三用fft做谱分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、（1）x(n)={10<=n<=5构造dft函数计算的10点和20点的0其他dft并画出图形。function[xk]=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=[xnzeros(1,4)];xn2=[xnzeros(1,14)];[xk1]=dft(xn1,N1);[xk2]=dft(xn2,

2、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('fftN=8');n=0:15;x=[ones(1,4),zeros(1,12)];subplo

3、t(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('fftN=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));gridon;ylabel('fftN=8');N=16;n=0:N-1;q=n/N;x=cos(pi.*n./4

4、);y=fft(x,N);subplot(2,2,2);stem(q,abs(y));gridon;ylabel('fftN=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));gridon;ylabel('fftN=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));gridon;y

5、label('fftN=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));gridon;ylabel('fftN=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/f

6、s);y=fft(x,N);subplot(3,1,2);stem(q,abs(y));gridon;ylabel('fftN=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));gridon;ylabel('fftN=128-128');