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

数字信号处理综合实验 一、实验目的 1．掌握MATLAB的程序设计方法； 2．掌握数字信号处理的基本理论和基本方法； 3．掌握语音信号的采集与处理方法； 4．掌握用MATLAB设计FIR和IIRMATLAB对信号进行分析和处理的方法． 二、实验原理 1.采样定理：在进行模拟/数字信号的转换过程中，当采样频率fs.max大于信号中最高频率fmax的2倍时，即：fs.max=2fmax,则采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定理又称奈奎斯特定理。 2.数字滤波器设计的基本原理：数字滤波器是由常数乘法器、加法器和延迟单元组成的一种算法或装置，是应用最广泛的线性时不变系统。它对输入序列进行处理，去除其频谱中不需要的成分，保留有效成分。 （1）对于IIR滤波器，其逼近问题就是寻找滤波器的各项系数，使其系统函数逼近一个所要求的特性。先设计一个合适的模拟滤波器，然后变换成满足约定指标的数字滤波器。用双线形变换法设计IIR数字滤波器的过程：a.将设计性能指标中的关键频率点进行“预畸”b.利用“预畸”得到的频率点设计一个模拟滤波器。c.双线形变换，确定系统函数 含两个频率分量：，用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；=64;两序列有效长度一样，分辨率一样，由于FFT点数不同，导致结果不同。N=1664时，对序列进行截断，产生频谱泄露，当N=128=2N64时，无频谱泄漏。 实验结果： (1) az=[1.0000 -8.6797 35.5999 -90.8275 159.8643 -203.5543 191.1691 -132.3641 66.2101 -22.7930 4.8632 -0.4877] bz=[5.2944e-008 5.8238e-007 2.9119e-006 8.7357e-006 1.7471e-005 2.446e-005 2.446e-005 1.7471e-005 8.7357e-006 2.9119e-006 5.8238e-007 5.2944e-008] 可得IIR滤波器的阶数为12阶。 通带和阻带的数字边界频率: 中心频率为： 选择窗函数和。阻带衰减不小于50dB，因此选择海明窗()。滤波器的过渡带宽为，因此窗口长度为: 线性相位延迟常数为: 根据理想边界频率和线性相位延迟常数，求理想单位脉冲响应: 窗函数与理想单位脉冲响应相乘，得到线性相位FIR低通滤波器的单位脉冲响应: 可得FIR滤波器的阶数为132阶 (2) az=[1.0000 8.6797 35.5999 90.8275 159.8643 203.5543 191.1691 132.3641 66.2101 22.7930 4.8632 0.4877] bz=[5.2944e-00