表2－1 几种典型波形的傅立叶变换表 名称 波形函数()f t 波形图 频谱函数()F ω 频谱图 矩形脉冲 ,2 0, 2 E t t τ τ ?
名称 波形函数()f t 波形图 频谱函数()F ω 频谱图 梯形脉冲 1 111 110, 22,222, 222, 222t E t t E t E t t τττττττττττττ? ≥?? ???+-<<- ??-?? ?? ?-<?
附录A 拉普拉斯变换及反变换 1.拉氏变换的基本性质 附表A-1 拉氏变换的基本性质 1()([n n k f t dt s s -+= +∑?个
2．常用函数的拉氏变换和z变换表 附表A-2 常用函数的拉氏变换和z变换表
3． 用查表法进行拉氏反变换 用查表法进行拉氏反变换的关键在于将变换式进行部分分式展开，然后逐项查表进行反变换。设)(s F 是s 的有理真分式，即 11 10 111)()()(a s a s a s a b s b s b s b s A s B s F n n n n m m m m ++++++++==---- （m n >） 式中，系数n n a a a a ,,...,,110-和011,, ,,m m b b b b -都是实常数；n m ,是正整数。按代数定理 可将)(s F 展开为部分分式。分以下两种情况讨论。 （1）0)(=s A 无重根：这时，F(s)可展开为n 个简单的部分分式之和的形式，即 ∑=-=-++-++-+-=n i i i n n i i s s c s s c s s c s s c s s c s F 122 11)( （F-1） 式中，n s s s ,,,21 是特征方程A(s)＝0的根；i c 为待定常数，称为()F s 在i s 处的留数，可按下列两式计算：lim()()i i i s s c s s F s →=- （F-2） 或 i s s i s A s B c ='= )() ( （F-3）
式中，)(s A '为)(s A 对s 的一阶导数。根据拉氏变换的性质，从式（F-1）可求得原函数为 []??????-==∑=--n i i i s s c L s F L t f 11 1 )()(＝1 i n s t i i c e =∑ (F-4) （2）0)(=s A 有重根：设0)(=s A 有r 重根1s ，F(s)可写为 ()) ()()() (11n r r s s s s s s s B s F ---= + = n n i i r r r r r r s s c s s c s s c s s c s s c s s c -++-++-+-++-+-++-- 11 111 111)()()( 式中，1s 为F(s)的r 重根，1+r s ，…，n s 为F(s)的n r -个单根；其中，1+r c ，…，n c 仍按式(F-2)或式(F-3)计算，r c ，1-r c ，…，1c 则按下式计算： )()(lim 11 s F s s c r s s r -=→ 11lim [()()]i r r s s d c s s F s ds -→=- )()(lim !11)() (1s F s s ds d j c r j j s s j r -=→- (F-5) )()(lim )!1(11)1() 1(11s F s s ds d r c r r r s s --=--→ 原函数)(t f 为 [])()(1 s F L t f -= ??????-++-++-+-++-+-=++---n n i i r r r r r r s s c s s c s s c s s c s s c s s c L 11 111 1111)()() ( t s n r i i t s r r r r i e c e c t c t r c t r c ∑+=---+?? ????+++-+-=112211 1 )!2()!1( （F-6）
实验八 利用FFT 实现快速卷积 一、 实验目的 （1） 通过这一实验，加深理解FFT 在实现数字滤波（或快速卷积）中的重要作用，更好的利用FFT 进行数字信号处理。 （2） 进一步掌握循环卷积和线性卷积两者之间的关系。 二、 实验原理与方法 数字滤波器根据系统的单位脉冲响应h(n)是有限长还是无限长可分为有限长单位脉冲响应（Finite Impulse Response ）系统（简记为FIR 系统）和无限长单位脉冲响应（Infinite Impulse Response ）系统（简记为IIR 系统）。 对于FIR 滤波器来说，除了可以通过数字网络来实现外，也可以通过FFT 的变换来实现。 一个信号序列x(n)通过FIR 滤波器时，其输出应该是x(n)与h(n)的卷积： ∑+∞ -∞ =-= =m m n h m x n h n x n y )()()(*)()( 或 ∑+∞ -∞ =-= =m m n x m h n x n h n y ) ()()(*)()( 当h(n)是一个有限长序列，即h(n)是FIR 滤波器，且10-≤≤N n 时 ∑-=-=1 0) ()()(N m m n x m h n y 在数字网络（见图6.1）类的FIR 滤波器中，普遍使用的横截型结构（见下图6.2 图6.1 滤波器的数字网络实现方法 图6.2 FIR 滤波器横截型结构 y(n) y(n) -1-1-1-1
应用FFT 实现数字滤波器实际上就是用FFT 来快速计算有限长度列间的线性卷积。 粗略地说，这种方法就是先将输入信号x(n)通过FFT 变换为它的频谱采样 值X(k)，然后再和FIR 滤波器的频响采样值H(k)相乘，H(k)可事先存放在存储器中，最后再将乘积H(k)X(k)通过快速傅里叶变换（简称IFFT ）还原为时域序列，即得到输出y(n)如图6.3所示。 图6.3 数字滤波器的快速傅里叶变换实现方法 现以FFT 求有限长序列间的卷积及求有限长度列与较长序列间的卷积为例来讨论FFT 的快速卷积方法。 （1） 序列)(n x 和)(n h 的列长差不多。设)(n x 的列长为1N ，)(n h 的列长为2N ，要求 )()(n x n y =N ∑-=-==1 ) ()()(*)()(N r r n h r x n h n x n h 用FFT 完成这一卷积的具体步骤如下： i. 为使两有限长序列的线性卷积可用其循环卷积代替而不发生混叠，必须选择循环卷积长度121-+≥N N N ，若采用基2-FFT 完成卷积运 算，要求m N 2=（m 为整数）。 ii. 用补零方法使)(n x ，)(n h 变成列长为N 的序列。 ?? ?-≤≤-≤≤=10 10)()(11N n N N n n x n x ?? ?-≤≤-≤≤=10 1 0)()(22N n N N n n h n h iii. 用FFT 计算)(),(n h n x 的N 点离散傅里叶变换 )()(k X n x FFT ??→? )()(k H n h FFT ??→? iv. 做)(k X 和)(k H 乘积，)()()(k H k X k Y ?= v. 用FFT 计算)(k Y 的离散傅里叶反变换得 y(n)
实验二快速傅立叶变换（FFT ）及其应用 一、实验目的 1.在理论学习的基础上，通过本实验，加深对FFT 的理解，熟悉FFT 子程序。 2.熟悉应用FFT 对典型信号进行频谱分析的方法 3.了解应用FFT 进行信号频谱分析过程中可能出现的问题以便在实际中正确应用FFT 。 二、实验原理 在各种信号序列中，有限长序列信号处理占有很重要地位，对有限长序列，我们可以使用离散Fouier 变换(DFT)。这一变换不但可以很好的反映序列的频谱特性，而且易于用快速 算法在计算机上实现，当序列x(n)的长度为N 时，它的DFT 定义为： 1 0()()N kn N n X k x n W -==∑， 2n j N N W e -=反换为：10 1()()N kn N k x n X k W N --==∑有限长序列的DFT 是其Z 变换在单位圆上的 等距采样，或者是序列Fourier 变换的等距采样，因此可以用于序列的谱分析。 FFT 并不是与DFT 不同的另一种变换，而是为了减少DFT 运算次数的一种快速算法。它是对变换式进行一次次分解，使其成为若干小点数的组合，从而减少运算量。常用的FFT 是以2为基数的，其长度 N=2L ，它的效率高，程序简单使用非常方便，当要变换的序列长度不等于2的整数次方时，为了使用以2为基数的FFT ，可以用末位补零的方法，使其长度延长至2的整数次方。 在运用DFT 进行频谱分析的过程中可能产生几种问题： (1) 混叠 序列的频谱时被采样信号的周期延拓，当采样速率不满足Nyquist 定理时，就会发生频谱混叠，使得采样后的信号序列频谱不能真实的反映原信号的频谱。 避免混叠现象的唯一方法是保证采样速率足够高，使频谱混叠现象不致出现，即在确定采样频率之前，必须对频谱的性质有所了解，在一般情况下，为了保证高于折叠频率的分量不会出现，在采样前，先用低通模拟滤波器对信号进行滤波。 (2) 泄漏 实际中我们往往用截短的序列来近似很长的甚至是无限长的序列，这样可以使用较短的DFT 来对信号进行频谱分析，这种截短等价于给原信号序列乘以一个矩形窗函数，也相当于在频域将信号的频谱和矩形窗函数的频谱卷积，所得的频谱是原序列频谱的扩展。 泄漏不能与混叠完全分开，因为泄漏导致频谱的扩展，从而造成混叠。为了减少泄漏的影响，可以选择适当的窗函数使频谱的扩散减至最小。 DFT 是对单位圆上Z 变换的均匀采样，所以它不可能将频谱视为一个连续函数，就一定意义上看，用DFT 来观察频谱就好像通过一个栅栏来观看一个图景一样，只能在离散点上看到真实的频谱，这样就有可能发生一些频谱的峰点或谷点被“尖桩的栅栏”所拦住，不能别我们观察到。 减小栅栏效应的一个方法就是借助于在原序列的末端填补一些零值，从而变动DFT 的点数，这一方法实际上是人为地改变了对真实频谱采样的点数和位置，相当于搬动了每一根“尖桩栅栏”的位置，从而使得频谱的峰点或谷点暴露出来。 用FFT 可以实现两个序列的圆周卷积。在一定的条件下，可以使圆周卷积等于线性卷积。一般情况，设两个序列的长度分别为N1和N2，要使圆周卷积等于线性卷积的充要条件是
实验一 快速傅里叶变换之报告 一 、实验目的 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(A,N); h = (abs(Y)); h=h/max(h(1:N)); for n=1:N; string1=strcat('X(',num2str(n-1), ')=',num2str(h(n))); disp(string1); f=(n/T)/N; end stem([0:N-1]/N/T,h); xlabel('?μ?ê/HZ'); ylabel('??·ùX￡¨ejw￡?'); title('·ù?μì?D?'); 上述代码命令中，将FFT 变换后的数字变量K ，在画图时转换成频域中的频率f 。这主 要是根据数字频率与模拟域频率之间的关系： T Ω=ω 其中ω、Ω分别为数字和模拟域中的频率，且N k πω2= f π2=Ω 于是有： NT k f = 运算结果： X(1)=1 X(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)=1 b）信号频率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)=1 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)=1 X(25)= X(26)= X(27)= X(28)= X(29)= X(30)= X(31)=
快速傅立叶变换(FFT)算法实验 摘要：FFT（Fast Fourier Transformation），即为快速傅里叶变换，是离散傅里叶变换的快速算法，它是根据离散傅里叶变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。这种算法大大减少了变换中的运算量，使得其在数字信号处理中有了广泛的运用。本实验主要要求掌握在CCS环境下用窗函数法设计FFT快速傅里叶的原理和方法；并且熟悉FFT快速傅里叶特性；以及通过本次试验了解各种窗函数对快速傅里叶特性的影响等。 引言： 快速傅里叶变换FFT是离散傅里叶变换DFT的一种快速算法。起初DFT的计算在数字信号处理中就非常有用，但由于计算量太大，即使采用计算机也很难对问题进行实时处理，所以并没有得到真正的运用。1965年J.W.库利和T.W.图基提出快速傅里叶变换，采用这种算法能使计算机计算离散傅里叶变换所需要的乘法次数大为减少,特别是被变换的抽样点数N越多,FFT算法计算量的节省就越显著。从此，对快速傅里叶变换（FFT）算法的研究便不断深入，数字信号处理这门新兴学科也随FFT的出现和发展而迅速发展。根据对序列分解与选取方法的不同而产生了FFT的多种算法，基本算法是基2DIT和基2DIF。FFT 的出现，使信号分析从时域分析向频域分析成为可能，极大地推动了信号分析在各领域的实际应用。FFT在离散傅里叶反变换、线性卷积和线性相关等方面也有重要应用。
一、 实验原理： FFT 并不是一种新的变换，它是离散傅立叶变换（DFT ）的一种快速算法。由于我们在计算DFT 时一次复数乘法需用四次实数乘法和二次实数加法；一次复数加法则需二次实数加法。每运算一个X （k ）需要4N 次复数乘法及2N+2（N-1）=2（2N-1）次实数加法。所以整个DFT 运算总共需要4N^2次实数乘法和N*2(2N-1)=2N(2N-1)次实数加法。如此一来，计算时乘法次数和加法次数都是和N^2成正比的，当N 很大时，运算量是可观的，因而需要改进对DFT 的算法减少运算速度。 根据傅立叶变换的对称性和周期性，我们可以将DFT 运算中有些项合并。我们先设序列长度为N=2^L ，L 为整数。将N=2^L 的序列x(n)(n=0,1,……，N-1)，按N 的奇偶分成两组，也就是说我们将一个N 点的DFT 分解成两个N/2点的DFT ，他们又从新组合成一个如下式所表达的N 点DFT ： ∑∑∑∑∑-=+-=-=++ = + =-≤≤= =1 )12(1 20 2为奇 为偶 10 )12()2()()(10, )()]([)(N r k r N N r rk N n nk N n nk N N n nk N W r x W r x W n x W n x N k W n x n x DFT k X
信号与系统的基本思想：把复杂的信号用简单的信号表示，再进行研究。 怎么样来分解信号？任何信号可以用Delta 函数的移位加权和表示。只有系统是线性时不变系统，才可以用单位冲激函数处理，主要讨论各个单位冲激函数移位加权的响应的叠加能得到总的响应。 线性系统（齐次性，叠加定理） 时不变系统 对一个系统输入单位冲激函数，得到的响应为h(t).表征线性时不变系统的非常重要的东西，只要知道了系统对单位冲击函数的响应，就知道了它对任何信号的响应，因为任何信号都可以表示为单位冲激函数的移位加权和。 例如：d(t)__h(t) 那么a*d(t-t0)__a*h(t-t0) -()= ()(t-)d f t f τδττ∝∝? 的响应为-y()=()(-)t f h t d τττ∝ ∝ ? 记为y(t)=f(t)*h(t),称为f(t)和h(t)的卷积 总结为两点：对于现行时不变系统，任何信号可以用单位冲激信号的移位加权和表示，任何信号的响应可以用输入函数和单位冲激函数响应的卷积来表示 连续时间信号和系统的频域分析 时域分析的重点是把信号分解为单位冲激函数的移位加权和，只讨论系统对单位冲激函数的响应。而频域的分析是把信号分解为各种不同频率的正弦函数的加权和，只讨论系统对sinwt 的响应。都是把信号分解为大量单一信号的组合。
周期函数可以展开为傅里叶级数，将矩形脉冲展开成傅里叶级数，得到傅里叶级数的系数 n A sin F = T x x τ 其中0=2 nw x τ。 取样函数sin ()=x S a x 。产生一种震荡，0点的值最大，然后渐渐衰减直至0 第一：对于傅里叶级数的系数，n 是离散的，所以频谱也是离散状的每条谱线都出现在基波频率的整数倍上，其包络是取样函数。 第二：谱线的间距是0w .。零点是0=2nw x τ，02w =T π是谱的基波频率。如果τ不变，T 增大，那么0w 减小，当T 非常大的时候，0w 非常小，谱线近似连续，越来越密，幅度越来越小。 傅里叶变换：非周期函数 正变换：--F jw)= ()iwt f t e dt ∝ ∝?（ 反变换：-1()=()2jnwt f t F jw e dw π ∝∝ ? 常用函数的傅里叶变换（典型非周期信号的频谱）
中南大学 DSP技术实验报告 实验名称：快速傅立叶变换(FFT)算法实验专业班级：信息0602 学生姓名：张倩曦（学号：24） 指导老师：陈宁 完成日期： 2009年12月2日 中南大学·信息科学与工程学院
快速傅立叶变换(FFT)算法实验一．实验目的 1．掌握用窗函数法设计FFT 快速傅里叶的原理和方法； 2．熟悉FFT 快速傅里叶特性； 3．了解各种窗函数对快速傅里叶特性的影响。 二．实验设备 PC 兼容机一台，操作系统为Windows2000(或Windows98，WindowsXP，以下默认为Windows2000)，安装Code Composer Studio 软件。 三．实验原理 1．FFT 的原理和参数生成公式： 公式（1）FFT 运算公式 FFT 并不是一种新的变换，它是离散傅立叶变换（DFT）的一种快速算法。由于我们在计算DFT 时一次复数乘法需用四次实数乘法和二次实数加法；一次复数加法则需二次实数加法。 每运算一个X（k）需要4N 次复数乘法及2N+2（N-1）=2（2N-1）次实数加法。所以整个DFT运算总共需要4N^2 次实数乘法和N*2(2N-1)=2N(2N-1)次实数加法。如此一来，计算时乘法次数和加法次数都是和N^2 成正比的，当N 很大时，运算量是可观的，因而需要改进对DFT 的算法减少运算速度。 根据傅立叶变换的对称性和周期性，我们可以将DFT 运算中有些项合并。我们先设序列长度为N=2^L，L 为整数。将N=2^L 的序列x(n)(n=0,1,……，N-1)，
按N的奇偶分成两组，也就是说我们将一个N 点的DFT 分解成两个N/2 点的DFT，他们又重新组合成一个如下式所表达的N 点DFT： 一般来说，输入被假定为连续的。当输入为纯粹的实数的时候，我们就可以利用左右对称的特性更好的计算DFT。 我们称这样的RFFT 优化算法是包装算法：首先2N 点实数的连续输入称为“进包”。其次N 点的FFT 被连续运行。最后作为结果产生的N 点的合成输出是“打开”成为最初的与DFT 相符合的2N 点输入。使用这一思想，我们可以划分FFT 的大小，它有一半花费在包装输入O（N）的操作和打开输出上。这样的RFFT 算法和一般的FFT 算法同样迅速，计算速度几乎都达到了两次DFT的连续输入。下列一部分将描述更多的在TMS320C55x 上算法和运行的细节。 5．程序流程图：
附录A 拉普拉斯变换及反变换 419
420
421 3． 用查表法进行拉氏反变换 用查表法进行拉氏反变换的关键在于将变换式进行部分分式展开，然后逐项查表进行反变换。设)(s F 是s 的有理真分式 1110 111)()()(a s a s a s a b s b s b s b s A s B s F n n n n m m m m ++++++++= =---- （m n >） 式中系数n n a a a a ,,...,,110-，m m b b b b ,,,110- 都是实常数；n m ,是正整数。按代数定理可将)(s F 展开为部分分式。分以下两种情况讨论。 ① 0)(=s A 无重根 这时，F(s)可展开为n 个简单的部分分式之和的形式。 ∑=-=-++-++-+-=n i i i n n i i s s c s s c s s c s s c s s c s F 122 11)( （F-1） 式中，n s s s ,,,21 是特征方程A(s)＝0的根。i c 为待定常数，称为F(s)在i s 处的留数，可按下式计算： )()(lim s F s s c i s s i i -=→ （F-2） 或 i s s i s A s B c ='= )() ( （F-3） 式中，)(s A '为)(s A 对s 的一阶导数。根据拉氏变换的性质，从式（F-1）可求得原函数 []??????-==∑=--n i i i s s c L s F L t f 11 1 )()(＝t s n i i i e c -=∑1 (F-4) ② 0)(=s A 有重根 设0)(=s A 有r 重根1s ，F(s)可写为 ()) ()()() (11n r r s s s s s s s B s F ---= + = n n i i r r r r r r s s c s s c s s c s s c s s c s s c -++-++-+-++-+-++-- 11 111 111)()()( 式中，1s 为F(s)的r 重根，1+r s ，…, n s 为F(s)的n-r 个单根；
快速傅里叶变换实验报告 班级： 姓名： 学号:
快速傅里叶变换 一．实验目的 1.在理论学习的基础上，通过本实验加深对快速傅立叶变换的理解； 2.熟悉并掌握按时间抽取FFT 算法的程序； 3.了解应用FFT 进行信号频谱分析过程中可能出现的问题，例如混淆、泄漏、栅栏效应等，以便在实际中正确应用FFT 。 二．实验内容 1.仔细分析教材第六章‘时间抽取法FFT ’的算法结构，编制出相应的用FFT 进行信号分析的C 语言（或MATLAB 语言）程序； 2.用FFT 程序分析正弦信号 ()sin(2)[()(*)],(0)1y t f t u t u t N T t u π=---∞<<+∞=设 分别在以下情况进行分析并讨论所得的结果： a ） 信号频率f ＝50Hz ，采样点数N=32,采样间隔T=0.000625s b ） 信号频率f ＝50Hz ，采样点数N=32,采样间隔T=0.005s c ） 信号频率f ＝50Hz ，采样点数N=32,采样间隔T=0.0046875s d ） 信号频率f ＝50Hz ，采样点数N=32,采样间隔T=0.004s e ） 信号频率 f ＝50Hz ，采样点数N=64,采样间隔T=0.000625s f ） 信号频率f ＝250Hz ，采样点数N=32,采样间隔T=0.005s g ） 将c ） 信号后补32个0，做64点FFT 三．实验要求 1.记录下实验内容中各种情况下的X (k)值，做出频谱图并深入讨论结果，说明参数的变化对信号频谱产生哪些影响。频谱只做模特性，模的最大值＝1，全部归一化；
2.打印出用C 语言（或MATLAB 语言）编写的FFT 源程序，并且在每一小段处加上详细的注释说明； 3.用C 语言（或MATLAB 语言）编写FFT 程序时，要求采用人机界面形式： N , T , f 变量均由键盘输入，补零或不补零要求设置一开关来选择。 四．实验分析 对于本实验进行快速傅里叶变换，依次需要对信号进行采样，补零（要求补零时），码位倒置，蝶形运算，归一化处理并作图。 此外，本实验要求采用人机界面形式，N,T,F 变量由键盘输入，补零或不补零设置一开关来选择。 1.采样 本实验进行FFT 运算，给出的是正弦信号，需要先对信号进行采样，得到有限 长序列()n x ， N n ...... 2,1,0= Matlab 实现： t=0:T:T*(N-1); x=sin(2*pi*f*t); 2.补零 根据实验要求确定补零与否，可以用if 语句做判断，若为1，再输入补零个数， 并将补的零放到采样得到的序列的后面组成新的序列，此时新的序列的元素个数等于原采样点个数加上补零个数，并将新的序列个数赋值给N 。 Matlab 实现： a=input('是否增加零点? 是请输入1 否请输入0\n'); if (a) ZeroNum=input('请输入增加零点的个数:\n'); else ZeroNum=0; end if (a) x=[x zeros(1, ZeroNum)];%%指令zeros(a,b)生成a 行b 列全0矩阵，在单行矩阵x 后补充0 end N=N+ZeroNum; 3.码位倒置 本实验做FFT 变换的级数为M ，N M 2log =
快速傅里叶变换实验报告 快速傅里叶变换实验报告 机械34班刘攀 2019010558 一、基本信号（函数）的FFT变换 1. x(t)=sin(ω0t+)+sin2ω0t+cos3ω0t 6 1) 采样频率fs=8f0，截断长度N=16； 取ω0=2πrad/s，则f0=1Hz，fs=8Hz，频率分辨率?f=?f=fs=0.5Hz。 Nπ最高频率fc=3f0=3Hz，fs>2fc，故满足采样定理，不会发生混叠现象。截断长度T=2T0，整周期截取，不会发生栅栏效应。理论上有一定的泄漏，但在整周期截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图如下： 幅值误差?A=0，相位误差??=0。 2) 采样频率fs=8f0，截断长度N=32； 取ω0=2πrad/s，则f0=1Hz，fs=8Hz，频率分辨率?f=?f=fs=0.25Hz。 N最高频率fc=3f0=3Hz，fs>2fc，故满足采样定理，不会发生混叠现象。截断长度T=4T0，整周期截取，不会发生栅栏效应。理论上有一定的泄漏，但在整周期截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图如下： 幅值误差?A=0，相位误差??=0。 2. x(t)=sin(ω0t+π 6)+sin11ω0t 1) 采样频率fs=8f0，截断长度N=16； 取ω0=2πrad/s，则f0=1Hz，fs=8Hz，频率分辨率?f=?f=fs=0.5Hz。 N最高频率 fc=11f0=11Hz，fs 漏，但在整周期截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图：
由上图可以看出，并未体现出11f0的成分，说明波形出现混叠失真。为了消除混叠 现象，应加大采样频率，使之大于等于 22Hz。 f0处的幅值误差?A=0，11f0处由于出现 了混叠现象，幅值误差没有意义；相位误差??=0。 2) 采样频率fs=32f0，截断长度N=32； 取ω0=2πrad/s，则f0=1Hz，fs=32Hz，频率分辨率?f=?f=fs=1Hz。 N最高频率 fc=11f0=11Hz，fs>2fc，故满足采样定理，不会发生混叠现象。 漏，但在整周期截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图： 该频谱图体现出了f0和11f0的成分，说明未失真，且幅值均为1,。幅值误差?A=0，相位误差??=0。 3. x(t)=0t 1) 采样频率fs=8f0，截断长度N=16； 取ω0=2πrad/s，则f0=1Hz，fs=8Hz，频率分辨率?f=?f=fs=0.5Hz。 N最高频率f cf 0Hz，fs>2fc，故满足采样定理，不会发生混叠现象。 频谱图： 在忽略旁瓣信号的情况下，可近似认为： x(t)≈0.9098cos(3ω0t+56.9520?) 故幅值误差?A=0.9096-1=-0.0904，相位误差??=56.9520?。 2) 采样频率fs=32f0，截断长度N=32； 取ω0=2πrad/s，则f0=1Hz，fs=32Hz，频率分辨率?f=?f=fs=1Hz。N最高频率f cf 0Hz，fs>2fc，故满足采样定理，不会发生混叠现象。 频谱图： 在忽略旁瓣信号的情况下，可近似认为：
时域信号 弧频率表示的 傅里叶变换 注释 1 | 线性 2时域平移 3频域平移, 变换2的频域对应 \ 4 如果值较大，则会收缩 到原点附近，而会扩 散并变得扁平. 当| a | 趋向无 穷时，成为Delta函数。 5 傅里叶变换的二元性性质。通过 交换时域变量和频域变量 得到. 6 / 傅里叶变换的微分性质 7变换6的频域对应
8 表示和的卷积—这 就是卷积定理 - 9 矩形脉冲和归一化的sinc函数 10变换10的频域对应。矩形函数是理想的低通滤波器，sinc函数是这类滤波器对反因果冲击的响应。 11- tri是三角形函数 12变换12的频域对应 13高斯函数exp( ? αt2) 的傅里叶变换是他本身. 只有当Re(α) > 0时，这是可积的。 ￥14 15 16》 a>0
18δ(ω) 代表狄拉克δ函数分布. 这个变换展示了狄拉克δ函数的重要性：该函数是常函数的傅立叶变换 【 19 变换23的频域对应20由变换3和24得到. 21` 由变换1和25得到，应用了欧拉公 式: cos(at) = (e iat + e?iat) / 2. 22由变换1和25得到 23这里, n是一个自然数. δ(n)(ω) 是狄拉克δ函数分布的n阶微分。这个变换是根据变换7和24得到的。将此变换与1结合使用，我们可以变换所有多项式。 / 24此处sgn(ω)为符号函数；注意此变换与变换7和24是一致的. 25变换29的推广. 17变换本身就是一个公式
26【 变换29的频域对应. 27此处u(t)是单位阶跃函数; 此变换根据变换1和31得到. 28u(t)是单位阶跃函数，且a > 0. 34狄拉克梳状函数——有助于解释或理解从连续到离散时间的转变.
常用拉普拉斯变换总结 1、指数函数 00)(≥
? ? ∞ -∞ -∞ ----==0 d d ][t s e s e t t te t L st st st 2 01d 1s t e s st == ?∞- 6、正弦函数 0sin 0 )(≥
快速傅里叶变换实验报告 机械34班 攀 2013010558 一、 基本信号（函数）的FFT 变换 1. 000()sin()sin 2cos36 x t t t t πωωω=+++ 1) 采样频率08s f f =，截断长度N=16； 取02ωπ=rad/s ，则0f =1Hz ，s f =8Hz ，频率分辨率f ?=s f f N ?==0.5Hz 。 最高频率c f =30f =3Hz ，s f >2c f ，故满足采样定理，不会发生混叠现象。 截断长度02T T =，整周期截取，不会发生栅栏效应。理论上有一定的泄漏，但在整周期 截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图如下：
幅值误差0A ?=，相位误差0??=。 2) 采样频率08s f f =，截断长度N=32； 取02ωπ=rad/s ，则0f =1Hz ，s f =8Hz ，频率分辨率f ?=s f f N ?==0.25Hz 。 最高频率c f =30f =3Hz ，s f >2c f ，故满足采样定理，不会发生混叠现象。 截断长度04T T =，整周期截取，不会发生栅栏效应。理论上有一定的泄漏，但在整周期 截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图如下：
幅值误差0A ?=，相位误差0??=。 2. 00()sin()sin116 x t t t πωω=++ 1) 采样频率08s f f =，截断长度N=16； 取02ωπ=rad/s ，则0f =1Hz ，s f =8Hz ，频率分辨率f ?=s f f N ?==0.5Hz 。 最高频率c f =110f =11Hz ，s f <2c f ，故不满足采样定理，会发生混叠现象。 截断长度02T T =，整周期截取，不会发生栅栏效应。理论上有一定的泄漏，但在整周期 截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图：
快速傅里叶变换实验
————————————————————————————————作者：————————————————————————————————日期： ?
实验七快速傅里叶变换实验 20１1０１0５４1?机１4 林志杭 一、实验目的 1.加深对几个特殊概念的理解：“采样”……“混叠”;“窗函数”（截断）……“泄漏”；“非整周期截取”……“栅栏”。 ２．加深理解如何才能避免“混叠”，减少“泄漏”，防止“栅栏”的方法和措施以及估计这些因素对频谱的影响。 3．对利用通用微型计算机及相应的FFＴ软件，实现频谱分析有一个初步的了解。 二、实验原理 为了实现信号的数字化处理，利用计算机进行频谱分析――计算信号的频谱。由于计算机只能进行有限的离散计算(即ＤFT)，因此就要对连续的模拟信号进行采样和截断。而这两个处理过程可能引起信号频谱的畸变，从而使DFT的计算结果与信号的实际频谱有误差。有时由于采样和截断的处理不当，使计算出来的频谱完全失真。因此在时域处理信号时要格外小心。 时域采样频率过低，将引起频域的“混叠”。为了避免产生“混叠”,要求时域采样时必须满足采样定理,即:采样频率fｓ必须大于信号中最高频率fc的2倍（ｆｓ＞2fc)。因此在信号数字处理中，为避免混叠，依不同的信号选择合适的采样频率将是十分重要的。 频域的“泄漏”是由时域的截断引起的。时域的截断使频域中本来集中的能量向它的邻域扩散(如由一个δ(f）变成一个siｎc(f），而泄漏的旁瓣将影响其它谱线的数值。时域截断还会引起“栅栏效应”,对周期信号而言，它是由于截断长度不等于周期信号的周期的整数倍而引起的。因此避免“栅栏”效应的办法就是整周期截断。 综上所述,在信号数字化处理中应十分注意以下几点： 1．为了避免“混叠”,要求在采样时必须满足采样定理。 为了减少“泄漏”，应适当增加截断长度和选择合适的窗 对信号进行整周期截取,则能消除“栅栏数应”。 增加截断长度，则可提高频率分辨率。 三、预习内容 熟悉Mａtlab语言、函数和使用方法；利用Ｍatｌab所提供的ＦFT函数编写程序。 四、实验内容及步骤 调通所编写的程序,对下列信号〔函数〕进行离散FFT变换,根据题目的要求……FFT变换点数〔截断长度〕及采样频率,计算各点的傅里叶变换值,画出频谱图,对典型的谱线标出其幅值及相角。 （-)内容: 1． t t t t x 3 cos 2 sin ) 6 sin( )(ω ω π ω+ + + = 代码： N=ｉnput(＇N='); n=input（'ｎ=')；ｔ=1：1：Ｎ；
实验七 快速傅立叶变换（FFT ）实验 一 实验目的 1． 熟悉CCS 集成开发环境； 2． 了解FFT 的算法原理和基本性质； 3． 熟悉DSP 中cmd 文件的作用及对它的修改； 4． 学习用FFT 对连续信号和时域信号进行频谱分析的方法； 5． 利用DSPLIB 中现有的库函数； 6． 了解DSP 处理FFT 算法的特殊寻址方式； 7． 熟悉对FFT 的调试方法。 二 实验内容 本实验要求使用FFT 变换对一个时域信号进行频谱分析，同时进行IFFT 。这里用到时域信号可以是来源于信号发生器输入到CODEC 输入端，也可以是通过其他工具计算获取的数据表。本实验使用Matlab 语言实现对FFT 算法的仿真，然后将结果和DSP 分析的结果进行比较，其中原始数据也直接来自Matlab 。 三 实验原理 一个N 点序列][k x 的DFT ][m X ,以及IDFT 分别定义为： 1,,1,0,][][1 0-==∑-=N m W k x m X km N N k 1,,1,0,][1 ][1 -== --=∑ N k W m X N k x km N N m 如果利用上式直接计算DFT,对于每一个固定的m,需要计算N 次复数乘法，N-1次加法,对于N 个不同的m,共需计算N 的2次方复数乘法,N*(N-1)次复数加法.显然,随着N 的增加,运算量将急剧增加, 快速傅里叶算法有效提高计算速度（本例使用基2 FFT 快速算法）,利用FFT 算法只需(N/2)logN 次运算。 四 知识要点 . 1、 CMD 文件的功能及编写 2、 一种特殊的寻址方式：间接寻址 间接寻址是按照存放在某个辅助寄存器的16位地址寻址的。C54x 的8个辅助寄存器 （AR0—AR7）都可以用来寻址64K 字数据存储空间中的任何一个存储单元。 3、 TMS320C54x DSPLIB 中关于FFT 变换的一些函数的调用（SPRA480B.pdf ） 利用DSPLIB 库时，在主程序中要包含头文件：54xdsp.lib 4、 FFT 在CCS 集成开发环境下的相关头文件 #include
CEＮTＲＡL SOUTH UNＩVEＲSIＴY 数字信号处理实验报告 题目快速傅里叶变换 学生姓名 学院物理与电子学院 专业班级电子信息科学与技术110５班 学号140４11０72
实验五快速傅里叶变换 一、实验仪器 PC机一台、JQ—ＳＯＰC开发系统实验箱及辅助软件(DSP Ｂuiｌdeｒ、Matｌab/Siｍulink、QuaｒtｕｓII、Modｅｌsim)。 二、实验目得 1、了解快速傅里叶变换得基本结构组成、 2、学习使用DＳP Bｕiｌder设计FFT。 三、实验原理 １、ＦＦT得原理: 快速傅里叶变换(FFT)就是离散傅里叶变换(DＦＴ)得一种高效运算方法,它大大简化了DFT得运算过程,使运算时间缩短几个数量级、FFＴ算法可以分为按时间抽取(ＤIＴ)与按频率抽取(ＤIＦ)两类,输入也可分为实数与复数两种情况。八点时间抽取基—2FFＴ算法信号流图如图１示: 图１8点基-2 DIT-ＦFＴ信号流图 四、实验步骤 1、将桌面得my_fft_８。mdｌ拷贝到“D:＼Pｒｏgram Filｅｓ\ＭATＬAB７１\worｋ”(ＭATLAB安装目录下得wｏrk文件夹)处,并双击打开。
图5－１快速傅里叶变换系统图 图5—2 快速傅里叶变换子系统1图
图5—３快速傅里叶变换子系统2图
图5—3 快速傅里叶变换子系统3图 2、点击工具栏即可开始系统级sｉmｕlink仿真,以验证该模型得正确性、在仿真进行过程中分别将三个输入控制开关打到000、０01、０10、0１1、100以选择五组输入数据进行FFＴ运算。 (1)当开关打到０00时选择第一组数据｛２.0,２.0,4。0,7。0,3.0,5。0,5、0,8.0｝,其运算结果应为３6、-2、４1+3.84i、—4+8i、0、42１９＋１.８４4i、—8、０。４102-１、84ｉ、-4-8i、－２、4２２－3.８４4ｉ。 (2)当开关打到0０１时选择第二组数据｛1。1,５。0,1０.5,１5、3,20。２,25。7,30、6,４0。1｝,其运算结果应该为148、5、－16、1+52、35i、－１９.8+２4、7i、-2２、０2+12。2５ｉ、—23.7、-22、1-12。15ｉ、－19.8－24。7i、－16.９—５2、45i、 (3)当开关打到０1０时选择第三组数据{5.6,5、0,１5、8,20.4,25、２,35、1,40、２,4５。０},其运算结果结果应该为192。3、-2３.3９＋63、１9i、—２5。2+２5、3i、-１5.６9+14。49i、—１８、7、－15。８1—14、３9ｉ、－２５。２-25。3i、－23。5１-63、29ｉ。 (4)当开关打到０11选择第四组数据{10、２,１5。3,18。1,20。3,24.2,30、０,３５.2,42、3},