實(shí)驗(yàn)一快速傅里葉變換

實(shí)驗(yàn)一快速傅里葉變換之報(bào)告一、實(shí)驗(yàn)?zāi)康?、在理論學(xué)習(xí)的基礎(chǔ)上，通過(guò)本實(shí)驗(yàn)加深對(duì)快速傅立葉變換的理解；2、熟悉并掌握按時(shí)間抽取FFT算法的程序；3、了解應(yīng)用FFT進(jìn)行信號(hào)頻譜分析過(guò)程中可能出現(xiàn)的問(wèn)題，例如混淆、泄漏、柵欄效應(yīng)等，以便在實(shí)際中正確應(yīng)用FFT。二實(shí)驗(yàn)內(nèi)容信號(hào)頻率F＝50Hz，采樣點(diǎn)數(shù)N=32,采樣間隔T=matlab程序代碼為：F=50;T=;N=32;n=0:N-1;t=n*T;A=sin(2*pi*F*t);figure;Y=fft(A,N);h=(abs(Y));h=h/max(h(1:N));forn=1:N;string1=strcat('X(',num2str(n-1),')=',num2str(h(n)));disp(string1);f=(n/T)/N;endstem([0:N-1]/N/T,h);xlabel('?μ?ê/HZ');ylabel('??·ùX￡¨ejw￡?');title('·ù?μì?D?');上述代碼命令中，將FFT變換后的數(shù)字變量K，在畫圖時(shí)轉(zhuǎn)換成頻域中的頻率f。這主要是根據(jù)數(shù)字頻率與模擬域頻率之間的關(guān)系：其中、分別為數(shù)字和模擬域中的頻率，且于是有：運(yùn)算結(jié)果：X(1)=1X(2)=X(3)=X(4)=X(5)=X(6)=X(7)=X(8)=X(9)=X(10)=X(11)=X(12)=X(13)=X(14)=X(15)=X(16)=X(17)=X(18)=X(19)=X(20)=X(21)=X(22)=X(23)=X(24)=X(25)=X(26)=X(27)=X(28)=X(29)=X(30)=X(31)=1信號(hào)頻率F＝50Hz，采樣點(diǎn)數(shù)N=32,采樣間隔T=同理可將a）中F、N、T，參數(shù)改成要求值（以下均是如此），即可得，X(0)=X(1)=X(2)=X(3)=X(4)=X(5)=X(6)=X(7)=X(8)=1X(9)=X(10)=X(11)=X(12)=X(13)=X(14)=X(15)=X(16)=X(17)=X(18)=X(19)=X(20)=X(21)=X(22)=X(23)=X(24)=1X(25)=X(26)=X(27)=X(28)=X(29)=X(30)=X(31)=信號(hào)頻率F＝50Hz，采樣點(diǎn)數(shù)N=32,采樣間隔T=X(0)=X(1)=X(2)=X(3)=X(4)=X(5)=X(6)=X(7)=1X(8)=X(9)=X(10)=X(11)=X(12)=X(13)=X(14)=X(15)=X(16)=X(17)=X(18)=X(19)=X(20)=X(21)=X(22)=X(23)=X(24)=X(25)=1X(26)=X(27)=X(28)=X(29)=X(30)=X(31)=此時(shí)采樣總時(shí)間t=*32=，信號(hào)周期T1=1/50=，則t/T1=，不為整數(shù)，即采樣與原信號(hào)沒(méi)有相同的基頻，因此事不能運(yùn)用FFT的，其結(jié)果也必將是錯(cuò)誤的。信號(hào)頻率F＝50Hz，采樣點(diǎn)數(shù)N=32,采樣間隔T=X(0)=X(1)=X(2)=X(3)=X(4)=X(5)=X(6)=1X(7)=X(8)=X(9)=X(10)=X(11)=X(12)=X(13)=X(14)=X(15)=X(16)=X(17)=X(18)=X(19)=X(20)=X(21)=X(22)=X(23)=X(24)=X(25)=X(26)=1X(27)=X(28)=X(29)=X(30)=X(31)=同c）一樣，此時(shí)采樣總時(shí)間t=*32=，信號(hào)周期T1=1/50=，則t/T1=，不為整數(shù)，即采樣與原信號(hào)沒(méi)有相同的基頻，因此事不能運(yùn)用FFT的，其結(jié)果也必將是錯(cuò)誤的。信號(hào)頻率F＝50Hz，采樣點(diǎn)數(shù)N=64,采樣間隔T=X(0)=X(1)=X(2)=1X(3)=X(4)=X(5)=X(6)=X(7)=X(8)=X(9)=X(10)=X(11)=X(12)=X(13)=X(14)=X(15)=X(16)=X(17)=X(18)=X(19)=X(20)=X(21)=X(22)=X(23)=X(24)=X(25)=X(26)=X(27)=X(28)=X(29)=X(30)=X(31)=X(32)=X(33)=X(34)=X(35)=X(36)=X(37)=X(38)=X(39)=X(40)=X(41)=X(42)=X(43)=X(44)=X(45)=X(46)=X(47)=X(48)=X(49)=X(50)=X(51)=X(52)=X(53)=X(54)=X(55)=X(56)=X(57)=X(58)=X(59)=X(60)=X(61)=X(62)=1X(63)=滿足奈奎斯采樣定律，結(jié)果正確。信號(hào)頻率F＝250Hz，采樣點(diǎn)數(shù)N=32,采樣間隔T=X(0)=X(1)=X(2)=X(3)=X(4)=X(5)=X(6)=X(7)=X(8)=1X(9)=X(10)=X(11)=X(12)=X(13)=X(14)=X(15)=X(16)=X(17)=X(18)=X(19)=X(20)=X(21)=X(22)=X(23)=X(24)=1X(25)=X(26)=X(27)=X(28)=X(29)=X(30)=X(31)=顯然，由于采樣頻率小于2倍的最高頻率，不滿足奈奎斯采樣定理，發(fā)生頻譜混疊，因此其FFT不能正確的反應(yīng)其真實(shí)的頻譜。將c）信號(hào)后補(bǔ)32個(gè)0，做64點(diǎn)FFTMatlab程序代碼：T=;F1=50;N=64;t=[0:32-1]*T;t2=[tzeros(1,32)];x=sin(F1*2*pi*t2);figure;Y=fft(A,N);h=(abs(Y));h=h/max(h(1:N));forn=1:N;string1=strcat('X(',num2str(n-1),')=',num2str(h(n)));disp(string1);f=(n/T)/N;endstem([0:N-1]/N/T,h);xlabel('?μ?ê/HZ');ylabel('??·ùX￡¨ejw￡?');title('·ù?μì?D?');運(yùn)行結(jié)果：X(0)=X(1)=X(2)=X(3)=X(4)=X(5)=X(6)=X(7)=X(8)=X(9)=X(10)=X(11)=X(12)=X(13)=X(14)=X(15)=X(16)=1X(17)=X(18)=X(19)=X(20)=X(21)=X(22)=X(23)=X(24)=X(25)=X(26)=X(27)=X(28)=X(29)=X(30)=X(31)=X(32)=X(33)=X(34)=X(35)=X(36)=X(37)=X(38)=X(39)=X(40)=X(41)=X(42)=X(43)=X(44)=X(45)=X(46)=X(47)=X(48)=1X(49)=X(50)=X(51)=X(52)=X(53)=X(54)=X(55)=X(56)=X(57)=X(58)=X(59)=X(60)=X(61)=X(62)=X(63)=由于在時(shí)間域內(nèi)信號(hào)加零，致使振幅譜中出現(xiàn)很多其他