语音信号滤波去噪——使用双线性变换法设计的切比雪夫II型滤波器
课程设计任务书
语音信号滤波去噪——使用脉冲响应不变法设计的巴特沃斯滤波器
摘要本课程设计主要运用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析，根据结果和学过的理论得出合理的结论。
关键词课程设计；滤波去噪；巴特沃斯滤波器；脉冲响应不变法；MATLAB
1 引言
本课程设计主要利用麦克风采集一段8000Hz，8k的单声道语音信号，并绘制波形观察其频谱，再用MATLAB利用脉冲响应不变法设计一个巴特沃斯滤波器，将该语音信号进行滤波去噪处理。
1.1 课程设计目的
《数字信号处理》课程设计是在学生完成数字信号处理和MATLAB的结合后的基本实验以后开设的。本课程设计的目的是为了让学生综合数字信号处理和MATLAB并实现一个较为完整的小型滤波系统。这一点与验证性的基本实验有本质性的区别。开设课程设计环节的主要目的是通过系统设计、软件仿真、程序安排与调试、写实习报告等步骤，使学生初步掌握工程设计的具体步骤和方法，提高分析问题和解决问题的能力，提高实际应用水平。
1.2课程设计的要求
（1）滤波器指标必须符合工程设计。
（2）设计完后应检查其频率响应曲线是否满足指标。
（3）处理结果和分析结论应该一致，而且应符合理论。
（4）独立完成课程设计并按要求编写课程设计报告。
2 设计原理
用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析。
2.1 IIR滤波器
I IR滤波器设计方法有间接法和直接法，间接法是借助于模拟滤波器的设计方法进行的。其设计步骤是：先设计过渡模拟滤波器得到系统函数H（s），然后将H（s）按某种方法转换成数字滤波器的系统函数H(z)。FIR滤波器比鞥采用间接法，常用的方法有窗函数法、频率采样发和切比雪夫等波纹逼近法。对于线性相位滤波器，经常采用FIR 滤波器。
对于数字高通、带通滤波器的设计，通用方法为双线性变换法。可以借助于模拟滤波器的频率转换设计一个所需类型的过渡模拟滤波器，再经过双线性变换将其转换策划那个所需的数字滤波器。具体设计步骤如下：
（1）确定所需类型数字滤波器的技术指标。
（2）将所需类型数字滤波器的边界频率转换成相应的模拟滤波器的边界频率，转换公式为Ω=2/T tan(0.5ω)
(3)将相应类型的模拟滤波器技术指标转换成模拟低通滤波器技术指标。
（4）设计模拟低通滤波器。
（5）通过频率变换将模拟低通转换成相应类型的过渡模拟滤波器。
（6）采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。
我们知道，脉冲响应不变法的主要缺点是会产生频谱混叠现象，使数字滤波器的频响偏离模拟滤波器的频响特性。为了克服之一缺点，可以采用双线性变换法。
下面我们总结一下利用模拟滤波器设计IIR数字低通滤波器的步骤：
（1）确定数字低通滤波器的技术指标：通带边界频率、通带最大衰减，阻带截止频率、阻带最小衰减。
（2）将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。
（3）按照模拟低通滤波器的技术指标设计及过渡模拟低通滤波器。
（4）用双线性变换法，模拟滤波器系统函数转换成数字低通滤波器系统函数。
设IIR滤波器的输入序列为，则IIR滤波器的输入序列和输出序列之间的关系可用下列方程式表示：
(2-1)
其中和是滤波器的系数，且中至少有一个非零。与上述差分方程对应的传递
函数为：
(2-2) 由传递函数可以发现无限长单位冲激响应滤波器有如下特点：
（1）单位冲激响应是无限长的。
（2）系统传递函数在有限z平面上有极点存在。
（3）结构上存在着输出到输入的反馈，也就是结构上是递归型的。
设是某个模拟滤波器的频率响应，那么低通滤波器在幅度平方响应上的技术指标给出为：
（2-3）
（2-4）式中是通带波纹参数，是通带截止频率以rad/s（弧度每秒）计，A是阻带衰减参数，以及是阻带截止频率以rad/s计。
（2-5）
（2-6）
参数和A是分别与以dB计的参数和有关的，这些关系是
（2-7）
（2-8）波纹和的绝对标尺是通过下式与和A有关的：
利用幅度平方响应给出的模拟滤波器要求（2-1）和（2-2）式不包含任何相位信息。现在，为了求s域的系统函数，考虑
那么有
（2-12）或者有
（2-9）
因此，幅度平方函数的零点和极点相对于轴是以镜像对称方式分布的【2】。
接下来，我们根据语音信号的特点给出有关滤波器的技术指标：
低通滤波器的性能指标：
fp=1200Hz，fs=1200Hz，ft=8000Hz, As=20db ,Ap=1dB
2.2 巴特沃斯滤波器
MATLAB信号处理工具箱函数buttp buttor butter是巴特沃斯滤波器设计函数，其有5种调用格式，本课程设计中用到的是[N,wc]=butter(N,wc,Rp,As,’s’),该格式用于计算巴特沃斯模拟滤波器的阶数N和3dB截止频率wc。
函数butter,cheby1和ellip设计IIR滤波器时都是默认的双线性变换法，所以在设计滤波器时只需要代入相应的实现函数即可。下面我们设计的IIR数字滤波器的主要程序：
wp=2*pi*Fp/Ft;
ws=2*pi*Fs/Ft;
fp=2*Ft*tan(wp/2);
fs=2*Ft*tan(ws/2);
[n,wn]=buttord(wp,ws,Ap,As,'s');
[b,a]=butter(n,wn,'s');
[num,den]=bilinear(b,a,1);
[h,w]=freqz(num,den);
3.设计步骤
3.1设计流程图
语音信号滤波去噪——使用脉冲不变响应法设计的巴特沃斯滤波器的设计流程如图2所示：
图1 脉冲响应不变法巴特沃斯滤波器对语音信号去噪流程图
3.2语言信号的采集
图3是基于PC机的语音信号采集过程:通过计算机录音系统录制1秒：
图2 语音采集过程
图3 录音过程截图3.3语音信号的频谱分析
3.4滤波器设计
将数字滤波器的设计指标设为通带截止频率fb=1100HZ,阻带频率fc=1200HZ,通带波纹Ap=1dB，阻带波纹As=20dB,要求确定H(z)。设计程序如下：
Fp=1200;
Ft=8000;
As=20;
Ap=1;
wp=2*pi*Fp/Ft;
ws=2*pi*Fs/Ft;
fp=2*Ft*tan(wp/2);
fs=2*Ft*tan(ws/2);
[n,wn]=buttord(wp,ws,Ap,As,'s');
[b,a]=butter(n,wn,'s');
[num,den]=bilinear(b,a,1);
[h,w]=freqz(num,den);
Buttor波形如下图5：
图5 IIR低通滤波器
3.5 信号滤波过程
[y,fs,nbits]=wavread('3333.wav');%载入波形N=length(y)%计算输入波形长度
Y=fft(y,N);%进行快速傅里叶变换
sound(y);%播放声音
figure（4）；%时域图
figure（5）；plot（abs（Y））；%滤波前频谱图
Fp=1200;%阻带截止频率
Fs=1100;%通带截止频率
Ft=8000;%采集频率
As=20;%通带波纹Ap=1
Ap=1;%阻带波纹As=20
wp=2*pi*Fp/Ft;
ws=2*pi*Fs/Ft;
fp=2*Ft*tan(wp/2);
fs=2*Ft*tan(ws/2);
[n,wn]=buttord(wp,ws,Ap,As,'s'); %求低通滤波器的阶数和截止频率[b,a]=butter(n,wn,'s'); %求S域的频率响应的参数
[num,den]=bilinear(b,a,1); %双线性变换实现S域到Z域的变换[h,w]=freqz(num,den); %根据参数求出频率响应
figure(3)
plot(w*8000*0.5/pi,abs(h));
z=filter(num,den,y);
sound(z);
m=z; %求滤波后的信号
figure(1)
subplot(2,2,3);
plot(abs(m),'r');
title('滤波后信号的频谱');
grid;
subplot(2,2,4);
plot(z,'b');
title('滤波后的信号波形');
grid;
subplot(2,2,2);
plot(y,'b');
title('滤波前信号的波形');
grid;
subplot(2,2,1);
plot(abs(Y),'r');
title('滤波前信号的频谱');
grid;
figure(2);
p=angle(m);
q=angle(Y);
subplot(2,1,1);plot(q,'b');title('滤波前相位');
grid;
subplot(2,1,2);plot(p,'b');title('滤波后相位');
grid;
所得波形为：
图6.1滤波前后的时域对比波形和滤波前后频域频谱图
图6.2滤波前后相位比较图
3.5结果分析
由图6.1中滤波前后波形比较可看出，经过滤波后的波形比原波形的振幅有所减小，去除了很多由于噪声所产生的干扰；从滤波前后的频谱比较可以看出经过滤波后除了原本的声音外，中间由于噪声产生的频谱波形已经滤除；由图6.2滤波前后相位比较图可看出由于经过滤波，相位变得稀疏；经过MATLAB仿真，听滤波前后的声音，可以听出有明显的滤波效果。因此利用脉冲响应不变法设计的巴特沃斯滤波器已经达到了设计的要求。
4.出现的问题及解决方法
在这次的课程设计中我们设计的滤波器对语音信号的滤噪处理或多或少会因为我们的疏忽或者实践能力不够还是会遇到问题的，当遇到问题时我们不但运用自己学到的知识还得到同学和老师的帮助将其解决，不断完善和修改。设计中出现的问题及解决方法如下：
（1）设计的滤波器达不到要求，因为在采集语音信号时，录音机的音频格式应该是PCM 8000kHz，8位，单声道。若由非人发出的声音可能会被滤波器滤去，且声音文件最好在一秒左右。
（2）在运用Matlab设计滤波器时，当编辑完前面两条程序时无法放出声音，后来发现我们应当把采集的语音信号wav文件放到Matlab的work文件夹中。
（3）还要在滤波器性能曲线的wc处画一根竖线，这样更方便看出结果，其中wc 处线的确定还需计算出wb/pi的值。
（4）所有的时间波形横坐标都要化为时间，滤波前后频谱的横坐标应是频率，这样在观察通带截止频率和阻带截止频率时更加精确，误差较小。
（5）要区分清楚ft、fs、fp所代表的意思。
5.结束语
这次课程设计，给我留下了很深的印象。虽然只是短暂的一周，但在这期间，却让我受益匪浅。
通过这次课程设计，使我对语音信号有了全面的认识，对数字信号处理的知识又有了深刻的理解，在之前数字信号与处理的学习以及完成课后的作业的过程中，已经使用过MATLAB，对其有了一些基础的了解和认识。本次经过一周的课程设计，让我们更熟悉了数字信号处理的基本知识和MATLAB的m语言，让我们把课上的理论知识运用到实际中去，让我们更近一步地巩固了课堂上所学的理论知识，并能很好地理解与掌握数字信号处理中的基本概念、基本原理、基本分析方法。
经过这次的DSP课程设计，让我有机会将自己学到的理论知识运用到实际中，提高了自己的动手能力和思维能力。在课程设计中发现自己的不足，所以在今后的学习和生活中我们要更加努力，学习好我们的专业知识并要能运用到实际。
总的来说，通过这次设计使我了解了MATLAB的使用方法，学会分析滤波器的优劣和性能，提高了分析和动手实践能力。同时我相信，进一步加强对MATLAB的学习与研究对我今后的学习将会起到很大的帮助！
参考文献
[1] 程佩青数字信号处理清华大学出版社 2007年2月
[2] 从玉良数字信号处理原理及matlab实现电子工业出版社 2005 年7月
[3] 刘波 matlab信号处理电子工业出版社2006年1月
IIR滤波器语音去噪处理
课程设计(论文) 题目基于IIR数字滤波器的有噪语音信号的 处理 课程设计（论文）任务书学院：电气工程学院 题目：基于IIR数字滤波器的有噪语音 信号的处理
起止时间：2016年10月25日至16年11月20日 学生姓名： 专业班级： 指导教师： 教研室主任： 院长： 2016年11 月20 日
摘要：滤波器设计在数字信号处理中占有极其重要的地位，IIR滤波器是滤波器设计的重要组成部分。课题基于MATLAB有噪音语音信号处理的设计与实现，综合运用数字信号处理的理论知识对加噪声语音信号进行时域，频域分析和滤波。通过理论推导得出相应结论，再利用MATLAB 作为编程工具进行计算机实现。在设计实现的过程中，用巴特沃斯、切比雪夫和双线性变法设计IIR数字滤波器，并利用MATLAB 作为辅助工具完成设计中的计算与图形的绘制。通过对对所设计滤波器的仿真和频率特性分析，可知利用MATLAB信号处理工具处理工具箱可以有效快捷地设计IIR数字滤波器，结果的各项性能指标均达到指定要求。 关键词：MATLAB；IIR滤波器；有噪音语音信号
Abstract: filter design plays an important role in digital signal processing, IIR filter is an important part of the filter design. Research based on MATLAB design and implementation of the noise processing of speech signal, the theory of knowledge of the integrated use of digital signal processing in time domain of speech signal plus noise, frequency domain analysis and filtering. Through theoretical derivation corresponding conclusion using MATLAB as a programming tool for computer implementation. In the design process, with Butterworth, Chebyshev and bilinear method of IIR digital filter design, and use MATLAB as a tool Complete the drawing calculation and graphic design. Through the simulation and frequency characteristic analysis on the design of filter, MATLAB signal processing tools processing toolbox can effectively and quickly design IIR digital filter based on the results of the performance indicators to meet the specified requirements. Keywords: MATLAB; IIR filter; noisy speech signal
切比雪夫1型数字低通滤波器
目录 1. 数字滤波器的设计任务及要求 (2) 2. 数字滤波器的设计及仿真 (2) 2.1数字滤波器的设计 (3) 2.2数字滤波器的性能分析 (3) 3. 数字滤波器的实现结构对其性能影响的分析 (8) 3.1数字滤波器的实现结构一及其幅频响应 (10) 3.2数字滤波器的实现结构二及其幅频响应 (12) 3.3 数字滤波器的实现结构对其性能影响的小结 (12) 4. 数字滤波器的参数字长对其性能影响的分析 (13) 4.1数字滤波器的实现结构一参数字长及幅频响应特性变化 4.2数字滤波器的实现结构二参数字长及幅频响应特性变化 4.3 数字滤波器的参数字长对其性能影响的小结 (16) 5. 结论及体会 (16) 5.1 滤波器设计、分析结论 (16) 5.2 我的体会 (16) 5.3 展望 (16)
1.数字滤波器的设计任务及要求 1. 设计说明 每位同学抽签得到一个四位数，由该四位数索引下表确定待设计数字滤波器的类型及其设计方法，然后用指定的设计方法完成滤波器设计。 要求：滤波器的设计指标： 低通： (1)通带截止频率πrad (id) pc 32 ln = ω (2)过渡带宽度πrad ) (i d 160 10log tz ≤?ω (3)滚降dB αroll 60= 其中，i d — 抽签得到那个四位数（学号的最末四位数），本设计中i d =0201。 2. 滤波器的初始设计通过手工计算完成； 3. 在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种以上合适的滤波器 结构进行分析）； 4. 在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响； 5. 以上各项要有理论分析和推导、原程序以及表示计算结果的图表； 6. 课程设计结束时提交设计说明书。 2.数字滤波器的设计及仿真 2.1数字滤波器（编号0201）的设计 数字滤波器是数字信号处理的重要工具之一，它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序，而数字滤波器处理精度高、体积小、稳定、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。 本次课程设计使用MATLAB 信号处理箱和运用切比雪夫法设计数字滤波器，将手工计算一个切比雪夫I 型的IIR 的低通模拟滤波器的系统函数，并在MATLAB 的FDATool 设计工具分析其性能指标。
滤波器语音信号去噪讲解
******************* 实践教学 ******************* 兰州理工大学 计算机与通信学院 2013年春季学期 信号处理课程设计 题目：基于MATLAB的FIR滤波器语音信号去噪专业班级：通信工程（1）班 姓名：王兴栋 学号：10250114 指导教师：陈海燕 成绩：
摘要 语音信号在数字信号处理中占有极其重要的地位，因此选择通过对语音信号的研究来巩固和掌握数字信号处理的基本能力十分具有代表性。对数字信号处理离不开滤波器，因此滤波器的设计在信号处理中占有极其重要的地位。而MATLAB 软件工具箱提供了对各种数字滤波器的设计。本论文“在MATLAB平台上实现对语音信号的去噪研究与仿真”综合运用了数字信号处理的各种基本知识，进而对不带噪语音信号进行谱分析以及带噪语音信号进行谱分析和滤波处理。通过理论推导得出相应的结论，再通过利用MATLAB作为编程工具来进行计算机实现比价已验证推导出来的结论。在设计过程中，通过设计FIR数字滤波器和IIR数字滤波器来完成滤波处理。在设计过程中，运用了MATLAB对整个设计中的图形的绘制和一些数据的计算以及仿真。 关键字滤波器；MATLAB；仿真；滤波
前言 语音是语言的声学表现，是人类交流信息最自然、最有效、最方便的手段。随着社会文化的进步和科学技术的发展，人类开始进入了信息化时代，用现代手段研究语音处理技术，使人们能更加有效地产生、传输、存储、和获取语音信息，这对于促进社会的发展具有十分重要的意义，因此，语音信号处理正越来越受到人们的关注和广泛的研究。 语音信号是信息技术处理中最重要的一门科学，是人类社会几步的标志。那么什么是语音？语音是人类特有的功能，也是人类获取外界信息的重要工具，也是人与人交流必不可少的重要手段。那么什么又是信号？那信号是什么呢？信号是传递信息的函数。离散时间信号——序列——可以用图形来表示。 语音信号处理是一门用研究数字信号处理研究信号的科学。它是一新兴的信息科学，同时又是综合多个学科领域的一门交叉科学。语音在我们的日常生活中随时可见，也随处可见，语音很大程度上可以影响我们的生活。所以研究语音信号无论是在科学领域上还是日常生活中都有其广泛而重要的意义。 本论文主要介绍的是的语音信号的简单处理。本论文针对以上问题，运用数字信号学基本原理实现语音信号的处理，在matlab7.0环境下综合运用信号提取，幅频变换以及傅里叶变换、滤波等技术来进行语音信号处理。我所做的工作就是在matlab7.0软件上编写一个处理语音信号的程序，能对语音信号进行采集，并对其进行各种处理，达到简单语音信号处理的目的。 对语音信号的研究，本论文采用了设计两种滤波器的基本研究方法来达到研究语音信号去噪的目的，最终结合图像以及对语音信号的回放，通过对比，得出结论。
切比雪夫Ⅰ型滤波器和切比雪夫Ⅱ型滤波器的对比
设计流程图如下： 设计思想： 首先设计一个源信号和一个混合信号，通过其频谱对比得出最大和最小通带，最大和最小阻带;然后再根据得到的参数来设计切比雪夫滤波器，最后通过切比雪夫Ⅰ型滤波器和切比雪夫Ⅱ型滤波器的对比来得出那种效果好。 切比雪夫滤波器设计原理：
切比雪夫滤波器的振幅具有等波纹特性，它有两种形式：1)振幅特性在通带内是等波纹的、在阻带内是单调的切比雪夫I 型滤波器；2)振幅特性在通带内是单调的、在阻带内是等波纹的切比雪夫II型滤波器，采用何种形式的切比雪夫滤波器取决于实际用途. 切比雪夫滤波器的设计方法就是将逼近精确度均匀分布在整个通带内，或者均匀分布在整个阻带内，或者均匀分布在两者之内，这样就可以使滤波器阶数大大降低。 切比雪夫I型滤波器平方幅度响应函数表示为: 2 ) (Ωj G=[1+2εC2 N (Ω)]2/1- 其中ε<1(正数)，它与通带波纹有关,ε越大，波纹也越大；C N (Ω) 是切比雪夫多项式，它被定义为： C N (Ω)=cos(Narccos(Ω)),Ω≤1, C N (Ω)=cosh(Narcosh(Ω)),Ω>1. 而切比雪夫II型滤波器平方幅度响应函数表示为： ) (Ωj G2={1+2ε{ C2 N (Ω)/[2 N (Ω/c Ω)]2}}1- 其中ε<1(正数)，表示波纹变化情况；c Ω为截止频率；N为滤波器的阶次,也 是C N ( N Ω Ω/)的阶次。 源信号编码及其图形： t=-1:0.01:1 y=(cos(2*pi*10*t)+cos(2*pi*40*t)); N=length(y); fx=fft(y); df=100/N; n=0:N/2; f=n*df; subplot(2,1,1); plot(f,abs(fx(n+1))*2/N); grid; title('源波形频谱')
语音信号滤波去噪——使用双线性变换法设计的切比雪夫II型滤波器
课程设计任务书
语音信号滤波去噪——使用脉冲响应不变法设计的巴特沃斯滤波器 摘要本课程设计主要运用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析，根据结果和学过的理论得出合理的结论。 关键词课程设计；滤波去噪；巴特沃斯滤波器；脉冲响应不变法；MATLAB 1 引言 本课程设计主要利用麦克风采集一段8000Hz，8k的单声道语音信号，并绘制波形观察其频谱，再用MATLAB利用脉冲响应不变法设计一个巴特沃斯滤波器，将该语音信号进行滤波去噪处理。 1.1 课程设计目的 《数字信号处理》课程设计是在学生完成数字信号处理和MATLAB的结合后的基本实验以后开设的。本课程设计的目的是为了让学生综合数字信号处理和MATLAB并实现一个较为完整的小型滤波系统。这一点与验证性的基本实验有本质性的区别。开设课程设计环节的主要目的是通过系统设计、软件仿真、程序安排与调试、写实习报告等步骤，使学生初步掌握工程设计的具体步骤和方法，提高分析问题和解决问题的能力，提高实际应用水平。 1.2课程设计的要求 （1）滤波器指标必须符合工程设计。 （2）设计完后应检查其频率响应曲线是否满足指标。
（3）处理结果和分析结论应该一致，而且应符合理论。 （4）独立完成课程设计并按要求编写课程设计报告。 2 设计原理 用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析。 2.1 IIR滤波器 I IR滤波器设计方法有间接法和直接法，间接法是借助于模拟滤波器的设计方法进行的。其设计步骤是：先设计过渡模拟滤波器得到系统函数H（s），然后将H（s）按某种方法转换成数字滤波器的系统函数H(z)。FIR滤波器比鞥采用间接法，常用的方法有窗函数法、频率采样发和切比雪夫等波纹逼近法。对于线性相位滤波器，经常采用FIR 滤波器。 对于数字高通、带通滤波器的设计，通用方法为双线性变换法。可以借助于模拟滤波器的频率转换设计一个所需类型的过渡模拟滤波器，再经过双线性变换将其转换策划那个所需的数字滤波器。具体设计步骤如下： （1）确定所需类型数字滤波器的技术指标。 （2）将所需类型数字滤波器的边界频率转换成相应的模拟滤波器的边界频率，转换公式为Ω=2/T tan(0.5ω) (3)将相应类型的模拟滤波器技术指标转换成模拟低通滤波器技术指标。 （4）设计模拟低通滤波器。 （5）通过频率变换将模拟低通转换成相应类型的过渡模拟滤波器。 （6）采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。 我们知道，脉冲响应不变法的主要缺点是会产生频谱混叠现象，使数字滤波器的频响偏离模拟滤波器的频响特性。为了克服之一缺点，可以采用双线性变换法。 下面我们总结一下利用模拟滤波器设计IIR数字低通滤波器的步骤： （1）确定数字低通滤波器的技术指标：通带边界频率、通带最大衰减，阻带截止频率、阻带最小衰减。
语音信号滤波去噪
一、设计的目的和意义 数字滤波器和快速傅立叶变换(FFT)等是语音信号数字处理的理论和技术基础，是20世纪60年代形成的一系列数字信号处理的理论和算法。在数字信号处理中,滤波器的设计占有极其重要的地位。而其中，FIR数字滤波器和IIR数字滤波器是重要组成部分。Matlab具有功能强大、简单易学、编程效率高等特点,深受广大科技工作者的喜爱。特别是Matlab中还具有信号分析工具箱,所以对于使用者，不需要具备很强的编程能力,就可以方便地进行信号分析、处理和设计。利用Matlab中的信号处理工具箱，可以快速有效的设计各种数字滤波器。本论文基于Matlab语音信号处理的设计与实现，综合运用数字信号处理的相关理论知识，对加噪声语音信号进行时域、频域分析并滤波。而后通过理论推导得出相应结论，再利用Matlab作为编程工具进行计算机实现工作。 本次课程设计的课题为《基于DSP的语音信号滤波去噪》，运用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析，根据结果和学过的理论得出合理的结论。 二、设计原理： 2.1 巴特沃斯滤波器 巴特沃斯滤波器是电子滤波器的一种。巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑。巴特沃斯滤波器的特性是通频带内的频率响应曲线最大限度平坦，没有起伏，而在组频带则逐渐下降为零。在振幅的对数对角频率的波得图上，从某一边界角频率开始，振幅随着角频率的增加而逐步减少，趋向负无穷大。 其振幅平方函数具有如2-1式： 式中，N为整数，称为滤波器的阶数，N越大，通带和阻带的近似 性越好，过渡带也越陡。如下图2.1所示：
信号系统处理 语音信号滤波处理
数字信号处理实验报告 实验名称：基于MATLAB对语音信号 进行分析及滤波处理 院系：物联网工程学院 班级：电子信息工程1101 姓名：
一、实验目的 综合计算运用数字信号处理的理论知识进行频谱分析和滤波器设计，通过理论推导得出相应的结论，培养发现问题、分析问题和解决问题的能力。并利用MATLAB作为工具进行实现，从而复习巩固课堂所学的理论知识，提高对所学知识的综合应用能力，并从实践上初步实现对数字信号的处理。此外，还系统的学习和实现对语音信号处理的整体过程，从语音信号的采集到分析、处理、频谱分析、显示和储存。 二、实验要求 1.分析原始语音信号的时域特性和频谱特性。 2.设计一个IIR滤波器，用该滤波器对语音信号进行滤波处理，分析滤波后信号的时域特性和频谱特性。 3.设计一个FIR滤波器，用该滤波器对语音信号进行滤波处理，分析滤波后信号的时域特性和频谱特性。 三、实验原理 1.采样定理：在进行模拟/数字信号的转换过程中，当采样频率fs.max大于信号中最高频率fmax的2倍时(fs.max>=2fmax)，采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定理又称奈奎斯特定理。 采样频率越高，即采样的时间间隔越短，对声音波形的表示越精确。 2.时域信号的FFT分析：信号的频谱分析就是计算机信号的傅里叶变换。连续信号与系统的傅里叶分析显然不便于用计算机进行计算，使其应用受到限制。而FFT是一种时域和频域均离散化的变换，适合数值运算，成为用计算机分析离散信号和系统的有力工具。对连续信号和系统，可以通过时域采样，应用DFT进行近似谱分析。 3.IIR数字滤波器设计原理利用双线性变换设计IIR滤波器（巴特沃斯数字低通滤波器的设计），首先要设计出满足指标要求的模拟滤波器的传递函数Ha(s),然后由Ha(s)通过双线性变换可得所要设计的IIR滤波器的系统函数H(z)。如果给定的指标为数字滤波器的指标，则首先要转换成模拟滤波器的技术指标，这里主要是边界频率Wp和Ws的转换，对ap和as指标不作变化。边界频率的转换关系为∩=2/T tan(w/2)。接着，按照模拟低通滤波器的技术指标根据相应设计公式
切比雪夫Ⅱ型低通滤波器
一、设计一个切比雪夫Ⅱ型低通滤波器 wp=0.2*pi; %通带边界频率； ws=0.4*pi; %阻带截止频率； rp=1; %通带最大衰减； rs=80; %阻带最小衰减； Fs=1000 %假设抽样脉冲1000hz [N,Wn]=cheb2ord(wp,ws,rp,rs,'s'); %Chebyshev II型滤波器参数计算（模拟域）； [Z,P,K]=cheby2(N,rs,Wn,'s'); %构造Chebyshev II型滤波器（零极点模型）；[H,W]=zp2tf(Z,P,K); %将零极点模型转化成传递函数的模型； figure(1); freqs(H,W); %在Figure1上显示滤波器的幅频响应及相频响应；[P,Q]=freqs(H,W); %返回滤波器的冲击响应的复数形式； figure(2); plot(Q*Fs/(2*pi),abs(P));grid; %在Figure2上显示幅频特性曲线； xlabel('频率/Hz'); ylabel('幅值');
二、设计一个高通Chebyshow型数字滤波器 wp=100;ws=80;Fs=300;rp=1;rs=45; %数字滤波器的各项指标； WP=100*2*pi; %把数字滤波器的频率特征转换成模拟滤波器的频率特征； WS=300*2*pi; [N,Wn]=cheb2ord(WP,WS,rp,rs,'s'); %Chebyshev II型滤波器参数计算（模拟域）； [Z,P,K]=cheb2ap(N,rs); %创建Chebyshev滤波器原型； [A,B,C,D]=zp2ss(Z,P,K); %表达式从零极点增益形式转换成状态方程形式； [AA,BB,CC,DD]=lp2hp(A,B,C,D,Wn); %实现低通到高通滤波器类型的转换；[a,b,c,d]=bilinear(AA,BB,CC,DD,Fs); %采用双线性变换法，从模拟高通到数字高通； [P,Q]=ss2tf(a,b,c,d); %表达式从状态方程形形式转换成传输函数形式；figure(1); freqz(P,Q); %绘出频率响应； [H,W]=freqz(P,Q); figure(2); plot(W*Fs/(2*pi),abs(H));grid; xlabel('频率/Hz'); ylabel('幅值');
语音信号处理实验报告
通信与信息工程学院 信息处理综合实验报告 班级：电子信息工程1502班 指导教师： 设计时间：2018/10/22-2018/11/23 评语： 通信与信息工程学院 二〇一八年 实验题目：语音信号分析与处理 一、实验内容 1. 设计内容 利用MATLAB对采集的原始语音信号及加入人为干扰后的信号进行频谱分析，使用窗函数法设计滤波器滤除噪声、并恢复信号。 2．设计任务与要求 1. 基本部分
（1）录制语音信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。 （2）对所录制的语音信号加入干扰噪声，并对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。 （3）分别利用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman 窗几种函数设计数字滤波器滤除噪声，并画出各种函数所设计的滤波器的频率响应。 （4）画出使用几种滤波器滤波后信号时域波形和频谱，对滤波前后的信号、几种滤波器滤波后的信号进行对比，分析信号处理前后及使用不同滤波器的变化；回放语音信号。 2. 提高部分 （5）录制一段音乐信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。 （6）利用MATLAB产生一个不同于以上频段的信号；画出信号频谱图。 （7）将上述两段信号叠加，并加入干扰噪声，尝试多次逐渐加大噪声功率，对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。 （8）选用一种合适的窗函数设计数字滤波器，画出滤波后音乐信号时域波形和频谱，对滤波前后的信号进行对比，回放音乐信号。 二、实验原理 1.设计原理分析 本设计主要是对语音信号的时频进行分析，并对语音信号加噪后设计滤波器对其进行滤波处理，对语音信号加噪声前后的频谱进行比较分析，对合成语音信号滤波前后进行频谱的分析比较。 首先用PC机WINDOWS下的录音机录制一段语音信号，并保存入MATLAB软件的根目录下，再运行MATLAB仿真软件把录制好的语音信号用audioread函数加载入MATLAB仿真软件的工作环境中，输入命令对语音信号进行时域，频谱变换。 对该段合成的语音信号，分别用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman窗几种函数在MATLAB中设计滤波器对其进行滤波处理，滤波后用命令可以绘制出其频谱图，回放语音信号。对原始语音信号、合成的语音信号和经过滤波器处理的语音信号进行频谱的比较分析。 2.语音信号的时域频域分析 在Matlab软件平台下可以利用函数audioread对语音信号进行采样,得到了声音数据变量y,同时把y的采样频率Fs=44100Hz放进了MATALB的工作空间。
MATLAB对语音信号加随机噪声及去噪程序
%对语言信号做原始的时域波形分析和频谱分析[y,fs,bits]=wavread('C:\Documentsand?Settings\Administrator\桌面\cuocuo.wav'); %??sound(y,fs)??????%回放语音信号 n=length(y)??%选取变换的点数? y_p=fft(y,n);??????%对n点进行傅里叶变换到频域 f=fs*(0:n/2-1)/n;???%对应点的频率 figure(1) subplot(2,1,1); plot(y);????????????????????%语音信号的时域波形图 title('原始语音信号采样后时域波形'); xlabel('时间轴') ylabel('幅值A') subplot(2,1,2); plot(f,abs(y_p(1:n/2)));?????%语音信号的频谱图 title('原始语音信号采样后频谱图'); xlabel('频率Hz'); ylabel('频率幅值'); %对音频信号产生噪声 ??L=length(y)????????%计算音频信号的长度 ??noise=0.1*randn(L,2);??%产生等长度的随机噪声信号(这里的噪声的大小取决于随机函数的幅度倍数） ??y_z=y+noise;????????%将两个信号叠加成一个新的信号——加噪声处理??? ??%sound(y_z,fs) %对加噪后的语音信号进行分析 n=length(y);??%选取变换的点数? y_zp=fft(y_z,n);??????%对n点进行傅里叶变换到频域 f=fs*(0:n/2-1)/n;???%对应点的频率 figure(2) subplot(2,1,1); plot(y_z);????????????????????%加噪语音信号的时域波形图 title('加噪语音信号时域波形'); xlabel('时间轴') ylabel('幅值A') subplot(2,1,2); plot(f,abs(y_zp(1:n/2)));?????%加噪语音信号的频谱图 title('加噪语音信号频谱图'); xlabel('频率Hz'); ylabel('频率幅值');
根据MATLAB的加噪语音信号的滤波
《计算机仿真技术》 基于MATLAB的加噪语音信号的滤波学生姓名： 专业：电子信息工程 班级： 学号： 指导教师： 完成时间：2017年12月
一．滤波器的简述 在MATLAB环境下IIR数字滤波器和FIR数字滤波器的设计方法即实现方法，并进行图形用户界面设计，以显示所介绍迷你滤波器的设计特性。 在无线脉冲响应（IIR）数字滤波器设计中，先进行模拟滤波器的设计，然后进行模拟-数字滤波器转换，即采用脉冲响应不变法及双线性Z变化法设计数字滤波器，最后进行滤波器的频带转换。在有限脉冲响应（FIR）数字滤波器设计中，讨论了FIR线性相位滤波的特点和用窗口函数设计FIR数字滤波器两个问题。两类滤波器整个过程都是按照理论分析、编程设计、集体实现的步骤进行的。为方便分析直观者直观、形象、方便的分析滤波器的特性，创新的设计出图形用户界面---滤波器分析系统。整个系统分为两个界面，其内容主要包括四个部分：System（系统）、Analysis（分析）、Tool（工具）、Help（帮助）。 数字滤波在DSP中占有重要地位。数字滤波器按实现的网络结构或者从单位脉冲响应，分为IIR（无限脉冲响应）和FIR（有限脉冲响应）滤波器。如果IRR 滤波器和FIR滤波器具有相同的性能，那么通常IIR滤波器可以用较低的阶数获得高的选择性，执行速度更快，所有的存储单元更少，所以既经济又高效。二．设计要求 1.在matlab平台上录制一段语音信号； 2.完成语音信号的谱分析； 3.对语音信号进行加噪以及加噪后信号的谱分析； 4.选择合适的滤波器进行滤波，确定相关指标； 5.实现滤波过程，显示滤波后的结果，并进行谱分析。 三．实验内容与步骤 1、语音信号的录入
基于切比雪夫I型的高通滤波器设计Matlab
设计题目基于切比雪夫I型的数字高通滤波器的设计 设计要求 设计一个9阶切比雪夫I型高通滤波器,通带纹波为10dB,下边界频率为400 / rad s ,并绘出其幅频响应曲线 设计过程1．系统设计方案 1.1 Matlab的简介和主要功能： 简介：MATLAB 是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。使用 MATLAB，您可以较使用传统的编程语言（如 C、C++ 和 Fortran）更快地解决技术计算问题。 MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 MATLAB 函数集）扩展了 MATLAB 环境，以解决这些应用领域内特定类型的问题。 MATLAB 提供了很多用于记录和分享工作成果的功能。可以将您的 MATLAB 代码与其他语言和应用程序集成，来分发您的 MATLAB 算法和应用。 主要功能:1.此高级语言可用于技术计算 2.此开发环境可对代码、文件和数据进行管理 3.交互式工具可以按迭代的方式探查、设计及求解问题 4.数学函数可用于线性代数、统计、傅立叶分析、筛选、优化以及数值 积分等 5.二维和三维图形函数可用于可视化数据 6.各种工具可用于构建自定义的图形用户界面 7.各种函数可将基于 MATLAB 的算法与外部应用程序和语言（如 C、 C++、Fortran、Java、COM 以及 Microsoft Excel）集成 1.2 开发算法和应用程序 开发算法和应用程序 MATLAB 提供了一种高级语言和开发工具，使您可以迅速地开发并分析算法和应用程序。
切比雪夫低通滤波器讲解
课程设计 课程名称：数字信号处理 题目编号： 0202 题目名称：切比雪夫Ⅱ型IIR低通滤波器 专业名称：电子信息工程 班级：电子1204班 学号： 20124470411 学生姓名：刘春阳 任课教师：黄国玉 2015年09月30日
课程设计任务书
目录 1. 数字滤波器的设计任务及要求（编号202） (2) 2. 数字滤波器的设计及仿真 (3) 2.1数字滤波器（编号202）的设计 (3) 2.2数字滤波器（编号202）的性能分析 (6) 3. 数字滤波器的实现结构对其性能影响的分析 (7) 3.1数字滤波器的实现结构一（直接型）及其幅频响应 (8) 3.2数字滤波器的实现结构二（级联型）及其幅频响应 (10) 3.3 数字滤波器的实现结构对其性能影响的小结 (10) 4. 数字滤波器的参数字长对其性能影响的分析 (11) 4.1数字滤波器的实现结构一（直接型）参数字长及幅频响应特性变化 (12) 4.2数字滤波器的实现结构二（级联型）参数字长及幅频响应特性变化 (14) 4.3 数字滤波器的参数字长对其性能影响的小结 (16) 5. 结论及体会 (16) 5.1 滤波器设计、分析结论 (16) 5.2 我的体会 (16) 5.3 展望 (17)
1.数字滤波器的设计任务及要求（0202） 每位同学抽签得到一个四位数，由该四位数索引下表一确定待设计数字滤波器的类型及其设计方法， 然后用指定的设计方法完成滤波器设计。 要求： （1）滤波器设计指标：通带截止频 pc ln ()32 d rad i πω= , 过渡带宽度 10 tz () 160 log d rad i πω?≤ ，滚降roll 60dB α=; 其中， id —抽签得到那个四位数（题目编号) （2）滤波器的初始设计通过手工计算完成； （3）在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种 以上合适的滤波器结构进行分析）； （4）在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响； （5） 以上各项要有理论分析和推导、原程序以及表示计算结果的图表； （6）课程设计结束时提交设计说明书。 2.数字滤波器的设计及仿真 2.1 数字滤波器（编号202）的设计 随着信息和数字时代的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。在现代通信系统中，由于信号中经常混有各种复杂的成分，因此很多信号的处理都是基于滤波器而进行的。所以数字滤波器在数字信号处理中起着举足轻重的作用。而数字滤波器的设计都要以模拟滤波器为基础的，这是因为模拟滤波器的设计方法都已发展的相当成熟，且有典型的模拟滤波器供我们选择。如巴特沃斯滤波器、切比雪夫滤波器等。本次课程设计将手工计算一个切比雪夫II 型的IIR 的低通模拟滤波器的系统函数，并在MATLAB 的FDATool 设计工具分析其性能指标。
切比雪夫低通滤波器设计
摘要 随着信息与数字时代的到来,数字信号处理已成为当今一门极其重要的学科与技术领域。在现代通信系统中,由于信号中经常混有各种复杂成分,因此很多信号的处理都就是基于滤波器而进行的。所以,数字滤波器在数字信号处理中起着举足轻重的作用。而数字滤波器的设计都要以模拟滤波器为基础的,这就是因为模拟滤波器的理论与设计方方法都已发展的相当成熟, 且有典型的模拟滤波器供我们选择。,如巴特沃思滤波器、切比雪夫滤波器等。 本次课程设计将运用MATLAB设计一个基于切比雪夫低通滤波器,并出所设计滤波器的幅度及幅度衰减特性。 关键词:模拟低通滤波切比雪夫
1课题描述 数字滤波器就是数字信号处理的重要工具之一,它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序,而数字滤波器处理精度高、体积小、稳定、重量轻、灵活、不存在阻抗匹配问题,可以实现模拟滤波器无法实现的特殊功能。故本课题使用MATLAB 信号处理箱与运用切比雪夫法设计数字低通滤波器。 2设计原理 2、1切比雪夫滤波器介绍 在巴特沃兹滤波器中,幅度响应在通带与阻带内都就是单调的。因此,若滤波器的技术要求就是用最大通带与阻带的逼近误差来给出的话,那么,在靠近通带低频端与阻带截止频率以上的部分都会超出技术指标。一种比较有效的途径就是使逼近误差均匀地分布于通带或阻带内,或同时在通带与阻带内都均匀分布,这样往往可以降低所要求的滤波器阶次。通过选择一种具有等波纹特性而不就是单调特性的逼近方法可以实现这一点。切比雪夫型滤波器就具有这种性质:其频率响应的幅度既可以在通带中就是等波纹的,而在阻带中就是单调的(称为I 型切比雪夫滤波器),也可以在通带中就是单调的,而在阻带中就是等波纹的(称为II 型切比雪夫滤波器)。I 型切比雪夫滤波器的幅度平方函数就是 2|)(|Ωj H C =)/(11 22c N C ΩΩ+ε (2、1) 式中为N 阶切比雪夫多项式,定义为 )cos cos()(1x N x C N -= (2、2) 从定义切比雪夫多项式可以直接得出由 )(x C N 与)(1x C N -求)(1x C N +的递推公式。将三角恒等式代入 (2、2)式,得 )(1x C N +=2x -)(x C N ) (1x C N - (2、3) 从 (2、2)式我们注意到,当0
语音信号的处理与滤波
广西工学院 数字信号处理课程设计设计题目：语音信号的处理与滤波 系别： 学号： 姓名： 班级： 指导教师： 完成日期：
目录 1. 摘要 (3) 2.MATLAB简介 (3) 3.设计目的 (4) 4.设计内容 (4) 5.设计原理 (4) 6.设计步骤 (5) 7.总结与分析 (10) 8.参考资料 (10)
摘要 本课题分析了数字信号处理课程的重要性及特点，可以帮助理解与掌握课程中的基本概念、基本原理、基本分析方法；并利用MATLAB对语音信号进行分析和处理，要求采集语音信号后，在MATLAB软件平台进行频谱分析。用设计的数字滤波器对语音信号进行了滤波 MATLAB简介 MATLAB 是一种对技术计算高性能的语言。它集成了计算，可视化和编程于一个易用的环境中，在此环境下，问题和解答都表达为我们熟悉的数学符号。典型的应用有： ?数学和计算 ?算法开发 ?建模，模拟和原形化 ?数据分析，探索和可视化 ?科学与工程制图 ?应用开发，包括图形用户界面的建立 MATLAB是一个交互式的系统，其基本数据元素是无须定义维数的数组。这让你能解决很多技术计算的问题，尤其是那些要用到矩阵和向量表达式的问题。而要花的时间则只是用一种标量非交互语言(例如C或Fortran)写一个程序的时间的一小部分。 . 名称“MATLAB”代表matrix laboratory（矩阵实验室）。MATLAB最初是编写来提供给对由LINPACK和EINPACK工程开发的矩阵软件简易访问的。今天，MATLA B使用由LAPACK和ARPACK工程开发的软件，这些工程共同表现了矩阵计算的软件中的技术发展。 MATLAB已经与许多用户输入一同发展了多年。在大学环境中，它是很多数学类、工程和科学类的初等和高等课程的标准指导工具。在工业上，MATLAB是高产研究、开发和分析所选择的工具。 MATLAB以一系列称为工具箱的应用指定解答为特征。对多数用户十分重要的是，工具箱使你能学习和应用专门的技术。工具箱是是MATLAB函数（M-文件）的全面的综合，这些文件把MATLAB的环境扩展到解决特殊类型问题上。具有可用工具箱的领域有：信号处理，控制系统神经网络，模糊逻辑，小波分析，模拟等等。
巴特沃斯、切比雪夫滤波器的仿真与实现
华北科技学院课程设计任务书 2013 — 2014 学年第二学期 电子信息工程学院（系、部）通信工程专业 B111 班级课程名称：移动通信 设计题目：巴特沃斯、切比雪夫滤波器的仿真与实现完成期限：自16 周至 18 周共 3 周
目录 1．前言 (3) 1.1 MATLAB (3) 1.2 滤波器的概念 (5) 1.2.1滤波器的原理 (6) 1.2.2理想滤波器与实际滤波器 (6) 1.2.3 滤波器的分类 (7) 2.设计目的 (9) 3.设计原理 (9) 3.1.模拟低通滤波器的设计指标及逼近方法 (9) 3.2.巴特沃斯低通滤波器的设计方法 (10) 3.3.切比雪夫滤波器的设计方法 (14) 4.详细设计与系统分析 (21) 4.1程序设计 (21) 4.1.1巴特沃斯滤波器 (21) 4.1.2切比雪肤滤波器 (23) 4.2同一滤波器不同参数的比较 (25) 4.2.1巴特沃斯滤波器 (25) 4.2.2切比雪夫滤波器 (27) 4.3不同滤波器同一阶数的比较 (30) 4.3.1低通滤波器 (30) 4.3.2高通滤波器 (30) 4.3.3带通滤波器 (31) 4.3.4带阻滤波器 (31) 5.心得体会 (32) 6.参考文献 (32)
摘要：利用MATLAB设计滤波器,可以按照设计要求非常方便地调整设计参数,极大地减轻了设计的工作量,有利于滤波器设计的最优化。MATLAB因其强大的数据处理功能被广泛应用于工程计算,其丰富的工具箱为工程计算提供了便利,利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器,设计简单方便。本文介绍了在MATLAB R2009a 环境下滤波器设计的方法和步骤。关键词：滤波器，MATLAB 1．前言 1.1 MATLAB MATLAB是美国MathWorks公司开发的一种功能极其强大的高技术计算语言和内容极其丰富的软件库，集数值计算、矩阵运算和信号处理与显示于一身。该软件最初是由美国教授Cleve Moler 创立的。1980年前后，他在教线性代数课程时，发现用其他高级语言编程时极不方便，便构思开发了MATLAB，即矩阵实验室（Matrix Laboratory）。该软件利用了当时代表数值线性代数领域最高水平的EISPACK和LINPACK两大软件包，并且利用Fortran 语言编写了最初的一套交互式软件系统，MATLAB的最初版本便由此产生了。 最初的MATLAB由于语言单一，只能进行矩阵的运算，绘图也只能用原始的描点法，内部函数只有几十个，因此功能十分简单。1984年该公司推出了第一个MATLAB的商业版，并用C语言作出了全部改写。现在的MATLAB程序是MathWorks公司用C语言开发的，第一版由steve Bangert主持开发编译解释程序，Steve Kleiman完
基于MATLAB的语音信号滤波处理
基于MATLAB的语音信号滤波处理 题目：基于MATLAB的语音信号滤波处理 课程：数字信号处理 学院：电气工程学院 班级： 学生： 指导教师： 二O一三年十二月
目录CONTENTS 摘要 一、引言 二、正文 1.设计要求 2.设计步骤 3.设计内容 4.简易GUI设计 三、结论 四、收获与心得 五、附录
一、引言 随着Matlab仿真技术的推广，我们可以在计算机上对声音信号进行处理，甚至是模拟。通过计算机作图，采样，我们可以更加直观的了解语音信号的性质，通过matlab编程，调用相关的函数，我们可以非常方便的对信号进行运算和处理。 二、正文 2.1 设计要求 在有噪音的环境中录制语音，并设计滤波器去除噪声。 2.2 设计步骤 1.分析原始信号，画出原始信号频谱图及时频图，确定滤波器类型及相关指标； 2.按照类型及指标要求设计出滤波器，画出滤波器幅度和相位响应，分析该滤波器是否符合要求； 3.用所设计的滤波器对原始信号进行滤波处理，画出滤波后信号的频谱图及时频图； 4.对滤波前的信号进行分析比对，评估所设计滤波器性能。 2.3 设计内容 1.原始信号分析
分析信号的谱图可知，噪音在1650HZ和3300HZ附近的能量较高，而人声的能量基本位于1000HZ以下。因此，可以设计低通滤波器对信号进行去噪处理。 2.IIR滤波器设计 用双线性变换法分别设计了巴特沃斯低通滤波器和椭圆低通滤波器和带阻滤波器： ①巴特沃斯滤波器 fp=800;fs=1300;rs=35;rp=0.5; 程序代码如下： fp=800;fs=1300;rs=35;rp=0.5;Fs=44100; wp=2*Fs*tan(2*pi*fp/(2*Fs));ws=2*Fs*tan(2*pi*fs/(2*Fs)); [n,wn]=buttord(wp,ws,rp,rs,'s'); [b,a]=butter(n,wn,'s'); [num,den]=bilinear(b,a,Fs); [h,w]=freqz(num,den,512,Fs);
语音信号去噪
语音信号去噪 摘要：在现代各种通信系统中，由于自然界中的各种各样的复杂噪声不免会掺杂在其中，数 字信号处理这门经典学科恰好能够解决这个问题，其中最通用的方法就是利用滤波器来滤除 这些杂波噪声，FIR数字滤波器就是滤波器设计的基本部分。本论文研究的主要内容就是基 于Matlab软件仿真设计一个数字滤波器，将掺杂在语音信号中的高频噪音消除，在此将分 析消除高频噪音前后语音信号的时域及频域特性，对比分析即可验证滤波前后特性差别。在 本课题中，将利用简单的窗函数法来设计FIR数字滤波器，通过Matlab仿真说明所设计滤 波器的正确性。通过这次毕业设计，将会进一步理解语音信号原理分析及滤波处理，为更好 的设计滤波器打好基础。 关键词：Matlab；窗函数法；FIR数字滤波器 目录 1 引言 (2) 1.1 课题研究现状 (2) 1.2 课题研究目的 (2) 1.3 课题研究内容 (2) 1.4 MATLAB软件设计平台简介 (3) 2 原始语音信号采集与处理 (3) 2.1 课题设计步骤及流程图 (3) 2.2 语音信号处理 (4) 2.2.1 语音信号的采集 (4) 2.2.2 语音信号的时域频谱分析 (5) 2.2.3 语音信号加噪与频谱分析 (7) 3 FIR数字滤波器的设计 (9) 3.1 数字滤波器基本概念 (9) 3.2 常用窗函数介绍 (10) 3.3 FIR数字滤波器概述 (10) 3.4 FIR滤波器的窗函数设计 (10) 3.5 滤波器的编程实现 (13) 3.6 用滤波器对加噪语音信号进行滤波 (14) 3.7 回放语音信号 (17) 4 结论 (18) 致谢 (19) 参考文献 (20)