西安交大数字信号处理实验详解.docx

数字信号处理实验——快速傅里叶变换学院：班级：姓名： 学号：快速傅里叶变换一、实验目的在理论学习的基础上，通过本实验加深对快速傅立叶变换的理解；熟悉并掌握按时间抽取FFT算法的程序；了解应用FFT进行信号频谱分析过程中可能出现的问题，例如混淆、泄漏、栅栏效应等，以便在实际中正确应用FFT。二、实验要求记录下实验内容中各种情况下的X (k)值，做出频谱图并深入讨论结果，说明参数的变化对信号频谱产生哪些影响。频谱只做模特性，模的最大值＝1，全部归一化；打印出用C语言（或MATLAB 语言）编写的FFT源程序，并且在每一小段处加上详细的注释说明；用C语言（或MATLAB 语言）编写FFT程序时，要求采用人机界面形式：N , T , f 变量均由键盘输入，补零或不补零要求设置一开关来选择。三、实验内容仔细分析教材第六章‘时间抽取法FFT ’的算法结构，编制出相应的用FFT进行信号分析的C语言（或MATLAB 语言）程序；用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四、程序编写FFT运算：先进行采样，后对采样序列进行码位倒置，再进行蝶形算法。1.采样用for循环来将采样点逐个存入数组x(j)中。for j=0:1:N-1x(j+1)=sin(2*pi*f*j*T); end2.补零采用insert语句对是否补零来进行判断，如果补零，就将补零数加到采样点数，且补零的那些点在数组中都存为0。insert=input(是否补0？[0-否；1-是]:); if insert==1 ZERO=input(补零数：); for j=N:1:N+ZERO-1x(j+1)=0; endN=N+ZERO; end3.码位倒置M为fft运算的级数，在这里被用为表示数组中数据下表的二进制代码位数上限。用for循环来实现对数组中的每个数据进行处理。其中dec2bin（）函数将十进制代码转换为二进制代码，fliplr（）将二进制的矩阵进行左右对称的翻转，即实现了码位倒置，bin2dec（）在讲二进制代码转换为十进制代码。M=log2(N); for t=1:1:N s=dec2bin(t-1,M); s=fliplr(s); s=bin2dec(s); A(s+1)=x(t); end4.fft蝶形运算 对fft运算采用了三个for循环来实现。第一个循环实现fft每一级的运算，所以它的范围应该为1—M,循环一次做一级的fft运算。第二个循环实现每一级运算中的分组，所以其分组与所处的运算级数有关，例如第一级有16的组，第二级有8个组，所以其上限为2^(L-1)-1。第三个循环实现每一个组中的fft运算，但每个组中进行的fft运算次数与分组有关，所以运算时对要进行蝶形运算的采样点的下标要用k和级数来控制，而对于系数W与分组和级数有关。for L=1:1:M %级数for J=0:1:(2^(L-1)-1) %分组for k=(J+1):2^L:N %组内fft运算 T=A(k)+A(k+2^(L-1))*exp((-i*2*pi*J*2^(M-L))/N);A(k+2^(L-1))=A(k)-A(k+2^(L-1))*exp((-i*2*pi*J*2^(M-L))/N);A(k)=T; endendend5.归一化对采样点的幅值去绝对值，再用每一个绝对值除以其中的最大值，这样就可以得到归一化后的图形x=abs(A); y=max(x); X=x/y;6.画图用for循环来实现采样点的幅值与采样点的对应绘图。Axis（）规定图形的纵轴范围（0,1）和横轴范围（0，N）。for j=1:1:Nstem(j-1