![课程设计(论文)-基于窗函数法的FIR数字高通滤波器设计.doc_第1页](http://file.renrendoc.com/FileRoot1/2019-11/16/bc412659-8f3f-4e83-beab-028e41cf8167/bc412659-8f3f-4e83-beab-028e41cf81671.gif)
![课程设计(论文)-基于窗函数法的FIR数字高通滤波器设计.doc_第2页](http://file.renrendoc.com/FileRoot1/2019-11/16/bc412659-8f3f-4e83-beab-028e41cf8167/bc412659-8f3f-4e83-beab-028e41cf81672.gif)
![课程设计(论文)-基于窗函数法的FIR数字高通滤波器设计.doc_第3页](http://file.renrendoc.com/FileRoot1/2019-11/16/bc412659-8f3f-4e83-beab-028e41cf8167/bc412659-8f3f-4e83-beab-028e41cf81673.gif)
![课程设计(论文)-基于窗函数法的FIR数字高通滤波器设计.doc_第4页](http://file.renrendoc.com/FileRoot1/2019-11/16/bc412659-8f3f-4e83-beab-028e41cf8167/bc412659-8f3f-4e83-beab-028e41cf81674.gif)
![课程设计(论文)-基于窗函数法的FIR数字高通滤波器设计.doc_第5页](http://file.renrendoc.com/FileRoot1/2019-11/16/bc412659-8f3f-4e83-beab-028e41cf8167/bc412659-8f3f-4e83-beab-028e41cf81675.gif)
已阅读5页,还剩15页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
郑州轻工业学院课程设计说明书题目:基于窗函数法的FIR数字高通滤波器设计姓名:院(系):专业班级:学号:指导教师:成绩:时间:2015 年6 月29日至2015年7月3日郑州轻工业学院课程设计任务书题目基于窗函数法的FIR数字高通滤波器设计专业、班级电子信息工程学号姓名主要内容、基本要求、主要参考资料等:主要内容:利用MATLAB软件读取一段含有噪声的 .wav格式的语音信号,然后基于FFT对该信号进行频谱分析;基于含噪语音信号的频谱确定滤波器的参数,利用窗函数法设计一个FIR数字高通滤波器,并利用所设计的滤波器对信号进行滤波处理。比较滤波前后语音信号的时域波形及频谱,分析滤波前后的语音变化。基本要求:1、基于含噪语音信号的频谱确定滤波器的参数;2、分别采用矩形窗、汉明窗和布莱克曼窗设计FIR数字高通滤波器;3、掌握利用wavread函数读取、播放 .wav格式语音信号的方法;4、对语音信号进行滤波,绘制滤波前后信号的时域波形及频谱;5、回放语音信号,分析滤波前后的语音变化。主要参考资料:1、从玉良.数字信号处理原理及其MATLAB实现M.北京:电子工业出版社.2009.72、胡广书.数字信号处理理论、算法与实现M.北京:清华大学出版社.2003,8完成期限:2010.6.222010.6.26指导教师签名:课程负责人签名:2013年6月22日目录摘要11 数字滤波器简介22 FIR数字滤波器设计22.1 FIR数字滤波器的原理22.2设计工具22.3用窗函数法设计FIR数字滤波器32.3.1常用窗函数32.4 FIR数字滤波器的一般设计步骤:53 MATLAB仿真设计64 主要程序75 仿真结果图116 总结13参考资料14基于窗函数法的FIR数字高通滤波器设计摘要:数字滤波器在图像处理、数据传输等场合具有广泛应用,其设计是信号处理的核心问题之一,软件实现数字滤波优势体现在滤波器参数的改变伴随滤波器性能的改变.阐述了数字滤波器的设计方法,讨论了线性相位的条件和幅度特性,并以窗函数在MATLAB软件中实现了滤波器的仿真设计。仿真表明,借助计算机软件设计数字滤波器的过程简便易行,滤波性能指标均能达到指定的要求!对数字滤波器的推广应用有一定的促进作用。关键词:数字滤波器;MATLAB;窗函数;高通滤波1 数字滤波器简介数字信号处理(Digital Signal Processing,简称DSP)是一门涉及许多学科而又广泛应用于许多领域的新兴学科。20世纪60年代以来,随着计算机和信息技术的飞速发展,数字信号处理技术应运而生并得到迅速的发展。在过去的二十多年时间里,数字信号处理已经在通信等领域得到极为广泛的应用。数字滤波器的种类很多,从功能上我们可以将其分为低通、高通、带通和带阻等滤波器(输入输出都为数字的滤波器称之为数字滤波器)数字滤波器是通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。(数字滤波器从单位脉冲响应分类,可以分为无限脉冲响应(IIR)滤波器和有限脉冲响应(FIR)滤波器(由数字信号处理的一般理论可知,IIR滤波器的特征是具有无限持续时间的冲激响应,而 FIR 滤波器的冲激响应只能持续一定的时间。2 FIR数字滤波器设计2.1 FIR数字滤波器的原理FIR滤波器的传递函数为:可得FIR滤波器的系统系数差分方程,当FIR滤波器的系数满足下列中心对称条件时:b(n)=b(N-1-n)或b(n)=-b(N一1一n) (22)线性相位滤波器的相位滞后和延迟在整个频带上是相等且不变的。对于一个N阶的线性相位FIR滤波器,群延迟为常数,即滤波后的信号简单地延迟常数个时间步长。这一特性使通带频率内信号通过滤波器后仍保持原有波形形状而无相位失真。2.2设计工具在MATLAB中有三种方法可以实现设计FIR数字滤波器:程序设计法、FDATool设计法和SPTool设计法,程序设计法具有设计简便、可移植性强等特点,故在本设计中采用的是程序设计法。而在程序设计方法中,又有很多方法可以实现FIR 滤波器的设计,如窗函数设计法、频率采样设计法和最优化设计法等,这些方法都比较快捷方便,其中最常用的是窗函数设计法。2.3用窗函数法设计FIR数字滤波器设所希望设计的滤波器传输函数为H(ei),h(n)是与其对应的单位脉冲响应,因此由已知的H(ei),求出h(n),经过Z变换可得到滤波器的系统函数。但一般情况下,H(ei)是逐段稳定的,在边界频率处有不连续点,因而h(n)是无限时宽的,且是非因果序列(但是从实现的角度来说,我们希望得到一个长度为 Z 的线性相位滤波器,因此只能通过对h(n)进行加窗得出。窗函数设计滤波器的基本思想,就是根据给定的滤波器技术指标,选择滤波器的阶数N和合适的窗函数w(n)。即用一个有限长度的窗口函数序列w(n)来截取一个无限长的序列h(n)获得一个有限长序列h(n),即H(n)=h(n)*w(n),这样我们用一个有限长的序列H(n)去代替h(n),肯定会引起误差,表现在频域就是通常所说的吉布斯效应该效应引起通带内和阻带内的波动性,尤其使阻带的衰减小,影响滤波器的性能。为了减小吉布斯效应,从原理上来说可以通过加大N,但实验表明,加大N虽然可以使H(w)过渡带变窄,但对带内波动以及阻带衰减并没有多大的影响。因此,我们希望能找到不同的窗函数,使其主瓣包含更多的能量,同时旁瓣幅度减小,从而使通带%阻带波动减小,加大阻带衰减。工程实际中常用的窗函数有6种,即矩形窗、三角形窗、汉宁窗、哈明窗、布莱克曼窗和凯泽窗,不同窗函数的过渡带和阻带衰减特性有所不同。2.3.1常用窗函数1汉宁(Hanning)窗汉宁窗函数的时域形式表示为:(3-1)频域形式为 (3-2)汉宁窗函数的最大旁瓣值比主瓣值低31dB,但是主瓣宽度比矩形窗函数的主瓣宽度增加了1倍,为8/N。hanning函数:生成汉宁窗调用方式:(1) w = hanning(n):输入参数n是窗函数的长度;输出参数w是由窗函数的值组成的n阶向量。注意:此函数不返回是零点的窗函数的首尾两个元素。(2) w = hanning(n,symmetric):与上面相类似。(3) w = hanning(n,periodic):此函数返回包括为零点的窗函数的首尾两个元素。图2.3.1汉宁窗及其频谱特性2汉明(Hamming)窗: 函数的时域形式可以表示为(3-3)频域形式为:(3-4)其中,为矩形窗函数的幅度频率特性函数。海明窗函数的最大旁瓣值比主瓣值低41dB,但它和汉宁窗函数的主瓣宽度是一样大的。Hamming函数:生成海明窗调用方式(1) w = hamming(n):输入参数n是窗函数的长度;输出参数w是由窗函数的值组成的n阶向量。(2) w = hamming(n,sflag):参数sflag用来控制窗函数首尾的两个元素值;其取值为symmetric或periodic;默认值为symmetric。3布莱克曼窗函数的时域形式可以表示为(3-5)它的频域特性为(3-6)其中,为矩形窗函数的幅度频率特性函数。布莱克曼窗函数的最大旁瓣值比主瓣值低57dB,但是主瓣宽度是矩形窗函数的主瓣宽度的3倍,为12/N。Blackman函数:生成海明窗调用方式(1) w = blackman (n):输入参数n是窗函数的长度;输出参数w是由窗函数的值组成的n阶向量。(2) w = blackman (n,sflag):参数sflag用来控制窗函数首尾的两个元素值;其取值为symmetric或periodic;默认值为symmetric。图2.3.2 布莱克曼窗及其频谱特性2.4 FIR数字滤波器的一般设计步骤:(1) 按照任务要求确定滤波器的指标;(2) 用因果稳定的离散线性时不变系统的系统函数逼近这一性能指标要求;(3) 利用有限精度算法实现该系统函数,包括选择运算结构、适当的字长和有效位处理的方法;(4) 进行实际技术实现,常采用软件、硬件或者软硬件结合的方法实现。3MATLAB仿真设计本文的设计利用MATLAB的FDAtool采用窗函数法来设计语音高通数字滤波器。在MATLAB command窗口下,执行fdatool命令,将会出现如下窗口步骤如下:选择一段时间较短的单声道音频,用计算机windows附件中的录音机将其他格式的语音转换成。Wav格式并命名为gaot保存到MATLAB的work文件夹中,也可以自行创建work文件夹。用MATLAB中的wavread()函数将.wav语音转换成数组格式数据,用wavwrite()函数将数组格式数据转换成.wav格式语音,利用MATLAB的FDAtool工具设计所需FIR滤波器系数h(n)。由于本文实验所用音频录制时的抽样频率为44.1kH,所以在设计h(n)系数时应采用44.1kH,其他参数设置为,Response选择Highpass,Design选择FIR的window通过频率Fp=44.5kH,截止频率Fs=20kH.将设计好的高通滤波器系数导出到工作空间,命名为high。4 主要程序x1,fs1,bits1=wavread(F:MATLAB课程设计workSHE.wav); %soundsc(x1,fs1);x2,fs2,bits2=wavread(F:MATLAB课程设计workHighpass.wav);%soundsc(x2,fs2);figure(1)subplot(2,1,1);plot(x1) %做原始语音信号以44.1k采样后的时域图形title(原始语音采样后时域信号); xlabel(时间轴 n);ylabel(幅值 A); subplot(2,1,2);plot(x2) %做加噪声语音信号以44.1k采样后的时域图形title(加噪声语音采样后时域信号); xlabel(时间轴 n); ylabel(幅值 A); pause;X1=fft(x1,4096);figure(2)subplot(2,1,1);f=fs1*(0:2047)/4096;plot(f,abs(X1(1:2048); title(原始信号频谱);X2=fft(x2,4096); subplot(2,1,2);plot(f,abs(X2(1:2048); title(加噪声信号频谱); pause;% fs3=44100; %设定采样频率% fp1=250;fp2=7500; %第一截止频率% fs1=180;fs2=7570; %第二截止频率 % As=60; %最小阻带衰减% Ws1=(fp1+fs1)/fs3;Ws2=(fp2+fs2)/fs3; %截止频率归一化处理% w=(fp1-fs1)/fs3; %求归一化过渡带% M=ceil(As-7.95)/(14.36*w); %计算所需滤波器的阶数% window=boxcar(M+1); %生成长度为M+1的矩形窗% %window=hamming(M+1); %生成长度为M+1的汉宁窗% %window=blackman(M+1); %生成长度为M+1的布莱克曼窗% b,a=fir1(M,Ws1,Ws2,window); %生成设计的fir滤波器% figure(3)% plot(b);% pause;% figure(4)% freqz(b); %绘制幅频和相频响应曲线%wp1=0.007*pi; ws1=0.017*pi;%tr_width=ws1-wp1;%求过渡带宽度%M=ceil(6.6*pi/tr_width)+1;%求得所需窗函数的长度%wc1=(ws1+wp1)/2; %求截止频率%hd=ideal_hp1(wc1,M);%求得理想带阻的冲击响应%w_ham=(hamming(M);%得到长度为M的汉宁窗%h=hd .* w_ham;%利用窗函数截短%figure(3) %freqz(hd); %理想高通滤波器的幅频响应频谱%figure(4) %freqz(h); %加窗之后的幅频响应频谱%pause;wp1=0.007*pi; ws1=0.017*pi;tr_width=ws1-wp1;%求过渡带宽度M=ceil(6.6*pi/tr_width)+1;%求得所需窗函数的长度wc1=(ws1+wp1)/2; %求截止频率hd=ideal_hp1(wc1,M);%求得理想带阻的冲击响应w_ham=(blackman(M);%得到长度为M的布莱克曼窗h=hd .* w_ham;%利用窗函数截短figure(3) freqz(hd); %理想高通滤波器的幅频响应频谱figure(4) freqz(h); %加窗之后的幅频响应频谱 pause;f2=filter(h,1,x2);F0=fft(f2,4096);figure(5)subplot(3,1,1);plot(f,abs(X1(1:2048); title(原始信号频谱);subplot(3,1,2);plot(f,abs(X2(1:2048); title(滤波前的频谱)xlabel(Hz);ylabel(fuzhi);subplot(3,1,3)plot(f,abs(F0(1:2048);title(滤波后的频谱)xlabel(Hz);ylabel(fuzhi);pause;figure(6)subplot(2,1,1);plot(x1) %做原始语音信号以44.1k采样后的时域图形title(原始语音采样后时域信号); xlabel(时间轴 n);ylabel(幅值 A); subplot(2,1,2);plot(f2) %做滤波之后语音信号以44.1k采样后的时域图形title(滤波之后语音采样后时域信号); xlabel(时间轴 n);ylabel(幅值 A); pause;sound(x1,44100); sound(f2,44100); %播放滤波后的语音信号functionhd=ideal_hp1(Wc,N) %compute the ideal highpassfiter unit pulse respondencehd(n) %wc: cutoff frequency %N: window
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 端午节周记4篇
- 学期学习计划范本集锦七篇
- 简易停车场施工合同范本(标准版)
- 个人土地租赁合同书范本(标准版)
- 牛皮纸腹膜合加工合同
- 馕坑维修合同
- 奶羊饲料开发票的合同
- 木工包工包料劳务合同
- 辽宁省辽阳市2025届高三生物二模考试试题含解析
- 广东省汕头市苏湾中学2024-2025学年初三毕业班第一次调研测试化学试题含解析
- 服务器硬件运维巡检报告模板8页
- 设备维修申请单(新)
- 120万吨劣质煤洁净加工项目可行性研究报告
- 实习前班会(共7页)
- 硫酸泄漏专项应急预案
- 盾构施工监测方案
- TCCES 6003-2021 预制混凝土构件用金属预埋吊件
- 第一讲-视频拍摄入门(上)PPT优秀课件
- 乐器名称中英文对照
- 《四大名著》中高考常考的内容
- 冷冻式干燥机用户手册(2018.07版)(说明书)
评论
0/150
提交评论