基2FFT算法比较.docVIP

下载本文档

9
0
约 4页
2017-06-08 发布于重庆
举报
版权申诉

基2FFT算法比较.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文档。上传文档

基2FFT算法比较

脑与认知科学小班讨论61.推导N 16点的基2--FFT表达式，并按时间抽取流程图。 ‘ 2.编写N点的基2--按时间抽取的FFT算法（fft_m.m），并与系统自带的FFT函数、采用dft直接算法比较，画出不同点数所需要的时间对比图 main.m 。要求用matlab语言编写，提交到课程中心。 1.系统自带的fft函数 clear,close all, clc Fs 8000; T 1/ Fs; N 4096; t 0:N-1 *T; y 1*sin 2*pi*500*t +2*sin 2*pi*3000*t ; Y fft y ;plot abs Y ; v*dt*cosx; 2.基2function y fft_m x %disp [xr;xi] ; disp \n Np length x ; Mp ceil log2 Np ; if power 2,Mp Np x [x zeros 1, 2^Mp-Np ]; %error DòáDxμ?3¤?èó|??ê?2μ???êy′??Y ; end N length x ; M log2 N ; xr real x ; xi imag x ; %Np N/2+1; % calculate rotate factor for i 1:Np rcf i cos i-1 *2*pi/Np ; icf i sin i-1 *2*pi/Np ; end % for original location calculation % change order [xr,xi] daoxu xr,xi ; % xi daoxu xi ; % disp [xr;xi] ; disp \n for L 1:M B power 2,L-1 ; P1 power 2,M-L ; k1 power 2,L ; for J 0:B-1 P J*P1; for k J+1:k1:N rPartKB xr k+B *rcf P+1 +xi k+B *icf P+1 ; iPartKB xi k+B *rcf P+1 -xr k+B *icf P+1 ; xr k+B xr k -rPartKB; xi k+B xi k -iPartKB; xr k xr k +rPartKB; xi k xi k +iPartKB; end %fprintf %10.5g,xr end % if L 2 % disp [xr;xi] ; disp \n % % end end y xr+xi*j; y y 1:Np ; 3.对模拟信号y t [ sin 2*pi*1000t + sin 2*pi*3000t +随机噪声，采样率16kHz]和数字信号z t [ 1 1 1 1 0 0 1 1 1 0 0 1 0 1 0 1 0,比特率为1kbps]进行频谱分析。 clear,close all, clc Fs 16000; T 1/ Fs; N 4096; t 0:N-1 *T; y 1*sin 2*pi*500*t +2*sin 2*pi*3000*t ; Y fft_m y ; Z fft y ; figure plot abs Y ;