有关基于MATLAB实现对语音信号频谱毕业论文写作资料-论文写作网

简介:关于本文可作为相关专业信号语音论文写作研究的大学硕士与本科毕业论文信号语音论文开题报告范文和职称论文参考文献资料。

信号语音论文范文

摘要：语音信号处理技术是语音处理领域中新近发展起来的一个学科分支．MATLAB是一个数据分析和处理功能十分强大的工程实用软件,运用它来进行语音信号的采集、分析和处理相当便捷.文章介绍了在Matlab环境中如何采集语音信号和语音信号采集后的频谱分析处理,并通过实例分析了语音信号处理的Matlab.

关键字：MATLAB,语音信号；频谱分析；傅里叶变换

中图分类号：TN911.6

文献标识码：A

文章编号：1673-1131(2010)04-015-03

一、引言

随着计算机技术和信息技术的发展,语音交互已经成为人机交互的必要手段,而语音信号的采集和处理是人机交互的前提和基础.Matlab是美国Math Works公司推出的一种面向工程和科学计算的交互式计算软件,它以矩阵运算为基础．把计算、可视化、程序设计融合到了一个简单易用的交互式工作环境中.同时由于Matlab是一个数据分析和处理功能十分强大的工程实用软件．它的信号处理与分析工具箱为语音信号分析提供了十分丰富的功能函数,利用这些功能函数可以快捷而又方便地完成语音信号的处理和分析以及信号的可视化,使人机交互更加便捷.

二、语音信号的采集

Windows附件的娱乐中有个录音机程序,通过它可以驱动声片采集、播放和简单处理语音信号.为了清晰的阐明对语音信号进行频谱分析,首先使用Windows的录音机功能可以用麦克风直接录制一段女声“世界上最遥远的距离不是生与死,而是我就站在你面前你确不知道我爱你”,并保存为W论文范文格式文件,供MATLAB相关函数直接读取、写入或播放.本文以W论文范文格式语音信号作为分析处理的输入数据.

例：在MATLAB中使用w论文范文read函数读取录制得到的声音采样频率Fs为llkHz,且声音是8位单声道的声音文件"test.w论文范文”.并利用w论文范文play函数播放声音,可清晰地听到语音.然后绘制显示出语音数据的波形图（如图l所示）.

Clear

S等于w论文范文read(&,acute,C:\Documents and Settings\administrator\桌面 \test.w论文范文)；%读取语音数据

Plot(S),title&,acute,原始语音信号画波形图

W论文范文play(S),%播放语音

三、语音信号分析及MATLAB的实现

3.1语音信号的频谱分析过程

傅里叶频谱分析是语音信号频域分析中广泛采用的一种方法.语音波是一个非平稳过程,因此适用于周期、瞬变或平稳随机信号的标准傅里叶变换不能直接表示语音信号,而应该用短时傅里叶变换对话音信号的频谱进行分析,相应的频谱称为“短时谱”.

进行频谱分析时,在时域数据进行短时FFT处理之前都要进行加窗处理.在FFT处理之后,普通频谱分析可以进行频域上的滤波处理,从而使频谱更加平滑.最后IFFT处理观察恢复后的时域信号图形,频谱分析过程.

3.2短时傅里叶变换（频谱分析）

傅里叶频谱分析的基础是傅里叶变换,用傅里叶变换及其反变换可以求得傅里叶谱、自相关函数、功率谱．倒谱.由于语音信号的特性是随着时间缓慢变化的,由此引出语音信号的短时分析.如同在时域特征分析中用到的一样,这里的傅里叶频谱分析也采用相同的短时分析技术.

信号x(n)的短时傅里叶变换定义为

x.(P∥)等于∑x(m)w(n - m)e-wn

～∞

式中．w(n)为窗口函数.

可以从两个角度理解函数X,,(erv)的物理意义：第一种解释是,当n固定时,如n-n,则(ejW^)是将窗函数的起点移至no处截取信号x(n),再做傅里叶变换而得到的一个频谱函数.这是直接将频率轴方向来理解的.另一种解释是从时间轴方向来理解,当频率固定时,例如珊等于q．则(P)可以看作是信号经过一个中心频率为的带通滤波器产生的输出.这是因为窗口函数∞（通常具有低通频率响应,而指数Pp对语音信号x(n有调制的作用,使频谱产生移位,即将x(n)频谱中对应于频率的分量平移到零频. 3.3短时傅里叶反变换（倒谱分析）傅里叶变换建立了信号从时域到频域的变换桥梁．而傅里叶反变换则建立了信号从频域到时域的变换桥梁这两个域之间的变换为一对一映射关系.我们知道,x．(ejm)可以看作是加窗后函数的傅里叶变换,为了实现反变换,将X.(e&,acute,“)进行频率采样,即令mk等于2丌k//,则有

x.0“)等于∑b-等

崩等于叫∞

其中,L为频率采样点数.将J．(e)在时域n上每隔R个样本采样,则可令

Yr(eWk)等于X,R(ejwk),n等于rR,r等于l,2,等

用这r0^）求出其离散傅里叶反变换y,如下：

(eF")P jw—

+∞

而y等于∑

∞

可以证明．x(n)和y(n)之间只相差一个比例因子,它们的关系如下：

y(n)等于 x(n)Wp j0)

3.4 MATLAB处理语音信号

Matlab语言是一种数据分析和处理功能十分强大的计算机应用软件,它可以将声音文件变换为离散的数据文件,然后利用其强大的矩阵运算能力处理数据,如数字滤波、傅里叶变换、时域和频域分析、声音回放以及各种分析图的呈现等等.下面以语音信号的波形图,频谱图和倒谱图（如图3所示）分析为例来说明Matlab在语音信号处理中的具体实现方法. %分段fourier变换,每段长为512．每段加矩形窗,连续断之间重叠的长度100

S等于w论文范文read(&,acute,C:\Documents and Settings\administrator\桌面\test.w论文范文’)；

subplot(3,l,l),plot(S),title(’波形图’)

w论文范文play(S),%播放语音

n等于size(S),

d等于100,%连续段之间的重叠长度

L等于512,

k等于L-d,

t等于fix(n/k),

frame_length等于512,

%采用滑动窗进行短时fft变换

r等于(rectwin(frame_length》,

for a-l:t

nl等于(L—d)+(a—1)+l；

n2等于(L-d)*(a-l)+framelength,

y等于S(nl:n2),

论文范文等于fft(y&,acute,.*r&,acute,,512)；

Xl(nl:n2)等于论文范文(l:frame_length),

end

FFTSIZE等于8000, %FFT长度为8000点,每点对应于1Hz

Y等于zeros(FFTSIZE,l),

Y等于20*loglO(abs(Xl,

subplot(3,1,2),plot( Y(1:400 0,title(‘频谱图’)；%画前4000点

%IFFT变换

Tl等于length(XI),

tl等于fix(Tl/L),

for i等于l:tl

z等于ifft(Xl“i—1)+L+1:1*L),51 2),

z-real(z./r&,acute,),

y“i-1)+k+1:i+k)等于z(1:k),

end

subplot(3,1,3),plot(y),title(倒谱图)

四、结束语

本文主要是基于高级面向对象开发语言MATLAB的基本特征,以及MATLAB强大的工具箱功能,实现信号的预处理,信号的时域处理和信号的频域处理.其中实现信号的时域处理过程中,设计了加窗处理；在信号频域处理过程中,实现对低于20kHz语音信号进行频谱分析.同时．通过上述实例说明,基于MATLAB可以很好的达到对语音信号的频谱的分析处理.

参考文献

[1]姚天任,数字语音处理[M]

武汉：华中科技大学出版社,1991

[2]刘卫国．MATLAB程序设计教程[M]

北京：中国水利水电出版社,2005

[3]张登奇,蒋敏．基于MATLAB的音频信号处理技术实现[J] 电脑知识与技术,2009 (2)