FFT实验报告.docxVIP

快速傅里叶变换实验 实验名称 快速傅里叶变换一．实验目的在理论学习的基础上，通过本实验加深对快速傅立叶变换的理解；熟悉并掌握按时间抽取FFT算法的程序；了解应用FFT进行信号频谱分析过程中可能出现的问题，例如混淆、泄漏、栅栏效应等，以便在实际中正确应用FFT。二.实验要求记录下实验内容中各种情况下的X (k)值，做出频谱图并深入讨论结果，说明参数的变化对信号频谱产生哪些影响。频谱只做模特性，模的最大值＝1，全部归一化；打印出用C语言（或MATLAB 语言）编写的FFT源程序，并且在每一小段处加上详细的注释说明；用C语言（或MATLAB 语言）编写FFT程序时，要求采用人机界面形式：N , T , f 变量均由键盘输入，补零或不补零要求设置一开关来选择。三.实验内容1.仔细分析教材第六章‘时间抽取法FFT ’的算法结构，编制出相应的用FFT进行信号分析的C语言（或MATLAB 语言）程序；2.用FFT程序计算有限长度正弦信号 分别在以下情况下所得的DFT结果并进行分析和讨论：信号频率f＝50Hz，采样点数N=32,采样间隔T=0.000625s信号频率f＝50Hz，采样点数N=32,采样间隔T=0.005s信号频率f＝50Hz，采样点数N=32,采样间隔T=0.0046875s信号频率f＝50Hz，采样点数N=32,采样间隔T=0.004s信号频率f＝50Hz，采样点数N=64,采样间隔T=0.000625s信号频率f＝250Hz，采样点数N=32,采样间隔T=0.005s 将c） 信号后补32个0，做64点FFT（1）实验结果：信号的频率为50Hz，采样周期为0.000625s，则采样频率为1600Hz，在数字域，采样频率对应的是2*pi，由于sin信号的频谱为在pi/16和-pi/16上的冲激，所以时域采样后，在频域以2*pi为周期进行周期延拖。加窗后，卷积一个sinc函数，由于窗的长度为32，所以所以其主瓣的宽度为pi/16，因为进行DFT要有32个点，所以第k个点对应的频谱为（k-1）*pi/16,由上面的sinc函数的主瓣宽度知，除了pi/16和31*pi/16（由于频域的周期延拖），其他的点处，sinc函数都为0，DFT的结果只有在第1个点和第31个点有大小相等的值，其它点为0.变化的是采样间隔变为0.005s，这样采样频率就变为200Hz，同样采样频率对应数字域的2*pi，所以sin信号的谱在数字域对应的频率为+0.5*pi和-0.5*pi,周期延拖后在1.5*pi处也会有，故在0到2*pi中有两个冲激，窗的宽度还是32，所以sinc函数的主瓣宽度和上题一样，因为进行DFT要有32个点，所以第k个点对应的频谱为（k-1）*pi/16,由上面的sinc函数的主瓣宽度知，除了pi/2和3*pi/2（由于频域的周期延拖），其他的点处，sinc函数都为0，所以DFT的结果只有在第8个点和第24个点有大小相等的值，其它点为0.采样周期为0.0046875s，所以采样频率为213.33Hz，不是整数，因为在数字域采样频率对应的是2*pi，所以可以知道sin函数在数字域里对应的是0.234*pi，和-0.234*pi，可见不是pi/16的整数倍，所以其卷积sinc函数后，由于sinc函数的主瓣仍然保持为pi/16，所以频率在k*pi/16上不为0而是两个sinc函数的叠加，由于sinc函数从主瓣向两边延伸时，幅度渐渐衰减，所以最终频谱为在7*pi/16,8*pi/16和24*pi/16,25*pi/16这4个距离sinc函数主瓣最近的频率处会有比较大的幅值，其他点依次衰减。 本题与c）题的区别在于最终做了一个64点的DFT，其它均没有变化，就是说，在对频谱采样前，本题的频谱和c）的是相同的。在c中采样的间隔是pi/16,而在本题中采样的间隔是pi/32。所以在2*n点上的值与c）对应的n点的值相同；在2*n+1点上除了在第15和第49点上有峰值外，其余点为0。（2）FFT的matlab代码：function y=myfft(x) %定义函数N=length(x); %输入数据组的长度n=log2(N);for m=1:N; %变址运算 k=dec2bin(m-1,n); k=fliplr(k); k=bin2dec(k); X(k+1)=x(m);endW=exp(-j*2*pi/N); %Wfor i=1:n %i层运算 for m=1:N/(2^i)