使用频率采样法设计的频率采样型FIR滤波器课程设计.docx

1、使用频率采样法设计的频率采样型FIR滤波器课程设计长沙理工大学数字信号处理课程设计报告学院计算机与通信工程专业班级学号通信工程学生姓名 指导教师课程成绩 完成日期课程设计成绩评定学 院计算机与通信工程专 业 通信工程班 级 学 号 学生姓名 指导教师 完成日期 2015年3月13日指导教师对学生在课程设计中的评价评分项目优良中及格不及格课程设计中的创造性成果学生掌握课程内容的程度课程设计完成情况课程设计动手能力文字表达学习态度规范要求课程设计论文的质量指导教师对课程设计的评定意见综合成绩 指导教师签字 2015 年3月15 日课程设计任务书计算机与通信工程 系 通信工程 专业课程名称数字信号处

2、理课程设计时间 20142015学年第二学期12周学生姓名指导老师题目语音信号滤波去噪一使用频率米样法设计的频率米样型 FIR滤波器主要内容：用麦克风采集一段语音信号，绘制波形并观察其频谱特点，加入一个带 外单频噪声，使用频率米样法设计 个满足指标的 FIR滤波器，对该含噪语音信号进行滤波去噪处理，画出频率采样型滤波器结构图。比较滤波前后的波 形和频谱并进行分析，根据结果和学过的理论得出合理的结论。 与不冋滤波方 法的同学比较各种滤波方法性能的优劣。要求：(1) 滤波器指标必须符合工程实际。(2) 设计完后应检查其频率响应曲线是否满足指标。(3) 处理结果和分析结论应该 致，而且应符合理论。(

3、4) 独立完成课程设计并按要求编写课程设计报告书。应当提交的文件：(1) 课程设计学年论文。(2) 课程设计附件(主要是源程序)。语音信号滤波去噪 使用频率采样法设计的频率采样型 FIR 滤波器学生姓名： 指导老师：摘 要 此次课程设计是使用频率采样法设计一个频率采样型滤波器， 对一段加噪的语 音信号进行滤波去噪处理，并且对滤波前后的时域波形和频谱分析滤波性能。此次课程 设计仿真平台为 MATLAB7.0 ，开发工具是 M 语言编程。首先利用 Windows 录音机录 制一段语音信号，并人为加入单频噪声， 然后对信号进行频谱分析以确定所加噪声频率， 再设计滤波器进行滤波去噪处理，最后比较滤波前

4、后的波形和频谱并进行分析，并用 VISIO 绘制滤波器的结构图。由分析结果可知，滤波后的语音信号与原始语音信号基本 一致，即设计的 FIR 滤波器能够去除信号中所加单频噪声，即达到了设计目的。关键词 MATLAB ；滤波去噪；频率采样法；频率采样型滤波器； VISIO1 引言滤波去噪 1是信号处理中一种最基本但十分重要的技术。利用滤波可以从复杂的信 号中提取所需的信号，抑制不需要的信号。滤波器就是这样一种可以在时域和频域对信 号进行滤波处理的系统。通常情况下，有用信号和干扰信号是在不同频段上的，于是通 过对滤波器的频率特性精心设计就能达到滤波的目的。 本课程设计是采用频率采样法设 计频率抽样型

5、滤波器，从而对语音信号滤波去噪。通过对比滤波前后的波形图及回放滤 波前后的语音信号，来判断滤波器对噪声信号确实有滤除作用。1.1 课程设计目的(1) 熟悉使用 MATLAB ；(2) 了解 FIR 滤波器原理及结构；(3) 利用所学数字信号处理相关知识用 MATLAB 设计一个 FIR 滤波器；(4) 提高自己动手的能力。(5) 对加噪的语音信号进行滤波去噪处理， 比较滤波前后的时域波形和频谱并进行分 析；1.2 课程设计的要求(1) 滤波器指标必须符合工程设计。(2) 设计完后应检查其频率响应曲线是否满足指标。(3) 处理结果和分析结论应该一致，而且应符合理论。(4) 独立完成课程设计并按要

6、求编写课程设计报告。1.3 设计平台本课程设计仿真平台为 MATLAB7.0 。 MATLAB 的名称源自 Matrix Laboratory， 1984 年由美国 Mathworks 公司推向市场。它是一种科学计算软件，专门以矩阵的形式 处理数据。 MATLAB 将高性能的数值计算和可视化集成在一起，并提供了大量的内置 函数，从而被广泛地应用于科学计算、控制系统、信息处理等领域的分析、仿真和设计 工作。1993 年 MathWorks 公司从加拿大滑铁卢大学购得 MAPLE 软件的使用权， 从而 以MAPLE为 引擎”开发了符号数学工具箱(Symbolic Math Toolbox)。2 设

7、计原理用 windows 自带的录音机录制一段语音信号， 绘制波形并观察其频谱， 给定相应技 术指标，用频率采样法设计的一个满足指标的频率采样型 FIR 滤波器，对该语音信号 进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析。2.1 FIR 滤波器的设计FIR(Fi nite Impulse Resp on se滤波器3:有限长单位冲激响应滤波器，又称为非递归 型滤波器，是数字信号处理系统中最基本的元件，它可以在保证任意幅频特性的同时具 有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。 因此， FIR 滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。这个

8、公式就为FIR滤波器提供了另外一种结构，这种结构由两部分级联组成其中级联的第一部分为梳状滤波器，其结构如下图所示:_Z -N图2-1梳状滤波器结构图第二部分由N个谐振器组成的谐振柜。它是由N个一阶网络并联组成，而这每一个一阶网络都是一个谐振器其结构如下图所示:H(k) Hk(z1-图2-2 一阶谐振器频率抽样型结构特点:2(1) 它的系数H(k)直接就是滤波器在Wk k处的频率响应。因此，控制滤波器的N频率响应是很直接的。(2) 结构有两个主要缺点：(a) 所有的相乘系数及H(k)都是复数，应将它们先化成二阶的实数，这样乘起来较复 杂，增加乘法次数，存储量。(b) 所有谐振器的极点都是在单位圆

9、上，由wNk决定考虑到系数量化的影响，当系数量化时，极点会移动，有些极点就不能被梳状滤波器的零点所抵消。(零点由延时单元决 定，不受量化的影响)系统就不稳定了。(3) 将一阶网络合并为二阶网络(a)第 k和第N-k个谐振器合并为一个实系数的二阶网络，因为 h(n)是实数，它的DFT也是圆周共轭对称的。(2.5)H(k) H*(N k) k 1,2,3丄 N 1因此，可以将第k和第N-k个谐振器合并为一个二阶网络Hk(z)H(k)1 rWNkz 1H (N k)1讥叫1H(k)1 rWNkz 1*H (k)1 rWNk*z 1 H(k)1 z 1Wn k* WNkr r2WNkWNk*z 210

10、k 1kz12rcos(2)N其中：0k 2ReH( k), 1k 2rReH( k)W,(2.6)(b)第k和第N-k个谐振器合并为一个二阶网络的极点在单位圆内， 而不是在单位圆上，因而从频率响应的几何解释可知，它相当于一个有限Q的谐振器。其谐振频率为:Wk#0*- 2 h詞 )、r | j |图2-3二阶网络结构图除了共轭复根外，还有实根。当N=偶数时，有一对实根，它们分别为k 0,k 两点2H(0 和 H n (z)1 rz 2当N=奇数时，只有一个实根z=r(k=0),即只有H(z)H (号)1 rz 1H0(z).(2.7)(c)修正频率抽样结构流图(N=偶数)H(z)图2-4修正频

11、率抽样结构流图(N=偶数)(1 rNz N) *严7N 1 rzh(N2)1 rz 1N/2 1 1 0k 1kz1 2 k 1 1 z 2rcos( k) N(2.8)修正频率抽样结构流图(N=奇数)H(z)2.3频率采样法设计思路：图2-5修正频率抽样结构流图(N=奇数)N N 1(1 r Z)NH(0)1 rz 1N/2 110k 1kZz 12rcos( k)N(2.9)这种设计方法是从频域进行设计的一种方法，首先给定一个希望逼近的频率响应H(k)Hd(ej )|k 0,1,2,.,N 1知道H(k)后,由IDFT定义，可以用这N个采样值H(k)来惟一确定有限长序列 即h(n)1 N

12、1-H(k)WNnkN k on 0,1,2,., N 1(2.10)h( n)，(2.11)(2.12)H(z) h(n)zn 0N 1H (ejw ) h(n )e jwnn 0(2.13)内插公式：H(ejN 1 2) H(k) kk 0 N(2.14)式中，(是内插函数(、sin( N / 2) j (n 1)/2) eN sin( /2)(2.15)四种线性相位的FIR滤波器如下表2-1所示。表2-1四种线性相位的FIR滤波器甲4)1型AT迦奇数就砺关于皿 血偶对称 厘型風)其于沪D、Qjt偶对荊”关于tu=t裔对称23型加对乂城Ml)，N炯数换血关于=0、JL亦奇对称V- 1 JT

13、口型城)关于F、，加奇对称，关于过=謝禺对称2 2逼近误差及其改进措施这样设计所得到的频响H(ej 3与要求的理想频响Hd(ej 3会有一定的差别(1)采样点上滤波器的实际频率响应是严格地和理想频率响应数值相等的。但是在采样点之间的频响则是由各采样点的加权内插函数的延伸叠加而成的 ,因而有一定的逼近误差， 误差大小取决于理想频率响应曲线形状。(3) 理想频率响应特性变化越平缓，则内插值越接近理想值，逼近误差越小。(4) 如果采样点之间的理想频率特性变化越陡， 则内插值与理想值的误差就越大，因而在理想频率特性的不连续点附近，就会产生肩峰和起伏。滤波器性能的改善：(1) 增加过渡带采样点，它可以大

14、大减少振荡，阻带衰减也可以得到进一步改善。一 般一点到二点的过渡带采样即可得到满意的结果。(2) 增加采样点密度，过渡带的宽度与采样点数 N成反比。但N值意味着或长度的增加，滤波器运算量必然增大 。3设计步骤3.1设计流程图图3-1流程图3.2录制语音信号使用电脑录音机录制一段语音信号，属性设置为 8000Hz, 8位单声道的音频格式。如图3.1-2所示：图3-2语音信号设置然后在MATLAB软件平台下，利用函数wavread对语音信号进行采样，源程序为：x,fs,bits=wavread(ht.wav)，记住采样频率和采样点数，matlab 实现得：fs=8000; bits=8。3.3语音

15、信号的频谱分析在MATLAB中编写加噪之后时域，频域对比图的源程序如下 ：x,fs,bits=wavread(D:ht.wav);%输入参数为文件的全路径和文件名，输出的第一个参数是每个样本的值， fs是生成该波形文件时的采样率，bits是波形文件每样本的编码位数。sou nd(x,fs,bits); %按指定的采样率和每样本编码位数回放N=length(x); %计算信号x的长度fn=2500; %单频噪声频率，此参数可改t=0:1/fs:(N-1)/fs; %计算时间范围，样本数除以采样频率x=x; y=x+0.1*sin(fn*2*pi*t); % 加噪sou nd(y,fs,bits)

16、; %应该可以明显听出有尖锐的单频啸叫声X=abs(fft(x); Y=abs(fft(y); %对原始信号和加噪信号进行fft变换，取幅度谱X=X(1:N/2); Y=Y(1:N/2); % 截取前半部分f=0:deltaf:fs/2-deltaf; %计算频谱频率范围运行结果如下图所示：由图3.3-1可以看出，在频域为2500Hz处加入了一个单频噪声，而加入噪声之后， 时域的波形出现了明显失真，通过听取原声音信号 x和加噪声音信号y,可以明显听到y语音信号中有一明显尖锐噪声。3.4滤波器设计设计一个带阻滤波器，将单频信号滤出去，源程序如下所示：fpd=2400;fsd=2450;fsu=2

17、550;fpu=2600;Rp=0.5;As=45; % 带阻滤波器设计指标fcd=(fpd+fsd)/2;fcu=(fpu+fsu)/2;df=mi n( (fsd-fpd),(fpu-fsu);%计算上下边带中心频率fcd fcu，和频率间隔dfwp=fcd/fs*2*pi;ws=fcu/fs*2*pi;dw=df/fs*2*pi; % 将 Hz 为单位的模拟频率换算为 rad为单位的数字频率M=41;alpha=(M-1)/2;% 设置抽样点数l=0:M-1;wl=2*pi*l/M;T1=0.6025,T2=0.127%过渡带的抽样值NP仁fix(M*fpd/8000);NS仁fix(M

18、*fsd/8000);NS2=ceil(M*fsu/8000);NP2=ceil(M*fpu/8000);%各边界频率的抽样值1Hrs=o nes(1,NP1-2),0.6025,0.127,zeros(1,NP2-NP1+1),0.127,0.6025,o nes(1,M-2*NP2-5),0.6025,0.127,zeros(1,NP2-NP1+1),0.127,0.6025,o nes(1,NP1-3);Hdr=1,1,0,0,1,1;%带阻滤波器的参数设置wdl仁 fpd*2/fs;wdl2=fsd*2/fs;wdl3=fsu*2/fs;wdl4=fpu*2/fs;%rad wdl=0

19、,wdl1,wdl2,wdl3,wdl4,1;k仁 0:floor(M-1)/2);k2=floor(M-1)/2)+1:M-1;an gH=-alpha*(2*pi)/M*k1,alpha*(2*pi)/M*(M-k2);H=Hrs.*exp(j*a ngH);h=real(ifft(H,M);db,mag,pha,grd,w=freqz_m(h,1);Hr,ww,a,L=hr_type2(h);程序运行结果如下图3.3-2所示：ffj F-igure 11 : -： : 4: | ； - - - -I- - _ju 20呦i弋.申器悟=响应0.60.5 1wpiSEdht VIqmi In

20、sert Tools Desktoo LVindow Help图3-4滤波器参数图由图3-4可以看出，滤波器的衰减大于设定值 As=45,满足性能指标，滤波器的衰减可由增加过度带宽来得到。上图中的抽样值 M二41,因此阻带的带宽比较宽。可以通过调大M值来使得阻带的带宽变窄，但这样会使得滤波器的阶数相应的增加，使得其结 构图标的复杂，为了是滤波器能达到相应的性能要求，又使得其结构变得简单，最终取M值为413.5信号滤波处理原程序如下所示:%进行滤波处理 z=filter(h,1,y); %用设计好的滤波器对y进行滤波Z=fft(z);Z=Z(1:le ngth(Z)/2); % 计算频谱取前一半

21、程序运行结果如下图3.3-3所示:图3-5滤波前后时域，频域对比图频谱取前一半由上图3-5可以看出，在滤波之后时域图能得到恢复，频域图中的单频噪声信号也的到滤除，说明了设计的滤波器能滤除加入的噪声信号，因此说达到了 设计的要求3.6结果分析开始通过分析决定设计一个带阻的滤波器来滤除加入的单频噪声， 根据噪声的频率来设计阻带的范围。在采集语音信号后，按照步骤用频率采样法设计频率采样型滤波器。 由图3-4可知,设计的滤波器达到要求。我们观察到图3-5滤波前后语音信号的波形对比图， 发现时域 波形中加干扰噪声后有明显的变化，不过经过滤波后几乎没有变化，说明设计的滤波器 达到要求。再通过听取原始语音信

22、号，加噪信号 y，和滤波之后的信号z。对比之后，发现滤波器确实滤除了噪声。从理想的角度考虑，该带阻滤波器的阻带带宽应该可以变得更窄， 让滤波效果更好， 但是这样的采样值会变得非常大，考虑到实际的情况，通过权衡之后，决定牺牲带宽来 使得滤波器的阶数降低，因此在上图 3-5中我们可以看出，在噪声频谱左右两边的信号 也被滤除了。3.6滤波器结构设计通过计算得出滤波器的系数如下图所示:2. QQOQL 0000-E97661. 0000乱 QOOD-o.7i0*397(L 000 Q-1.90681.00G02. 00400. MIS-0,9S831. 0000-1.79231. 00002. QOH

23、-0, 97370.9737L 0000-I.43S91.00002. QD0Q0.9534-fl,9534L OWO-1.1410K00O02. 0QC0-0.M7S0.A275L 0000-1.21251,00002, 00000,7 B9621. 0000-0,45541. ODDO2-0DCD-0, 8&960,95961. 0000-0.47601.00002- 00000. 81797* sin昇 ovoo-CL 3808I. OOOQ-直 771ETT151” 0000-0,07661.0000队 2540叭 72057+ 7205L 00000.22941. 00000. 00

24、00-0, 663Q. 6653.oooo0.S30Q1.0000o.oooaQ. 3162MML 0Q00e.aiti1. 0000o. oodo(L ontio0.4)251.0000i.osn1.00000.2540此9G1-Q+ 酣1, 0000i.33071.0000-0.4011Q” *0911. 0000LS4301.00002. 00000.-0.33801. 00*0L719!l.OOQO2. aaoo-0. 26500.咖0LQOOQI.S550EQOGQ2. 00000.1*04-0+ E9Q41. 0000L94741.00002-0000-0. 11470, tl47

25、1. 00001.991!1.00D01. 00000-0383-0.03*3L 000Q-L 00000图3-6滤波器系数根据以上求出的系数，利用 VISIO软件来绘制滤波器的结构图如下图所示:z1z1z154-1-0.6653 0.2540-0.33804z1z1-1z1-1.8550-0.7715 2.0000.3380 1.20502.0000-0.9737-0.8596 2.0000h = Z1 0.8596z11.9766 Z1 0.99710.1904 2.0000-0.1904 丄1.7923 Z1 0.9737-1-0.9275 2.00001.4410 z： 0.92750

26、.3808 Jz1 0.7715-0.2294 z=0.6653-1 z1-1.54301.9068-1Z1 -0.9883z11.63590.9534 2.0000I 片Z1 -0.9534-1z10.8962 2.00001.2125 z1-0.89620.8179 2.00000.6760-1-1-1z1-0.8179 Jz10.7205 1.2050-0.4091 O.2540 07 z1 0.4091 kz1-0.2650 2.000-1.7191 *z1 0.R650-1 z1-0.1147 2.0000応41十00383 . F -0.0383 *I -11.00002.0000

27、-1 9474 z1 0J147-1 z1z1图3-7滤波器的结构图y(n)4 出现的问题及解决办法 在这次课程设计当中，由于基础不扎实，出现了很多问题，即有 MATLAB 软件操作不 当，也有知识掌握程度不够出现的各种问题。1，在一开始录制语音信号时，没有修改参数值，导致频率抽样过高，后在老师的 提醒之下将频率修改为 8000hz。2，在调用音乐文件时没有将文件放在 MATLAB 的工作文件夹下面， 到时文件找不 到，还有程序中用到的各种函数都没有放入正确的位置。3，在绘制加噪前后频率对比图时，留白过多，对比不明显，通过 axis 函数对横纵坐标进行限定4，在使用 visio 软件绘制滤波器

28、的结构图时，没有合理规划整个页面，导致重新绘 制。5，在使用 visio 绘制结构图时，前面滤波器参数设置出现错误，所以滤波器的系数 发生了改变，只能重新修改参数。6，在设计滤波器的过程中，没有准确理解每一条指令代表的含义，导致程序前后 不对应，出现很多错误，通过 matlab 中的错误提示，准确找到错误的那一行代码，进行 修改。7，在设计滤波器的过程当中，滤波器的衰减小于开始所设置的值，通过牺牲过渡 带和调节过渡带的采样值，即 T1 和 T2 来使得衰减大于所设定值 As。8，频谱图没有进行去前半段的处理。 9，最后听取滤波后声音，没有将其保存，在请教了同学之后，学会了如何保存滤 波后的声音文件。10，在设计滤波器的过程当中，抽样点数过高，是的滤波器的阶数很高，导致绘图 困难。通过降低抽样的点数，来平衡滤波器阶数和阻带带宽之间的矛盾。5 结束语这是第三次进行课程设计了，在前面的课程设计当中我们学习到了如何使用 matlab,所以对于matlab软件的使用并没有那么陌生了，尽管如此， 在使用matlab的过 程当中还是出现了很多错误，比如说忘记添加函数文件，参数前后不对应