数字信号处理实验指导书2016.doc

实验一 信号、系统及系统响应 - 1 - 一 实验目的 - 1 - 二 实验原理 - 1 - 三 实验内容及步骤 - 3 - 四 思考题 - 4 - 五 实验报告要求 - 5 - 实验二 应用FFT对信号进行频谱分析 - 6 - 一 实验目的 - 6 - 二 实验原理与方法 - 6 - 三 实验内容及步骤 - 8 - 四 思考题 - 9 - 五 实验报告要求 - 9 - 实验三 用双线性变换法设计IIR滤波器 - 11 - 一 实验目的 - 11 - 二 实验原理与方法 - 11 - 三 实验内容及步骤 - 14 - 四 思考题 - 15 - 五 实验报告要求 - 15 - 实验四 用窗函数设计FIR滤波器 - 16 - 一 实验目的 - 16 - 二 实验原理和方法 - 16 - 三 实验内容及步骤 - 19 - 四 思考题 - 21 - 五 实验报告要求 - 21 - 附录 MATLAB 下的数字信号处理实现示例 - 22 - 1 信号、系统和系统响应 - 22 - 2 用FFT 进行信号的频谱分析 - 26 - 3 窗函数法设计FIR 滤波器 - 28 - 4 IIR 滤波器的实现 - 31 - 5 窗函数设计FIR滤波器基本示例 - 34 - 实验一 信号、系统及系统响应 一 实验目的 1．熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对采样定理的理解。 2．熟悉离散信号和系统的时域特性。 3．熟悉线性卷积的计算编程方法：利用卷积的方法，观察、分析系统响应的时域特性。 4．掌握序列傅氏变换的计算机实现方法，利用序列的傅氏变换对离散信号、系统及系统响应进行频域分析。 二 实验原理 1．连续时间信号的采样 采样是从连续时间信号到离散时间信号的过渡桥梁，对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生的变化以及信号内容不丢失的条件，而且有助于加深对拉氏变化、傅氏变换、z变换和序列傅氏变换之间关系的理解。 对一个连续时间信号进行理想采样的过程可以表示为该信号的一个周期冲激脉冲的乘积，即 (1-1) 其中是连续信号的理想采样，是周期冲激脉冲 (1-2) 它也可以用傅立叶级数表示为： (1-3) 其中T为采样周期，是采样角频率。设是连续时间信号的双边拉氏变换，即有： (1-4) 此时理想采样信号的拉氏变换为 (1-5) 作为拉氏变换的一种特例，信号理想采样的傅立叶变换 (1-6) 由式(1-5)和式(1-6)可知，信号理想采样后的频谱是原信号频谱的周期延拓，其延拓周期等于采样频率。根据Shannon取样定理，如果原信号是带限信号，且采样频率高于原信号最高频率分量的2倍，则采样以后不会发生频谱混淆现象。 在计算机处理时，不采用式(1-6)计算信号的频谱，而是利用序列的傅立叶变换计算信号的频谱，定义序列，根据Z变换的定义，可以得到序列x(n)的Z变换为： (1-7) 以代替上式中的z，就可以得到序列x(n)的傅立叶变换 (1-8) 式（1-6）和式（1-8）具有如下关系： (1-9) 由式（1-9）可知，在分析一个连续时间信号的频谱时，可以通过取样将有关的计算转化为序列傅立叶变换的计算。 2.有限长序列分析 一般来说，在计算机上不可能，也不必要处理连续的曲线，通常，我们只要观察、分析在某些频率点上的值。对于长度为N的有限长序列 (1-10) 一般只需要在之间均匀地取M个频率点，计算这些点上的序列傅立叶变换 (1-11) 其中，是一个复函数，它的模就是幅频特性曲线。 3．信号卷积 一个线性时不变离散系统的响应y(n)可以用它的单位冲激响应h(n)和输入信号x(n)的卷积来表示： (1-12) 根据傅立叶变换和Z变换的性质，与式（1-12）对应应该有 (1-13) (1-14) 式（1-12）告诉我们可以通过对两个序列的移位、相乘、累加计算信号响应；而式（1-14）告诉我们卷积运算也可以在频域上用乘积实现。 三 实验内容及步骤 1、分析理想采样信号序列的特性： 产生理想采样信号序列，对信号进行理想采样，可以得到一个理想的采样信号序列：，其中A为幅度因子，是衰减因子，是频率。T为采样周期。根据实验内容的需要，这些参量请设计为在程序运行过程中输入。 首先使。 （1）先选用采样频率为1000Hz，T=1/1000，观察所得理想采样信号的幅频特性，在折叠频率以内和给定的理想幅频特性无明显差异，并做记录； （2）改变采样频率为300Hz，T=1/300，观察所得到的频谱特性曲线的变化，并做记录； （3）进一步减小采样频率为200Hz，T=1/200，观察频谱“混淆”现象是否明显存在，说