- 1、本文档共11页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 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点傅里叶变
您可能关注的文档
- 基于MATLAB的图像分块方法.pdf
- 基于MATLAB的微分方程求解.pdf
- 基于MATLAB的维纳滤波在图像恢复中的应用.pdf
- 基于MATLAB的相平面分析方法的研究.pdf
- 基于MATLAB的杨氏双缝干涉实验仿真.pdf
- 基于MATLAB的语音信号采集和分析系统的可视化设计.pdf
- 基于MATLAB的语音信号去噪方法应用.pdf
- 基于MATLAB的圆度评定方法.pdf
- 基于matlab的正交振幅调制与解调.doc
- 基于Matlab的直流电机速度控制.pdf
- 地理-内蒙古锡林郭勒盟2023~2024学年第一学期高三年级全盟统考试题和答案.docx
- 历史|湖北元月调考暨湖北省部分市州元月高三期末联考历史试卷及答案.docx
- 历史-内蒙古锡林郭勒盟2023~2024学年第一学期高三年级全盟统考试题和答案.docx
- 数学(理科)-内蒙古锡林郭勒盟2023~2024学年第一学期高三年级全盟统考试题和答案.docx
- 历史丨金太阳24~239C湖北省十堰市2024届高三上学期1月调研考试历史试卷及答案.docx
- 青岛某高层办公楼地下防水施工方案(sbs卷材).docx
- 土木工程实习报告表格.doc
- 工程流体力学教学课件ppt作者闻建龙工程流体力学习题答案部分.doc
- 围手术期管理制度(1)_0.doc
- 2023年公务员多省联考《申论》题(内蒙古旗县卷)及参考答案.doc
文档评论(0)