數(shù)字信號(hào)處理 實(shí)驗(yàn)報(bào)告四

實(shí)驗(yàn)四用FFT 做譜分析一、實(shí)驗(yàn)?zāi)康模?1、加深對(duì)DFT算法原理及基本性質(zhì)的理解。 2、熟悉FFT算法原理及子程序的應(yīng)用。3、學(xué)習(xí)FFT對(duì)連續(xù)信號(hào)和時(shí)域離散信號(hào)作譜分析的方法。了解可能出現(xiàn)的分析誤差和原因，以便在實(shí)際中正確應(yīng)用FFT。二、實(shí)驗(yàn)數(shù)據(jù)：>>x=[2,1,3,2,1,5,1];>>h=[1,2,-1,-3];>>N=length(x)+length(h)-1;>>n=0:N-1;>>x=[x,zeros(1,N-length(x))];>>h=[h,zeros(1,N-length(h))];>>X=fft(x);H=fft(h);>>Y=X.*H;y=ifft(Y);>>stem(n,x,'.');xlabel('n');ylabel('x(n)');M文件：mydft.m：functiony=mydft(x);N=length(x);n=0:N-1;k=n;nk=n'*k;WN=exp(-j*2*pi/N);WNnk=WN.^nk;Xk=x*WNnk;y=Xk;M文件：mydifft.m：functiony=mydifft(x)m=nextpow2(x);N=2^m;iflength(x)<N x=[x,zeros(1,N-length(x))];endnxd=bin2dec(fliplr(dec2bin([1:N]-1,m)))+1;y=x(nxd);formm=1:m Nmr=2^mm; u=1; WN=exp(i*2*pi/Nmr) forj=1:Nmr/2 fork=j:Nmr:N kp=k+Nmr/2; t=y(kp)*u; y(kp)=y(k)-t; y(k)=y(k)+t; end u=u*WN; endendM文件：Sy4_1.mn=[0:1:1023];t=n*0.001;x=20*sin(60*pi*t)+5*cos(120*pi*t);y1=fft(x);plot(y1);y2=mydft(x);plot(y2);y3=mydifft(x);plot(y3);>>Sy4_1記錄三種方法所耗費(fèi)的時(shí)間：1.k=input('k=');k=2；x=randn(1,2^k);tic,X=fft(x),toc;X=Columns1through3-1.6851-0.5579+1.9533i1.0707Column4-0.5579-1.9533iElapsedtimeis0.137266seconds.2.tic,X=mydft(x),toc;X=Columns1through3-1.6851-0.5579+1.9533i1.0707+0.0000iColumn4-0.5579-1.9533iElapsedtimeis1.354649seconds.3.tic,X=mydifft(x),toc;WN=-1.0000+0.0000iWN=0.0000+1.0000iX=Columns1through3-1.6851-0.5579-1.9533i1.0707Column4-0.5579+1.9533iElapsedtimeis1.0