学习matlab滤波器设计.pdfVIP

一、课程设计目的： 1. 进一步理解数字信号处理的基本概念、基本理论和基本方法； 2. 熟悉在 Windows环境下语音信号采集的方法； 3. 学会用 MATLAB软件对信号进行分析和处理； 4. 综合运用数字信号处理理论知识，掌握用 MATLAB软件设计 FIR 和 IIR 数 字滤波器的方法； 5. 提高依据所学知识及查阅的课外资料来分析问题解决问题的能力。 二、课程设计内容： 1. 语音信号的采集 利用 windows 下的录音机录制一段自己的话音， 时间控制在 1 秒左右； 并对语音信号进 行采样，理解采样频率、采样位数等概念。 2. 语音信号的频谱分析 利用函数 fft 对采样后语音信号进行快速傅里叶变换，得到信号的频谱特性。 3. 设计数字滤波器 采用窗函数法和双线性变换法设计下列要求的三种滤波器，根据语音信号的特点给出 有关滤波器的性能指标： 1）低通滤波器性能指标， fp=1000Hz ，fc=1200Hz ， As=100dB ，Ap=1dB； 2 ）高通滤波器性能指标， fc=4800Hz ，fp=5000Hz ，As=100dB，Ap=1dB； 3）带通滤波器性能指标， fp 1=1200Hz，fp 2=3000Hz，fc 1=1000Hz，fc 2 =3200Hz，As=100dB， Ap=1dB。 4. 对语音信号进行滤波 比较用两种方法设计的各滤波器的性能，然后用性能好的滤波器分别对采集的语音信 号进行滤波；并比较滤波前后语音信号的波形及频谱，分析信号的变化。 5. 回放语音信号，感觉滤波前后的声音变化。 三、实验原理 （一）基于双线性 Z 变换法的 IIR 数字滤波器设计 由于 的频率映射关系是根据 推导的，所以使 j Ω 轴每隔 2 π/Ts 便映射到 单位圆上一周， 利用冲激响应不变法设计数字滤波器时可能会导致上述的频域混叠现象。 为 了克服这一问题，需要找到由 s 平面到 z 平面的另外的映射关系，这种关系应保证： 1) s 平面的整个 j Ω 轴仅映射为 z 平面单位圆上的一周； 2) 若 G(s) 是稳定的，由 G(s) 映射得到的 H(z) 也应该是稳定的； 3) 这种映射是可逆的，既能由 G(s) 得到 H(z) ，也能由 H(z) 得到 G(s) ； 4) 如果 G(j0)=1 ，那么 。 双线性 Z 变换满足以上 4 个条件的映射关系，其变换公式为 双线性 Z 变换的基本思路是： 首先将整个 s 平面压缩到一条从– π/Ts 到 π/Ts 的带宽 为 2 π/Ts 的横带里， 然后通过标准的变换关系 将横带变换成整个 z 平面上去， 这样 就得到 s 平面与 z 平面间的一一对应的单值关系。 在 MATLAB中，双线性 Z 变换可以通过 bilinear 函数实现，其调用格式为： [Bz ，Az] ＝bilinear(B ，A，Fs) ； 其中 B，A 为模拟滤波器传递函数 G(s) 的分子分母多项式的系数向量，而 Bz，Az 为数 字滤波器的传递函数 H(z) 的分子分母多项式的系数向量。 （二）窗函数法设计 FIR 数字滤波器 窗函数法就是设计 FIR 数字滤波器的最简单的方法。 它在设计 FIR 数字滤波器中有很重 要的作用， 正确地选择窗函数可以提高设计数字滤波器的性能， 或者在满足设计要求的情况 下，减小 FIR 数字滤波器的阶次。 常用的窗函数有以下几种： 矩形窗 (Rect