吹管音乐滤波去噪:基于汉宁创的fir滤波器

上传人:aa****6 文档编号:38188242 上传时间:2018-04-28 格式:DOC 页数:21 大小:982KB
返回 下载 相关 举报
吹管音乐滤波去噪:基于汉宁创的fir滤波器_第1页
第1页 / 共21页
吹管音乐滤波去噪:基于汉宁创的fir滤波器_第2页
第2页 / 共21页
吹管音乐滤波去噪:基于汉宁创的fir滤波器_第3页
第3页 / 共21页
吹管音乐滤波去噪:基于汉宁创的fir滤波器_第4页
第4页 / 共21页
吹管音乐滤波去噪:基于汉宁创的fir滤波器_第5页
第5页 / 共21页
点击查看更多>>
资源描述

《吹管音乐滤波去噪:基于汉宁创的fir滤波器》由会员分享,可在线阅读,更多相关《吹管音乐滤波去噪:基于汉宁创的fir滤波器(21页珍藏版)》请在金锄头文库上搜索。

1、 吹管音乐滤波去噪-使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 吹管音乐滤波去噪吹管音乐滤波去噪 基于汉宁窗的基于汉宁窗的 FIR 滤波器滤波器学生姓名:学生姓名: 指导老师:指导老师: 摘 要 从网站上下载一段吹管乐器演奏音乐,利用 CE 软件对音乐进行编辑。绘制波形并观察其频谱特点,加入一个带外单频噪声,用汉宁窗设计一个满足指标的 FIR 滤波器,对该含噪音乐信号进行滤波去噪处理,比较滤波前后波形和频谱并进行分析,根据结果和学过的理论得出合理结论。与不同信源相同滤波方法的同学比较各种信源的特点,与相同信源不同滤波方法的同学比较各种滤波方法性能优劣。关键词 滤波去噪;FIR 滤波器

2、;汉宁窗;MATLAB1 引引 言言本课程设计主要针对一段吹管音乐信号,在进行加噪后,利用窗函数设计法选择汉宁窗设计的 FIR 滤波器,对加噪后的吹管音乐信号进行滤波去噪处理,并对前后时域波形和频域波形进行对比分析的程序设计。1.1 课程设计目的课程设计目的本次课设中的主要目的是让学生在熟悉 Matlab 语言环境,掌握其语言编程规则的前提下,利用汉宁窗设计一个符合要求的 FIR 滤波器来实现音乐信号的滤波去噪,并绘制滤波前后的时域波形和频谱图。根据图形分析判断滤波器设计的正确性。通过本次课设,我们能够学会如何综合运用课堂上学会的理论知识,增强自己的动能力与联系实际的能力,为以后的工作奠定基础

3、。1.2 课程设计的要求课程设计的要求(1)滤波器指标必须符合工程实际。(2)设计完后应检查其频率响应曲线是否满足指标。(3)处理结果和分析结论应该一致,而且应符合理论。吹管音乐滤波去噪-使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 (4)独立完成课程设计并按要求编写课程设计报告书。1.3 课程设计平台课程设计平台课程设计的主要设计平台是 MATLAB 7.0。MATLAB 的名称源自 Matrix Laboratory ,它是美国 MathWorks 公司生产的一个为科学和工程计算专门设计的交互式大型软件,是一个可以完成各种精确计算和数据处理的、可视化的、强大的计算工具。它集图示和精

4、确计算于一身,在应用数学、物理、化工、机电工程、医药、金融和其他需要进行复杂数值计算的领域得到广泛应用。它不仅是一个在各类工程设计中便于使用的计算工具,而且也是一个在数学、数值分析和工程计算等课程教学中的优秀的教学工具,在世界各地的高等院校中十分流行,在各类工业应用中更有不俗的表现。MATLAB 可以在几乎所有的 PC 机和大型计算机上运行,适用于 Windows、UNIX 等各种系统平台。1MATLAB 软件包括五大通用功能:数值计算功能(Nemeric);符号运算功能(Symbolic);数据可视化功能(Graphic);数据图形文字统一处理功能(Notebook)和建模仿真可视化功能(S

5、imulink)。其中,符号运算功能的实现是通过请求 MAPLE 内核计算并将结果返回到 MATLAB 命令窗口。该软件有三大特点:一是功能强大;二是界面友善、语言自然;三是开放性强。2MATLAB 在信号与系统中的应用主要包括符号运算和数值计算仿真分析。由于信号与系统课程的许多内容都是基于公式演算,而 MATLAB 借助符号数学工具箱提供的符号运算功能能基本满足信号与系统课程的需求。例如,解微分方程、傅里叶正反变换、拉普拉斯正反变换、z 正反变换等。MATLAB 在信号与系统中的另一主要应用是数值计算与仿真分析,主要包括函数波形绘制、函数运算、冲激响应与阶跃响应仿真分析、信号的时域分析、信号

6、的频谱分析、系统的 S 域分析、零极点图绘制等内容。数值计算仿真分析可以帮助学生更深入理解信号与系统的理论知识,并为将来使用 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)()(NnnjjenheH一般将表示成如下形式:)(jeH(2-2))()()(j gjeHeH式中,是的实函数(可以去负值) 。与前面的表示形式,即)(gH相比, 与不同。与 不同。为了区别于幅频)()()(j gjeHeH)(gH)()(响应函数和相频响应函数,称为幅频特性函数,称为相频特性)(jeH)()(gH)(函数。第一类线

9、性相位 FIR 滤波器的相位特性函数是的严格线性函数:(2-3) 吹管音乐滤波去噪-使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 2.3 窗口设计法窗口设计法窗口设计法是一种通过截断和计权的方法使无限长非因果序列成为有限长脉冲响应序列的设计方法。通常在设计滤波器之前,应该先根据具体的工程应用确定滤波器的技术指标。在大多数实际应用中,数字滤波器常常被用来实现选频操作,所以指标的形式一般为在频域中以分贝值给出的相对幅度响应和相位响应。6窗口设计法基本步骤如下:(1)根据过渡带宽及阻带衰减要求,选择窗函数的类型并估计窗口长度 N。窗函数的类型可根据最小阻带衰减 AS 独立选择。(2)根据待求

10、滤波器的理想频率响应求出理想单位脉冲响应 hd(n)。(3)由性能指标确定窗函数 W(n)和长度 N。(4)求得实际滤波器的单位脉冲响应 h(n),h(n)即为所设计 FIR 滤波器系数向量 b(n)。( )( ).( )dh nhn W n常见的窗函数性能表如下表2-1所示。图2.1 常见窗函数性能表名称滤波器 过渡带宽最小阻带 衰减名称滤波器 过渡带宽最小阻带衰 减矩形1.8/M21dBPARZENWIN6.6/M56dB巴特利特6.1/M25dBFLATTOPWIN19.6/M108dB汉宁6.2/M44dBGAUSSWIN5.8/M60dB汉明6.6/M51dBBARTHANNWIN3

11、.6/M40dB布莱克曼11/M74dBBLACKMANHARRIS16.1/M109dBBOHMANWIN5.8/M51.5dBCHEBWIN15.2/M113dBNUTTALLWIN15.4/M108dBTUKEYWIN2.4/M22dB2.4 汉宁窗汉宁窗(Hanning window)吹管音乐滤波去噪-使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 汉宁窗函数是余弦平方函数,又称之为升余弦函数,它的时域形式可以表为:(2-( )0.5(1 cos(2)1kw kn8)其中 k=1,2,,k。它的频域幅度特性函数为:(2-1()222( )0.5( )0.25()()11NjRRR

12、WWWWeNN9)其中为矩形窗函数的幅度频率特性函数。汉宁窗函数的最大旁瓣值比主瓣值低( )RW31dB,但是主瓣宽度比矩形窗函数的主瓣宽度增加了一倍,为。汉宁窗函数的时域幅8 N度与频域幅度特性曲线的 MATLAB 实现的曲线图如图 2-1 所示。图 2-1 汉宁窗函数的时域幅度与频域幅度特性曲线3 设计步骤设计步骤3.1 设计流程图设计流程图本课程设计主要是从网站上下载一段吹管乐器演奏音乐,利用 CE 软件对音乐进行编辑。绘制波形并观察其频谱特点,加入一个带外单频噪声,用汉宁窗设计一个满足指标的 FIR 滤波器,对该含噪音乐信号进行滤波去噪处理,比较滤波前后波形和频谱并进行吹管音乐滤波去噪

13、-使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 分析,根据结果和学过的理论得出合理结论。程序的设计流程图如下图 3-1 所示。 开始下载一段吹管音乐信号,用 CE 软件编辑格式为 wav.加入单频干扰噪声对吹管音乐信号进行频谱分析,画 出干扰前后的时域和频域波形图利用汉宁窗设计 FIR 滤波器 对吹管音乐信号进行滤波比较滤波前后的时域波形和频谱图,并 回放音乐信号,验证是否达到去噪效果是否达到去 噪效果结束NOYES图 3-1 程序设计流程图 3.2 编辑语音信号编辑语音信号在网上下载一段音乐,再利用 CE 软件将其转换成单声道的.格式文件,再3mpwav将此.格式音乐控制在 10 秒

14、内,以减少设计中的误差。然后在 Matlab 软件平台下,wav利用函数 wavread 对语音信号进行采样,记住采样频率和采样点数。CE 软件操作界面如图 3-2 所示。吹管音乐滤波去噪-使用汉宁窗设计的 FIR 滤波器 第页 共 20 页 图 3-2 CE 软件操作界面3.3 语音加噪处理语音加噪处理采集完成后在信号中加入一个单频噪声,绘制原音乐信号和加噪后的音乐信号的时域和频域的波形图。首先,输入原始音乐信号并播放一次。调用程序如下:x,fs,bits=wavread(h:2013DSPPurpleBambooTune.wav); % 输入参数为文件的全路径和吹管音乐滤波去噪-使用汉宁窗

15、设计的 FIR 滤波器 第页 共 20 页 文件名,输出的第一个参数是每个样本的值,fs 是生成该波形文件时的采样率,bits 是波形文件每样本的编码位数。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); % 对原始信号和加噪信号进行 fft 变换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(时

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 学术论文 > 毕业论文

电脑版 |金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号