实验四应用快速傅里叶变换对信号进行频谱分析
用FFT对信号做频谱分析
用FFT对信号做频谱分析傅里叶变换（Fourier Transform）是一种将信号从时域转换到频域的数学方法，可用于信号的频谱分析。
通过傅里叶变换，我们可以将时域上的信号转换为频域上的频谱，帮助我们理解信号的频率组成以及各个频率分量的强弱。
频谱分析是对信号进行频率分析的过程，是了解信号在频域上的特性和频率成分的一种方法。
通过频谱分析，我们可以获得信号的频率分布情况，帮助我们了解信号的频率成分、频率峰值等信息。
在进行频谱分析时，常用的方法之一是采用快速傅里叶变换（FFT）。
FFT是一种高效的算法，能够快速计算离散傅里叶变换（DiscreteFourier Transform）。
下面将详细介绍FFT在频谱分析中的应用。
首先，我们需要将待分析的信号转换为数字信号，并对其进行采样，得到一个离散的信号序列。
然后，使用FFT算法对这个离散信号序列进行傅里叶变换，得到信号的频谱。
在进行FFT之前，需要进行一些预处理工作。
首先，需要将信号进行加窗处理，以减少泄露效应。
加窗可以选择矩形窗、汉宁窗、汉明窗等，不同的窗函数对应不同的性能和应用场景。
其次，需要对信号进行零填充，即在信号序列末尾添加零值，以增加频谱的分辨率。
零填充可以提高频谱的平滑度，使得频域上的分辨率更高。
接下来，我们使用FFT算法对经过加窗和零填充的信号序列进行傅里叶变换。
FFT算法将离散信号变换为离散频谱，得到信号的频率成分和强度。
FFT结果通常呈现为频率和振幅的二维图像，横轴表示频率，纵轴表示振幅。
通过观察频谱图像，我们可以得到一些关于信号的重要信息。
首先，我们可以观察到信号的频率成分，即信号在不同频率上的分布情况。
在频谱图像中，高峰表示信号在该频率上强度较高，低峰表示信号在该频率上强度较低。
其次，我们可以通过峰值的位置和强度来分析信号的主要频率和频率成分。
频谱图像上的峰值位置对应着信号的主要频率，峰值的高度对应着信号在该频率上的强度。
最后，我们还可以通过观察频谱图像的整体分布情况，来获取信号的频率范围和频率分布的特点。
应用快速傅里叶变换对信号进行频谱分析实验报告
应用快速傅里叶变换对信号进行频谱分析实验报告实验报告：快速傅里叶变换在信号频谱分析中的应用【引言】傅里叶分析是一种重要的信号处理方法，可将时域信号转换为频域信号，并且可以分解信号的频谱成分。
传统的傅里叶变换算法在计算复杂度方面较高，为了降低计算的复杂度，人们提出了快速傅里叶变换（FFT）算法。
本实验旨在通过应用快速傅里叶变换对信号进行频谱分析，研究信号的频谱特性。
【实验目的】1.了解傅里叶变换的基本原理，研究其在信号处理中的应用；2.学习快速傅里叶变换算法的原理和优点；3.通过实验操作，观察信号的频谱特性，分析实验结果。
【实验原理】1. 傅里叶变换（FT）：对于一个连续时间域信号x(t)，其傅里叶变换可表示为X(ω) = ∫[t=−∞,∞]x(t)e^(-jωt)dt，其中X(ω)表示频域上的信号分量，ω为角频率。
2.快速傅里叶变换（FFT）算法：FFT是一种离散时间域信号的频谱分析方法，具有较低的计算复杂度。
FFT算法使用了分治法的思想，将信号分解为较小的频谱分量，并通过递归计算得到完整的频谱图。
3.FFT算法的步骤：1)若信号长度为N，则将其分为两个长度为N/2的子信号；2)对子信号进行FFT变换；3)将两个子信号拼接起来，得到完整信号的频谱分量。
【实验步骤】1.准备实验材料和装置：计算机、FFT分析软件、信号发生器等；2.设置信号发生器的输出参数，例如频率、幅度等；3.连接信号发生器和计算机，打开FFT分析软件；4.在FFT软件中选择输入信号通道，设置采样参数等；5.开始实验，观察计算机屏幕上的频谱图；6.调整信号发生器的参数，重复第5步，记录实验结果；7.结束实验，关闭设备。
【实验结果与分析】我们选择了一个简单的正弦波信号作为输入信号，信号频率设置为100Hz，幅度设置为1V。
在进行频谱分析之前，我们通过示波器观察到一个明显的正弦波信号。
接下来，我们将信号输入到计算机上的FFT分析软件中，进行频谱分析。
实验四应用快速傅里叶变换对信号进行频谱分析
实验四应用快速傅里叶变换对信号进行频谱分析引言：频谱分析是信号处理领域中的重要技术之一，可以用于研究信号的频率特性和频域内的信号成分。
傅里叶变换是一种能将时域信号转换为频域信号的数学工具，通过将信号分解成一系列频率分量来分析信号。
快速傅里叶变换（FFT）是一种高效的计算傅里叶变换的方法，尤其适合实时信号处理。
实验目的：1.理解傅里叶变换在频谱分析中的应用；2.掌握使用FFT对信号进行频谱分析的方法；3.实现频谱分析并得出相应的频谱图。
实验器材和材料：1.信号源（例如信号发生器）；2.电脑或数字信号处理器（DSP）；3.音频线或数据线连接信号源和电脑或DSP。
实验步骤：1.确定实验所需信号源的类型和参数，例如正弦信号、方波信号或任意信号；2.连接信号源和电脑或DSP，确保信号源输出的信号能够被电脑或DSP接收；3. 在电脑或DSP上选择合适的软件或编程语言环境，例如MATLAB、Python或C；4.编写程序或命令以控制信号源产生相应的信号，并将信号输入到电脑或DSP中；5.读取信号，并使用FFT对信号进行傅里叶变换；6.分析得到的频谱数据，绘制频谱图；7.对得到的频谱图进行解读和分析。
实验注意事项：1.在选择信号源和连接电脑或DSP时，注意信号源的输出范围和电脑或DSP的输入范围，避免信号超出范围导致损坏设备；2.根据实际需要选择合适的采样率和采样点数，以保证能够对信号进行充分的频谱分析；3.在进行FFT计算时，注意选择适当的窗函数和重叠率，以克服频谱分析中的泄漏效应。
实验结果与讨论：通过对信号进行频谱分析，我们可以得到信号的频率特性和频域内的成分信息。
根据得到的频谱图，我们可以分析信号的主要频率分量、功率谱密度以及可能存在的干扰或噪声。
通过对频谱图的解读和分析，可以帮助我们理解信号的特征和变化规律，为后续的信号处理和应用提供有价值的信息。
结论：本实验通过应用快速傅里叶变换对信号进行频谱分析，从而得到信号在频域内的成分信息并绘制出频谱图。
应用FFT实现信号频谱分析
应用FFT实现信号频谱分析一、快速傅里叶变换（FFT）原理快速傅里叶变换是一种将时域信号转换为频域信号的算法，它通过将信号分解为不同频率的正弦波的和，来实现频谱分析。
FFT算法是一种高效的计算DFT（离散傅里叶变换）的方法，它的时间复杂度为O(nlogn)，在实际应用中得到广泛使用。
二、FFT算法FFT算法中最基本的思想是将DFT进行分解，将一个长度为N的信号分解成长度为N/2的两个互为逆序的子信号，然后对这两个子信号再进行类似的分解，直到分解成长度为1的信号。
在这一过程中，可以通过频谱折叠的性质，减少计算的复杂度，从而提高计算效率。
三、FFT实现在实际应用中，可以使用Matlab等软件来实现FFT算法。
以Matlab 为例，实现FFT可以分为以下几个步骤：1.读取信号并进行预处理，如去除直流分量、归一化等。
2. 对信号进行FFT变换，可以调用Matlab中的fft函数，得到频域信号。
3.计算频谱，可以通过对频域信号进行幅度谱计算，即取频域信号的模值。
4.可选地，可以对频谱进行平滑处理，以降低噪音干扰。
5.可选地，可以对频谱进行归一化处理，以便于分析和比较不同信号的频谱特性。
四、应用1.音频处理：通过分析音频信号的频谱，可以实现音频特性的提取，如频率、振幅、共振等。
2.图像处理：通过分析图像信号的频谱，可以实现图像特征的提取，如纹理、边缘等。
3.通信系统：通过分析信号的频谱，可以实现信号的调制解调、频谱分配等功能。
4.电力系统：通过分析电力信号的频谱，可以实现电力质量分析、故障检测等。
总结：应用FFT实现信号频谱分析是一种高效的信号处理方法，通过将时域信号转换为频域信号，可以实现对信号频谱特性的提取和分析。
在实际应用中，我们可以利用FFT算法和相应的软件工具，对信号进行频谱分析，以便于进一步的研究和应用。
FFT算法分析实验实验报告
FFT算法分析实验实验报告一、实验目的快速傅里叶变换（Fast Fourier Transform，FFT）是数字信号处理中一种非常重要的算法。
本次实验的目的在于深入理解 FFT 算法的基本原理、性能特点，并通过实际编程实现和实验数据分析，掌握 FFT 算法在频谱分析中的应用。
二、实验原理FFT 算法是离散傅里叶变换（Discrete Fourier Transform，DFT）的快速计算方法。
DFT 的定义为：对于长度为 N 的序列 x(n)，其 DFT 为X(k) ＝∑n=0 到 N-1 x(n) e^（j 2π k n ／ N) ，其中 j 为虚数单位。
FFT 算法基于分治法的思想，将 N 点 DFT 分解为多个较小规模的DFT，从而大大减少了计算量。
常见的 FFT 算法有基 2 算法、基 4 算法等。
三、实验环境本次实验使用的编程语言为 Python，主要依赖 numpy 库来实现 FFT 计算和相关的数据处理。
四、实验步骤1、生成测试信号首先，生成一个包含不同频率成分的正弦波叠加信号，例如100Hz、200Hz 和 300Hz 的正弦波。
设定采样频率为 1000Hz，采样时间为 1 秒，以获取足够的采样点进行分析。
2、进行 FFT 计算使用 numpy 库中的 fft 函数对生成的测试信号进行 FFT 变换。
3、频谱分析计算 FFT 结果的幅度谱和相位谱。
通过幅度谱确定信号中各个频率成分的强度。
4、误差分析与理论上的频率成分进行对比，计算误差。
五、实验结果与分析1、幅度谱分析观察到在 100Hz、200Hz 和 300Hz 附近出现明显的峰值，对应于生成信号中的频率成分。
峰值的大小反映了相应频率成分的强度。
2、相位谱分析相位谱显示了各个频率成分的相位信息。
3、误差分析计算得到的频率与理论值相比，存在一定的误差，但在可接受范围内。
误差主要来源于采样过程中的量化误差以及 FFT 算法本身的近似处理。
应用快速傅里叶变换对信号进行频谱分析实验报告
应用快速傅里叶变换对信号进行频谱分析2.1 实验目的1、通过本实验，进一步加深对DFT 算法原理和基本性质的理解，熟悉FFT 算法原理和FFT 子程序应用2、掌握应用FFT 对信号进行频谱分析的方法。
3、通过本次实验进一步掌握频域采样定理。
4、了解应用FFT 进行信号频谱分析过程中可能出现的问题，以便在实际中正。
确应用FFT 。
2.2实验原理与方法对于有限长序列我们可以使用离散傅里叶变换（DFT ）。
这一变换不但可以好地反映序列的频域特性，而且易于用快速傅里叶变换在计算机上实现当序列x(n)的长度为N 时，它的离散傅里叶变换为：10()[()]()N knN n X k DFT x n x n W -===∑其中(2/)j N N W e π-=，它的反变换定义为：11()[()]()N kn Nk x n IDFT X k X k WN--===∑比较Z 变换公式，令k N z W -=则10()|()[()]k NN nkN z W n X z x n W DFT x n --====∑因此有()()|k Nz W X k X z -==。
所以，X(k)是x(n)的Z 变换在单位圆上的等距采样，或者说是序列傅里叶变换的等距采样。
DFT 是对序列傅里叶变换的等距采样，因此可以用于对序列的频谱分析。
在运用DFT 进行频谱分析的过程中有可能产生三种误差： 1、混叠现象序列的频谱是原模拟信号频谱的周期延拓，周期为2/T π。
因此，当采样频率小于两倍信号的最大频率时，经过采样就会发生频谱混叠，使采样后的信号序列频谱不能真实反映原信号的频谱。
2、泄漏现象实际中信号序列往往很长，常用截短的序列来近似它们，这样可以用较短的DFT 对信号进行频谱分析，这种截短等价于给原信号序列乘以一个矩形函数。
这样得到的频谱会将原频谱扩展开。
3、栅栏效应DFT 是对单位圆上Z 变换的均匀采样，所以它不可能将频谱视为一个连续函数。
实验三用FFT对信号作频谱分析_实验报告
实验三用FFT对信号作频谱分析_实验报告一、实验目的1.学习使用FFT（快速傅里叶变换）对信号进行频谱分析；2.掌握频谱分析的基本原理和方法；3.熟悉使用MATLAB进行频谱分析的操作。
二、实验原理FFT是一种基于傅里叶变换的算法，可以将时域信号转换为频域信号，并将信号的频谱特征展示出来。
在频谱分析中，我们通过分析信号的频谱可以获得信号的频率、幅值等信息，从而对信号的性质和特征进行研究。
对于一个连续信号，我们可以通过采样的方式将其转换为离散信号，再利用FFT算法对离散信号进行频谱分析。
FFT算法可以将信号从时域转换到频域，得到离散的频谱，其中包含了信号的频率分量以及对应的幅值。
MATLAB中提供了fft函数，可以方便地对信号进行FFT分析。
通过对信号进行FFT操作，可以得到信号的频谱图，并从中提取出感兴趣的频率信息。
三、实验步骤1.准备工作：（2）建立新的MATLAB脚本文件。
2.生成信号：在脚本中，我们可以通过定义一个信号的频率、幅值和时间长度来生成一个信号的波形。
例如，我们可以生成一个频率为1000Hz，幅值为1的正弦波信号，并设置信号的时间长度为1秒。
3.对信号进行FFT分析：调用MATLAB中的fft函数，对信号进行FFT分析。
通过设置采样频率和FFT长度，可以得到信号的频谱。
其中，采样频率是指在单位时间内连续采样的次数，FFT长度是指离散信号的样本点数。
4.绘制频谱图：调用MATLAB中的plot函数，并设置x轴为频率，y轴为幅值，可以绘制出信号的频谱图。
频谱图上横坐标表示信号的频率，纵坐标表示信号的幅值，通过观察可以得到信号的频率分布情况。
四、实验结果在实验过程中，我们生成了一个频率为1000Hz，幅值为1的正弦波信号，并对其进行FFT分析。
通过绘制频谱图，我们发现信号在1000Hz处有最大幅值，说明信号主要由这一频率成分组成。
五、实验总结本实验通过使用FFT对信号进行频谱分析，我们可以方便地从信号的波形中提取出频率分量的信息，并绘制出频谱图进行观察。
实验一应用快速傅里叶变换对信号进行频谱分析
实验一应用快速傅里叶变换对信号进行频谱分析快速傅里叶变换（Fast Fourier Transform, FFT）是一种高效的算法，用于将时域信号转换为频域信号。
频谱分析是通过对信号进行傅里叶变换来研究信号的频率成分和频率分布的过程。
在实验中，我们将使用FFT算法来对一个信号进行频谱分析。
首先，我们需要了解一些基本概念。
信号的频谱表示了信号在不同频率下的能量分布。
频率表示了信号中发生变化的速度，而幅度则表示了信号在该频率下的强度。
通过对信号进行FFT变换，我们可以将信号从时域转换为频域，得到信号的频谱。
在实验中，我们将使用Python语言来实现信号的FFT变换和频谱分析。
首先，我们需要导入一些必要的库。
import numpy as npimport matplotlib.pyplot as plt我们将创建一个测试信号，然后使用FFT函数对其进行变换和分析。
#创建一个测试信号fs = 1000 # 采样率T = 1 / fs # 采样周期t = np.arange(0, 1, T) # 时间序列f1=10#第一个频率成分f2=100#第二个频率成分A1=2#第一个频率成分的幅度A2=0.5#第二个频率成分的幅度y = A1 * np.sin(2 * np.pi * f1 * t) + A2 * np.sin(2 * np.pi * f2 * t) # 合成信号接下来，我们使用FFT函数对信号进行变换，并绘制其频谱图。
#使用FFT对信号进行变换Y = np.fft.fft(y)#计算频谱N = len(Y) # 信号的长度freq = np.fft.fftfreq(N, T) # 计算频率轴powspec = np.abs(Y) ** 2 / N # 计算功率谱#绘制频谱图plt.figureplt.plot(freq, powspec)plt.xlabel('Frequency (Hz)')plt.ylabel('Power Spectrum')plt.title('Spectrum Analysis')plt.show在频谱图中，横轴表示频率，纵轴表示功率谱，即信号在不同频率下的能量分布。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
二、实验原理与方法
一个连续信号 x a (t ) 的频谱可以用它的傅立叶变换表示为
+∞
∫x
a
(t )e − jΩt dt
（2-1）
如果对该信号进行理想采样，可以得到采样序列
x(n) = x a (nT )
同样可以对该序列进行 z 变换，其中 T 为采样周期
（2-2）
X ( z ) = ∑ x( n) z − n
− ( n− p ) q ,0 ≤ n ≤ 15 x a (n) = e 0, else
2
（2）衰减正弦序列： xb (n) =
e −αn sin 2πfn,0 ≤ n ≤ 15 0, else
（3）三角波序列：
n + 1,0 ≤ n ≤ 3 xc (n) = 8 − n,4 ≤ n ≤ 7 0, else 4 − n,0 ≤ n ≤ 3 x d (n) = n − 3,4 ≤ n ≤ 7 0, else
如果您在阅读过程中发现疏漏和错误，请您尽快和编者取得联系 network@ cxh@
5
M
。它的运算效率高，程
序比较简单，使用也十分地方便。当需要进行变换的序列的长度不是 2 的整数次方的时候， 为了使用以 2 为基的 FFT，可以用末尾补零的方法，使其长度延长至 2 的整数次方。IFFT 一般可以通过 FFT 程序来完成，比较式（2-7）和（2-8） ，只要对 X(k)取共轭，进行 FFT 运 算，然后再取共轭，并乘以因子 1/N，就可以完成 IFFT。
如果您在阅读过程中发现疏漏和错误，请您尽快和编者取得联系 network@ cxh@ 4
中国科学技术大学电子工程与信息科学系
多媒体通信实验室 （Copyright 2000）
五、实验报告要求
1、在实验报告中简述实验目的和实验原理要点。 2、在实验报告中附上在实验过程中记录的各个信号序列的时域和幅频特性曲线，分析所得 到的结果图形，说明各个信号的参数变化对其时域和幅频特性的影响。 3、总结实验中的主要结论。 4、回答思考题。
X (e jω ) =
1 +∞ ω − 2πm ) Xa( j ∑ T T −∞
（2-6）
即序列的频谱是采样信号频谱的周期延拓。从式（2-6）可以看出，只要分析采样序列的频
1
如果您在阅读过程中发现疏漏和错误，请您尽快和编者取得联系 network@ cxh@
中国科学技术大学电子工程与信息科学系
多媒体通信实验室 （Copyright 2000）
实验二
应用 FFT 对信号进行频谱分析
一、 实验目的
1、在理论学习的基础上，通过本次实验，加深对快速傅里叶变换的理解，熟悉 FFT 算法及 其程序的编写。 2、熟悉应用 FFT 对典型信号进行频谱分析的方法。 3、了解应用 FFT 进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用 FFT。
三、实验内容及步骤
（一）编制实验用主程序及相应子程序 1、在实验之前，认真复习 DFT 和 FFT 有关的知识，阅读本实验原理与方法和实验附录部 分中和本实验有关的子程序，掌握子程序的原理并学习调用方法。 2、编制信号产生子程序及本实验的频掊分析主程序。实验中需要用到的基本信号包括：
（1）高斯序列：
中国科学技术大学电子工程与信息科学系
多媒体通信实验室 （Copyright 2000）
谱，就可以得到相应的连续信号的频谱。注意：这里的信号必须是带限信号，采样也必须满 足 Nyquist 定理。 在各种信号序列中， 有限长序列在数字信号处理中占有很重要的地位。 无限长的序列也 往往可以用有限长序列来逼近。对于有限长的序列我们可以使用离散傅立叶变换（DFT） ， 这一变换可以很好地反应序列的频域特性， 并且容易利用快速算法在计算机上实现当序列的 长度是 N 时，我们定义离散傅立叶变换为：
X ( z ) | z =W − k = ∑ x(n)W Nnk = DFT [ x(n)]
N
N −1 n =0
（2-9）
可以得到 X ( k ) = X ( z ) | z = W
−k N
=e
j
2π k N
, W N 是 z 平面单位圆上幅角为 ω =
−k
2π k 的点， N
就是将单位圆进行 N 等分以后第 k 个点。所以，X(k)是 z 变换在单位圆上的等距采样，或 者说是序列傅立叶变换的等距采样。 时域采样在满足 Nyquist 定理时， 就不会发生频谱混淆； 同样地，在频率域进行采样的时候，只要采样间隔足够小，也不会发生时域序列的混淆。 DFT 是对序列傅立叶变换的等距采样，因此可以用于序列的频谱分析。在运用 DFT 进 行频谱分析的时候可能有三种误差，分析如下： （1）混淆现象 从式（2-6）中可以看出，序列的频谱是采样信号频谱的周期延拓，周期是 2π/T，因此 （这里指的是实信号） 当采样速率不满足 Nyquist 定理， 即采样频率 f s = 1 / T 小于两倍的信号 频率时， 经过采样就会发生频谱混淆。 这导致采样后的信号序列频谱不能真实地反映原信号 的频谱。所以，在利用 DFT 分析连续信号频谱的时候，必须注意这一问题。避免混淆现象 的唯一方法是保证采样的速率足够高，使频谱交叠的现象不出现。这就告诉我们，在确定信 号的采样频率之前，需要对频谱的性质有所了解。在一般的情况下，为了保证高于折叠频率 的分量不会出现，在采样之前，先用低通模拟滤波器对信号进行滤波。 （2）泄漏现象 实际中的信号序列往往很长，甚至是无限长序列。为了方便，我们往往用截短的序列来 近似它们。这样可以使用较短的 DFT 来对信号进行频谱分析。这种截短等价于给原信号序 列乘以一个矩形窗函数。 而矩形窗函数的频谱不是有限带宽的， 从而它和原信号的频谱进行 卷积以后会扩展原信号的频谱。值得一提的是，泄漏是不能和混淆完全分离开的，因为泄露 导致频谱的扩展，从而造成混淆。为了减小泄漏的影响，可以选择适当的窗函数使频谱的扩 散减到最小。 （3）栅栏效应
（4）反三角序列：
如果您在阅读过程中发现疏漏和错误，请您尽快和编者取得联系 network@ cxh@
3
中国科学技术大学电子工程与信息科学系
多媒体通信实验室 （Copyright 2000）
（二）上机实验内容 1、观察高斯序列的时域和频域特性 ①固定信号 x a ( n) 中的参数 p=8，改变 q 的值，使 q 分别等于 2，4，8。观察它们的时 域和幅频特性，了解 q 取不同值的时候，对信号时域特性和幅频特性的影响。②固定 q=8， 改变 p，使 p 分别等于 8，13，14，观察参数 p 变化对信号序列时域及幅频特性的影响。注 意 p 等于多少时， 会发生明显的泄漏现象， 混淆现象是否也随之出现？记录实验中观察到的 现象，绘制相应的时域序列和幅频特性曲线。 2、观察衰减正弦序列的时域和幅频特性 ①令α=0.1 并且 f=0.0625，检查谱峰出现的位置是否正确，注意频谱的形状，绘制幅 频特性曲线。②改变 f=0.4375，再变化 f=0.5625，观察这两种情况下，频谱的形状和谱峰出 现的位置，有无混淆和泄漏现象发生？说明产生现象的原因。 3、观察三角波序列和反三角波序列的时域和幅频特性 观察两者的序列形状和频谱曲线有 ①用 8 点 FFT 分析信号 xc ( n) 和 x d ( n) 的幅频特性， 什么异同？（注意：这时候的 x d ( n) 可以看作是 xc ( n) 经过圆周移位以后得到的）绘制两者 的序列和幅频特性曲线。②在的 xc ( n) 和 x d ( n) 末尾补零，用 16 点 FFT 分析这两个信号的 幅频特性，观察幅频特性发生了什么变化？两个信号之间的 FFT 频谱还有没有相同之处？ 这些变化说明了什么？ （三）* MatLab 上机内容 1、在 MatLab 下重复（二）上机实验内容的所有要求，将 MatLab 的输出结果同自己程序的 输出结果进行比较。 2、将 xb ( n) 信号的长度 N 设为 63，用 MatLab 中 randn(1,N)函数产生一个噪声信号 w(n)， 计算将这个噪声信号叠加到 xb ( n) 上以后新信号 y ( n) = xb ( n) + w( n) 的频谱，观察发生 的变化并记录。 3、在步骤 2 的基础上，改变参数α和 f，观察在出现混淆现象和泄漏现象的时候有噪声的 y(n)信号的频谱有什么变化，是否明显？
如果您在阅读过程中发现疏漏和错误，请您尽快和编者取得联系 network@ cxh@ 2
中国科学技术大学电子工程与信息科学系
多媒体通信实验室 （Copyright 2000）
因为 DFT 是对单位圆上 z 变换的均匀采样，所以它不可能将频谱视为一个连续函数。这样 就产生了栅栏效应，从某种角度来看，用 DFT 来观看频谱就好像通过一个栅栏来观看一幅 景象，只能在离散点上看到真实的频谱。这样的话就会有一些频谱的峰点或谷点被“栅栏” 挡住，不能被我们观察到。减小栅栏效应的一个方法是在源序列的末端补一些零值，从而变 动 DFT 的点数。这种方法的实质是认为地改变了对真实频谱采样的点数和位置，相当于搬 动了“栅栏”的位置，从而使得原来被挡住的一些频谱的峰点或谷点显露出来。注意，这时 候每根谱线多对应的频率和原来的已经不相同了。 从上面的分析过程可以看出，DFT 可以用于信号的频谱分析，但必须注意可能产生的 误差，在应用过程中要尽可能减小和消除这些误差的影响。 快速傅立叶变换 FFT 并不是与 DFT 不相同的另一种变换， 而是为了减少 DFT 运算次数 的一种快速算法。它是对变换式（2-7）进行一次次的分解，使其成为若干小点数 DFT 的组 合，从而减小运算量。常用的 FFT 是以 2 为基数，其长度 N = 2
四、思考题
1、 实验中的信号序列 xc ( n) 和 x d ( n) ，在单位圆上的 z 变换频谱 X c (e jϖ ) 和 X d (e jϖ ) 会 相同吗？如果不同，你能说出哪一个低频分量更多一些吗？为什么？ 2、 对一个有限长序列进行离散傅里叶变换（DFT） ，等价于将该序列周期延拓后进行傅里 叶级数（DFS）展开。因为 DFS 也只是取其中一个周期来运算，所以 FFT 在一定条件 下也可以用以分析周期信号序列。 如果实正弦信号 sin (2πfn ) ， f = 0.1 ， 用 16 点的 FFT 来做 DFS 运算，得到的频谱是信号本身的真实谱吗？