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