弹拨音乐信号滤波去噪—三角形窗设计的FIR滤波器2.docVIP

学生姓名：夏霞 指导老师： 摘 要 本课程设计主要内容是通过使用三角形窗设计一个FIR滤波器以对弹拨音乐信号进行滤波去噪处理。本设计首先通过麦克风采集一段弹拨音乐信号，依据对该信号的频谱分析，给定相关指标。MATLAB7.0和TI C6713为平台，M语言编程，采用三角形窗设计满足指标的FIR滤波器，以该音乐信号进行滤波去噪处理。通过对比滤波前后的波形图，，该设计基本满足设计要求。关键词 课程设计；滤波去噪；FIR滤波器；三角窗函数；MATLAB 1 引 言 。 1.2课程设计的要求 （1）算法结果和分析结论应该一致，而且应符合理论（）处理结果和分析结论应该一致，而且应符合理论。 （）独立完成课程设计并按要求编写课程设计报告书。 1.3设计平台 2 设计原理2.1 滤波器原理 FIR(Finite Impulse Response)滤波器：有限长单位冲激响应滤波器，是数字信号处理系统中最基本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。 有限长单位冲激响应（FIR）滤波器有以下特点： (1) 系统的单位冲激响应h(n)在有限个n值处不为零； (2) 系统函数H(z)在|z|0处收敛，极点全部在z=0处（因果系统）； (3) 结构上主要是非递归结构，没有输出到输入的反馈，但有些结构中（例如频率抽样结构）也包含有反馈的递归部分。 FIR滤波器的系统函数用下式表示：。H(n)就是FIR滤波器的单位脉冲响应。FIR滤波器最重要的优点就是由于不存在系统极点，FIR滤波器是绝对稳定的系统。 2.2窗口设计法 窗函数设计法的基本思想是用FIRDF逼近洗完的滤波特性。设希望逼近的滤波器的频率响应函数为，其单位脉冲响应为表示。为了设计简单方便，通常选择为具有片段常数特性的理想滤波器。因此是无限长非因果序列，不能直接作为FIRDF的单位脉冲响应。窗函数设计法就是截取为有限长的一段因果序列，并用合适的窗函数进行加权作为FIRDF的单位脉冲响应h(n)。 用窗函数法设计FIRDF的具体设计步骤如下： （1）构造希望逼近的频率响应函数。 （2）求出。 （3）加窗得到FIRDF的单位脉冲响应h(n)。h(n)=w(n)。式中，w(n)称为窗函数，其长度为N。如果要求设计第一类线性相位FIRDF，则要求h(n)关于(N-1)/2点偶对称。而关于n=t点偶对称，所以要求t=(N-1)/2。同时要求w(n)关于(N-1)/2点偶对称。 2.3三角形窗 三角形窗：利用w=triang(n)的形式得到窗函数，其中n为窗函数的长度，而返回值w为一个n阶的向量，它的元素由窗函数的值组成，w=triang(N-2)等价于bartlett（N）。 三角形窗的旁瓣峰值为25db，主瓣峰值为6.1pi/M,其最小阻带衰减是25db。 2.4 硬件指标 TMS320C6713 是TI 公司推出的一款TMS320C6000系列的浮点DSP 芯片。其片内有8 个并行的处理单元,单字节字长为32 位,从而每周期可以执行8 条32 位指令;他具有强大的外设支持能力, 32 位外部存储器接口( EMIF) 可以很方便地和SRAM , EPROM , FLASH , SB2SRAM 和SDRAM 等同步和异步存储器或者512 MB 的外部存储空间连接 。 一个DSP 硬件系统可以分为最小硬件系统板和外围接口电路2 部分,DSP 最小板作为开发DSP 系统的基础,一般主要包括电源、复位电路、时钟电路、外部存储器总线接口电路、仿真器接口电路等部分,缺一不可。作为控制系统的最小板,需在其外围接入扩展板,以使系统能够实现相应功能,为此,最小板设计扩展板接口,实现DSP与扩展板及其他芯片通信的目的。 图2-1是TMS320C67xx CPU 结构框图，图2-2是其数据路径。 图2-1 TMS320C67xx CPU (DSP内核)的结构框图 图2-2 TMS320C67xx CPU (DSP内核)的数据路径 3设计步骤 3.1 采集音乐信号 在Matlab软件平台滤波器的性能指标fcd=1400；fcu=1950；Rp=1；As=40，图3-1为原始语音的m文件程序图，图3-2为加噪后的m文件程序图，图3-3为去噪的m文件程序图。 图3-1 原始语音的m文件程序图 图3-2 加噪后的m文件程序 图3-3 去噪后的m文件程序 3.2 CCS仿真代码 #include stdio.h #include c6x.h #include math.h #define pi 3.1415927 #define