基于matlab语音信号采集课程设计

数字信号处理课程设计报告课设题目语音信号的采集与处理学院信息科学与工程学院专业电子信息工程班级电子0903姓名学号指导教师张海云、李锦彬时间2012123120130111一引言介绍二MATLAB工具简介三课程设计任务31语音信号的采集32语音信号的频谱分析33设计数字滤波器和画出其频率响应给出各滤波器的性能指标34用滤波器对信号进行滤波35回放语音信号，分析滤波前后的语音变化四课程设计原理及方案41用窗函数法设计FIR滤波器42用双线性变换法设计IIR数字滤波器五语音信号的采集及频谱分析六设计滤波器及对语音信号进行滤波61FIR低通滤波器设计及信号滤波62FIR高通滤波器设计及信号滤波63FIR带通滤波器设计及信号滤波64IIR低通滤波器设计及信号滤波65IIR高通滤波器设计及信号滤波66IIR带通滤波器设计及信号滤波七课程设计心得体会八参考文献一引言随着信息时代和数字世界的到来，数字信号处理已成为今一门极其重要的学科和技术领域。数字信号处理在通信语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理应用中，数字滤波器十分重要并已获得广泛应用。二MATLAB工具简介1MATLAB是矩阵实验室（MATRIXLABORATORY）之意，现已发展成为适合多学科，多种工作平台的功能强大的大型软件，已经成为线性代数，自动控制理论，数理统计，数字信号处理，时间序列分析，动态系统仿真等高级课程的基本教学工具；2MATLAB的语言特点（1）语言简洁紧凑，使用方便灵活，库函数极其丰富。MATLAB程序书写形式自由，利用起丰富的库函数避开繁杂的子程序编程任务，压缩了一切不必要的编程工作。由于库函数都由本领域的专家编写，用户不必担心函数的可靠性。（2）运算符丰富。由于MATLAB是用C语言编写的，MATLAB提供了和C语言几乎一样多的运算符，灵活使用MATLAB的运算符将使程序变得极为简短。（3）MATLAB既具有结构化的控制语句（如FOR循环，WHILE循环，BREAK语句和IF语句），又有面向对象编程的特性。（4）程序限制不严格，程序设计自由度大。例如，在MATLAB里，用户无需对矩阵预定义就可使用。（5）程序的可移植性很好，基本上不做修改就可在各种型号的计算机和操作系统上运行。（6）MATLAB的图形功能强大。在FORTRAN和C语言里，绘图都很不容易，但在MATLAB里，数据的可视化非常简单。MATLAB还具有较强的编辑图形界面的能力。（7）MATLAB的缺点是，它和其他高级程序相比，程序的执行速度较慢。由于MATLAB的程序不用编译等预处理，也不生成可执行文件，程序为解释执行，所以速度较慢。（8）功能强大的工具箱是MATLAB的另一特色。MATLAB包含两个部分核心部分和各种可选的工具箱。核心部分中有数百个核心内部函数。其工具箱又分为两类功能性工具箱和学科性工具箱。功能性工具箱主要用来扩充其符号计算功能，图示建模仿真功能，文字处理功能以及与硬件实时交互功能。功能性工具箱用于多种学科。而学科性工具箱是专业性比较强的。（9）源程序的开放性。开放性也许是MATLAB最受人们欢迎的特点。除内部函数以外，所有MATLAB的核心文件和工具箱文件都是可读可改的源文件，用户可通过对源文件的修改以及加入自己的文件构成新的工具箱。三课程设计任务31语音信号的采集利用WINDOWS下的录音机，录制一段自己的话音，时间在3秒左右，然后在MATLAB软件平台下，利用函数WAVREAD对语音信号进行采样，记住采样频率和采样点数，采样率8000HZ，单通道，格式WAV。注如果没有办法录制自己语音，也可以用现成的语音也可以，文件格式WAV32语音信号的频谱分析在MATLAB中，可以利用函数FFT对信号进行快速傅立叶变换，得到信号的频谱特性，要求学生首先画出语音信号的时域波形，然后对语音信号进行频谱分析。33设计数字滤波器和画出其频率响应给出各滤波器的性能指标；给定滤波器的性能指标如下1低通滤波器的性能指标FB1000HZ,FC1200HZ,AS100DB,AP1DB2高通滤波器的性能指标FB2700HZ,FC3000HZ,AS100DB,AP1DB3带通滤波器的性能指标FB11200HZ,FB23000HZ,AS100DB,AP1DBFC11000HZ,FC23200HZ,AS100DB,AP1DB采用窗函数法和双线性变换法设计上面要求的3种滤波器，并画出滤波器的频率响应；34用滤波器对信号进行滤波然后用自己设计的滤波器对采集到的信号进行滤波，画出滤波后信号的时域波形及频谱，并对滤波前后的信号进行对比，分析信号的变化；35回放语音信号，分析滤波前后的语音变化；四滤波器设计原理及设计方案41用窗函数法设计FIR滤波器根据过渡带宽及阻带衰减要求，选择窗函数的类型并估计窗口长度N（或阶数MN1），窗函数类型可根据最小阻带衰减AS独立选择，因为窗口长度N对最小阻带衰减AS没有影响，在确定窗函数类型以后，可根据过渡带宽小于给定指标确定所拟用的窗函数的窗口长度N，设待求滤波器的过渡带宽为W，它与窗口长度N近似成反比，窗函数类型确定后，其计算公式也确定了，不过这些公式是近似的，得出的窗口长度还要在计算中逐步修正，原则是在保证阻带衰减满足要求的情况下，尽量选择较小的N，在N和窗函数类型确定后，即可调用MATLAB中的窗函数求出窗函数WD（N）。根据待求滤波器的理想频率响应求出理想单位脉冲响应HDN，如果给出待求滤波器频率应为HD，则理想的单位脉冲响应可以用下面的傅里叶反变换式求出在一般情况下，HDN是不能用封闭公式表示的，需要采用数值方法表示；从W0到W2采样N点，采用离散傅里叶反变换（IDFT）即可求出。用窗函数WDN将HDN截断，并进行加权处理，得到如果要求线性相位特性，则HN还必须满足12JNDDHNHEDHNN根据上式中的正、负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。要根据所设计的滤波特性正确选择其中一类。例如，要设计线性相位低通特性可选择HNHN1N一类，而不能选HNHN1N一类。验算技术指标是否满足要求，为了计算数字滤波器在频域中的特性，可调用FREQZ子程序，如果不满足要求，可根据具体情况，调整窗函数类型或长度，直到满足要求为止。42用双线性变换法设计IIR数字滤波器脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到/T/T之间，再用ZEST转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的/T/T一条横带里；第二步再通过标准变换关系ZES1T将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，映射关系如图1所示。图1双线性变换的映射关系为了将S平面的整个虚轴J压缩到S1平面J1轴上的/T到/T段上，可以通过以下的正切变换实现（1）式中,T仍是采样间隔。当1由/T经过0变化到/T时，由经过0变化到，也即映射了整个J轴。将式（1）写成1HNN将此关系解析延拓到整个S平面和S1平面，令JS，J1S1，则得再将S1平面通过以下标准变换关系映射到Z平面ZES1T从而得到S平面和Z平面的单值映射关系为23式（2）与式（3）是S平面与Z平面之间的单值映射关系，这种变换都是两个线性函数之比，因此称为双线性变换式（1）与式（2）的双线性变换符合映射变换应满足的两点要求。首先,把ZEJ，可得4即S平面的虚轴映射到Z平面的单位圆。其次，将SJ代入式（4），得因此由此看出，当0时，|Z|1。也就是说，S平面的左半平面映射到Z平面的单位圆内，S平面的右半平面映射到Z平面的单位圆外，S平面的虚轴映射到Z平面的单位圆上。因此，稳定的模拟滤波器经双线性变换后所得的数字滤波器也一定是稳定的。双线性变换法优缺点双线性变换法与脉冲响应不变法相比，其主要的优点是避免了频率响应的混叠现象。这是因为S平面与Z平面是单值的一一对应关系。S平面整个J轴单值地对应于Z平面单位圆一周，即频率轴是单值变换关系。这个关系如式（4）所示，重写如下上式表明，S平面上与Z平面的成非线性的正切关系，如图2所示。由图2看出，在零频率附近，模拟角频率与数字频率之间的变换关系接近于线性关系；但当进一步增加时，增长得越来越慢，最后当时，终止在折叠频率处，因而双线性变换就不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象，从而消除了频率混叠现象。图2双线性变换法的频率变换关系但是双线性变换的这个特点是靠频率的严重非线性关系而得到的，如式（4）及图2所示。由于这种频率之间的非线性变换关系，就产生了新的问题。首先，一个线性相位的模拟滤波器经双线性变换后得到非线性相位的数字滤波器，不再保持原有的线性相位了；其次，这种非线性关系要求模拟滤波器的幅频响应必须是分段常数型的，即某一频率段的幅频响应近似等于某一常数（这正是一般典型的低通、高通、带通、带阻型滤波器的响应特性），不然变换所产生的数字滤波器幅频响应相对于原模拟滤波器的幅频响应会有畸变，如图3所示。图3双线性变换法幅度和相位特性的非线性映射对于分段常数的滤波器，双线性变换后，仍得到幅频特性为分段常数的滤波器，但是各个分段边缘的临界频率点产生了畸变，这种频率的畸变，可以通过频率的预畸来加以校正。也就是将临界模拟频率事先加以畸变，然后经变换后正好映射到所需要的数字频率上。五语音信号的采集及频谱分析用WINDOWS下的录音机，用单声道录制一段音乐或声音，时间在5S内。然后MATLAB软件平台下，利用函数WAVREAD对语音信号进行采样，记住采样频率和采样点数。对语音信号进行快速傅立叶变换，在一个窗口同时画出信号的时域波形图和频谱图，分析语音信号的频谱特点Y,FS,NBITSWAVREADMUSICWAV加载语音信号利用函数WAVREAD对语音信号进行采样WAVPLAYY,FS回放语音信号NLENGTHY语音信号长度YFFTY,N用傅里叶变换求频域SUBPLOT2,1,1PLOTYTITLE图1原始信号时域波形图SUBPLOT2,1,2PLOTABSYTITLE图2原始信号频谱图六设计滤波器及对语音信号进行滤波61FIR低通滤波器设计及信号滤波X1,FS,BITSWAVREADMUSICWAVDERTA_FSFS/LENGTHX1设置频谱的间隔，分辨率，这里保证了X轴的点数必须和Y轴点数一致FSFSFP11000FS11200AS1100WP12PIFP1/FSWS12PIFS1/FSBF1WS1WP1WC1WP1WS1/2M1CEILAS1795/2286BF11按凯泽窗计算滤波器阶数N1M11BETA101102AS187WINDOWKAISERN1,BETA1求凯泽窗窗函数B1FIR1M1,WC1/PI,WINDOWWC1/PI为归一化,窗函数法设计函数FIGURE2FREQZB1,1,512H,WFREQZB,A,N,1中B和A分别为离散系统的系统函数分子、分母多项式的系数向量，返回量H则包含了离散系统频响在0PI范围内N个频率等分点的值（其中N为正整数），W则包含了范围内N个频率等分点。调用默认的N时，其值是512。TITLEFIR低通滤波器的频率响应X1_LOWFILTERB1,1,X1对信号进行低通滤波,YFILTERB,A,X，输入X为滤波前序列，Y为滤波结果序列，B/A提供滤波器系数，B为分子，A为分母SOUNDX1_LOW,FS,BITSFIGURE3SUBPLOT2,1,1PLOTX1_LOWTITLE信号经过FIR低通滤波器时域SUBPLOT2,1,2PLOTFS/2DERTA_FSFS/2DERTA_FS,ABSFFTSHIFTFFTX1_LOWTITLE信号经过FIR低通滤波器（频域）62FIR高通滤波器设计及信号滤波X1,FS,BITSWAVREADMUSICWAVDERTA_FSFS/LENGTHX1设置频谱的间隔，分辨率，这里保证了X轴的点数必须和Y轴点数一致FSFSAS2100FP23000FS22700WP22PIFP2/FSWS22PIFS2/FSBF2WP2WS2WC2WP2WS2/2M2CEILAS2795/2286BF21按凯泽窗计算滤波器阶数N2M21BETA201102AS287WINDOWKAISERN2,BETA2求凯泽窗窗函数B2FIR1M2,WC2/PI,HIGH,WINDOWFIGURE4FREQZB2,1,512数字滤波器频率响应TITLEFIR高通滤波器的频率响应X1_HIGHFILTERB2,1,X1对信号进行高通滤波SOUNDX1_HIGH,FS,BITSFIGURE5SUBPLOT211PLOTX1_HIGHTITLE信号经过FIR高通滤波器时域SUBPLOT212PLOTFS/2DERTA_FSFS/2DERTA_FS,ABSFFTSHIFTFFTX1_HIGHTITLE信号经过FIR高通滤波器（频域）63FIR带通滤波器设计及信号滤波FSFSAS3100FP31200,3000FS31000,3200WP32PIFP3/FSWS32PIFS3/FSBF3WP31WS31WC3WP3BF3/2M3CEILAS3795/2286BF31按凯泽窗计算滤波器阶数N3M31BETA301102AS387WINDOWKAISERN3,BETA3求凯泽窗窗函数B3FIR1M3,WC3/PI,BANDPASS,WINDOW带通滤波器FIGURE6FREQZB3,1,512数字滤波器频率响应TITLEFIR带通滤波器的频率响应X1_DAITONGFILTERB3,1,X1对信号进行带通滤波SOUNDX1_DAITONG,FS,BITSFIGURE7SUBPLOT211PLOTX1_DAITONGTITLE信号经过FIR带通滤波器时域SUBPLOT212PLOTFS/2DERTA_FSFS/2DERTA_FS,ABSFFTSHIFTFFTX1_DAITONGTITLE信号经过FIR带通滤波器（频域）64IIR低通滤波器设计及信号滤波X1,FS,BITSWAVREADMUSICWAVDERTA_FSFS/LENGTHX1设置频谱的间隔，分辨率，这里保证了X轴的点数必须和Y轴点数一致FSFSFP1I1000FS1I1200WP1I2PIFP1I/FSWS1I2PIFS1I/FSRP1I1RS1I100TS1/FSWP1I2/TSTANWP1I/2WS1I2/TSTANWS1I/2按频率转换公式进行转换，预畸变N1I,WN1ICHEB1ORDWP1I,WS1I,RP1I,RS1I,S计算模拟滤波器的最小阶数B1I,A1ICHEBY1N1I,RP1I,WN1I,S设计模拟原型滤波器BZ1I,AZ1IBILINEARB1I,A1I,FS运用双线性变换法得到数字滤波器传递函数FIGURE10FREQZBZ1I,AZ1I,512,FSTITLE切比雪夫1型低通滤波器的频率响应X1_LOWIFILTERBZ1I,AZ1I,X1对信号进行低通滤波,YFILTERB,A,X，输入X为滤波前序列，Y为滤波结果序列，B/A提供滤波器系数，B为分子，A为分母SOUNDX1_LOWI,FS,BITSFIGURE11SUBPLOT211PLOTX1_LOWITITLE信号经过IIR低通滤波器时域SUBPLOT212PLOTFS/2DERTA_FSFS/2DERTA_FS,ABSFFTSHIFTFFTX1_LOWITITLE信号经过IIR低通滤波器（频域）65IIR高通滤波器设计及信号滤波X1,FS,BITSWAVREADMUSICWAVDERTA_FSFS/LENGTHX1设置频谱的间隔，分辨率，这里保证了X轴的点数必须和Y轴点数一致IIR高通按钮消息响应函数FSFSFP2I3000通带截止频率FS2I2700组带截止频率RP2I1RS2I100WP2I22PIFP2I/FSWS2I22PIFS2I/FSWP2ITANWP2I2/2WS2ITANWS2I2/2WPHI1归一化WSHIWP2I/WS2I高通频率转换成低通频率指标估计滤波器的阶数N2I,WN2ICHEB1ORDWPHI,WSHI,RP2I,RS2I,S设计滤波器B2I,A2ICHEBY1N2I,RP2I,WN2I,SNUMI2,DENI2LP2HPB2I,A2I,WP2IBZ2I,AZ2IBILINEARNUMI2,DENI2,05FIGURE12FREQZBZ2I,AZ2I,512,FSTITLE切比雪夫1型IIR高通滤波器的频率响应X1_HIGHIFILTERBZ2I,AZ2I,X1对信号进行高通滤波,YFILTERB,A,X，输入X为滤波前序列，Y为滤波结果序列，B/A提供滤波器系数，B为分子，A为分母SOUNDX1_HIGHI,FS,BITSFIGURE13SUBPLOT211PLOTX1_HIGHITITLE信号经过IIR高通滤波器时域SUBPLOT212PLOTFS/2DERTA_FSFS/2DERTA_FS,ABSFFTSHIFTFFTX1_HIGHITITLE信号经过IIR高通滤波器（频域）66IIR带通滤波器设计及信号滤波X1,FS,BITSWAVREADMUSICWAVDERTA_FSFS/LENGTHX1设置频谱的间隔，分辨率，这里保证了X轴的点数必须和Y轴点数一致FSFSFDP1I1200FDP2I3000通带截止频率FDS1I1000FDS2I3200组带截止频率RP3I1RS3I100WDP1I2PIFDP1I/FSWDP2I2PIFDP2I/FSWDS2I2PIFDS1I/FSWDS2I2PIFDS2I/FSWDP1ITANWDP1I/2WDP2ITANWDP2I/2WDS1ITANWDS2I/2WDS2ITANWDS2I/2按频率转换公式进行转换，预畸变BWDP2IWDP1IW3IWDP1IWDP2I/WDS1I调整截止频率参数。W0WDP1IWDP2I05WP3I1归一化，求低通原型的参数。WS3IWDP1IWDP2IWDS2I2/BWDS2IN3I,WN3IBUTTORDWP3I,WS3I,RP3I,RS3I,S计算模拟滤波器的最小阶数B3I,A3IBUTTERN3I,WN3I,S设计模拟原型滤波器NUMI3,DENI3LP2BPB3I,A3I,W0,BBZ3I,AZ3IBILINEARNUMI3,DENI3,05双线性变换H,WFREQZBZ3I,AZ3I,512FIGURE14SUBPLOT211PLOTW/PI,20LOG10ABSHGRIDYLABELGAINTITLE巴特沃兹IIR带通滤波器AXIS01805SUBPLOT212PLOTW/PI,ANGLEHGRIDXLABELW/PI,YLABELPHASEX1_DAITONGIFILTERBZ3I,AZ3I,X1对信号进行带通滤波,YFILTERB,A,X，输入X为滤波前序列，Y为滤波结果序列，B/A提供滤波器系数，B为分子，A为分母SOUNDX1_DAITONGI,FS,BITSFIGURE15SUBPLOT211PLOTX1_DAITONGITITLE信号经过IIR带通滤波器时域SUBPLOT212PLOTFS/2DERTA_FSFS/2DERTA_FS,ABSFFTSHIFTFFTX1_DAITONGITITLE信号经过IIR带通滤波器（频域）七课程设计心得体会通过这两周的课程设计，自己能够更好的使用MATLAB进行程序设计，对滤波器有了进一步的了解，初步学会设计简单的滤波器。刚开始的时候，先分析要做什么，怎么来做，如何来更好的完成此次的课程