数字信号的 FFT 分析

《数字信号的 FFT 分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、DSP实验报告实验名称：实验二数字信号的FFT分析1．实验要求数字信号处理的一个重要分支就是信号分析，而信号分析的基本工具是离散傅立叶变换。利用傅立叶变换和级数所形成的频谱分析技术作为处理连续信号的重要工具已经应用得很久了,1956年库力（Cooley）和图基（Tukey）所发展的近似频谱的快速算法为频谱分析的数字信号的谱分析铺平了道路。因此，DFT（FFT）得到广泛应用。本次实验设计了两个内容：(1)已知信号这里，N=25，Q=0.9+j0.3。可以推导出，首先根据这个式子计算X(k)的理论值，然后计算输入序列x

2、(n)的32个值，再利用基2时间抽选的FFT算法，计算x(n)的DFTX(k)，与X(k)的理论值比较（要求计算结果最少6位有效数字）。(2)假设信号x(n)由下述信号组成：这个信号有两根主谱线0.3pi和0.302pi靠的非常近，而另一根谱线0.45pi的幅度很小，请选择合适的长度N和窗函数，用DFT分析其频谱，得到清楚的三根谱线。通过本次实验，应该掌握：(a)用傅立叶变换进行信号分析时基本参数的选择。(b)经过离散时间傅立叶变换（DTFT）和有限长度离散傅立叶变换（DFT）后信号频谱上的区别，前者DTFT时间域

3、是离散信号，频率域还是连续的，而DFT在两个域中都是离散的。(c)离散傅立叶变换的基本原理、特性，以及经典的快速算法（基2时间抽选法），体会快速算法的效率。(d)获得一个高密度频谱和高分辨率频谱的概念和方法，建立频率分辨率和时间分辨率的概念，为将来进一步进行时频分析（例如小波）的学习和研究打下基础。(e)建立DFT从整体上可看成是由窄带相邻滤波器组成的滤波器组的概念，此概念的一个典型应用是数字音频压缩中的分析滤波器，例如DVDAC3和MPEGAudio。2.程序分析第一部分：1.25点DFT变换Q=0.9+j*0.

4、3;N=25;n=0:1:N-1;k=0:1:N-1;WN=exp(-j*2*pi/N);X=(1-Q.^N)./(1-Q*WN.^k);X=vpa(X,6)magX=abs(X(1:1:25));%abs为幅值k=0:1:24;w=2*pi*k/25;stem(w/pi,magX);%stem函数作图title('25点DFT');xlabel('2π内幅度特性');得到X(k)的理论值为：X=[1.83992+2.88851*i,10.0793+7.6341*i,0.751948-5.86055*i,0.231

5、478-2.56282*i,0.28018-1.62973*i,0.333453-1.18569*i,0.374736-0.920169*i,0.407012-0.73897*i,0.433275-0.603825*i,0.455542-0.496182*i,0.475155-0.405862*i,0.493047-0.326701*i,0.50991-0.254627*i,0.526307-0.18669*i,0.542748-0.120528*i,0.559746-0.0540151*i,0.577892+0.

6、0150219*i,0.597934+0.0890957*i,0.620917+0.171456*i,0.64843+0.266736*i,0.683099+0.382148*i,0.729676+0.529983*i,0.797838+0.733482*i,0.910733+1.04306*i,1.13973+1.59277*i]2.计算输入序列x(n)的32个值n=0:1:N-1;x=[Q.^n,zeros(1,32-N)];%zeros为零矩阵x=vpa(x,6)序列的后几位数添为0，得到结果为x=[1.0,

7、0.9+0.3*i,0.72+0.54*i,0.486+0.702*i,0.2268+0.7776*i,-0.02916+0.76788*i,-0.256608+0.682344*i,-0.43565+0.537127*i,-0.553224+0.352719*i,-0.603717+0.15148*i,-0.588789-0.0447828*i,-0.516476-0.216941*i,-0.399746-0.35019*i,-0.254714-0.435095*i,-0.0987144-0.467999*i,0

8、.0515569-0.450814*i,0.181645-0.390265*i,0.28056-0.296745*i,0.341528-0.182903*i,0.362246-0.0621539*i,0.344667+0.0527352*i,0.29438+0.150862*i,0.219684+0.22409*i,0.130488+0.267586*i,0