实验2FIR滤波器设计切比雪夫一致逼近法.docVIP

实时数字信号处理 FIR数字滤波器设计 姓名： 专业： 学号： 日期： 目录 第1章FIR数字滤波器设计原理 TOC \o 1-5 \h \z HYPERLINK \l bookmark1 \o Current Document \h FIR数字滤波器设计方法 3 HYPERLINK \l bookmark2 \o Current Document \h FIR数字滤波器的切比雪夫逼近法设计原理 3 HYPERLINK \l bookmark3 \o Current Document \h 策2章FIR数字滤波器设计 4 HYPERLINK \l bookmark4 \o Current Document \h 2.1语音文件滤波处理的需求分析 4 HYPERLINK \l bookmark5 \o Current Document \h FIR陷波滤波器参数设计 4 第3章FIR数字滤波器滤波实验与结果分析 5 HYPERLINK \l bookmark6 \o Current Document \h 3.1 MATLAB数字滤波实验 5 3.1.1滤波实验流程设计 5 3.1.2语音文件噪声添加 5 3.1.3语音文件滤波处理 6 HYPERLINK \l bookmark7 \o Current Document \h 3.2语音文件滤波实验结果 6 3.2.1滤波实验结果 6 3.2.2滤波实验结果分析 8 FIR数字滤波器设计方法 有限冲击响应数字滤波器(FIR，Finite Impulse Response)是-?种全零点的系 统。一个离散时间系统H(z)=B(z)/A(z)，在分母多项式A(z)的系数全为0，即口丁 变为全为零点的FIR系统。 目前FTR数字滤波器的设计方法主要是建立在对理想滤波器特性做某种近 似的基础上，总共冇三种方法：窗函数法、频率抽样法及最佳一致逼近法。 FIR数字滤波器的切比雪夫逼近法设计原理 最佳一致逼近法，是在所需要的区lXI［a,b］内，使误差函数E(x)=|p(x)-f(x)|较 均匀一致，并且通过合理选择p(x),使E(x)的最大值En达到最小。切比雪夫逼 近理论解决了 pOQ的存在性、唯一性及如何构造的问题。 切比雪夫最佳一致逼近的基本思想是，对于给定区W［a, b］上的连续函数 f(x),在所冇n次多项式的集合屮，寻找一个多项式，使其在［a, b］上对f(x)的偏 差和其他一切属于集合的多项式对f(x)的偏差相比是最小的。 切比雪夫逼近理论指出，这样的多项式是存在的，并且是唯一的。构造这种 最佳一致逼近多项式的方法为交错点组定理。 En=max|p(x)-f(x)| E(x)=p(x)-f(x) p(x)是f(x)最佳一致逼近多项式的充要条件是E(x)在［a，b］上至少存在n+2 个交错点，使得E(x：)=En或者-En，同时要满足E(Xi)=-E(xi+1)。 该n+2个点为交错点组，也为E(x)的极值点。 利用切比雪夫逼近理论设计F1R数字滤波器时，满足最佳一致逼近的滤波器 又具有等纹波性质。存在通带纹波以及阻带纹波。 在实际求解过程中，利用数值分析中的Remez算法，靠一次次迭代来求得 一组交错点组。 第2章FIR数字滤波器设计 2.1语音文件滤波处理的需求分析 语音文件在传输过程中会混入一些噪声，凶此本实验的0的是模拟消除语音 文件噪声。实验过程中的噪声由人为添加，调用滤波程序，观察实验前后的语音 文件频谱。 FIR陷波滤波器参数设计 实验屮人为添加的噪声分为两种：一种为100Hz的低频噪声；另一种为 3000Hz的高频噪声。采样频率为8000Hz。 根据FIR陷波滤波器的设计方法，该滤波器的滤波范围应该包含3个通带以 及2个阻带。两个阻带分别在100Hz与3000Hz处。 由比例换算公式可以得出： 低频阻带：100/8000=0.0125 高频阻带：3000/8000=0.375 因此在0一0.5的范围之pj设置阻带与通带如下： 通带：0.0—0.01 阻带：0.012—0.014 通带：0.016—0.372 阻带：0.374—0.376 通带：0.378—0.5 设置权值： 通带：8 阻带：1 滤波器阶数设置为512阶。 3.1 MATLAB数字滤波实验 3.1.1滤波实验流程设计 (1)读入语音信号，作DFT绘制频谱阁； ⑵加入低频100Hz噪声与高频3000Hz噪声 使用FIR陷波滤波器处理语音文件 处理后的语音文件作DFT绘制频谱阁 观察分析结果 3.1.2语音文件噪声添加 添加低频100Hz与高频3000Hz止弦噪声的程序如下图3.1与3.2所示。 p = var (u) : %功車 snbr = 1/(10^(snb/10)):