下载文档原格式
语音信号的采集和频谱分析：[y,fs,bits]=wavread('voice'); %读取音频信息（双声道，16位，频率44100Hz）sound(y,fs,bits); %回放该音频Y=fft(y,4096); %进行傅立叶变换subplot(211);plot(y);title('声音信号的波形');subplot(212)plot(abs(Y));title('声音信号的频谱');窗函数设计低通滤波器：fp=1000;fc=1200;as=100;ap=1;fs=22000;wp=2*fp/fs;wc=2*fc/fs;N=ceil((as-7.95)/(14.36*(wc-wp)/2))+1;beta=0.1102*(as-8.7);window=Kaiser(N+1,beta);b=fir1(N,wc,window);freqz(b,1,512,fs);结果：滤波：[y,fs,bits]=wavread('voice');d=filter(b,a,y);D=fft(d);subplot(211)plot(d);title('滤波后的声音波形')subplot(212)plot(abs(D))title('滤波后的声音频谱')回放：sound(d,fs,bits)与滤波之前相比，噪音明显降低了许多。
过零率的计算要用下面的代码：zcr = zeros(size(y,1)1);delta= 0.02;for i=1:size(y,1)x=y(i,:);for j=1;length(x)-1if x(j)*x(j+1)<0 &abs(x(j)-x(j+1))>deltazcr(i)=zcr(i)+1;endendend其中设置了门限delta=0.02。
这是个经验值，可以进行细微的调整。
【毕业设计】语音信号的采集与分析河南农业大学本科生毕业论文（设计）任务书论文（设计）题目语音信号的采集与分析学院专业班级学号姓名2009年月日语音信号的采集与分析作者：123 指导老师：456摘要语音信号的采集与分析技术是一门涉及面很广的交叉科学，它的应用和发展与语音学、声音测量学、电子测量技术以及数字信号处理等学科紧密联系。
其中语音采集和分析仪器的小型化、智能化、数字化以及多功能化的发展越来越快，分析速度较以往也有了大幅度的高。
本文简要介绍了语音信号采集与分析的发展史以及语音信号的特征、采集与分析方法，并通过PC机录制自己的一段声音，运用Matlab进行仿真分析，最后加入噪声进行滤波处理，比较滤波前后的变化。
关键词：语音信号，采集与分析，MatlabAudio signal acquisition and analysisAuthor：zhuyousong Teacher guidance：lifuqiangAbstractSpeech signal acquisition and analysis techniques are a wide range of cross-scientific，Its application and development of voice study, sound measurement study, electronic measuring technology, and digital signal processing disciplines, such as close contact。
Collection and analysis of voice one of the small-scale equipment, intelligence, digital and multi-functional development of more and more quickly, faster than the previous analysis has been substantially high。
南昌工程学院《语音信号的采集与分析》课程设计题目语音信号的采集与分析课程名称语音信号处理系院信息工程学院专业通信工程班级 10通信工程2班学生姓名刘敏学号 2010103362设计地点电子信息楼指导教师邹宝娟设计起止时间：2013年12月9日至2013年12月20日目录一、需求分析 (4)1.1选题背景及意义 (4)1.2设计要求 (4)二、系统总体设计 (4)2.1 系统设计思路 (4)2.2 功能结构图及功能说明 (4)2.3 工作原理 (6)三、系统详细设计 (6)3.1 语音信号的matlab仿真的数据分析 (6)3.2 程序代码分析 (12)四、调试与维护 (14)4.1 调试过程的问题与维护 (14)五、结束语 (15)六、参考文献 (16)七、指导教师评阅（手写） (17)一、需求分析1.1选题背景及意义该设计主要是介绍语音信号的采集与分析方法，通过PC机录制自己的一段声音，运用Matlab提供的函数进行仿真分析，并画出采样后语音信号的时域波形和频谱图，对所采集的语音信号加入干扰随机高斯噪声，对加入噪声的信号进行播放，并进行时域和频谱分析；对比加噪前后的时域图和频谱图，分析讨论采用什么样的滤波器进行滤除噪声。
1.2设计要求(1)通过PC机录制自己的一段声音“南昌工程学院刘敏”；(2)运用MATLAB中信号处理相关的函数对语音信号进行时域、频域上的分析，如短时能量，短时平均过零率，语谱图等；(3)运用MATLAB对语音信号进行综合与分析，包括语音信号的调制，叠加，和滤波等。
二、系统总体设计2.1 系统设计思路系统的整体设计思路包括语音信号的录制，语音信号的采集，语音信号的分析，其中语音信号的分析又包括了语音信号的时域分析和频域分析，语音信号的加噪处理和滤噪设计分析。
2.2 功能结构图及功能说明实际工作中，我们可以利用windows自带的录音机录制语音文件，声卡可以完成语音波形的A/D转换，获得WAVE文件，为后续的处理储备原材料。
实验二：语音信号的频域分析实验目的：以MATLAB 为工具，研究语音信号的频域特性，以及这些特性在《语音信号处理》中的应用情况。
实验要求：利用所给语音数据，分析语音的频谱、语谱图、基音频率、共振峰等频域参数。
要求会求取这些参数，并举例说明这些参数在语音信号处理中的应用。
实验内容：1、 语音信号的频谱分析1.1加载“ma1_1”语音数据。
基于DFT 变换，画出其中一帧数据（采样频率为8kHz ，帧长为37.5ms ，每帧有300个样点）的频域波形（对数幅度谱）。
load ma1_1;x = ma1_1 (4161:4460); plot (x)N = 1024; k = - N/2:N/2-1;X = fftshift (fft (x.*hann (length (x)),N));plot (k,20*log10 (abs(X))), axis ([0 fix(N/2) -inf inf ])已知该帧信号的时域波形如图（a ）所示，相应的10阶LPC 谱如图（b ）所示。
问题1：这帧语音是清音还是浊音？基于DFT 求出的对数幅度谱和相应的LPC 谱相比，两者有什么联系和区别？问题2：根据这帧基于DFT 的对数幅度谱，如何估计出共振峰频率和基音周期？问题3：时域对语音信号进行加窗，反映在频域，其窗谱对基于DFT 的对数幅度谱有何影响？如何估计出窗谱的主瓣宽度？1.2对于浊音语音，可以利用其频谱)(ωX 具有丰富的谐波分量的特点，求出其谐波乘积谱：∏==R r r X HPSx 1)()(ωω式中，R 一般取为5。
在谐波乘积谱中，基频分量变得很大，更易于估计基音周期。
1.3加载“vowels.mat”语音数据，分别画出一帧/i/和一帧/u/（采样频率为10kHz，帧长为30ms，每帧有300个样点）的基于DFT的对数幅度谱。
其Matlab代码如下：load vowelsx = vowels.i_1(2001:2300);N = 1024; k= -N/2:N/2-1;X = fftshift (fft (x.*hann (length(x)),N));plot (k,20*log10(abs(X))), axis([0 fix(N/2) 0 100])x = vowels.u_1(2001:2300);N= 1024; k = -N/2:N/2-1;X = fftshift (fft (x.*hann(length(x)),N));plot (k,20*log10(abs(X))), axis([0 fix(N/2) 0 100])1.4画出一帧清音语音的基于DFT的对数幅度谱。
语音信号采样和频谱分析一．实验目的（1）掌握傅里叶变换的物理意义，深刻理解傅里叶变换的内涵；（2）了解MATLAB 对声音信号的处理指令；（3）了解计算机存储信号的方式及语音信号的特点；（4）加深对采样定理的理解；（5）加深学生对信号分析工程应用的理解，拓展学生在信号分析领域的综合应用能力。
二．实验内容本实验利用MATLAB 指令录制一段语音信号，观察其时域波形并进行傅里叶变换，观察其频域的频谱。
根据该信号的频谱构成，选择三种不同的采样频率重新录制该语音信号，并试听回放效果，进行比较，以验证采样定理，并了解MATLAB 对声音信号的处理指令，加深对采样定理的理解。
关键词：傅里叶变换 信号采样三、实验原理语音信号是一种连续变化的模拟信号，而计算机只能处理和记录二进制的数字信号，因此，由自然音而得的音频信号必须用计算机的声音编辑工具，先进行语音采样，然后利用了计算机上的A/D 转换器，将模拟的声音信号变成离散的量化了的数字信号量化和编码，变成二进制数据后才能送到计算机进行再编辑和存储。
语音信号输出时，量化了的数字信号又通过D/A 转换器，把保存起来的数字数据恢复成原来的模拟的语音信号。
（1）应用MATLAB 进行声音的录制 （2）应用MATLAB 进行声音的播放 （3）语音信号的频谱分析 。
傅里叶变换建立了信号频谱的概念。
所谓傅里叶分析即分析信号的频谱（频率构成）、频带宽度等。
对语音信号的分析也不例外，也必须采用傅里叶变换这一工具。
对于连续时间信号)(t f ，其傅里叶变换)(ωF 为：⎰∞∞--=dt e t f F t j ωω)()(四、实验任务（1）应用MATLAB 进行声音的录制在MATLAB 命令窗口中键入“y=wavrecord(8000,8000,1)”，并按回车键，此时刻以后的1（8000/8000）秒时段内的声音信号将以y 为文件名，以数字声音信号.wav 格式存储在MATLAB 的工作空间里。
1 绪论1.1 课题的背景与意义通过语音传递倍息是人类最重要、最有效、最常用和最方便的交换信息的形式。
语言是人类持有的功能．声音是人类常用的工具，是相互传递信息的最主要的手段。
因此，语音信号是人们构成思想疏通和感情交流的最主要的途径。
并且，由于语言和语音与人的智力活动密切相关，与社会文化和进步紧密相连，所以它具有最大的信息容量和最高的智能水平。
现在，人类已开始进入了信息化时代，用现代手段研究语音信号，使人们能更加有效地产生、传输、存储、获取和应用语音信息，这对于促进社会的发展具有十分重要的意义。
让计算机能听懂人类的语言，是人类自计算机诞生以来梦寐以求的想法。
随着计算机越来越向便携化方向发展，随着计算环境的日趋复杂化，人们越来越迫切要求摆脱键盘的束缚而代之以语音输人这样便于使用的、自然的、人性化的输人方式。
作为高科鼓应用领域的研究热点，语音信号采集与分析从理论的研究到产品的开发已经走过了几十个春秋并且取得了长足的进步。
它正在直接与办公、交通、金融、公安、商业、旅游等行业的语音咨询与管理．工业生产部门的语声控制，电话、电信系统的自动拨号、辅助控制与查询以及医疗卫生和福利事业的生活支援系统等各种实际应用领域相接轨，并且有望成为下一代操作系统和应用程序的用户界面。
可见，语音信号采集与分析的研究将是一项极具市场价值和挑战性的工作。
我们今天进行这一领域的研究与开拓就是要让语音信号处理技术走人人们的日常生活当中，并不断朝更高目标而努力。
语音信号采集与分析之所以能够那样长期地、深深地吸引广大科学工作者去不断地对其进行研究和探讨，除了它的实用性之外，另一个重要原因是，它始终与当时信息科学中最活跃的前沿学科保持密切的联系．并且一起发展。
语音信号采集与分析是以语音语言学和数字信号处理为基础而形成的一门涉及面很广的综合性学科，与心理、生理学、计算机科学、通信与信息科学以及模式识别和人工智能等学科都有着非常密切的关系。
对语音信号采集与分析的研究一直是数字信号处理技术发展的重要推动力量。
语音信号的采集及预处理1.语音信号的录音、读入、放音等[x,fs,nbit]=wavread('D:\2.wav'); %fs=10000,nbit=16y=soundview('D:\2.wav')2.语音信号的分帧程序：[x,fs,nbit]=wavread('D:\2.wav');len=256;inc=128;y=enframe(x,len,inc);figure;subplot(2,1,1),plot(x)subplot(2,1,2),plot(y)3.语音信号加窗：程序：N=256;w = window('rectangle',N);w1 = window('hamming',N);w2 = window('hanning',N);wvtool(w,w1,w2)4.预加重程序：[x,fs,nbit]=wavread('D:\2.wav');len=256;inc=128;y=enframe(x,len,inc);z=filter([1-0.9375],1,y)figure(2)subplot(2,1,1),plot(y)subplot(2,1,2),plot(z)语音信号的时域分析1.语音信号的录音、读入、放音等：利用函数wavread对语音信号进行采样,记住采样频率和采样点数，给出以下语音的波形图（2.wav）。
[Y,FS,NBITS]= wavread('D:\2.wav')X= wavread('D:\2.wav')plot(X)2.短时能量分析：(1)首先对语音信号预加重；(2)对预加重后的语音信号进行分帧，帧长取N=256各样值点，帧移取128个样值点；(3)求短时能量。
org=wavread('D:\2.wav')wgt=filter([1 -0.9375],1,org)w1=enframe(wgt,256,128)amp=sum(abs(w1),2)plot(amp)3.短时过零率分析：求语音信号的短时过零率。
实验一 语音信号的采集及预处理一、实验目的在理论学习的基础上，进一步地理解和掌握语音信号预处理及短时加窗的意义及基于matlab 的实现方法。
二、实验原理1. 语音信号的录音、读入、放音等：练习matlab 中几个音频处理函数，利用函数wavread 对语音信号进行采样,记住采样频率和采样点数，给出以下语音的波形图（2.wav ）。
利用wavplay 或soundview 放音。
也可以利用wavrecord 自己录制一段语音，并进行以上操作(需要话筒)。
2. 语音信号的分帧：对语音信号进行分帧，可以利用voicebox 工具箱中的函数enframe 。
voicebox 工具箱是基于GNU 协议的自由软件，其中包含了很多语音信号相关的函数。
3． 语音信号的加窗：本步要求利用window 函数设计窗口长度为256(N=256)的矩形窗(rectwin)、汉明窗(hamming)及汉宁窗(hann))，利用wvtool 函数观察其时域波形图及频谱特性，比较得出结论。
观察整个信号加矩形窗及汉明窗后的波形，利用subplot 与reshape 函数将分帧后波形、加矩形窗波形及加汉明窗波形画在一张图上比较。
取出其中一帧，利用subplot 与reshape 函数将一帧语音的波形、加矩形窗波形及加汉明窗波形画在一张图上比较将得出结论。
4. 预加重：即语音信号通过一个一阶高通滤波器19375.01--z 。
三、实验步骤、实验程序、图形及结论1.语音信号的录音、读入、放音等程序：[x,fs,nbit]=wavread('D:\2.wav'); %fs=10000,nbit=16y=soundview('D:\2.wav')2.语音信号的分帧程序：[x,fs,nbit]=wavread('D:\2.wav');len=256;inc=128;y=enframe(x,len,inc);figure;subplot(2,1,1),plot(x)subplot(2,1,2),plot(y)3.语音信号加窗：程序：N=256;w = window('rectangle',N);w1 = window('hamming',N);w2 = window('hanning',N); wvtool(w,w1,w2)4.预加重程序：[x,fs,nbit]=wavread('D:\2.wav'); len=256;inc=128;y=enframe(x,len,inc);z=filter([1-0.9375],1,y)figure(2)subplot(2,1,1),plot(y)subplot(2,1,2),plot(z)四、思考题1.语音信号包括哪些预处理，作用分别是什么？2.不同窗口的优缺点，窗口长度如何选取？。
利用MATLAB软件对音频信号进行频谱分析与处理一、简介频谱分析是通过对信号的频率成分进行分析，它允许我们了解信号的特性，计算信号的能量分布，同时还可以用来定位造成干扰的频率组件，以及检测和分析信号的变化。
MATLAB是一种编程语言和科学计算软件，它可以非常便捷地实现对音频信号的频谱分析和处理。
二、实现方法1.导入音频信号在使用MATLAB进行频谱分析时，首先需要先将音频信号导入MATLAB环境中。
可以使用audioplayer和audioread函数来完成这一步骤，示例代码如下：[audioData, fs] = audioread(‘AudioFile.wav’);player = audioplayer(audioData, fs);play(player);其中audioData表示从wav文件中读取的音频数据，fs表示采样率，player表示存储audioData和fs的audioplayer实例，play函数可以播放音频文件。
2.信号预处理针对所记录的音频信号，需要进行一些基本的信号处理操作，包括去噪、均衡、时域平均等。
去噪可以用MATLAB内置的函数完成，例如：audioData_NoiseRemoved = denoise(audioData,‘meanspectrum’);均衡是指将频谱的一些区域调整到更好的水平，可以用equalizer函数实现：audioData_Equalized = equalizer(audioData, ‘bandwidth’, 0.2);时域平均则可以使用conv函数实现：audioData_Meaned = conv(audioData, [1/N 1/N ... 1/N]);3.频谱分析频谱分析的主要工作是计算信号的谱密度，也就是每一个频率分量的能量。
《信号与系统》大作业语音信号的采集与频谱分析——基于Matlab的语音信号处理学生姓名：学号：专业班级：电子工程学院卓越班指导老师：2015年6月22日摘要本设计用苹果手机自带的录音设备采集了原始语音，并导入了电脑转成wav格式，然后用MATLAB和Adobe audition对其进行时域分析。
接着利用傅里叶变换进行了频域分析，绘制频谱图，再录制一段加上歌曲的伴奏的语音与原唱进行了对比分析，得出了我与歌星在频域上的差别。
本设计给信号加了两种噪声并通过观察加噪后的频谱和试听回放效果比较加噪前后的差别，最后，设计了FIR数字低通滤波器和带通滤波器，分析滤波前后的频谱。
再次试听回放效果，得出结论。
关键词：语音、FFT、频谱图、噪声、滤波器AbstractThis design is based on the general function of Matlab and Adobe edition to deal with Audio signals. The original signals are collected by iPhone’s built-in recording equipment.First，I compare the file generated by myself with that of thesame song sang by a famous singer.The emphasis is generally laid on analysing the difference in frequncy domain,but time domain will be included too.After that,two noise signals are added to the original signal respectively and let them pass a filter to analyse it.In the two process mentioned before,I make comparison between the before and after frequency domain.Sampling Theorem is the base of my design.It is by sampling we can get discrete signals from the original one and draw the image in time domain.Also,fast fourier transform is employed(FFT)to get the signals in frequency domain.The ayalysis of frequency domain is the highlight of this design.Through this design,I can deepen my comprehension of principles of audio signals and I have learnt how to deal with it.Through met with much hindrance,I improved my skills finally.Keywords： audio signal、TTT、noise、filter1 绪论1.1课题的研究意义语音信号处理属于信息科学的一个重要分支，它是研究用数字信号处理技术对语音信号进行处理的一门新兴学科，同时又是综合性的多学科领域和涉及面很广的交叉学科，因此我们进行语言信号处理具有时代的意义。
学会运用MATLAB的信号处理功能，采集语音信号，并对语音信号进行滤波及变换处理，观察其时域和频域特性，加深对信号处理理论的理解，这为今后熟练使用MATLAB进行系统的分析仿真和设计奠定基础。
1.2 课题的目的与要求利用Matlab和Adobe audition对语音信号进行数字信号处理和分析，要求采集语音信号后，在软件平台进行频谱分析；并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器滤除噪声，再进行频谱分析和再播放来感受噪声滤除的效果。
2.理论依据采样频率、采样位数的概念，采样定理；时域信号的FFT分析；数字滤波器设计原理和方法，各种不同类型滤波器的性能比较。
2.1. 时域信号的FFT分析FFT即为快速傅氏变换，是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。
它利用“一分为二”的思想不断进行，点数越多，运算量的节约就越大，这就是FFT的优越性。
在MATLAB 的信号处理工具箱中函数FFT 和IFFT 用于快速傅立叶变换和逆变换。
函数FFT 用于序列快速傅立叶变换，其调用格式为y=fft(x)，其中，x 是序列，y 是序列的FFT ，x 可以为一向量或矩阵，若x 为一向量，y 是x 的FFT 且和x 相同长度；若x 为一矩阵，则y 是对矩阵的每一列向量进行FFT 。
如果x 长度是2的幂次方，函数fft 执行高速基－2FFT 算法，否则fft 执行一种混合基的离散傅立叶变换算法，计算速度较慢。
函数FFT 的另一种调用格式为y=fft(x,N)，式中，x ，y 意义同前，N 为正整数。
函数执行N 点的FFT ，若x 为向量且长度小于N ，则函数将x 补零至长度N ；若向量x 的长度大于N ，则函数截短x 使之长度为N ；若x 为矩阵，按相同方法对x 进行处理。
2.1.2采样定理采样定理，又称香农采样定律、奈奎斯特采样定律，是信息论，特别是通讯与信号处理学科中的一个重要基本结论。
采样是将一个信号（即时间或空间上的连续函数）转换成一个数值序列（即时间或空间上的离散函数），这个数值序列即是一个数字信号。
采样定理是将连续的模拟信号转变为离散的数字信号的方法。
对于一个时域连续的信号f(t) 利用开关函数s(t)进行采样，或者冲激函数序列进行采样（()T t δ）。
两种取样都具有一定的取样间隔，同时要求f(t)为频带带限信号，否则取样后将发生频谱混合重叠，破坏原信号。
采样定理公式如下：其中Ts 为采样间隔，设f(t)的频谱区间在(,)m m ωω-之间，即信号频带带限，需满足的采样条件为：2s m f f ≥，其中fs 为采样频率， 这样便完成了对于连续信号的采样形成离散的数字信号方便处理。
2.1.2 数字滤波器()()()()2s s s n t f t f nT t nT Sa ωδ∞=-∞=-*∑2m ωπ数字滤波器是一个输入输出都是离散时间信号的特定功能装置。
应用数字滤波器处理模拟信号（对应模拟频率）时，首先须对输入模拟信号进行限带、抽样和模数转换。
数字滤波器输入信号的数字频率（2π*f/fs,f为模拟信号的频率，fs为采样频率，注意区别于模拟频率),按照奈奎斯特取样定理，要使抽样信号的频谱不产生重叠，应小于折叠频率（ws/2=π)，其频率响应具有以2π为间隔的周期重复特性，且以折叠频率即ω=π点对称。
为得到模拟信号，数字滤波器处理的输出数字信号须经数模转换。
数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。
数字滤波器有低通、高通、带通、带阻和全通等类型。
它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。
应用最广的是线性、时不变数字滤波器，以及FIR滤波器。
从单位脉冲响应的角度，可以把数字滤波器分为：IIR滤波器（无限长单位冲激响应滤波器）和FIR滤波器（有限长单位冲激响应滤波器）3.设计思想3.1语音信号的采集此次我选择的歌曲是赵薇的《情深深雨蒙蒙》，用苹果手机自带的录音功能（Voice Memos）录音，和windows自带的录音机程序相比，优点是方便，而且可以忽略电脑运行本身的噪声。
然后将其导入到电脑，转换格式存为wave格式，并对从网上下载原唱的语音用同样的方式转换格式，并保存到正确的目录，以便以后能在matlab中打开。
同时手机自带的录音功能还能完成简单的音频剪辑操作，减掉多余的不需要的部分。
3.2 语音处理工具的选择① MatlabMatlab是大家都熟悉的一种科学计算软件，在信号与系统中运用很广泛，优点是可编程性强，且具有很好的可移植性和可扩展性，而且和专业信号与系统关联更大。
但缺点是数据都以矩阵形式进行处理，对我来说不够直观。
② Adobe auditionAdobe audition是一款专门处理音频的软件，可提供先进的编辑、混音、控制和效果处理功能，优点是很方便，比如采样频率可以直接选择，处理效果很直观，缺点是和信号与系统相关性较小，可编程性差。
这次设计的频谱分析主要是采用了matlab进行分析，但对原始语音频谱的分析以及由采样频率不同带来的变化中我同时采用了matlab和adobeaudition来分析。
3.3噪声函数的选择①rand函数产生随机噪声，具体语句将在下面的列表中列出②高斯白噪声Gaussian White Noise此噪声在《随机信号分析课程》中已学过，此次向原始信号中加入此噪声，加强理解。
3.4 滤波器的选择设计中选择的是FIR有限长单位冲激响应滤波器，相比于IIR滤波器，它有如下的优点：⑴ FIR滤波器，由于冲激响应有限长，可以用快速傅里叶变换，这样运算速度可以快很多；⑵ FIR滤波器可以做到严格的线性相位，而若要IIR得到严格的线性相位，则要大大增加滤波器的阶数，⑶ FIR滤波器主要采用非递归结构，因而从理论上以及时性从实际的有限精度的运算中，都是稳定的。
有限精度运算误差也较小，IIR滤波器必须采用递归的结构，极点必须在Z平面单位圆内，才能稳定，这种结构，运算中的四舍五入处理，有时会引起寄生振荡。
但是IIR也有相同技术指标下阶数较少，存储单元少，运算次数少的有点，此处不再加以赘述。
两种滤波器：①低通巴特沃斯滤波器巴特沃斯滤波器的特点为的通频带频率响应最为平滑，阻带缓慢衰减。
用窗函数法设计低通滤波器，具体指标为fp=3000Hz fc=4000Hz Ap=1dB As=100dB②高通滤波器具体指标为：fp=35ooHz fc=4000Hz Ap=1dB As=10dB3.5 Matlab函数的选择（1）[y,fs,bits]=wavread(‘filename’)wavread函数可读入语音信号y是采样获得的离散时域信号，fs是返回的采样频率， bits是采样位数，filename是音频信号的文件名。
（2）wavwrite(y,fs,bits,’filename’)wavwrite函数，改变采样频率，可以将22.05KHz或者11.025KHz的信号写回并观察这样的采样频率会产生什么样的结果。
