离散傅里叶变换.doc

实验一离散傅立叶变换一、实验目的 1. 在理论学习的基础上,通过本实验,加深对离散傅立叶变换( DFT )的理解。 2. 以正弦信号为例,学习和掌握利用离散傅立叶变换( DFT )分析信号频谱的方法。 3. 了解应用 DFT 进行信号频谱分析过程中可能出现的问题,以便在实际中正确应用 DFT 。 4. 学习和熟悉在计算机上作图显示和研究信号的离散频谱二、实验原理与方法在各种信号序列中,有限长序列的数字信号处理占有很重要地位,对有限长序列,我们可以使用离散傅立叶变换( DFT ), 这一变换不但可以很好地反映序列的频谱特性, 而且易于用快速算法在计算机上实现,当序列 x(n) 的长度为 N 时,它的 DFT 定义为: X(k )= 反变换为: x(n )= 有限长序列的 DFT 是其 Z 变换在单位圆上的等距采样,或者说是序列傅立叶变换的等距采样, 因此可以用于序列的谱分析。 DFT 要求被分析序列的长度为有限长,对无限长和很长的序列则需要截取。截取的结果使得信号的频谱分析出现频率泄漏。本实验分析的周期性信号,对周期信号,当截取的序列长度是信号周期的整数倍的时候, DFT 分析的频率离散点恰好采在周期信号的谐波分量和其它零值频率点上, 能正确反映信号的频谱。当截取的序列长度不为信号周期的整数倍时, DFT 分析的结果就会有频率泄漏的反映。三、实验内容 1. 设计说明编制程序用 DFT 分析正弦信号的频谱。具体包括以下几部分(1 )编程产生正弦抽样输入信号 x(t) = sin( Ω t)= sin(2 πf t) x(n) = x(nT) = sin(2 πfn T) 正弦信号的频率 f 、取样点数 N 和取样间隔 T 通过人机对话的方式由键盘输入,以便调整。(2 )编程实现 DFT 的运算 DFT 的计算公式为(3 )计算 DFT 的幅度(4 )作图画出 DFT 的幅度谱线。 2 .实验步骤将编制的程序输入计算机,调试运行正确后,进行以下实验。(1 )分别在下列几种情况下,输入正弦信号的频率 f 、采样点数 N 和采样周期 T ,观察和研究 DFT 的结果,并记录谱线的大体形状。(a)f= 50 Hz,N= 32,T= s (b)f= 50 Hz,N= 32,T= s (c)f= 50 Hz,N= 32,T= s (d)f= 150 Hz,N= 32,T= s (2 )修改程序,将输入信号改变成 50 Hz和 150 Hz 两个正弦信号分量的迭加,并令 50 Hz分量的幅度为 2, 150 Hz 分量的幅度为 1 。即 x(n) =2 sin(2 π 50n T)+ sin(2 π 150 n T) 取N= 32 ,分别在 T= s、T= s和T= s 情况下进行 DFT 分析,观察和记录结果。(3 )有余力的同学还可以将输入信号改成 50Hz 的周期性方波,再对其进行 DFT 分析。四、 C 语言参考程序/*用 DFT 分析正弦信号频谱的 C 语言程序*/ #include<> #include<> #include<> #include<> #define PI f