数字信号处理综合实验..docVIP

数字信号处理综合实验 一、实验目的 1．掌握MATLAB的程序设计方法； 2．掌握数字信号处理的基本理论和基本方法； 3．掌握语音信号的采集与处理方法； 4．掌握用MATLAB设计FIR和IIR数字滤波器的方法； 5．掌握用MATLAB对信号进行分析和处理的方法． 二、实验原理 1.采样定理：在进行模拟/数字信号的转换过程中，当采样频率fs.max大于信号中最高频率fmax的2倍时，即：fs.max=2fmax,则采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定理又称奈奎斯特定理。 2.数字滤波器设计的基本原理：数字滤波器是由常数乘法器、加法器和延迟单元组成的一种算法或装置，是应用最广泛的线性时不变系统。它对输入序列进行处理，去除其频谱中不需要的成分，保留有效成分。 （1）对于IIR滤波器，其逼近问题就是寻找滤波器的各项系数，使其系统函数逼近一个所要求的特性。先设计一个合适的模拟滤波器，然后变换成满足约定指标的数字滤波器。 用双线形变换法设计IIR数字滤波器的过程： a.将设计性能指标中的关键频率点进行“预畸”? b.利用“预畸”得到的频率点设计一个模拟滤波器。? c.双线形变换，确定系统函数 FIR滤波器的优点是有N个零点和N个极点，但其全部极点都在z平面的原点，FIR系统总是稳定的；在满足一定的对称条件下，可以实现严格的线性相位。缺点是由于没有非零极点，要获得较好的过渡带特性，必须以较高的滤波器阶数作为代价。 IIR滤波器设计可以通过模拟滤波器设计和各种变换法实现，而FIR滤波器不能利用模拟滤波器的设计理论，因为FIR滤波器的系统函数是的多项式，并非有理分式，无法找到与其对应的原型模拟滤波器。 3.频谱现象分析： （1）混叠：采样序列的频谱是被采样信号频谱的周期延拓，当采样频率不满足奈奎斯特采样定理的时候，就会发生混叠，使得刺痒后的序列信号的频谱不能真实的反映原采样信号的频谱。 （2）泄露：根据理论分析，一个时间的信号其频带宽度为无限，一个时间无限的信号其频带宽度则为有限。因此对一个时间有限的信号，应用DFT进行分析，频谱混叠难以避免。对一个时间无限的信号虽然频带有限，但在实际运算中，时间总是取有限值，在将信号截断的过程中，出现了分散的扩展谱线的现象，称之为频谱泄露或功率泄露。 （3）栅栏效应：DFT是对单位圆上Z变换的均匀采样，所以它不可能将频谱视为一个连续函数，就在一定意义上看，用DFT来观察频谱就好象通过一个栅栏来观看一个景象一样，只能在离散点上看到真实的频谱，这样就有可能发生一些频谱的峰点和谷点被“尖桩的栅栏”所挡住，不能被我们观察到。 三、实验内容 1、连续信号的频率分析 一个连续信号含两个频率分量：，用Hz的采样器对该信号进行采样，对采样后的序列分别用长度为16和128的FFT观察其频谱，你观察到什么现象？请解释原因。 2、数字滤波器的设计 分别用双线性变换法和窗口设计法设计IIR和FIR数字滤波器，包括低通、高通和带通三种类型。fs=8000HZ; 三种滤波器的性能指标如下： （1）低通滤波器：通带和阻带边界频率分别为1000Hz和1200Hz，通带波动为1dB，阻带最小衰减为50dB； （2）高通滤波器：通带和阻带边界频率分别为3000Hz和2800Hz，通带波动为1dB，阻带最小衰减为50dB； （3）带通滤波器：通带边界频率分别为1200Hz和3000Hz，阻带边界频率分别为1000Hz和3200Hz，通带波动为1dB，阻带最小衰减为50dB。 请给出每种滤波器的阶数、系统函数和对数幅频响应。 3、语音信号的采集与处理 （1）用Windows自带的录音机或其他录音软件录制一段语音。录制格式如下：采样频率8000Hz，8位，单声道。录音内容为你的姓名和学号，时长控制在5秒左右，说话不要有过多的停顿。最后格式为.wav; （2）对语音信号进行频谱分析，请画出时域波形和对数幅频特性。 （3）分别用“2、数字滤波器的设计”中的三种FIR滤波器对语音信号进行滤波，比较滤波前后语音信号的时域波形与幅频特性。并用MATLAB的sound函数回放语音，试听滤波前后语音的变化。请在实验结果与分析中解释实验现象。 （4）分别用“2、数字滤波器的设计”中的IIR低通滤波器和线性相位FIR低通滤波器对语音信号进行滤波，比较两种滤波器滤波后的时域波形与幅频特性。并用MATLAB的sound函数回放语音，试听两种滤波器滤波后语音的差异。请在实验结果与分析中解释实验现象。 四、实验结果与分析 1.实验结果： 实验结果分析：采样后序列：，N1=8; N2=64；所以N=64;两序列有效长度一样，分辨率一样，由于FFT点数不同，导致结果不同。N=166