基于MATLAB的FIR数字滤波器设计总结.docVIP

数字信号处理课程设计报告 课设名称：基于MATLAB的FIR数字滤波器设计 学 院：信息工程 专 业：通信工程 班 级：2012159 学 号：201215910 姓 名：张孟承 辅导老师：李亚萍 陈黎霞 时 间：2015.12.21——2015.12.30基于MATLAB的FIR数字滤波器设计 摘要： 论文针对窗函数法、频率采样法和最优化设计法三种设计方法，采用MATLAB 进行FIR滤波器的设计与仿真，并比较了三种方法的特点。结果表明，在同样的设计指标下，利用等波纹切比雪夫逼近法则的设计可以获得最佳的频率特性和衰耗特性，具有通带和阻带平坦，过渡带窄等优点。 关键词：FIR滤波器；MATLAB;频率采样法； 一、设计目的 本次课程设计的目的是利用MATLAB对语音信号进行数字信号处理和分析，要求学生对语音信号进行频谱分析，并对采集语音信号加入干扰，对加入信号后的频谱进行频谱分析，设计合适的滤波器滤除噪声，恢复原信号，学会GUI界面设 计方法，掌握基本操作原理。 二．设计要求 1.采集一段语音信号，并对其频谱进行分析。 2.设计滤波器，可以选用高通，低通或带通滤波器，滤除噪声，分析滤波后信号的时域或频域特征。 3.播放三段语音对比频谱，分析其异同。 三．FIR 数字滤波器设计 FIR滤波器设计的任务是选择有限长度的，使传输函数满足一定的幅度特性和线性相位要求。由于FIR 滤波器很容易实现严格的线性相位，所以FIR 数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。 设计过程一般包括以下三个基本问题： （1） 根据实际要求确定数字滤波器性能指标； （2） 用一个因果稳定的系统函数去逼近这个理想性能指标； （3） 用一个有限精度的运算去实现这个传输函数。 频率采样法： 频率采样法是从频域出发,根据频域采样定理，对给定的理想滤波器的频率响应加以等间隔的抽样 ，得到。 而在各采样点间的频率响应则是其的加权内插函数延伸叠加的结果。但对于一个无限长的序列，用频率采样法必然有一定的逼近误差,误差的大小取决于理想频响曲线的形状, 理想频响特性变换越平缓, 则内插函数值越接近理想值,误差越小。为了提高逼近的质量，可以通过在频率相应的过渡带内插入比较连续的采样点，扩展过渡带使其比较连续，从而使得通带和阻带之间变换比较缓慢，以达到减少逼近误差的目的。 选取w∈[0,2π]内N个采样点的约束条件为： （3-3） （1）增大阻带衰减三种方法： 1）加宽过渡带宽，以牺牲过渡带换取阻带衰减的增加。 2）过渡带的优化设计 利用线性最优化的方法确定过渡带采样点的值，得到要求的滤波器的最佳逼近（而不是盲目地设定一个过渡带值）。 3）增大N。如果要进一步增加阻带衰减，但又不增加过渡带宽，可增加采样点数N。代价是滤波器阶数增加，运算量增加。 直接从频域进行设计，物理概念清楚，直观方便；适合于窄带滤波器设计，这时频率响应只有少数几个非零值，但是截止频率难以控制。 典型应用：用一串窄带滤波器组成多卜勒雷达接收机，覆盖不同的频段，多卜勒频偏可反映被测目标的运动速度。 四．GUI设计效果图 图4.GUI效果初始界面图 五．课程设计的具体实现 1、 信号的采集 要求录制一段语音信号，并对录制的信号进行采样。在MATLAB中可以使用wavread（）函数完成。 2、 语音信号的读入与打开 在MATLAB中wavread();用于读取语音，采样值放在向量y中，fs表示采样频率(Hz)，bits表示采样位数。sound()；用于语音的播放，用于对声音的回放。 3、 语音信号的FFT变换 FFT即为快速傅氏变换，是离散傅氏变换的快速算法，在MATLAB的信号处理工具箱中函数FFT用于快速傅立叶变换。其调用格式为y=fft(x)，其中，x是序列，y是序列的FFT，x可以为一向量或矩阵，若x为一向量，y是x的FFT且和x相同长度；若x为一矩阵，则y是对矩阵的每一列向量进行FFT。函数FFT的另一种调用格式为y=fft(x,N)，式中，x，y意义同前，N为正整数。函数执行N点的FFT，若x为向量且长度小于N，则函数将x补零至长度N；若向量x的长度大于N，则函数截短x使之长度为N；若x 为矩阵，按相同方法对x进行处理。对语音信号进行FFT变换就是为了得到它的频域的图形，便于从图中观察出信号的幅度等特性。 4、含噪信号的合成 在MATL