2022年窗函数法设计方案低滤波器 .pdf
《2022年窗函数法设计方案低滤波器 .pdf》由会员分享,可在线阅读,更多相关《2022年窗函数法设计方案低滤波器 .pdf(17页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、摘要此次课程设计主要是要采集一段语音信号,并用MATLAB 软件绘制出语音信号波形并对语音信号进行截短、预处理等操作,观察其波形的变化并与原始语音信号波形加以对比。然后对该语音信号加入高斯白噪声,观察加噪后的波形及其频谱图,观察噪声对语音信号的影响。最后根据给定的相应技术指标,用汉宁窗设计一个满足指标的FIR 低通滤波器,对该语音信号进行滤波去噪处理,最后对前后时域和频域的波形图进行对比分析,从波形可以看出噪声被完全滤除,达到了语音不失真的效果。同时在课设过程中,通过与同组的其他同学交流,比较各种滤波方法性能的优劣并找到相对的最佳滤波方法。在此次课程设计中,以Windows XP系统为操作平台
2、 . 本次设计在 MATLAB 环境下,用窗函数法设计FIR 滤波器。通过了解不同的窗函数方法及性能设计 FIR 低通滤波器,并对对所设计的滤波器进行分析比较,得出各种方法设计的滤波器的优缺点,从而正确的选择FIR 数字滤波器的窗函数及设计方法。关键词 :FIR 滤波器, MATLAB ,窗函数,汉宁窗目录前言 1 一 MATLAB 程序设计的基本方法2 二数字信号处理的基本理论与方法2 2.1 设计理论依据 2 三利用 MATLAB 采集语言信号并分析4 3.1 语音的录入与打开4 3.2 录制语音信号过程4 3.3 时域信号的 FFT分析 5 四用 MATLAB 环境采用窗函数法设计数字滤
3、波器6 4.1 数字滤波器概念及原理6 4.2 数字滤波器设计的过程7 4.3 低通 FIR 滤波器基本原理 7 4.4 利用窗函数法设计线性相位FIR 数字滤波器 7 4.5 利用HANNING窗设计低通滤波器 9 五设计过程 10 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页,共 17 页5.1 流程图 10 5.3 信号频谱分析 12 5.4 信号的脉冲响应 12 设计总结 13 参考文献 14 附录错误!未定义书签。致谢 14 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页,共
4、 17 页前言数字化是控制系统的重要发展方向,而数字信号处理已在通信、语音、图像、自动控制、雷达、军事、航空航天等领域广泛应用。数字信号处理方法通常涉及变换、滤波、频谱分析、编码解码等处理。数字滤波是重要环节,它能满足滤波器对幅度和相位特性的严格要求,克服模拟滤波器所无法解决的电压和温度漂移以及噪声等问题。而有限冲激响应FIR 滤波器在设计任意幅频特性的同时能够保证严格的线性相位特性。利用MATLAB 工具软件的辅助设计,使得 FIR 滤波器具有快速、灵活、适用性强,硬件资源耗费少等特点。FIR 滤波器是最常用的组件之一,它完成信号预调、频带选择和滤波等功能。FIR 滤波器在截止频率的边沿陡峭
5、性能虽然不及IIR 滤波器,但是,考虑到FIR 滤波器严格的线性相位特性和不像IIR 滤波器存在稳定性的问题,FIR 滤波器能够在数字信号处理领域得到广泛的应用。 FIR 是有限冲激响应(Finite Impulse Response )的简称。由线性系统理论可知,在某种适度条件下,输入到线性系统的一个冲击完全可以表征系统。当我们处理有限的离散数据时,线形系统的响应(包括对冲击的响应)也是有限的。若线性系统仅是一个空间滤波器,则通过简单地观察它对冲击的响应,我们就可以完全确定该滤波器。通过这种方式确定的滤波器称为有限冲击响应(FIR)滤波器。 FIR 滤波器是在数字信号处理(DSP )中经常使
6、用的两种基本的滤波器之一。FIR 滤波器具有严格的相位特性,对于信号处理和数据传输是很重要的。目前FIR 滤波器的设计方法主要有三种:窗函数法、频率取样法和切比雪夫等波纹逼近的最优化设计方法。常用的是窗函数法和切比雪夫等波纹逼近的最优化设计方法。用窗函数设FIR 滤波器的基本思路:从时域出发设计h(n) 逼近理想hd(n) 。设理想滤波器的单位响应在时域表达为hd(n), 则 Hd(n) 一般是无限长的,且是非因果的,不能直接作为FIR 滤波器的单位脉冲响应。要想得到一个因果的有限长的滤波器单位抽样响应h(n) 最直接的方法是先将hd(n) 往右平移,再进行截断,即截取为有限长因果序列:h(n
7、)=hd(n)w(n), 并 用 合 适 的窗 函 数 进 行 加 权 作为FIR 滤波 器 的 单 位 脉 冲 响应 。MATLAB 设计FIR 滤波器有多种方法和对应的函数。窗函数设计法不仅在数字滤波器的设计中占有重要的地位,同时可以用于功率谱的估计,从根本上讲,使用窗函数的目的就是消除由无限序列的截短而引起的Gibbs 现象所带来的影响。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页,共 17 页一 MATLAB程序设计的基本方法MATLAB 是矩阵实验室( Matrix Laboratory)的简称,是美国MathWorks公司出品
8、的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境,主要包括MATLAB 和 Simulink两大部分。 MATLAB 和 Mathematica 、Maple 并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。它集图示和精确计算于一身,在应用数学、物理、化工、机电工程、医药、金融和其他需要进行复杂数值计算的领域得到广泛应用。它不仅是一个在各类工程设计中便于使用的计算工具,而且也是一个在数学、数值分析和工程计算等课程教案中的优秀的教案工具,在世界各地的高等院校中十分流行,在各类工业应用中更有不俗的表现。MATLAB 可以在几乎所有的P
9、C机和大型计算机上运行,适用于Windows 、UNIX等各种系统平台。总的来说,该软件有三大特点。一是功能强大。具有数值计算和符号计算、计算结果和编程可视化、数学和文字统一处理、离线和在线计算等功能;二是界面友善、语言自然。MATLAB 以复数处理作为计算单元,指令表达与标准教科书的数学表达式相近;三是开放性强。当学好 MATLAB 的同时,会更好的帮助自己去就解决一些难题,而且MATLAB 拥有非常好的发展前途,对我们未来的帮助也是不可限量的。二 数字信号处理的基本理论与方法2.1 设计理论依据精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4
10、 页,共 17 页图 2.1 原理示意图2.1.1 采样定理在进行模拟 / 数字信号的转换过程中,当采样频率fs.max 大于信号中,最高频率fmax的 2 倍时,即: fs.max=2fmax, 则采样之后的数字信号完整地保留了原始信号中的信息,一般实际应用中保证采样频率为信号最高频率的510 倍;采样定理又称奈奎斯特定理: 理想低通信道的最高大码元传输速率=2W*log2 N ( 其中 W是理想低通信道的带宽,N 是电平强度 ) 2.1.2 采样频率采样频率(也称为采样速度或者采样率)定义了每秒从连续信号中提取并组成离散信号的采样个数,它用赫兹(Hz)来表示。采样频率的倒数是采样周期或者叫
11、作采样时间,它是采样之间的时间间隔。采样频率只能用于周期性采样的采样器,对于非周期性采样的采样器没有规则限制。采样频率的常用的表示符号是 f_s 。2.1.3 采样位数与采样频率采样位数可以理解为声卡处理声音的解读度。这个数值越大,解读度就越高,录制和回放的声音就越真实。电脑中的声音文件是用数字0 和 1 来表示的,所以在电脑上录音的本质就是把模拟声音信号转换成数字信号。反之,在播放时则是把数字信号还原成模拟声音信号输出。采样频率是指录音设备在一秒钟内对声音信号的采样次数,采样频率越高声音的还原就越真实越自然。采样位数和采样频率对于音频接口来说是最为重要的两个指标,也是选择音频接口的两精选学习
12、资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页,共 17 页个重要标准。无论采样频率如何,理论上来说采样的位数决定了音频数据最大的力度范围。每增加一个采样位数相当于力度范围增加了6dB。采样位数越多则捕捉到的信号越精确。采样率越高,计算机摄取的图片越多,对于原始音频的还原也越加精确。三 利用 MATLAB 采集语言信号并分析3.1 语音的录入与打开利用 MATLAB 中的 wavread 命令来读入(采集)语音信号,将它赋值给某一向量。再将该向量看作一个普通的信号,对其进行FFT变换实现频谱分析,再依据实际情况对它进行滤波。对于波形图与频谱图(包括
13、滤波前后的对比图)都可以用MATLAB 画出。我们还可以通过 sound 命令来对语音信号进行回放,以便在听觉上来感受声音的变化。y,fs,bits=wavread(Blip,N1 N2),用于读取语音,采样值放在向量y 中,fs 表示采样频率 (Hz) ,bits表示采样位数。 N1 N2 表示读取从 N1点到 N2点的值(若只有一个N的点则表示读取前N点的采样值)。 sound(x,fs,bits)。 用于对声音的回放。向量y 则就代表了一个信号(也即一个复杂的“函数表达式”)也就是说可以像处理一个信号表达式一样处理这个声音信号。3.2 录制语音信号过程打开 PC 机上的酷狗音乐,用鼠标右
14、击播放列表中的任意一首歌曲,在弹出的下拉列表中单击工具,在工具栏选择制作铃声,然后显示如图3.1 所示的界面,点击添加歌曲,选择一首喜欢的歌,设置起点和终点,并将格式转换为WAV 格式,最后进行保存且与所写的程序对应,语音录音完成。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页,共 17 页图 3.1 录制语音信号图3.3 时域信号的 FFT分析函数的一种调用格式为:y=fft(x)(3.1 )其中, x 是序列, y 是序列的 FFT ,x 可以为一向量或矩阵,若x 为一向量, y 是 x 的FFT 。且和 x 相同长度。若 x 为一矩阵
15、,则 y 是对矩阵的每一列向量进行FFT 。如果 x 长度是 2 的幂次方,函数 fft执行高速基 2FFT算法;否则 fft执行一种混合基的离散傅立叶变换算法,计算速度较慢。函数 FFT的另一种调用格式为:y=fft(x,N)(3.2 )式中, x,y 意义同前, N为正整数。函数执行N点的 FFT 。若 x 为向量且长度小于N,则函数将 x 补零至长度 N 。若向量 x 的长度大于 N ,则函数截短 x 使之长度为 N。若 x 为矩阵,按相同方法对x 进行处理。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页,共 17 页经函数 fft求
16、得的序列 y 一般是复序列,通常要求其幅值和相位。MATLAB 提供求复数的幅值和相位函数: abs,angle ,这些函数一般和 FFT同时使用。四 用 MATLAB 环境采用窗函数法设计数字滤波器4.1 数字滤波器概念及原理4.1.1 数字滤波器概念数字滤波器是一种用来过滤时间离散信号的数字系统,通过对抽样数据进行数学处理来达到频域滤波的目的。根据其单位冲激响应函数的时域特性可分为两类:无限冲激响应(IIR )滤波器和有限冲激响应(FIR)滤波器。与IIR滤波器相比, FIR 的实现是非递归的,总是稳定的;更重要的是,FIR 滤波器在满足幅频响应要求的同时,可以获得严格的线性相位特性。因此
17、,它在高保真的信号处理,如数字音频、图像处理、数据传输、生物医学等领域得到广泛应用。4.1.2 数字滤波器原理数字滤波器的作用是利用离散时间系统的特性对输入信号波形( 或频谱 ) 进行加工处理,或者说利用数字方法按预定的要求对信号进行变换。数字滤波器可以理解为是一个计算程序或算法,将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列,并在转化过程中,使信号按预定的形式变化。数字滤波器有多种分类,根据数字滤波器冲激响应的时域特征,可将数字滤波器分为两种,即无限长冲激响应( IIR )滤波器和有限长冲激响应(FIR)滤波器。从性能上来说,IIR滤波器传输函数的极点可位于单位圆内的任何地方,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2022年窗函数法设计方案低滤波器 2022 函数 设计方案 滤波器
限制150内