语音信号的FIR滤波器处理课程设计.docx
《语音信号的FIR滤波器处理课程设计.docx》由会员分享,可在线阅读,更多相关《语音信号的FIR滤波器处理课程设计.docx(17页珍藏版)》请在冰豆网上搜索。
![语音信号的FIR滤波器处理课程设计.docx](https://file1.bdocx.com/fileroot1/2023-1/27/0e5029b2-cd21-4655-b791-1961787ab4f0/0e5029b2-cd21-4655-b791-1961787ab4f01.gif)
语音信号的FIR滤波器处理课程设计
DSP课程设计报告
语音信号的FIR滤波器处理
******
班级:
12电信3
学号:
************
*******
日期:
2015.06.02~2014.06.19
华南农业大学电子工程学院
摘要
随着信息与数字技术的发展,数字信号处理已经成为当今极其重要而学科与技术领域之一。
它在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。
在数字信号处理的基本方法中,通常会涉及到变换、滤波、频谱分析、调制解调和编码解码等处理。
其中滤波是应用非常广泛的一个环节,数字滤波器的理论和相关设计也一直都是人们研究的重点之一。
FIR滤波器的是非递归的,稳定性好,精度高;更重要的是,FIR滤波器在满足幅频响应要求的同时,可以获得严格的线性相位特征。
因此,它在高保真的信号处理,如数字音频、图像处理、数据传输和生物医学等领域得到广泛应用。
CSS集成开发环境使用CCS内置的软件仿真simulator对程序进行编译,调试和运行,主要用于检测目标程序运行的正确性和连贯性,并能通过仿真器与目标板连接,在目标板上实时观察效果。
在本次设计中,我们选择的课题是基于DSP的语言信号的FIR滤波处理。
首先利用MATLAB进行了仿真,得到滤波前后的时域波形和频谱。
然后通过调用MATLAB的分析工具FDATOOL,根据仿真结果导出了滤波器的相关参数,将原始信号数据和滤波器参数输入CCS进行DSP编程。
最后在DSP中实现了FIR低通滤波,并通过CCS的频谱分析功能查看了最终DSP的滤波效果。
关键词:
语音信号处理FIR滤波MATLABCCS
1.前言4
2.设计需求4
3.滤波器设计原理介绍4
3.1数字滤波器设计原理4
3.2FIR滤波器的基本结构5
3.3滤波器的特点5
3.4窗函数的介绍6
4.FIR滤波器的设计7
5.FIR滤波器的MATLAB设计8
6.FIR滤波器的DSP实现10
6.1CCS设计与仿真10
6.2CCS仿真结果13
6.调试问题及解决办法16
7.设计感想16
参考文献17
8.附录17
1前言
通过语音传递信息是人类最重要、最有效、最常用、最方便的交换信息的方法,因此,如何去除混杂在有用信号中的噪声并获得低损甚至无损的有用语音信号成为语音信号处理(即用数字信号处理技术和语言学知识对语音信号进行处理并提取有用信息)至关重要的问题。
但是因为语音识别技术对信号噪声非常敏感,系统本身和器件对系统的处理结果都有着不可忽视的干扰,因此很难实现效果较好的语音信号滤波处理。
而DSP因为运算速度快,片上资源丰富和能够实现复杂的线性和非线性算法等特性,在语音信号处理技术方面有明显的优势。
在数字信号处理中,滤波占有极其重要的地位。
数字滤波是语音信号处理、图像处理、模式识别、频谱分析等应用的基本处理算法。
用DSP芯片实现数字滤波除了具有稳定性好、精确度高、不受环境影响等优点外,还具有灵活性好等特点。
本文介绍了数字滤波器的设计基础及用窗函数法设计FIR滤波器的方法,运用MATHLAB语言实现了低通滤波器的设计并用CCS2.0进行观察效果。
2设计需求
1)滤波器的设计的原理介绍
2)FIR滤波器的设计
3)FIR滤波器的DSP实现
4)滤波器为低通滤波器,并用用语音信号通过滤波器后查看滤波效果是否达到设计要求,进行FIR滤波器的检验
根据老师的的意见,我们自己拟定了滤波器的指标:
Fs=8000hz,Wp=1500Hz,Ws=2000Hz,通带波纹为0.01,阻带波纹为0.1,N=37。
3滤波器设计原理介绍
3.1数字滤波器的设计原理
数字滤波器的设计问题就是寻找一组系数ai和bi,使得其性能在某种意义上逼近所要求的特性。
如果在s平面上去逼近,就得到模拟滤波器,如果在z平面上去逼近,则得到数字滤波器。
数字滤波是将输入的信号序列,按规定的算法进行处理,从而得到所期望的输出序列。
一个线性位移不变系统的输出序列y(n)和输入序列x(n)之间的关系,
应满足常系数线性差分方程:
(2.1.1)
x(n)为输入序列,y(n)为输出序列,ai、bi为滤波器系数,N为滤波器的阶数。
3.2FIR滤波器的基本结构
在式2.1.1中,若所有的ai均为0,则得FIR滤波器的差分方程:
(2.2.1)
对式(2.2.1)进行z变换,可得FIR滤波器的传递函数:
(2.2.2)
由此可得到FIR滤波器的结构如图1所示。
FIR滤波器的单位冲击响应h(n)是一个有限长序列。
若h(n)为实数,且满足偶对称或奇对称的条件,即h(n)=h(N-1-n)或h(n)=-h(N-1-n),则FIR滤波器具有线性相位特性。
3.3滤波器的特点
优点:
(1)很容易获得严格的线性相位,避免被处理的信号产生相位失真,这一特点在宽频带信号处理、阵列信号处理、数据传输等系统中非常重要;
(2)可以在幅度特性随意设计的同时,保证精确、严格的线性相位,还可得到多带幅频特性;
(3)极点全部在原点(永远稳定),无稳定性问题;
(4)任何一个非因果的有限长序列,总可以通过一定的延时,转变为因果序列,所以因果性总是满足;
(5)无反馈运算,运算误差小。
缺点:
(1)因为无极点,要获得好的过渡带特性,需以较高的阶数为代价;
(2)无法利用模拟滤波器的设计结果,一般无解析设计公式,要借助计算机辅助设计程序完成。
3.4窗函数简介
数字信号处理的主要数学工具是傅里叶变换。
而傅里叶变换是研究整个时间域和频率域的关系。
不过,当运用计算机实现工程测试信号处理时,不可能对无限长的信号进行测量和运算,而是取其有限的时间片段进行分析。
做法是从信号中截取一个时间片段,然后用观察的信号时间片段进行周期延拓处理,得到虚拟的无限长的信号,然后就可以对信号进行傅里叶变换、相关分析等数学处理。
无线长的信号被截断以后,其频谱发生了畸变,原来集中在f(0)处的能量被分散到两个较宽的频带中去了(这种现象称之为频谱能量泄漏)。
为了减少频谱能量泄漏,可采用不同的截取函数对信号进行截短,截断函数称为窗函数,简称为窗。
信号截短以后产生的能量泄漏现象是必然的,因为窗函数w(t)是一个频带无限的函数,所以即使原信号x(t)是有限带宽信号,而在截短以后也必然成为无限带宽的函数,即信号在频域的能量与分布被扩展了。
又从采样定理可知,无论采样频率多高,只要信号一经截短,就不可避免地引起混叠,因此信号截短必然导致一些误差。
泄漏与窗函数频谱的两侧旁瓣有关,如果两侧瓣的高度趋于零,而使能量相对集中在主瓣,就可以较为接近于真实的频谱,为此,在时间域中可采用不同的窗函数来截短信号。
4FIR滤波器的设计
4.1FIR滤波器设计
FIR滤波器的设计问题在于寻求一系统函数,使其应逼近滤波器要求的理想频率响应。
4.2窗函数设计的基本方法
从时域出发,设计h(n)逼近理想hd(n),设理想滤波器的单位脉冲响应为hd(n),则有
所求得的一般是无限长的,且是非因果的。
要想得到一个因果的有限长的滤波器h(n),最直接的方法是截断,或者说用一个窗口函数对进行加窗处理,即。
因此,应选择合适的窗函数。
4.3滤波器的
算法实现
FIR滤波器的输出表达式为
式中,为滤波器系数;x(n)表示滤波器在n时刻的输入;y(n)为n时刻的输出。
它的基本算法是一种乘法-累加运算,即不断地输入样本x(n),经过延时后,再进行乘法-累加,最后输出滤波结果y(n)。
1)线性缓冲区法
线性缓冲区法又称延迟线法。
其特点:
·对于N级的FIR滤波器,在数据存储器中开辟一个N单元的缓冲区(滑窗),用来存放最新的N个输入样本;
·从最老样本开始取数,每取一个样本后,将此样本向下移位;
·读完最后一个样本后,输入最新样本存入缓冲区的顶部。
2)循环缓冲区法
本次设计的FIR滤波器所采用的就是循环缓冲区法。
循环缓冲区法的特点如下:
·对于N级FIR滤波器,在数据存储器中开辟一个N单元的缓冲区(滑窗),用来存放最新的N个输入样本;
·从最新样本开始取数;
·读完最后一个样本(最老样本)后,输入最新样本来代替最老样本,而其他数据位置不变;
·用片内BK(循环缓冲区长度)寄存器对缓冲区进行间接寻址,使循环缓冲区地址首尾相邻。
5FIR滤波器的matlab实现
MATLAB是一种功能强、效率高、便于进行科学和工程计算的交互式软件包,它集数值分析、矩阵运算、信号处理和图形显示于一体,为用户提供了方便、友好的界面环境。
MATLAB中的工具箱(Toolbox)包含了许多实用程序。
它提供了多种FIR滤波器设计方法。
fir1函数用来设计标准频率响应的基于窗函数的FIR滤波器,可实现加窗线性相位FIR数字滤波器的设计。
具体语法如下:
b=fir1(n,Wn)
b=fir1(n,Wn,‘ftype’)
b=fir1(n,Wn,Window)
b=fir1(n,Wn,‘ftype’,Window)
其中n为滤波器的阶数;Wn为滤波器的截止频率;ftype为用来决定滤波器的类型,当ftype=high时,可设计高通滤波器;当ftype=stop时,可设计带阻滤波器。
Window为用来指定滤波器采用的窗函数类型,Window参数可采用的窗口函数有:
Boxcar,Hanning,Bartlett,Blackman,Kasier和chebwin等,其默认时为Hamming窗,从而得到滤波器的系数。
fir2函数用来设计有任意频率响应的各种加窗FIR滤波器。
具体语法如下:
b=fir2(n,f,m)
b=fir2(n,f,m,Window)
b=fir2(n,f,m,npt)
b=fir2(n,f,m,npt,Window)
b=fir2(n,f,m,npt,lap)
b=fir2(n,f,m,nptt,lap,Window)
其中n为滤波器的阶数;f为频率点矢量;m为幅度点矢量;Window用来指定所使用的窗函数类型,默认值为汉明(Hamming)窗;npt用来指定fir2函数对频率响应进行内插的点数;lap用来指定fir2函数在重复频率点附近插入的区域大小,从而得到滤波器的系数。
具体的代码见附录,输入之后可以导出滤波器的参数。
CCS的输入语音信号的辅助产生(MATLAB)
1、理想情况:
FIR滤波器的输入数据存储在indata.dat文件中,共300个点,由频率为1000Hz和2500Hz的两个余弦波叠加而成
2、为了验证设计的可行性,自身录制一个语音片段,主要信号为低频成分,其中夹杂着高频成分的噪声。
在MATLAB中产生一个indata.dat文件供CCS输入,在matlab中输入的代码为:
x=wavread('D:
\matlab\work\XINZAO.wav');%语音文件的路径
y=fft(x,1024);
plot(abs(y(1:
1000)));%采样100个点
x1=round(1024*x);
fid=fopen('indata.dat','w');%生成一个dat文件
fprintf(fid,'1651100E21000\n');
fprintf(fid,'%d\n',x1);%输出
6FIR滤波器的DSP实现
CCS是一种针对TMS320系列DSP的集成开发环境,在windows操作系统下,采用图形接口界面,提供环境配置、源文件编辑、程序调试、跟踪和分析等工具,提供了配置、建立、调试、跟踪和分析程序的工具,从而完成编辑、编译、链接、调试和数据分析等工作。
6.1CCS设计
图15CCS汇编程序流程图
与仿真
6.1.1程序分析
FIR滤波器的DSP实现主要有线性缓冲区法和循环缓冲区法两种。
线性缓冲区法又称延迟线法。
其特点:
(1)对于N级的FIR滤波器,在数据存储器中开辟一个N单元的缓冲区(滑窗),用来存放最新的N个输入样本;
(2)从最老样本开始取数,每取一个样本后,将此样本向下移位;
(3)读完最后一个样本后,输入最新样本存入缓冲区的顶部。
循环缓冲区法的特点如下:
(1)对于N级FIR滤波器,在数据存储器中开辟一个N单元的缓冲区(滑窗),用来存放最新的N个输入样本;
(2)从最新样本开始取数;
(3)读完最后一个样本(最老样本)后,输入最新样本来代替最老样本,而其他数据位置不变;
(4)用片内BK(循环缓冲区长度)寄存器对缓冲区进行间接寻址,使循环缓冲区地址首尾相邻。
本次设计的FIR滤波器所采用的就是循环缓冲区法。
对DSP进行初始化,并设置好存储单元。
为数字滤波作准备,将预先设计好的有N个抽头的FIR数字滤波器的冲激响应序列h(n)中的N个数值放入存储单元bn中,从indata数据存储段中读入抽样值,将抽样值逐一放入存储单元xn中。
主要循环程序如下:
STM#-1,AR0
RPTZA,#N;将ACC累加器清零,循环N次
MAC*AR2+0%,*AR3+0%,A
STHA,*(yn);保存输出
MVKD*(yn),*AR5+
将累加器清零,将存储单元bn与xn值对应相乘并累加到ACC中,循环N次将ACC中的值保存输出。
再从indata数据存储段中读入一个新的数据,代替原来的最老样本MVDD*AR4+,*AR2+0%,重复循环直到数据全部处理完。
最后将处理完的数据输入到outdata数据存储段中。
6.1.2数据输入
MATLAB辅助DSP实现FIR滤波器滤波,其总体过程是为DSP中编写处理程序提供滤波器系数。
在MATLAB中利用滤波器设计、分析工具(FDATOOL),根据指定的滤波器性能快速设计一个FIR,然后把滤波器系数以头文件形式导入CCS中,头文件中包含MATLAB辅助DSP实现FIR数字滤波器的阶数和系数数组,在MATLAB中调试、运行DSP程序并显示、分析处理后的数据。
打开CCSIDE中file的data加载初始化data数据,在Address中输入数据存放的起始地址indata,并设置数据栈长度。
6.1.3运行程序仿真
运行程序,查看输入输出波形,修改相应参数进行调试。
在View的Graph中单击Time/frequency出现graphpropertydialog框。
将显示类型,图形名称,起始地址,抽样点数,数据类型等分别进行设置,输出各种波形。
输出滤波前的信号时域波形。
首先,设置如下图所示,输出滤波前的信号时域波形。
图2滤波前语音信号的时域GraphPropertyDialog设置参数图
将DisplayType改为FFTMagnitude就可以看到滤波前信号的频谱,将StartAddress改为outdata就可以看到滤波后信号的时域波形和频谱。
6.2CCS仿真结果
图3滤波前语音信号的时域波形图
图4滤波前信号频谱波形图
6.2.2FIR滤波器滤波后的信号
图5滤波后信号时域波形图
图6滤波后信号频谱图
图3波前的信号,由高频信号和原信号叠加而成的波形较为杂乱,加噪后的信号幅值较大。
图5波后的信号,从图上可以看到信号幅度较低,波形较为清晰。
对比图4图6可以看到加躁后信号幅值较大,信号的能量集中在高频段,滤波后,高频成分大部分被滤去,信号的能量主要集中在低频段。
7调试问题及解决方法
1)实验时,未做初始化设置直接进入,没有装入gel文件及相关参数,导致输入代码编译链接后不能载入程序。
解决方法:
初始化设置软件仿真cpu和硬件仿真cpu,并设置其中所用的gel文件,保存并退出启动dsp选择忽略进入编辑编译初始化界面,由软件仿真cpu进入编辑编译工具。
2)在解决问题1后,载入程序,载入数据,并设置载入数据的初始地址和长度,运行程序,发现程序一直卡死在某个地方,无法显示正确的时域图和频域图。
解决方法:
换了一台电脑运行,发现问题已被解决,可能原因是自己电脑安装时装少了一些文件。
8设计感想
刚开始面对大作业的时候,我们都还是觉得很生疏,不知从哪里入手,尽管做过实验。
对于自己选的题目,具体需要做什么,能不能做出来,我们并不清楚。
但当时间一天天过去后,我们开始认真起来了。
从复习课本,再到复习实验,相互探讨编程思路,上网搜索资料,调试实验等,我们都在努力的自学着。
首先,是学习上的巩固。
一方面,通过复习课本和实验,对DSP的基础知识又有了很大的巩固。
其次,通过对用Mathlab实现FIR滤波器的设计,熟悉了matlab软件的一些相关的窗口函数以及相关功能的调用,如怎样实现窗函数的调用,怎样实现滤波,等等。
第二,对FIR的滤波器性能和作用有了更深一层的了解。
FIR滤波器的应用十分广泛,当今许多信号处理系统和图像处理系统等都要求信号具有线性相位特性。
在这方面,FIR滤波器有十分独特的优点,运用mathlab语言,我们能够很容易的设计出具有严格线性相位的FIR滤波系统,以及比较容易的实现。
此外,我们还知道了FIR滤波器是永远稳定的,因为FIR滤波器的冲激响应是有限长序列,其系统函数为一个多项式,它所含的的极点多为原点。
最后,也感谢徐梅宣老师的指导,跑了她办公室两次咨询情况,都很耐心解答,真给了一些写报告或者软件应用的建议。
参考文献
[1]邹彦,DSP原理及应用(修订版),北京:
电子工业出版社,2012.
[2]王宏.MATLAB6.5及其在信号处理中的应用.[M].北京:
清华大学出版社,2004
[3]高西全,丁玉美,阔永红.数字信号处理---原理、实现及应用.[M].北京:
电子工业出版社,2006
[4]胡广书编著,数字信号处理理论、算法与实现.北京:
清华大学出版社,2005;
附录
1.主程序fir.asm:
.title"FIR.ASM"
.mmregs
.defstart
.bssyn,1
xn.usect"xn",37
hn.usect"hn",37
indata.usect"indata",300
outdata.usect"outdata",256
.data
fir_coff:
.word-10*1024/10000,-36*1024/10000
.word0*1024/10000,65*1024/10000
.word32*1024/10000,-88*1024/10000
.word-93*1024/10000,90*1024/10000
.word184*1024/10000,-47*1024/10000
.word-297*1024/10000,-71*1024/10000
.word417*1024/10000,316*1024/10000
.word-524*1024/10000,-848*1024/10000
.word598*1024/10000,3108*1024/10000
.word4375*1024/10000
.word3108*1024/10000,598*1024/10000
.word-848*1024/10000,-524*1024/10000
.word316*1024/10000,417*1024/10000
.word-71*1024/10000,-297*1024/10000
.word-47*1024/10000,184*1024/10000
.word90*1024/10000,-93*1024/10000
.word-88*1024/10000,32*1024/10000
.word65*1024/10000,0*1024/10000
.word-36*1024/10000,-10*1024/10000
.text
start:
ssbxfrct;使用小数运算
stm#hn,ar1;系数首地址
rpt#36;将系数移入循
mvpdfir_coff,*ar1+;环缓冲区
stm#-1,ar0
stm#outdata,ar5
stm#indata,ar4
stm#xn+36,ar2
rpt#36;将输入数据移入
mvdd*ar4+,*ar2+0%;循环缓冲区
stm#xn+36,ar2;指向第一个输入
stm#hn+36,ar3
stm#37,bk;缓冲区大小37
stm#255,brc;块重复256次
stm#-1,ar0
rptbloop-1;块程序重复大小
rptza,#36;计算一个输出
mac*ar2+0%,*ar3+0%,a
stha,*(yn);保存输出
mvkd*(yn),*ar5+
mvdd*ar4+,*ar2+0%;读进一个输入
loop:
nop
bloop
.end
2.链接命令文件fir.cmd:
MEMORY
{PAGE0:
PARAM:
org=00E0hlen=1000h
PAGE1:
SPRAMorg=0060hlen=0020h
DARAM:
org=0080hlen=1380h
}
SECTIONS
{
.text:
>PARAMPAGE0
.data:
>PARAMPAGE0
.bss:
>SPRAMPAGE1
xn:
align(64){}>DARAMPAGE1
hn:
align(64){}>DARAMPAGE1
indata:
>DARAMPAGE1
outdata:
>DARAMPAGE1
}