基于MATLAB的语音信号设计,基于matlab的语音信号设计.doc

基于matlab的语音信号设计.doc

数字信号处理课程设计题目基于MATLAB的语音信号设计姓名院系电子信息工程系专业通信工程班级092学号910705241指导教师留黎钦成绩指导教师日期李永棠基于MATLAB的语言信号设计22012年6月基于MATLAB的语音信号设计(电子信息工程学系指导教师留黎钦)1课程设计目的1了解语音信号的产生、采集，能绘制语音信号的频率响应曲线及频谱图；2学会用MATLAB对语音信号进行分析和处理；3掌握用滤波器去除语音信号噪声的方法，观察去噪前后的语音信号。2课程设计基本要求1学会MATLAB的使用，掌握MATLAB的程序设计方法。2掌握在WINDOWS环境下语音信号采集的方法。3掌握数字信号处理的基本概念、基本理论和基本方法。4掌握MATLAB设计FIR和IIR数字滤波器的方法。5学会用MATLAB对信号进行分析和处理。3课程设计原理数字滤波器根据其冲激响应函数的时域特性，可分为2种，即无限长冲激响应(IIR)滤波器和有限长冲激响应(FIR)滤波器。FIR和IIR的滤波原理都是进行卷积，就是对数入信号进行某种计算。A利用模拟滤波器设计IIR数字滤波器方法(1)根据所给出的数字滤波器性能指标计算出相应的模拟滤波器的设计指标。(2)根据得出的滤波器性能指标设计出相应的模拟滤波器的系统函数HS。(3)根据得出的模拟滤波器的系统函数HS，经某种变换得到对该模拟滤波器相应的。BFIR滤波器的窗函数设计法FIR滤波器的设计方法有许多种，如窗函数设计法、频率采样设计法和最优化设计法等。窗函数设计法的基本原理是用一定宽度窗函数截取无限脉冲响应序列获得有限长的脉冲响应序列，主要设计步骤为(1)通过傅里叶逆变换获得理想滤波器的单位脉冲响应HDN。(2)由性能指标确定窗函数WN和窗口长度N。(3)求得实际滤波器的单位脉冲响应HN，HN即为所设计FIR滤波器系数向量BN。C在MATLAB中本次设计涉及到的几个重要函数(1)采样Y,FS,NBITWAVREAD，返回采样值放在向量Y中，FS为采样频率(HZ)，NBIT为采样位数。(2)快速傅里叶变换YFFT(Y，N),其中X为有限长序列，N为序列Y的长度，Y为序列Y的DFT。(3)数字滤波器的频率响应H,WFREQZB,A,N,返回数字滤波器的N点复频率响应，B和A是滤波器系数的分子和分母向量；H是复频率响应，W是频率点，N默认为512。(模拟滤波器的频率响应用函李永棠基于MATLAB的语音信号设计4数FREQS)。4课程设计内容基本要求(1)录制一段自己的语音信号，并对录制的信号进行采样；(2)画出采样后语音信号的时域波形和频谱图；(3)给定滤波器的性能指标，采用窗函数法或双线性变换设计FIR滤波器或设计IIR滤波器，并画出滤波器的频率响应；(4)然后用自己设计的滤波器对采集的语音信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；(5)回放语音信号。扩展部分(1)设计GUI界面进行对该语音处理系统的界面化管理。(2)实现MATLAB界面中的直接录音功能5详细设计过程与调试分析51语音信号的采集录制原始语音我们利用WINDOWS的录音机，录制了一段自己的声音“同志们，我们一起吃饭吧”，单声道，时间在3S左右。接着保存为WAV格式(存在MATLAB软件的WORK文件夹下，命名为“CHIFAN”)。52原始语音信号的分析利用函数WAVREAD对语音信号进行采样，并读取语音的相关信息，并通过函数SOUND听到采样后自己所录的一段声音，并画出它的时域和频域的图形。程序如下FILENAME CHIFANWAV Y1,FS,NBITSWAVREADFILENAMEFPRINTF 语音“S“的信息\N ,FILENAMEFPRINTF 语言长度G秒\N ,LENGTHY1/FSFPRINTF 取样频率G取样点/秒\N ,FSFPRINTF 解析度G位元/取样点\N ,NBITSSOUNDY1,FSMLENGTHY1TIME01/FSM1/FS时间轴的向量YF1FFTY1FIGURESUBPLOT2,2,1PLOTTIME,Y1AXIS0,4,15,15TITLE 原信号的时域图 XLABEL 时间 YLABEL 幅度 SUBPLOT2,2,2PLOT0FS/MFSM1/M,ABSYF1AXIS0,25000,0,1200TITLE 原始语音信号频谱 XLABEL 频率 YLABEL 幅度 运行结果李永棠基于MATLAB的语音信号设计5图1原始信号分析53对语音信号进行加噪处理利用函数AWGN给语音信号加上高斯白噪声，调用SIN函数来单音正弦信号。通过函数SOUND来听到加完噪声后的声音，并画出它的时域和频域的图形。程序如下FILENAME CHIFANWAV Y1,FS,NBITSWAVREADFILENAMEFPRINTF 语音“S“的信息\N ,FILENAMEFPRINTF 语言长度G秒\N ,LENGTHY1/FSFPRINTF 取样频率G取样点/秒\N ,FSFPRINTF 解析度G位元/取样点\N ,NBITSMLENGTHY1TIME01/FSM1/FS时间轴的向量GAWGNY1,20给语音信号加上高斯白噪声SOUNDG,FS回放加噪信号GFFTG对加噪后的语音信号进行DFT变换N05SIN2PI7000TIME 噪声为7KHZ的正弦信号Y2Y1NWAVPLAYY2,FS, ASYNC 播放加噪声后的语音信号YF2FFTY2FIGURESUBPLOT2,2,1PLOTTIME,GAXIS0,4,15,15TITLE 加高斯噪声后语音信号时域图 XLABEL 时间 YLABEL 幅度 SUBPLOT2,2,2PLOT0FS/MFSM1/M,ABSGAXIS0,25000,0,1200TITLE 加高斯噪声后语音信号频谱图 XLABEL 频率 YLABEL 幅度 SUBPLOT2,2,3PLOTTIME,Y2AXIS0,4,15,15TITLE 加正弦信号的时域图 XLABEL 时间 YLABEL 幅度 SUBPLOT2,2,4PLOT0FS/MFSM1/M,ABSYF2AXIS0,25000,0,1200TITLE 加正弦信号的频谱图 XLABEL 频率 YLABEL 幅度 运行结果李永棠基于MATLAB的语音信号设计6图2原始信号加噪信号的时域和频谱图54设计数字滤波器和画出其频谱响应541设计个巴特沃斯低通滤波器设定通带截止频率(FP1000)、阻带截止频率(FS3000)、通带波纹系数(RP05)、阻带波纹系数(RS40)、抽样频率(FC22050)。程序如下巴特沃斯低通滤波器FP1000FS3000RP05RS40FC22050WP2FP/FCWS2FS/FC将模拟域转化成数字域。N,WCBUTTORDWP,WS,RP,RS估算巴特沃思滤波器的阶数N和3DB截止频率WC。B,ABUTTERN,WC求传输函数的分子和分母多项式的系数。FIGUREFREQZB,A,256,FC求滤波器的频率响应。TITLE 巴特沃斯低通滤波器的幅频和相频特性图 运行结果图3巴特沃斯滤波器的性能李永棠基于MATLAB的语音信号设计7542用KAI