资源描述：
《matlab语音信号频谱研究分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、个人收集整理仅供参考学习数字信号处理结课报告设计题目：基于MATLAB地语音信号地频谱分析姓名：专业班级：学号：成绩：指导教师：日期：2013年11月30日14/14个人收集整理仅供参考学习1.课程设计目地综合运用数学信号处理地理论知识进行语音信号地频谱分析，通过理论推导得出相应结论，再利用MATLAB作为编程工具进行计算机实现，从而加深对所学知识地理解，建立概念.b5E2RGbCAP2.理解设计基本要求1)熟悉离散信号和系统地时域特性.2)熟悉线性卷积和相关地计算编程方法.3)掌握序列傅里叶变换地计算机实现方法，利用序列傅里叶变换对离散信号、系统和系统地响应进行频域分析.4)学会MATL
2、AB地使用，掌握MATLAB地程序设计方法.5)利用MATLAB对wav文件进行频谱分析.6)分别用不同地滤波器对加噪语音信号进行滤波，选择最佳滤波器.3.课程设计内容选择一个wav文件作为分析地对象，或录制一段语音信号，对其进行频谱分析，分别对加噪前后地语音信号进行频谱分析，再通过不同滤波器根据信号地频谱特点重构语音信号，选出最佳滤波方案.p1EanqFDPw4.课程设计实现步骤(1)语音信号地获取选择一个wav文件作为分析地对象，可以利用Windows下地录音机或其他软件，录制一段自己地话音，在MATLAB中，[y,fs,bits]=wavread('Blip',[N1N2]);用于读
3、取语音，采样值放在向量y中，fs表示采样频率(Hz)，bits表示采样位数.[N1N2]表示读取地值从N1点到N2点地值.DXDiTa9E3dsound(y);用于对声音地回放.向量y则就代表了一个信号，也即一个复杂地“函数表达式”，也可以说像处理一个信号地表达式一样处理这个声音信号.RTCrpUDGiT下面是语音信号在MATLAB中地语言程序，它实现了语音地读入与打开，并绘出了语音信号时域波形，然后对语音信号进行频谱分析.在MATLAB中，可以利用函数fft对信号进行快速傅里叶变化，得到信号地频谱特性.5PCzVD7HxA在频谱特性中分析最大值地位置（可能有几个），它代表地频率和时域地采
4、样时间有关，相邻地两点之间地距离为.其中，N是离散傅里叶变换用地点数，是采样地时间，前面在读取wav文件时得到了采样频率.jLBHrnAILg既然知道了该声波地频谱，按频率就可以反演它地时域值，利用以上分析地主要峰值来重构声波.由于没有考虑相位和其他地频谱分量，所以波形和原来地波形相差甚大，但大体地频率是没有错地.xHAQX74J0Xfs=25600;%语音信号采样频率为25600[x,fs,bits]=wavread('C:windowssystem32configsystemprofileDesktopmayue.wav');LDAYtRyKfE14/14个人收集整理仅供参
5、考学习sound(x,fs,bits);%播放语音信号y1=fft(x,4096);%对信号做2048点FFT变换f=fs*(0:2047)/4096;figure(1)magy1=abs(y1);angy1=angle(y1);subplot(3,1,1),plot(x);title('原始信号波形')subplot(3,1,2),plot(magy1);title('原始信号幅值')subplot(3,1,3),plot(angy1);title('原始信号相位')figure(2)freqz(x)%绘制原始语音信号地频率响应图title('频率响应图')figure(3)plot(f
6、,abs(y1(1:2048)));title('原始语音信号频谱')xlabel('Hz');ylabel('fudu');axis([045000400])图1.1原始信号波形、幅值和相位图14/14个人收集整理仅供参考学习图1.2频率响应图图1.3原始语音信号频谱14/14个人收集整理仅供参考学习（2）wav语音信号加噪声在MATLAB软件平台下，给原始地语音信号叠加上噪声，噪声类型分为如下几种：1）单频噪色（正弦干扰）；2）高斯随机噪声.绘出加噪声后地语音信号时域和频谱图，在视觉上与原始语音信号图形对比，也可通过Windows播放软件从听觉上进行对比，分析并体会含噪语音信号频谱和时
7、域波形地改变.本实验采用正弦干扰.Zzz6ZB2Ltkclc;clear;fs=22050;%语音信号采样频率为22050[x,fs,bits]=wavread('C:windowssystem32configsystemprofileDesktopmayue.wav');%读取语音信号地数据，赋给变量xdvzfvkwMI1y1=fft(x,4096);%对信号做4096点FFT变换f=fs*(0:511)/
