完整word版基于MATLAB的语音信号的采集与处理详解.docx
- 文档编号:4936222
- 上传时间:2022-12-11
- 格式:DOCX
- 页数:11
- 大小:202.14KB
完整word版基于MATLAB的语音信号的采集与处理详解.docx
《完整word版基于MATLAB的语音信号的采集与处理详解.docx》由会员分享,可在线阅读,更多相关《完整word版基于MATLAB的语音信号的采集与处理详解.docx(11页珍藏版)》请在冰豆网上搜索。
完整word版基于MATLAB的语音信号的采集与处理详解
数字信号处理
课程设计
题目:
基于MATLAB的语音信号的采集与处理
学院:
皖西学院
专业:
通信工程
班级:
通信1001班
学号:
20100134612010013494
姓名:
刘敏纵大庆
指导教师:
何富贵
摘要:
本次课程设计题目为<〈基于MATLAB的语音信号的采集与处理>〉.首先我们利用计算机上的录音软件获得语音信号,然后利用MATLAB对语音信号进行分析和处理,采集语音信号后,利用MATLAB软件平台进行频谱分析;并对所采集的语音信号加入干扰噪声,对加入噪声的信号进行频谱分析,设计合适的滤波器滤除噪声,恢复原信号!
1。
背景…………………………………………………………………1
2.设计目的……………………………………………………………2
3.设计原理……………………………………………………………2
4。
设计过程……………………………………………………………3
5。
实验代码及结果………………………………………………………4
5。
1语音信号的采集……………………………………………………4
5。
2语音信号加噪与频谱分析…………………………………………………7
5.3巴特沃斯滤波器的设计…………………………………………………9
5。
4比较滤波前后语音信号波形及频谱………………………………………10
6。
收获与体会……………………………………………………………12
参考文献…………………………………………………………………13
1.引言
数字信号处理是利用计算机或专用处理设备,以数值计算的方法对信号进行采集、抽样、变换、综合、估值与识别等加工处理,借以达到提取信息和便于应用的目的。
它在语音、雷达、图像、系统控制、通信、航空航天、生物医学等众多领域都获得了极其广泛的应用。
具有灵活、精确、抗干扰强、度快等优点。
数字滤波器,是数字信号处理中及其重要的一部分.随着信息时代和数字技术的发展,受到人们越来越多的重视。
数字滤波器可以通过数值运算实现滤波,所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活不存在阻抗匹配问题,可以实现模拟滤波器无法实现的特殊功能。
数字滤波器种类很多,根据其实现的网络结构或者其冲激响应函数的时域特性,可分为两种,即有限冲激响应(FIR,FiniteImpulseResponse)滤波器和无限冲激响应(IIR,InfiniteImpulseResponse)滤波器。
FIR滤波器结构上主要是非递归结构,没有输出到输入的反馈,系统函数H(z)在z=∞处收敛,极点全部在z=0处(因果系统),因而只能用较高的阶数达到高的选择性。
FIR数字滤波器的幅频特性精度较之于IIR数字滤波器低,但是线性相位,就是不同频率分量的信号经过FIR滤波器后他们的时间差不变,这是很好的性质。
FIR数字滤波器是有限单位脉冲响应有利于对数字信号的处理,便于编程,用于计算的时延也小,这对实时的信号处理很重要。
FIR滤波器因具有系统稳定,易实现相位控制,允许设计多通带(或多阻带)滤波器等优点收到人们的青睐[1]。
IIR滤波器采用递归型结构,即结构上带有反馈环路。
IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成,可以组合成直接型、正准型、级联型、并联型四种结构形式,都具有反馈回路.同时,IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果,如巴特沃斯滤波器等。
2.设计目的
本次课程设计的课题为《基于MATLAB的语音信号采集与处理》,学会运用MATLAB的信号处理功能,采集语音信号,并对语音信号进行滤波及变换处理,观察其时域和频域特性,加深对信号处理理论的理解,并为今后熟。
MATLAB在信号与系统中的应用主要包括符号运算和数值计算仿真分析.由于信号与系统课程的许多内容都是基于公式演算,而MATLAB借助符号数学工具箱提供的符号运算功能,能基本满足信号与系统课程的需求。
例如解微分方程、傅里叶正反变换、拉普拉斯正反变换和z正反变换等。
MATLAB在信号与系统中的另一主要应用是数值计算与仿真分析,主要包括函数波形绘制、函数运算、冲击响应与阶跃响应仿真分析、信号的时域分析、信号的频谱分析、系统的S域分析和零极点图绘制等内容。
数值计算仿真分析可以帮助学生更深入地理解理论知识,并为将来使用MATLAB进行信号处理领域的各种分析和实际应用打下基础.
3。
设计原理
利用MATLAB对语音信号进行分析和处理,采集语音信号后,利用MATLAB软件平台进行频谱分析;并对所采集的语音信号加入干扰噪声,对加入噪声的信号进行频谱分析,设计合适的滤波器滤除噪声,恢复原信号。
语音信号的“短时谱"对于非平稳信号,它是非周期的,频谱随时间连续变化,因此由傅里叶变换得到的频谱无法获知其在各个时刻的频谱特性。
如果利用加窗的方法从语音流中取出其中一个短断,再进行傅里叶变换,就可以得到该语音的短时谱。
MATLAB在信号与系统中的应用主要包括符号运算和数值计算仿真分析。
由于信号与系统课程的许多内容都是基于公式演算,而MATLAB借助符号数学工具箱提供的符号运算功能,能基本满足信号与系统课程的需求。
例如解微分方程、傅里叶正反变换、拉普拉斯正反变换和z正反变换等。
MATLAB在信号与系统中的另一主要应用是数值计算与仿真分析,主要包括函数波形绘制、函数运算、冲击响应与阶跃响应仿真分析、信号的时域分析、信号的频谱分析、系统的S域分析和零极点图绘制等内容.数值计算仿真分析可以帮助学生更深入地理解理论知识,并为将来使用MATLAB进行信号处理领域的各种分析和实际应用打下基础[3].
滤波器的设计可以通过软件或设计专用的硬件两种方式来实现。
随着MATLAB软件及信号处理工具箱的不断完善,MATLAB很快成为应用学科等领域不可或缺的基础软件。
它可以快速有效地实现数字滤波器的设计、分析和仿真,极大地减轻了工作量,有利于滤波器设计的最优化。
利用MATLAB中的随机函数产生噪声加入到语音信号中,模仿语音信号被污染,并对其进行频谱分析;设计巴特沃斯滤波器,并对被噪声污染的语音信号进行滤波,对滤波前后信号进行时域分析。
4.设计过程
基于声卡进行数字信号的采集.将话筒插入计算机的语音输入插口上,启动录音机。
按下录音按钮,对话筒说话,说完后停止录音。
要保存文件时,利用了计算机上的A/D转换器,把模拟的声音信号变成了离散的量化了的数字信号,放音时,它又通过D/A转换器,把保存的数字数据恢复为原来的模拟的声音信号。
在Matlab软件平台下可以利用函数wavread对语音信号进行采样,得到了声音数据变量x1,同时把x1的采样频率fs=22050Hz和数据位Nbits=16Bit放进了MATALB的工作空间。
然后通过freqz函数绘制原始语音信号的频率响应图。
接着对语音信号进行频谱分析,在Matlab中可以利用函数fft对信号进行快速傅里叶变换,得到信号的频谱图。
语音信号加噪与频谱分析在Matlab中人为设计一个固定频率5500Hz的噪声干扰信号,噪声信号通常为随机序列,在本设计中用正弦序列代替,干扰信号构建命令函数为d=[Au*sin(2*pi*5500*t)],给出的干扰信号为一个正弦信号,针对上面的语音信号,采集了其中一段。
再对噪音信号进行频谱变换得到其频谱图。
最后通过函数将原信号与噪声信号叠加在一起构成加噪的语音信号,采集他的加噪信号,接着对其进行fourier变换,得到加噪信号的频谱图,最后设计一个高通滤波器,对加噪的信号进行滤波处理,得到的信号与原始信号进行比较!
5.实验代码及结果
5.1语音信号的采集
源程序为:
fs=22050;%语音信号采样频率为22050
x1=wavread('D:
\MATLAB\bin\zx。
wav’);%读取语音信号的数据,赋给变量x1
sound(x1,22050);%播放语音信号
y1=fft(x1,1024);%对信号做1024点FFT变换
f=fs*(0:
511)/1024;%将0到511,步长为1的序列的值与fs相乘并除以1024的值,赋值给f
Plot(x1);%做原始语音信号的时域图形
图5—1:
语音信号的采集
通过freqz函数绘制原始语音信号的频率响应图5-2。
fs=22050;%语音信号采样频率为22050
x1=wavread('D:
\MATLAB\bin\zx。
wav');%读取语音信号的数据,赋给变量x1
sound(x1,22050);%播放语音信号
y1=fft(x1,1024);%对信号做1024点FFT变换
f=fs*(0:
511)/1024;%将0到511,步长为1的序列的值与fs相乘并除以1024的值,赋值给f
Freqz(x1);%绘制原始语音信号的频率响应图
图5-2:
原始语音信号的频率采样
然后对语音信号进行频谱分析,在Matlab中可以利用函数fft对信号行快速傅里叶变换,得到信号的频谱图5—3.
fs=22050;%语音信号采样频率为22050
x1=wavread('D:
\MATLAB\bin\zx。
wav');%读取语音信号的数据,赋给变量x1
sound(x1,22050);%播放语音信号
y1=fft(x1,1024);%对信号做1024点FFT变换
f=fs*(0:
511)/1024;%将0到511,步长为1的序列的值与fs相乘并除以1024的值,赋值给f
subplot(2,1,1);%创建两行一列绘图区间的第1个绘图区间
plot(abs(y1(1:
512)));%做原始语音信号的FFT频谱图
title('原始语音信号FFT频谱’);
subplot(2,1,2);
plot(f,abs(y1(1:
512)));%abs是绝对值,plot是直角坐标下线性刻度曲线
title('原始语音信号频谱’);
图5-3:
原始信号的FFT变换
5.2语音信号加噪与频谱分析
fs=22050;%语音信号采样频率为22050
x1=wavread(’D:
\MATLAB\bin\zx。
wav');%读取语音信号的数据,赋给变量x1
f=fs*(0:
511)/1024;%将0到511,步长为1的序列的值与fs相乘并除以1024的值,赋值给f
t=0:
1/fs:
(length(x1)-1)/fs;%将0到x1的长度减1后的值除以fs的值,且步长为1/fs的值,的序列的值,赋予t
Au=0.005;%噪声幅值
d=[Au*sin(2*pi*5500*t)]';%所加的噪声是正弦信号
x2=x1+d;%将正弦信号噪声加在语音信号上
sound(x2,22050);%播放语音信号
y1=fft(x1,1024);%对信号y1做1024点FFT变换
y2=fft(x2,1024);%对信号y2做1024点FFT变换
figure(4);%创建图形窗1
plot(t,x2);%做加噪后的信号时域
title(’加噪后的信号');
图5—4:
加噪声后的语音信号
5—3:
巴特沃斯滤波器的设计:
wp=0.25*pi;%通带截止频率
ws=0。
3*pi;%阻带截止频率
Rp=1;%通带最大衰减(db)
Rs=15;%阻带最大衰减(db)
Fs=22050;
Ts=1/Fs;
wp1=2/Ts*tan(wp/2);%将模拟指标转换成数字指标
ws1=2/Ts*tan(ws/2);%将模拟指标转换成数字指标
[N,Wn]=buttord(wp1,ws1,Rp,Rs,’s’);%选择滤波器的最小阶数(估算得到Butterworth低通滤波器的最小阶数N和3dB截止频率Wc)
[Z,P,K]=buttap(N);%创建Butterworth低通滤波器原型
[Bap,Aap]=zp2tf(Z,P,K);%将零极点增益转换为普遍分子,分母
[b,a]=lp2lp(Bap,Aap,Wn);%将普遍的分子和分母转换为以Wn为截止频率
[bz,az]=bilinear(b,a,Fs);%用双线性变换法实现模拟滤波器到数字滤波器的转换(模拟转换为数字)
[H,W]=freqz(bz,az);%求频率响应
figure(6)
plot(W*Fs/(2*pi),abs(H))%绘制Butterworth低通滤波器频率响应曲线
grid%添加图格
xlabel(’频率/Hz')
ylabel(’频率响应幅度’)
title('Butterworth’)
图5-5:
巴特沃斯滤波器
5-4比较滤波前后语音信号波形及频谱
wp=0.25*pi;%通带截止频率
ws=0。
3*pi;%阻带截止频率
Rp=1;%通带最大衰减(db)
Rs=15;%阻带最大衰减(db)
Fs=22050;
Ts=1/Fs;
wp1=2/Ts*tan(wp/2);%将模拟指标转换成数字指标
ws1=2/Ts*tan(ws/2);%将模拟指标转换成数字指标
[N,Wn]=buttord(wp1,ws1,Rp,Rs,'s');%选择滤波器的最小阶数(估算得到Butterworth低通滤波器的最小阶数N和3dB截止频率Wc)
[Z,P,K]=buttap(N);%创建Butterworth低通滤波器原型
[Bap,Aap]=zp2tf(Z,P,K);%将零极点增益转换为普遍分子,分母
[b,a]=lp2lp(Bap,Aap,Wn);%将普遍的分子和分母转换为以Wn为截止频率
[bz,az]=bilinear(b,a,Fs);%用双线性变换法实现模拟滤波器到数字滤波器的转换(模拟转换为数字)
[H,W]=freqz(bz,az);%求频率响应
figure(6)
plot(W*Fs/(2*pi),abs(H))%绘制Butterworth低通滤波器频率响应曲线
grid%添加图格
xlabel(’频率/Hz’)
ylabel(’频率响应幅度')
title('Butterworth')
f1=filter(bz,az,x2);%利用上面已做好的滤波器(已bz和az为特征的滤波器)对x2信号进行滤波,赋值给f1
figure(7)
subplot(2,1,1)%创建两行一列绘图区间的第1绘图区间
plot(t,x2);%做加噪后的信号时域图形
title('滤波前的时域波形’);
subplot(2,1,2)%创建两行一列绘图区间的第2绘图区间
plot(t,f1);
title('滤波后的时域波形’);
图5—6:
滤波前后图形比较
通过对比分析可知,滤波后的输出波形和原始语音加噪声信号的图形发生了一些变化,在3000-6000Hz之间的信号消失,出现0—1000和7000-8000之间的信号.滤波后的输出波形明显在对应时间幅度比原语音加噪声信号的要小,而且滤波的效果也与滤波器的选择有关,可以看出滤波器的性能差异以及参数的调节会对滤波器产生一定的影响,并且通过回放可以发现滤波前后的声音有变化。
低通滤波后,已很接近原来的声音,人耳几乎辨别不出。
从频谱图中我们还可以看出声音的能量信号主要集中在低频部分,说明高频语音信号被滤出,滤波器达到既定要求。
信噪比为0。
62。
第6章收获与体会
通过用巴特沃斯滤波器对原始信号进行滤波,对滤波前后的波形分析对比可知,总体可以恢复原始信号波形,效果很明显,稍有点误差,在以后的设计中需要改进。
滤波后的输出波形明显在对应时间幅度比原语音加噪声信号的要小,而且滤波的效果也与滤波器的选择有关,可以看出滤波器的性能差异以及参数的调节会对滤波器产生一定的影响,并且通过回放可以发现滤波前后的声音有变化.低通滤波后,已很接近原来的声音,人耳几乎辨别不出。
通过为期两周的MATLAB课程设计,我对MATLAB这个仿真软件有了更进一步的认识和了解。
在这两周时间里,我通过自己摸索,查阅资料,并且在指导老师朱老师的指导下完成了:
语音信号的采集及分析;给原始信号加上一个高频噪声;设计一个滤波器,滤除高频噪声;并最终将课程设计报告总结完毕.
在整个设计过程中我懂得了许多东西,也培养了独立思考和设计的能力,树立了对知识应用的信心,相信会对今后的学习工作和生活有非常大的帮助,并且提高了自己的动手实践操作能力,使自己充分体会到了在设计过程中的成功喜悦。
虽然这个设计做的不怎么好,但是在设计过程中所学到的东西是这次课程设计的最大收获和财富,使我终身受益。
在没有做课程设计以前,觉得课程设计只是对知识的单纯总结,但是通过这次课程设计发现自己的看法有点太片面,课程设计不仅是对前面所学知识的一种检验,也是对自己能力的一种提高,通过这次课程设计使自己明白了原来的那点知识是非常欠缺的,要学习的东西还很多,通过这次课程设计,明白学习是一个长期积累的过程,在以后的工作和生活中都应该不断的学习,努力提高自己的知识和综合素质。
参考文献
[1]:
高西全丁玉美数字信号处理第三版西安:
西安电子科技大学出版社,2008
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 完整 word 基于 MATLAB 语音 信号 采集 处理 详解