基于某MATLAB地加噪语音信号地滤波

基于某MATLAB地加噪语音信号地滤波

ID:47025401

大小:576.01 KB

页数:12页

时间:2019-06-28

基于某MATLAB地加噪语音信号地滤波_第1页
基于某MATLAB地加噪语音信号地滤波_第2页
基于某MATLAB地加噪语音信号地滤波_第3页
基于某MATLAB地加噪语音信号地滤波_第4页
基于某MATLAB地加噪语音信号地滤波_第5页
资源描述:

《基于某MATLAB地加噪语音信号地滤波》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、实用文档《计算机仿真技术》基于MATLAB的加噪语音信号的滤波学生姓名:专业:电子信息工程班级:学号:指导教师:完成时间:2017年12月标准文案实用文档一.滤波器的简述在MATLAB环境下IIR数字滤波器和FIR数字滤波器的设计方法即实现方法,并进行图形用户界面设计,以显示所介绍迷你滤波器的设计特性。在无线脉冲响应(IIR)数字滤波器设计中,先进行模拟滤波器的设计,然后进行模拟-数字滤波器转换,即采用脉冲响应不变法及双线性Z变化法设计数字滤波器,最后进行滤波器的频带转换。在有限脉冲响应(FIR)数字滤

2、波器设计中,讨论了FIR线性相位滤波的特点和用窗口函数设计FIR数字滤波器两个问题。两类滤波器整个过程都是按照理论分析、编程设计、集体实现的步骤进行的。为方便分析直观者直观、形象、方便的分析滤波器的特性,创新的设计出图形用户界面---滤波器分析系统。整个系统分为两个界面,其内容主要包括四个部分:System(系统)、Analysis(分析)、Tool(工具)、Help(帮助)。数字滤波在DSP中占有重要地位。数字滤波器按实现的网络结构或者从单位脉冲响应,分为IIR(无限脉冲响应)和FIR(有限脉冲响应)

3、滤波器。如果IRR滤波器和FIR滤波器具有相同的性能,那么通常IIR滤波器可以用较低的阶数获得高的选择性,执行速度更快,所有的存储单元更少,所以既经济又高效。二.设计要求1.在matlab平台上录制一段语音信号;2.完成语音信号的谱分析;3.对语音信号进行加噪以及加噪后信号的谱分析;4.选择合适的滤波器进行滤波,确定相关指标;5.实现滤波过程,显示滤波后的结果,并进行谱分析。标准文案实用文档三.实验内容与步骤1、语音信号的录入打开matlab平台,先使用R=audiorecorder(44100,16,

4、2)函数创建一个保存音频信息的对象,其中44100表示采样频率为44100Hz,16表示用16位存储,2表示两通道。再使用record(R)开始录音,对着电脑麦克风采集声音,录入的音频内容为“MATLAB的课程大作业”。stop(R)语句停止录入。然后将音频信息以数字矩阵存储。最后使用wavwrite函数保存这段音频。原始语音R=audiorecorder(44100,16,2);record(R);stop(R);myword=getaudiodata(R);plot(myspeech)wavwrit

5、e(myspeech,44100,16,'myword');程序图标准文案实用文档波形图:由输出的波形图可以看出,当采集“MATLAB的课程大作业”这句音频信息时,波动明显,但受到环境的影响其他地方也有轻微的波动。1、语音信号的谱分析wavread读取存储在电脑中的音频。suond函数回放语音信号,fft(y,n)则是对n点进行傅里叶变换,实现从时域到频域。然后用plot函数画出语音信号的时域波形和频域波形图。[y,fs,bits]=wavread('D:Matalbe123binsouds.wa

6、v');sound(y,fs);n=length(y);y_f=fft(y,n);f=fs*(0:n/2-1)/n;subplot211;plot(y);xlabel('时间s');ylabel('幅值');标准文案实用文档title('加噪前的时域波形');subplot212;plot(f,abs(y_f(1:n/2)));xlabel('频率Hz');ylabel('频率幅值');title('加噪前的频谱图');程序图:标准文案实用文档时域波形与频谱图如下:由时域波形图可以看出声音信号主要集中在

7、了2.5s~5.0s之间,在其他时间内只是有少量的波动。由频域波形图可以看出该声音信号的频率与人的发声频率相符合。1、语音加噪及谱分析对原始的语音信号进行加噪处理,使用randn函数产生一个与音频信号等长度的高斯随机噪声信号(噪声大小取决于随机函数的幅度倍数)。然后通过信号的叠加y_z=y+noise;,生成一个新的信号。再通过傅里叶变换得到信号在频域上的波形,最后画出时域与频域波形。L=length(y);noise=0.1*randn(L,2);y_z=y+noise;sound(y_z,fs);n

8、=length(y);y_zf=fft(y_z,n);f=fs*(0:n/2-1)/n;subplot211;plot(y_z);xlabel('时间s');标准文案实用文档ylabel('幅值');title('加噪后的时域波形');subplot212;plot(f,abs(y_zf(1:n/2)));xlabel('频率Hz');ylabel('频率幅值');title('加噪后的频谱图');程序图:波形图如下:当执行上述程序时,会听

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。