- 1、本文档共15页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
matlab语音信号采集及处理
电子科技大学
课 程 设 计 报 告
课程名称: 信号与系统
设计名称: 语音信号的采集和处理
姓 名: 肖燕平
学 号: 2011019090028
班 级: 通信九班
指导教师: 何春
起止日期: 2012.12.15-2012.12.20
基语音信号的采集和声音是由物体的振动产生,以声波的形式在介质中传播,介质主要可分为固 体, 液体以及气体。 声波振动内耳的听小骨, 这些振动被转化为微小的电子脑波, 它就是我们觉察到的声音。 内耳采用的原理与麦克风捕获声波或扬声器的发音一 样,它是移动的机械部分与气压波之间的关系。在国际标准中,人声的频率范围 是 300Hz~3400Hz,不同的人或乐器产生的声音频率不一致,通过对声音信号的 研究能够更好的处理声音信号的处理以及传输。Matlab 作为一款主要面对科学 计算、可视化以及交互式程序设计的高科技计算软件,能够很好的完成对声音信 号的分析和处理,快速的得出声音信号的时域图以及频域图。关键字:声音频率时域图频域图
1绪论
掌握语音信号采集的方法
掌握一种语音信号基音周期提取方法
了解Matlab的编程方法
本课题的设计要求方案
一:使用wavrecord录入自己的语音信号wavplay
函数进行播放。
二:使用plot再画出语音信号的时域波形频谱滤波器画出滤波后的信号时域、频域图
300Hz~3400Hz,于是给原始语音信号加入高频噪声,分析语音信号的特点
六:设计低通计滤波器将高频噪声滤除
2.1本课题的设计要求
1,使用wavrecord录入自己的语音信号wavplay函数进行播放。程序实现如下:
fs=8000; %采样频率为8000hz
y=wavrecord(3*fs,fs); %录入3s的声音
wavplay(y,fs); %播放已录入的声音信号
file=voice.mat; %创建一个数据文件,且赋给变量file
save(voice.mat,y); %将语音数据保存到数据文件voice.mat中
load voice.mat; %读出voicet.mat的语音数据
2,使用plot再画出该语音信号的时域波形
程序代码为:
subplot(211); %对图形窗口进行分割
plot(y);
xlabel(time n);
ylabel(amplitude);
title(原始语音信号);
再对原始波形进行用fft函数傅里叶变换后得到向量,使用plot画出其频谱。因为直接用fft得出的数据与频率不是对应的,fftshift可以将fft的直流分量移到频谱中心,即让正半轴部分和负半轴部分的图像分别关于各自的中心对称。fftshift对fft进行变换。
程序代码为:
a=fftshift(fft(y,l)/l); %对语音信号进行离散时间傅里叶变换,且将fft的直流分量移 到频谱中心
fd=fs/l;
fx=fd*(-l/2:l/2-1); %将横轴变为频率轴
fy=abs(a); %将纵轴变为频率幅度轴
subplot(212);
plot(fx,fy);
xlabel(频率Hz);
ylabel(幅度);
title(原始信号频谱);
3,设计切比雪夫滤波器利用cheb1ord计算滤波器阶数:[N, Wn] =cheb1ord(Wp, Ws, Rp, Rs);利用cheby1计算滤波器系数:[b,a] = cheby1(N,Wn)。
程序为:
fph=1000;
fsh=1200;
rp=2;
rs=15;
omegaph=fph/(fs/2); %归一化
omegash=fsh/(fs/2);
[N,wn]=cheb1ord(omegaph,omegash,rp,rs);
[BZ,AZ]=cheby1(N,rp,wn,low);
yl=filter(BZ,AZ,y);
subplot(211)
文档评论(0)