用fft对信号作频谱分析的程序.docVIP

下载本文档

12
0
约2.23千字
约 5页
2019-08-05 发布于江西
举报
版权申诉

用fft对信号作频谱分析的程序.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

实验（1） function shiyan1 x1=[1,1,1,1]; n1=0:3; n2=4:7; xa=n1+1; xb=8-n2; x2=[xa,xb]; xa1=4-n1; xb1=n2-3; x3=[xa1,xb1]; X18=fft(x1,8); k1=0:8-1; subplot(3,2,1); stem(k1,abs(X18),.); axis([0,8,0,1.2*max(abs(X18))]); title(8点DFT[x1]);ylabel(幅度); X116=fft(x1,16); k2=0:16-1; subplot(3,2,2); stem(k2,abs(X116),.); axis([0,16,0,1.2*max(abs(X116))]); title(16点DFT[x1]);ylabel(幅度); X28=fft(x2,8); subplot(3,2,3); stem(k1,abs(X28),.); axis([0,8,0,1.2*max(abs(X28))]); title(8点DFT[x2]);ylabel(幅度); X216=fft(x2,16); subplot(3,2,4); stem(k2,abs(X216),.); axis([0,16,0,1.2*max(abs(X216))]); title(16点DFT[x2]);ylabel(幅度); X38=fft(x3,8); subplot(3,2,5); stem(k1,abs(X38),.); axis([0,8,0,1.2*max(abs(X38))]); title(8点DFT[x3]);ylabel(幅度); X316=fft(x3,16); subplot(3,2,6); stem(k2,abs(X316),.); axis([0,16,0,1.2*max(abs(X316))]); title(16点DFT[x3]);ylabel(幅度); 实验（2） function shiyan2 n1=0:8-1; n2=0:16-1; x48=cos(pi/4*n1); x58=cos(pi/4*n1)+cos(pi/8*n1); x416=cos(pi/4*n2); x516=cos(pi/4*n2)+cos(pi/8*n2); X48=fft(x48,8); X416=fft(x416,16); X58=fft(x58,8); X516=fft(x516,16); k1=0:8-1; k2=0:16-1; subplot(2,2,1); stem(k1,abs(X48),.); axis([0,8,0,5]); title(8点DFT[x4]);ylabel(幅度); subplot(2,2,2); stem(k2,abs(X416),.); axis([0,16,0,10]); title(16点DFT[x4]);ylabel(幅度); subplot(2,2,3); stem(k1,abs(X58),.); axis([0,8,0,10]); title(8点DFT[x5]);ylabel(幅度); subplot(2,2,4); stem(k2,abs(X516),.); axis([0,16,0,10]); title(16点DFT[x5]);ylabel(幅度); 实验（3） function shiyan3 Fs=64;T=1/Fs; N1=16;N2=32;N3=64; n1=0:N1-1; n2=0:N2-1; n3=0:N3-1; x816=cos(8*pi*n1*T)+cos(16*pi*n1*T)+cos(20*pi*n1*T); x832=cos(8*pi*n2*T)+cos(16*pi*n2*T)+cos(20*pi*n2*T); x864=cos(8*pi*n3*T)+cos(16*pi*n3*T)+cos(20*pi*n3*T); X816=fft(x816,16); X816=fftshift(X816); X832=fft(x832,32); X832=fftshift(X832); X864=fft(x864,64); X864=fftshift(X864); Tp=N1*T;F=1/Tp; k=-N1/2:N1/2-1;fk=k*F; subplot(3,1,1);stem(fk,abs(X816),.); title(16点DFTx8);xlabel(f(Hz));ylabel(幅度); Tp=N2*T;F=1/Tp; k=-N2/2:N2/2-1;fk=k*F; subplot(3,1,2);stem(fk,abs(X832),.); title(32点DFTX8);xl