数字信号处理实验参考书new.docVIP

目 录 目录 1 实验一 信号、系统及系统响应 2 实验二 应用FFT对信号进行频谱分析 6 实验三 用双线性变换法设计IIR滤波器 11 实验四 用窗函数设计FIR滤波器 16 附录A C语言实现数字信号处理算法 21 附录B MATLAB的信号表示和处理 32 附录C MATLAB 下的数字信号处理实现示例 56 附录D MATLAB 下的数字滤波器设计示例 64 实验一 信号、系统及系统响应 一．实验目的 1．? 2．? 3．? 4．? 二．实验原理 （一）连续时间信号的采样 采样是从连续信号到离散时间信号的过渡桥梁，对采样过程的研究不仅可以了解采样前后信号时域何频域特性发生的变化以及信号内容不丢失的条件，而且有助于加深对拉氏变换、傅氏变换、Z变换和序列傅氏变换之间关系的理解。 对一个连续时间信号进行理想采样的过程可以表示为信号与一个周期冲激脉冲的乘积，即： （1-1） 其中，是连续信号的理想采样，是周期冲激脉冲 （1-2） 此时，采样信号的拉氏变换可以表示为：,其中是信号的拉氏变换 作为拉氏变换的一种特例，信号理想采样的傅立叶变换为： （1-3） 由（1-3）式可知，信号理想采样后的频谱式原来信号频谱的周期延拓，其延拓周期等于采样频率。根据香农定理，如果原信号是带限信号，且采样频率高于原信号最高频率的2倍，则采样后的离散序列不会发生频谱混叠现象。 在计算机处理时，不采用时计算信号的频谱，而是利用序列的傅立叶变换计算信号的频谱，定义序列： （1-4） 根据Z变换的定义，可以得到序列的Z变换为： （1-5） 以代替上式中的Z，就可以得到序列的傅立叶变换： （1-6） 式（1-3）和式（1-6）具有如下关系： （1-7） 由式可知，在分析一个连续时间信号的频谱时，可以通过取样将有关的计算转化为序列傅立叶变换的计算。 （二）有限长序列的分析 一般来说，在计算机上不可能、也不必要处理连续的曲线，通常，我们的做饭时只要观察、分析在某些频率点上的值。对于长度为N的有限长序列： （1-8） 一般只需要在之间均匀的取M个频率点，计算这些点上的序列傅立叶变换： （1-9） 其中，，。是一个复函数，它的模就是幅频特性曲线。 （三）信号卷积 一个线性时不变离散系统的响应可以用它的单位冲激响应和输入信号的卷积来表示： （1-6） 根据傅立叶变换和 Z变换的性质，应该有： （1-7） （1-8） 上两式告诉我们：可以通过对两个序列的移位、相乘和累加计算信号响应；卷积运算可以在频域通过乘积实现。 三．实验内容及步骤 (一) 编制实验用主程序及相应子程序 1．信号产生子程序，包括： （1）理想采样信号序列：对信号进行理想采样，可以得到一个理想的采样信号序列，其中为幅度因子，是衰减因子，是频率，为采样周期。根据实验内容的需要，这些参量请设计为在程序运行过程中输入。 （2）单位脉冲序列 （3）矩形序列 , 其中取：N=10 2．系统单位冲激响应序列产生子程序，本实验中用到两种FIR系统： （1） （2） 3．有限长序列线性卷积子程序，用于计算两个给定长度（分别为M和N）的序列的卷积，输出序列长度为L=M+N-1。 （二）上机实验内容 在编制以上各部分程序以后，编制主程序调用各个功能模块实现对信号、系统和系统响应的时域和频域分析，完成以下实验内容。 1．分析理想采样信号序列的特性 产生理想采样信号序列，使： （1）首先选用采样频率为1000Hz，T=1/1000，观察所得理想采样信号的幅频特性，在折叠频率以内和给定的理想幅频特性无明显差异，并作记录。 （2）改变采样频率为300Hz，T=1/300，观察所得理想采样信号的幅频特性曲线的变化，并作记录。 （3）进一步减小采样频率为200Hz，T=1/200，观察频谱混叠现象是否明显存在，说明原因，并记录此时的幅频特性曲线。 2．离散信号、系统和系统响应的分析 （1）观察信号和系统的时域和频域特性；利用线性卷积求信号通过系统以后的响应。比较系统响应和系统的时域和幅频特性。注意它们之间有无差异，绘出图形。 （2）观察信号和系统的时域和幅频特性，利用线性卷积求系统响应。判断响应序列图形及序列非零值长度是否与理论结果一致，说出一种定性判断响应序列图形正确与否的方法（提示：）。利用序列的傅立叶变换数值计算子程序求出，观察响应序列的幅频特性。定性判断结果是否正确。改变信号的矩形宽度，使N=5,重复以上动作，观察变化，记录改变参数前后的差异。 （3）将（2）中的信号换为，其中，，，，重复（2）中的实验各步；改变参数再重复（2）中的实验各步；改变参数重复（2）中的实验各步。观察参数的改变对信号与系统响应的时域和幅频特性的影响，绘制相应的图形。 3．卷积定律的验证。利