基于MATLAB的语音信号的采集与分析.doc

基于MATLAB的语音信号的采集与分析.doc

  1. 1、本文档共11页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
PAGE PAGE 1 基于MATLAB的语音信号的采集与分析 我们通过学习使用MATLAB仿真软件实现语音信号分析,加深对信号与系统这门课程所学习内容的理解,锻炼自学能力和动手能力。我们通过电脑的声卡采集声音信号,借助已有的知识和MATLAB对采集的声音信号进行时域波形和频域频谱的显示,研究男女声信号的差别,通过查找资料提取声音信号的基音频率,并通过大量测试确定门限值来自动判别男女声信号,最后对信号进行降采样处理并播放,重新绘制频谱图分析,验证抽样定理。 关键词:MATLAB仿真、频谱分析、基音频率、降采样、抽样定理。 1.音频信号的采集 我们所要分析的语音信号需要自行采集,所以信号分析的第一步就是采集音频信号。实现音频信号的采集最简单的办法就是通过电脑的声卡直接进行采集,这样采集到的音频信号虽然已经被转化成了数字量存储在电脑中,但通过查询我们了解到电脑录音所使用的采样频率是为44100Hz,完全保证了人类耳朵能听到所有声音频率分量的无失真采集,如果通过MATLAB软件采集还能够调节采样频率,所以能够完全满足我们实验的要求。 1.1使用MATLAB采集语音信号 通过上网查询,我们了解到MATLAB有自带的音频信号采集函数audiorecord(),通过它可以在程序运行时即时采集音频信号进行存储并处理,并且可以通过改变输入参数来改变采样频率,可以直接模拟降采样的过程,直观地验证抽样定理。 但鉴于我们需要重复进行试验和演示,即时采集信号显得繁琐且不必要,而且会增加我们非界面化编程的难度,所以我们放弃了这种方法。 1.2使用电脑录音机采集语音信号 通过电脑自带的录音机软件可以实现更简单的音频信号采集操作,虽然采样频率不可调节,但其固有的采样频率完全满足了我们对所采集信号的要求,可以通过MATLAB的降采样处理的到较低采样频率的信号。这样采集的音频信号会直接以文件的形式存储在电脑中,方便我们随时进行调用,方便分析与演示,所以我们决定采用这种方式实现语音信号的采集。 2.音频信号的图形显示 2.1音频信号的时域波形 通过电脑声卡采集后的声音信号在电脑中是以数字的形式存储的,用audioread()函数可以将音频文件读入MATLAB中,以矩阵的形式存储在变量中,其长度N即为采样点数,并且系统默认的采样信号为双声道,所以该矩阵是一个N×2的矩阵。通过观察,我们的录音文件的两列数据几乎完全相同,所以我去掉了一列的数据将其变为单声道来处理(通过实际实验,是否加入这个操作并不影响波形的显示,但在降采样分析中为保证降采样后信号长度进行插值处理时该操作保证了矩阵维度的匹配,是必要的)。 [x,fs]=audioread(file); % 读入声音文件 data=x(:,1); %取单声道 在时域中采样点数本身就代表了一个时间的概念,在量化过程中表现为矩阵的长度,而通过处理我们可以把采样点数转化为时间:t=N/fs,其中fs是采样频率。通过这样的转化我们可以把波形图的横轴以时间的形式表示出来。 n=0:length(x)-1; %建立一个信号等长的序列 time=n/fs; %建立时间序列,作为横坐标 用plot()函数进行作图,就得到了音频信号的时域波形,如图1所示。 figure(1); %图1:时域波形图 plot(time,data); %作图 title('音频信号波形图') %标题 xlabel('时间/s'); %标注横坐标 ylabel('幅值'); %标注纵坐标 grid on; %打开网格线 图1 音频信号时域波形图 2.2音频信号的频域频谱图 音频信号函数在进行傅里叶变换后可以在频域中作出频谱图,便于我们的分析。 MATLAB提供了一个快速傅里叶变换函数fft(),通过它我们可以把从文件中读取并存储的离散信号进行傅里叶变换得到其在频域中所对应矩阵。 我们知道,经过傅里叶变换我们得到的是复数,在MATLAB中是一个以复数为元素的矩阵,包含了幅度和相位的信息。开始时没有注意到这一点,我们直接用变换后的数据作图,得到了一个复平面中大量点的连线,没有参考意义。 随后通过查询我们了解到需要对傅里叶变换后的矩阵通过mag()函数进行取模处理得到幅度,才能作出频谱图。由于傅里叶变换的对称性,只有一半的数据是有效的,另一半的数据是对前一半的重复,所以我们得到了如下的代码,作出的频谱图如图2所示。 N=length(data); %取信号矩阵的长度 Y1=fft(data,N); %N点傅里叶变

文档评论(0)

anma + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档