1、 吹管音乐滤波去噪 -使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 吹管音乐滤波去噪基于汉宁窗的 FIR 滤波器学生姓名： 指导老师： 摘 要 从网站上下载一段吹管乐器演奏音乐，利用 CE 软件对音乐进行编辑。绘制波形并观察其频谱特点，加入一个带外单频噪声，用汉宁窗设计一个满足指标的 FIR 滤波器，对该含噪音乐信号进行滤波去噪处理，比较滤波前后波形和频谱并进行分析，根据结果和学过的理论得出合理结论。与不同信源相同滤波方法的同学比较各种信源的特点，与相同信源不同滤波方法的同学比较各种滤波方法性能优劣。关键词 滤波去噪；FIR 滤波器；汉宁窗；MATLAB1 引 言本课程设计主要针对一段
2、吹管音乐信号，在进行加噪后，利用窗函数设计法选择汉宁窗设计的 FIR 滤波器，对加噪后的吹管音乐信号进行滤波去噪处理，并对前后时域波形和频域波形进行对比分析的程序设计。1.1 课程设计目的本次课设中的主要目的是让学生在熟悉 Matlab 语言环境，掌握其语言编程规则的前提下，利用汉宁窗设计一个符合要求的 FIR 滤波器来实现音乐信号的滤波去噪，并绘制滤波前后的时域波形和频谱图。根据图形分析判断滤波器设计的正确性。通过本次课设，我们能够学会如何综合运用课堂上学会的理论知识，增强自己的动能力与联系实际的能力，为以后的工作奠定基础。1.2 课程设计的要求（1）滤波器指标必须符合工程实际。（2）设计完
3、后应检查其频率响应曲线是否满足指标。（3）处理结果和分析结论应该一致，而且应符合理论。吹管音乐滤波去噪 -使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 （4）独立完成课程设计并按要求编写课程设计报告书。1.3 课程设计平台课程设计的主要设计平台是 MATLAB 7.0。MATLAB 的名称源自 Matrix Laboratory ，它是美国 MathWorks 公司生产的一个为科学和工程计算专门设计的交互式大型软件，是一个可以完成各种精确计算和数据处理的、可视化的、强大的计算工具。它集图示和精确计算于一身，在应用数学、物理、化工、机电工程、医药、金融和其他需要进行复杂数值计算的领域得到
4、广泛应用。它不仅是一个在各类工程设计中便于使用的计算工具，而且也是一个在数学、数值分析和工程计算等课程教学中的优秀的教学工具，在世界各地的高等院校中十分流行，在各类工业应用中更有不俗的表现。MATLAB 可以在几乎所有的 PC 机和大型计算机上运行，适用于 Windows、UNIX 等各种系统平台。 1MATLAB 软件包括五大通用功能：数值计算功能（ Nemeric）；符号运算功能（Symbolic）；数据可视化功能（Graphic）；数据图形文字统一处理功能（Notebook ）和建模仿真可视化功能（Simulink）。其中，符号运算功能的实现是通过请求 MAPLE 内核计算并将结果返回到
5、 MATLAB 命令窗口。该软件有三大特点：一是功能强大；二是界面友善、语言自然；三是开放性强。 2MATLAB 在信号与系统中的应用主要包括符号运算和数值计算仿真分析。由于信号与系统课程的许多内容都是基于公式演算，而 MATLAB 借助符号数学工具箱提供的符号运算功能能基本满足信号与系统课程的需求。例如，解微分方程、傅里叶正反变换、拉普拉斯正反变换、z 正反变换等。 MATLAB 在信号与系统中的另一主要应用是数值计算与仿真分析，主要包括函数波形绘制、函数运算、冲激响应与阶跃响应仿真分析、信号的时域分析、信号的频谱分析、系统的 S 域分析、零极点图绘制等内容。数值计算仿真分析可以帮助学生更深
6、入理解信号与系统的理论知识，并为将来使用 MATLAB 进行信号处理领域的各种分析和实际应用打下基础 3。2 设计原理2.1 数字信号处理数字信号处理（Digital Signal Processing，简称 DSP）是一门涉及许多学科而又广泛吹管音乐滤波去噪 -使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 应用于许多领域的新兴学科 4。20 世纪 60 年代以来，随着计算机和信息技术的飞速发展，数字信号处理技术应运而生并得到迅速的发展。在过去的二十多年时间里，数字信号处理已经在通信等领域得到极为广泛的应用。数字信号处理是利用计算机或专用处理设备，以数字形式对信号进行采集、变换、滤波、
7、估值、增强、压缩、识别等处理，以得到符合人们需要的信号形式。 52.2 FIR 滤波器有限长单位脉冲响应数字滤波器（Finite Impulse Response Digital Filter,缩写 FIRDF）简称 FIR 滤波器，是数字信号处理系统中最基本的原件，其最大优点是可以实现线性相位滤波，可以在保证任意幅频特性的同时具有严格的线性相频特性，满足了在数字通信和图像传输与处理等应用场合对线性相位的要求。FIR 滤波器是全零点滤波器，硬件和软件实现结构简单，因而是十分稳定的系统。 6FIR 滤波器的设计方法主要分为两类：第一类是基于逼近理想滤波器器特性的方法包括窗函数法、频率采样法、和等
8、波纹最佳逼近法；第二类是最优设计法。本次课设采用的是第一类设计法中的窗函数法。设 FIR 滤波器的单位脉冲响应 的长度为 ，则其频率响应函数为)(nhN（2-1）10(Nnnjj eeH一般将 表示成如下形式：)(je（2-2）)()(jgj ee式中， 是 的实函数（可以去负值） 。与前面的表示形式，即)(gH相比， 与 不同。 与 不同。为了区别于幅频)(jj ee)(gH)()(响应函数 和相频响应函数 ，称 为幅频特性函数，称 为相频特性)(j g )(函数。第一类线性相位 FIR 滤波器的相位特性函数是 的严格线性函数：（2-3）吹管音乐滤波去噪 -使用汉宁窗设计的 FIR 滤波器
9、第页 共 20 页 2.3 窗口设计法窗口设计法是一种通过截断和计权的方法使无限长非因果序列成为有限长脉冲响应序列的设计方法。通常在设计滤波器之前，应该先根据具体的工程应用确定滤波器的技术指标。在大多数实际应用中，数字滤波器常常被用来实现选频操作，所以指标的形式一般为在频域中以分贝值给出的相对幅度响应和相位响应。 6窗口设计法基本步骤如下：（1）根据过渡带宽及阻带衰减要求，选择窗函数的类型并估计窗口长度 N。窗函数的类型可根据最小阻带衰减 AS 独立选择。（2）根据待求滤波器的理想频率响应求出理想单位脉冲响应 hd(n)。（3）由性能指标确定窗函数 W(n)和长度 N。（4）求得实际滤波器的单
10、位脉冲响应 h(n)，h(n) 即为所设计 FIR 滤波器系数向量 b(n)。()().dhnWn常见的窗函数性能表如下表2-1所示。图2.1 常见窗函数性能表名称 滤波器过渡带宽 最小阻带 衰减 名称 滤波器过渡带宽 最小阻带衰 减矩形 1.8/M 21dB PARZENWIN 6.6/M 56dB巴特利特 6.1/M 25dB FLATTOPWIN 19.6/M 108dB汉宁 6.2/M 44dB GAUSSWIN 5.8/M 60dB汉明 6.6/M 51dB BARTHANNWIN 3.6/M 40dB布莱克曼 11/M 74dB BLACKMANHARRIS 16.1/M 109d
11、BBOHMANWIN 5.8/M 51.5dB CHEBWIN 15.2/M 113dBNUTTALLWIN 15.4/M 108dB TUKEYWIN 2.4/M 22dB2.4 汉宁窗（Hanning window）吹管音乐滤波去噪 -使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 汉宁窗函数是余弦平方函数，又称之为升余弦函数，它的时域形式可以表为：（2-()0.5(1cos2)1kwkn8）其中 k=1,2,，k。它的频域幅度特性函数为：（2-1()22()0.5().()()1NjRRRWWeN 9）其中 为矩形窗函数的幅度频率特性函数。汉宁窗函数的最大旁瓣值比主瓣值低()R31
12、dB，但是主瓣宽度比矩形窗函数的主瓣宽度增加了一倍，为 。汉宁窗函数的时域幅8N度与频域幅度特性曲线的 MATLAB 实现的曲线图如图 2-1 所示。图 2-1 汉宁窗函数的时域幅度与频域幅度特性曲线3 设计步骤3.1 设计流程图本课程设计主要是从网站上下载一段吹管乐器演奏音乐，利用 CE 软件对音乐进行编辑。绘制波形并观察其频谱特点，加入一个带外单频噪声，用汉宁窗设计一个满足指标的 FIR 滤波器，对该含噪音乐信号进行滤波去噪处理，比较滤波前后波形和频谱并进行吹管音乐滤波去噪 -使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 分析，根据结果和学过的理论得出合理结论。程序的设计流程图如下
13、图 3-1 所示。 开始下载一段吹管音乐信号，用CE 软件编辑格式为 wav.加入单频干扰噪声对吹管音乐信号进行频谱分析，画出干扰前后的时域和频域波形图利用汉宁窗设计 FIR 滤波器对吹管音乐信号进行滤波比较滤波前后的时域波形和频谱图，并回放音乐信号，验证是否达到去噪效果是否达到去噪效果结束NOYES图 3-1 程序设计流程图 3.2 编辑语音信号在网上下载一段 音乐，再利用 CE 软件将其转换成单声道的 . 格式文件，再3mp wav将此. 格式音乐控制在 10 秒内，以减少设计中的误差。然后在 Matlab 软件平台下，wav利用函数 wavread 对语音信号进行采样，记住采样频率和采样
14、点数。CE 软件操作界面如图 3-2 所示。吹管音乐滤波去噪 -使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 图 3-2 CE 软件操作界面3.3 语音加噪处理采集完成后在信号中加入一个单频噪声，绘制原音乐信号和加噪后的音乐信号的时域和频域的波形图。首先，输入原始音乐信号并播放一次。调用程序如下：x,fs,bits=wavread(h:2013DSPPurpleBambooTune.wav); % 输入参数为文件的全路径和吹管音乐滤波去噪 -使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 文件名，输出的第一个参数是每个样本的值，fs 是生成该波形文件时的采样率， bits 是波形
15、文件每样本的编码位数。sound(x,fs,bits); % 按指定的采样率和每样本编码位数回放计算信号长度并加入噪声。调用程序如下：N=length(x); % 计算信号 x 的长度fn=2100; % 单频噪声频率t=0:1/fs:(N-1)/fs; % 计算时间范围，样本数除以采样频率x=x(:,1);y=x+sin(fn*2*pi*t); % 加入一个单频噪声sound(y,fs,bits); % 可以明显听出有尖锐的单频啸叫声绘制原始音乐信号和加入噪声后的音乐信号的时域和频谱波形图。调用程序如下：X=abs(fft(x); Y=abs(fft(y); % 对原始信号和加噪信号进行 f
16、ft 变换X=X(1:N/2); Y=Y(1:N/2); % 截取前半部分deltaf=fs/ N; % 计算频谱的谱线间隔f=0:deltaf:fs/2-deltaf; % 计算频谱频率范围figure(1);subplot(2,2,1);plot(t,x); % 布局为 2*2 的四个小图title(原始音乐信号);xlabel(时间(t);ylabel(幅度); %改变横纵坐标的范围axis(0,2,-1.5,1.5); %加上标题和横坐标名称grid on; % 加上网格subplot(2,2,2);plot(f,X);title(原始音乐信号频谱);xlabel(频率(f);ylab
17、el(幅度谱);axis(0,3000,0,3000);grid on; subplot(2,2,3);plot(t,y);title(加入干扰后的音乐信号);xlabel( 时间(t);ylabel(幅度);axis(0,2,-1.5,1.5); grid on;subplot(2,2,4);plot(f,Y);title(加入单频干扰后的音乐信号频谱);xlabel( 频率 (f);ylabel(幅度谱);axis(0,3000,0,3000);grid on;用绘图命令分别画出加噪前后的时域和频域波形，如下图 3-3 所示。吹管音乐滤波去噪 -使用汉宁窗设计的 FIR 滤波器 第页 共
18、20 页 图 3-3 吹管音乐信号加入单频噪声前后的时域与频谱波形图由上图可以看到，语音信号加入单频噪声后的时域波形比未加之前在幅度范围内有了明显的增加，在频谱方面可以看到除了在加了噪声之后的频谱图上的 2100Hz 出现一个明显的冲激信号外，其它地方均与未加时的原始吹管音乐信号频谱相同，这一现象表现在音乐播放时，可以听见一声尖锐的啸叫声。3.4 滤波器设计本次课程设计中主要应用汉宁窗设计出 FIR 滤波器。利用 Matlab 中的函数 freqz 画出各滤波器的频率响应，首先利用数字信号处理里面学过的知识，根据选定的参数，用汉宁窗函数法设计 FIR 数字滤波器，得到数字滤波器的参数 b，a。
19、其中 b 为系统函数的分子系数，a 为系统函数分母系数。再调用 freqz(b,a,512,fs)即可得到该滤波器的频率响应。主程序如下：fpd=1450;fsd=1650;fsu=2250;fpu=2350; % FIR 滤波器的上下截止频率Rp=1;As=37; % 带阻滤波器设计指标吹管音乐滤波去噪 -使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 fcd=(fpd+fsd)/2;fcu=(fpu+fsu)/2;df=min(fsd-fpd),(fpu-fsu); % 计算上下边带中心频率，和频率间隔wcd=fcd/fs*2*pi;wcu=fcu/fs*2*pi; dw=df/fs
20、*2*pi; % 将 Hz 为单位的模拟频率换算为 rad为单位的数字频率wsd=fsd/fs*2*pi; wsu=fsu/fs*2*pi; M=ceil(6.2*pi/dw)+1; % 计算汉宁窗设计该滤波器时需要的阶数n=0:M-1; % 定义时间范围w_ham=hanning(M); % 产生 M 阶的汉宁窗 hd_bs=ideal_lp(wcd,M)+ideal_lp(pi,M)-ideal_lp(wcu,M); % 调用自编函数计算理想带阻滤波器的脉冲响应h_bs=w_ham.*hd_bs; % 用窗口法计算实际滤波器脉冲响应db,mag,pha,grd,w=freqz_m(h_bs,1); % 调用自编函数计算滤波器的频率特性通过绘图工具可得出滤波器的波形图，如图 3-4 所示。图 3-4 FIR 滤波器的频率响应
