MATAB课程设计--语音信号的时域频域分析.doc

精品文档，知识共享！ 语音信号的采集 为了将原始模拟语音信号变为数字信号，必须经过采样和量化两个步骤，从而得到时间和幅度上均为离散的数字语音信号。取样是将时间上连续的语音信号离散化为一个样本序列。根据取样定理，当取样频率大于信号的两倍带宽时，取样过程不会丢失信息，且从取样信号中可以精确地重构原始信号波形。 经过预滤波和采样后，由A/D变换器变换为二进制数字码。这种防混叠滤波通常与模数转换器做在一个集成块内，因此目前来说，语音信号的数字化的质量还是有保证的。市面上购买到的普通声卡在这方面做的都很好，语音声波通过话筒输入到声卡后直接获得的是经过防混叠滤波、A/D变换、量化处理的离散的数字语音信号。 实际工作中，我们可以通过很多种方式录制语音信号，如手机、电脑、MP3、录音笔等，然后通过转换工具将音频格式转换为MATLAB可以识别的音频文件格式，这里的转换工具可以是千千静听或者其它的专用音频文件格式转换软件，可以获得任何我们想要的音频格式。本设计中我们利用windows自带的录音机录制语音文件，图1-1是基于PC机的语音信号录制过程，声卡可以完成语音波形的A/D转换，获得WAV文件，为后续的处理储备原材料。 wavA/D采样滤波麦克风声卡声音 wav A/D 采样 滤波 麦克风 声卡 Windows自带的录音机 Windows自带的录音机 图1-1 基于PC机的语音信号录制过程 本设计先用电脑自带的录音机录制自己的一段语音，保存在MATLAB软件的工作目录下，命名为“abc.wav”，然后经过MATLAB软件绘制出其时域波形图和频谱图。之后通过rand随机函数产生一个噪声，叠加到原信号中去，作为滤波器的输入。 2．语音信号的时、频域分析 2.1 语音信号的时域分析 语音信号的时域分析就是分析和提取语音信号的时域参数。进行语音分析时，最先接触到并且也是最直观的是它的时域波形。语音信号本身就是时域信号，因而时域分析是最早使用，也是应用最广泛的一种分析方法，这种方法直接利用语音信号的时域波形。时域分析通常用于最基本的参数分析及应用，如语音的分割、预处理、大分类等。 这种分析方法的特点是： = 1 \* GB3 ①表示语音信号比较直观、物理意义明确。 = 2 \* GB3 ②实现起来比较简单、运算量少。 = 3 \* GB3 ③可以得到语音的一些重要的参数。④只使用示波器等通用设备，使用较为简单等。 MATLAB数据采集箱中提供的进行语音信号分析的函数命令如下： wavread?:wavread?用于读取扩展名为“.wav”的声音文件。其调用形式为:?y?=?wavread?(‘filename’)?。其作用是读取wave?文件，将读取的采样数据送到y?中。 sound:音频信号是以向量的形式表示声音采样的。sound?函数用于将向量转换为声音，其调用形式为:sound?(y?,fs)?，作用是向扬声器送出向量y?中的音频信号(采样频率为fs)?。通过wavread和plot()函数即可显示语音信号的时域波形。图形如下图2-1上面的图所示。 2.2 语音信号的频域分析 语音信号的频域分析就是分析语音信号的频域持征。从广义上讲，语音信号的频域分析包括语音信号的频谱、功率谱、倒频谱、频谱包络分析等，而常用的频域分析方法有带通滤波器组法、傅里叶变换法、线件预测法等几种。本文介绍的是语音信号的傅里叶分析法。因为语音波是一个非平稳过程，因此适用于周期、瞬变或平稳随机信号的标准傅里叶变换不能用来直接表示语音信号，而应该用短时傅里叶变换对语音信号的频谱进行分析，相应的频谱称为“短时谱 ” 在MATLAB的信号处理工具箱中函数FFT和IFFT用于快速傅立叶变换和逆变换。函数FFT用于序列快速傅立叶变换，其调用格式为y=fft(x)，其中，x是序列，y是序列的FFT，x可以为一向量或矩阵，若x为一向量，y是x的FFT且和x相同长度；若x为一矩阵，则y是对矩阵的每一列向量进行FFT。函数FFT的另一种调用格式为y=fft(x,N)，式中，x，y意义同前，N为正整数。函数执行N点的FFT，若x为向量且长度小于N，则函数将x补零至长度N；若向量x的长度大于N，则函数截短x使之长度为N；若x 为矩阵，按相同方法对x进行处理。 利用上述函数即可画出语音信号的频谱图，如下图2-1下面的图所示。 3.对语音信号进行加噪处理 利用MATLAB中的随机函数(randn)产生噪声加入到语音信号中，模仿语音信号被污染，并对其进行频谱分析。程序见附录。 加入噪声后的时域图和频谱图如下： 图3-1 4. 数字滤波器的设计 4.1 滤波器的设计选择 设计数字滤波器的任务就是寻求一