实验一快速傅里叶变换.pdf

实验一快速傅里叶变换之报告、实验目的1、在理论学习的基础上通过本实验加深对快速傅立叶变换的理解；2、熟悉并掌握按时间抽取FFT算法的程序；3、了解应用FFT进行信号频谱分析过程中可能出现的问题例如混淆、泄漏、栅栏效应等以便在实际中正确应用FFT。二实验内容a)信号频率F=50Hz采样点数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(AN);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/Th);xlabel('?卩?e/HZ');•u£_£ylabel('??Xejw?'title(•u??D?'上述代码命令中将FFT变换后的数字变量在画图时转换成频域中的频率要是根据数字频率与模拟域频率之间的关系:其中分别为数字和模拟域中的频率且2f于是有:NT运算结果：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)=1b)信号频率F=50Hz采样点数N=32采样间隔T=同理可将a)中F、N、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)