《FIR数字滤波器的(海明)窗函数法设计》由会员分享，可在线阅读，更多相关《FIR数字滤波器的(海明)窗函数法设计（15页珍藏版）》请在金锄头文库上搜索。
1、课程设计说明书 NO.1FIR 数字滤波器的（海明）窗函数法设计1课程设计目的 （1）熟悉并掌握 MATLAB 中有关声音（wave ）录制、播放、存储和读取的函数。（2）加深对 FIR 数字滤波器设计的理解，并用窗函数法进行 FIR 数字滤波器的设计。（3）将设计出来的 FIR 数字滤波器利用 MATLAB 进行仿真。（4）对一段音频文件进行加入噪声处理，对带有噪声的文件进行滤波处理。2.设计方案论证2.1 Matlab 语言概述MATLAB 是一种以矩阵运算为基础的交互式程序语言，专门针对科学、工程计算及绘图的需求。随着版本的不断升级，内容不断扩充，功能更加强大，从而被广泛应用于仿真技术、
2、自动控制和数字信号处理领域。此高级语言可用于技术计算 此开发环境可对代码、文件和数据进行管理交互式工具可以按迭代的方式探查、设计及求解问题 数学函数可用于线性代数、统计、傅立叶分析、筛选、优化以及数积分等二维和三维图形函数可用于可视化数据 各种工具可用于构建自定义的图形用户界面 各种函数可将基于 MATLAB 的算法与外部应用程序和语言（如 C、C+、Fortran、Java、 COM 以及 Microsoft Excel）集成不支持大写输入，内核仅仅支持小写2.2 声音处理语音是人类获取信息的重要来源和利用信息的重要手段。语音信号处理是一门发展十分迅速，应用非常广泛的前沿交叉学科，同时又是一
3、门跨学科的综合性应用研究领域和新兴技术。声音是一种模拟信号，而计算机只能处理数字信息0 和 1。因此，首先要把模拟的声音信号变成计算机能够识别和处理的数字信号，课程设计说明书 NO.2这个过程称为数字化，也叫“模数转换” 。在计算机对数字化后的声音信号处理完后，得到的依然是数字信号。必须把数字声音信号转变成模拟声音信号，然后再输出到扬声器，这个过程称为“数模转换” 。2.3 数字滤波器的介绍数字滤波器(digital filter)是由数字乘法器、加法器和延时单元组成的一种装置。其功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。数字滤波器是一个离散时间系统（按预定的算法，
4、将输入离散时间信号转换为所要求的输出离散时间信号的特定功能装置） 。应用数字滤波器处理模拟信号时，首先须对输入模拟信号进行限带、抽样和模数转换。数字滤波器输入信号的抽样率应大于被处理信号带宽的两倍，其频率响应具有以抽样频率为间隔的周期重复特性，且以折叠频率即 12 抽样频率点呈镜像对称。为得到模拟信号，数字滤波器处理的输出数字信号须经数模转换、平滑。数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。2.4 FIR 滤波器基本结构FIR 滤波器的数学表达式为：式中：N 为 FIR 滤波器的抽头数；x(n)为第 n 时刻的输入样本；h(i)为 FIR滤波器第 i 级抽头系数。其
5、相应的 z 变换为：式中：z-i 为 N-1 阶多项式。在自适应处理、数据通信等领域中往往要求信号在传输过程中不能有明显的相位失真，FIR 滤波器可以做到线性相位满足此要求。F1R 滤波器实质上是一个分节的延迟线，把每一节的输出加权累加，得到滤波器的输出。对于 FIR 滤波器的单位脉冲响应 h(i)只要满足以下 2 个条件之一，则为线性相位滤波器。课程设计说明书 NO.3线性相位的 FIR 滤波器具有中心对称的特性，其对称中心在 N2 处。(2)由性能指标确定窗函数 w(n)和窗口长度 N(3)求得实际滤波器的单位脉冲响应 h(n)(4)检验滤波器性能。设计常用的窗函数有矩形窗、汉宁窗、海明窗
6、、凯撒窗等。其中：海明窗的旁瓣峰值小于主瓣峰值的 1，99963的能量集中在主瓣内且通过海明窗设计的 FIR 滤波器在较少的阶数下可以得到较小通带纹波，非常适合工程设计。2.5 利用 Windows 进行语音信号的采集利用 windows 下的录音机，按照（开始程序附件娱乐 录音机，文件 属性 立即转换 8000KHz，8 位，单声道）的顺序操作，录制一段自己的语音，录制时间为 5 秒，如图 1，图 2 所示，将自己录好的语音文件保存为“zf.wav ”。图 1 选择 windows 下的录音机课程设计说明书 NO.4图 2 用 8000Hz 采样录音2.6 语音信号的分析(1) 将上一步骤中
7、保存下来的语音信号文件“zf.wav”复制到计算机装有Matlab 软件的磁盘中相应 Matlab 目录中的“work”文件夹中:( C:Program FilesMATLAB71work)。(2) 双击桌面上 Matlab 软件的快捷图标 ，打开 Matlab 软件。(3) 在 Matlab 菜单栏中选择 “FilenewM-File”或是点击快捷按钮 ，打开 m 文件编辑器。(4) 在 m 文件编辑器中输入相应的指令将自己的语音信号导入 Matlab 工作台。(5) 画出原始语音信号 s 的波形，由于原始语音信号开始一段会是无用的语音信号，因此要截取掉，截取的一段语音信号为 1 至 1+f
8、s-1，即从 1 到7999，画出截取原始语音信号 s1 的波形，代码如下，波形如图 3，所示：close allclear allclcs, fs,bits=wavread(C:Program FilesMATLAB71workzf.wav); s1=s(1:8000);sound(s1,fs,bits);figure(1);subplot(211)课程设计说明书 NO.5plot(s)title(原始 语音信号)subplot(212)plot(s1)title(截短 语音信号);会得到下图：图 3 原始语音信号和截短语音信号(6) 对语音信号进行频谱分析，在 Matlab 中，利用函数
9、FFT()对信号进行快速傅里叶变换，得到信号的频谱特性，如图 4 所示wavwrite(s1,fs,s1.wav); %将被处理信号s1输出为语音文件 “s1.wav”S1=fft(s1);figure(2)subplot(311);plot(s1);title(截短 预处理语音信号)课程设计说明书 NO.6subplot(312)plot(abs(S1)title(预处 理语音信号频谱);subplot(313);k=0:4000;plot(k(1:4000)*1,abs(S1(1:4000);title(预处 理语音信号单边带频谱)图 4 截短预处理语音信号如图 4 所示，从右向左看，第一
10、个较大的波峰所在的频率即为 3db 截止频率，第二个波峰所对应的频率为通带截止频率 ，在图中可以读出 =610Hz；pf pf一般在 3db 截止频率右侧的波谷位置选择阻带截止频率 =750Hz。通带截止频stf率 、阻带截止频率 数值的确定，就可以确定滤波器的基本指标。pfstf课程设计说明书 NO.7图 4 里第二个图是信号 的 FFT 结果，即 是信号)(nx )()(nxDFTkX的实际频谱 采样，本设计中信号 的长度取的是)(nx )(DFTejXwL=8000 点，所以，图中， 的每两个相邻点之间的频率间隔大小，即频率分(k辨率： （Hz） ，所以，根据它的放大图，即图 4 中18
11、0信 号 的 长 度采 样 频 率f的第三个图的放大图，可以确定。HzffHzffL 76076,24240点点2.7 滤波器的设计2.7.1 滤波器的参数设定本次课设我设计的是一个线性 FIR 低通滤波器，利用的窗函数是 hamming 窗，如上（图 4）所述：所以： HzfzfHzf sstp 80,75,610通带截止频率为 (rad/sample)4791.62sppfw阻带截止频率为 (rad/sample)58.08ststf(rad/sample)19.47.059.pstw由于海明窗过渡带满足： N23.求得滤波器阶数 189Nstpstpc ff21217.0)(stpscf
12、fw942 1 N课程设计说明书 NO.8（1）给定所要求的频率响应函数 )(jdeH.|0( cjjd（2）求单位采样响应 )nhd )(jwddeHIDFTjn21)(Sac（3）海明窗 )(12os(46.05.)( nRNnw（4）滤波器的单位采样响应： )()(wnhd )(12cos46.05.)( nRNSanhc )()94cos(46.054.0)94(17.0 189 nRnnwSa c 2.7.2 滤波器的 MATLAB 仿真在 M 文件中继续编写代码，把计算出来的参数带入代码中。代码如下：%加噪声完成信号截取s, fs,bits=wavread(C:Program Fi
13、lesMATLAB71workzf.wav);s1=s(1:8000);sound(s1,fs,bits);figure(1);subplot(211)plot(s)title(原始 语音信号)subplot(212)plot(s1)title(截短 语音信号);wavwrite(s1,fs,s1.wav);S1=fft(s1);figure(2)subplot(311);课程设计说明书 NO.9plot(s1);title(截短 预处理语音信号)subplot(312)plot(abs(S1)title(预处 理语音信号频谱);subplot(313);k=0:4000;plot(k(1:4
14、000)*1,abs(S1(1:4000);title(预处 理语音信号单边带频谱)s2=awgn(s1,15); %完成加噪！15dbwavwrite(s2,fs,s2.wav);figure(3);subplot(211);plot(s2);title(加噪后 语音信号);subplot(212);S2=fft(s2);plot(abs(S1);title(加噪后信号频谱);figure(4)subplot(211);plot(s1);title(语 音信号);subplot(212);plot(s2);title(加噪后 语音信号);%滤波器完成相关参数配置wp=610*2*pi/8000;wst=750*2*pi/8000;wc=(wp+wst)/2;N=ceil(3.3*2*pi/(wst-wp)+1;r=(N-1)/2;hn1=fir1(N-1,wc/pi,low,hamming(N);%s3=conv(s2,hn1);wavwrite(s3,fs,s3.wav);S3=fft(s3);figure(5)freqz(hn1);title(滤 波器幅频特性与相频特性)figure(6)subplot(111)plot(hn1);title(滤 波器系统函数);课程设计说明书 NO.10figure(7)subplot(211)plot(s3)title(滤
