基于matlab的语音信号设计

《基于matlab的语音信号设计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、李永棠基于MATLAB的语言信号设计成绩指导教师：日期：《数字信号处理》课程设计题目：基于MATLAB的语音信号设计姓名：院系：电子信息工程系专业：通信工程班级：092学号：910705241指导教师：留黎钦2012年6月2李永棠基于MATLAB的语音信号设计基于MATLAB的语音信号设计（电子信息工程学系指导教师：留黎钦）1.课程设计目的1.了解语音信号的产生、采集，能绘制语音信号的频率响应曲线及频谱图；2.学会用MATLAB对语音信号进行分析和处理；3.掌握用滤波器去除语音信号噪声的方法，观察去噪

2、前后的语音信号。2.课程设计基本要求1.学会MATLAB的使用，掌握MATLAB的程序设计方法。2.掌握在Windows环境下语音信号采集的方法。3.掌握数字信号处理的基本概念、基本理论和基本方法。4.掌握MATLAB设计FIR和IIR数字滤波器的方法。5.学会用MATLAB对信号进行分析和处理。3.课程设计原理数字滤波器根据其冲激响应函数的时域特性，可分为2种，即无限长冲激响应（IIR）滤波器和有限长冲激响应（FIR）滤波器。FIR和IIR的滤波原理都是进行卷积，就是对数入信号进行某种计算。a.利用

3、模拟滤波器设计IIR数字滤波器方法（1）根据所给出的数字滤波器性能指标计算出相应的模拟滤波器的设计指标。（2）根据得出的滤波器性能指标设计出相应的模拟滤波器的系统函数H(S)。（3）根据得出的模拟滤波器的系统函数H(S)，经某种变换得到对该模拟滤波器相应的。b.FIR滤波器的窗函数设计法FIR滤波器的设计方法有许多种，如窗函数设计法、频率采样设计法和最优化设计法等。窗函数设计法的基本原理是用一定宽度窗函数截取无限脉冲响应序列获得有限长的脉冲响应序列，主要设计步骤为：（1）通过傅里叶逆变换获得理想滤波器

4、的单位脉冲响应hd(n)。（2）由性能指标确定窗函数W(n)和窗口长度N。（3）求得实际滤波器的单位脉冲响应h(n)，h(n)即为所设计FIR滤波器系数向量b(n)。c.在matlab中本次设计涉及到的几个重要函数（1）采样：[y,fs,nbit]=wavread，返回采样值放在向量y中，fs为采样频率（Hz），nbit为采样位数。（2）快速傅里叶变换：[Y]=fft（y，N）,其中x为有限长序列，N为序列y的长度，Y为序列y的DFT。（3）数字滤波器的频率响应：[h,w]=freqz(b,a,n),

5、返回数字滤波器的n点复频率响应，b和a是滤波器系数的分子和分母向量；h是复频率响应，w是频率点，n默认为512。（模拟滤波器的频率响应用函数freqs）。12李永棠基于MATLAB的语音信号设计4.课程设计内容基本要求：（1）录制一段自己的语音信号，并对录制的信号进行采样；（2）画出采样后语音信号的时域波形和频谱图；（3）给定滤波器的性能指标，采用窗函数法或双线性变换设计FIR滤波器或设计IIR滤波器，并画出滤波器的频率响应；（4）然后用自己设计的滤波器对采集的语音信号进行滤波，画出滤波后信号的时域波

6、形和频谱，并对滤波前后的信号进行对比，分析信号的变化；（5）回放语音信号。扩展部分：（1）设计GUI界面进行对该语音处理系统的界面化管理。（2）实现matlab界面中的直接录音功能5.详细设计过程与调试分析5.1语音信号的采集录制原始语音：我们利用Windows的录音机，录制了一段自己的声音(“同志们，我们一起吃饭吧”，单声道)，时间在3s左右。接着保存为WAV格式（存在matlab软件的work文件夹下，命名为“chifan”）。5.2原始语音信号的分析利用函数wavread对语音信号进行采样，并读

7、取语音的相关信息，并通过函数sound听到采样后自己所录的一段声音，并画出它的时域和频域的图形。12李永棠基于MATLAB的语音信号设计程序如下：fileName='chifan.wav';[y1,Fs,nbits]=wavread(fileName);fprintf('语音"%s"的信息：',fileName);fprintf('语言长度=%g秒',length(y1)/Fs);fprintf('取样频率=%g取样点/秒',Fs);fprintf('解析度=%g位元/取样点',nbi

8、ts);sound(y1,Fs);M=length(y1);time=0:1/Fs:(M-1)/Fs;%时间轴的向量yf1=fft(y1);figure;subplot(2,2,1);plot(time,y1);axis([0,4,-1.5,1.5]);title('原信号的时域图');xlabel('时间');ylabel('幅度');subplot(2,2,2);plot(0:Fs/(M):Fs*(M-1)/M,abs(yf1));axis([0,2