倒谱计算与分析..docVIP

河南师范大学计算机与信息技术学院 《视频语音处理技术》 倒谱计算与分析 学院名称： 计算机与信息工程学院 专业名称： 计算机科学与技术 年级班级： 姓 名： 学 号： 计算机与信息技术学院综合性、设计性实验报告 专业：计算机科学技术 年级/班级：2011级 2012—2013学年第一学期 课程名称 视频语音处理技术 指导教师 张新明 本组成员 学号姓名 1108114153 王彦秋 实验地点 计科楼324 实验时间 项目名称 倒谱计算与分析 实验类型 设计性 实验目的： 对语音信号进行同态分析可得到语音信号的倒谱参数。语音的倒谱是将语音的短时谱取对数后再进行IDFT得到的，所以浊音信号的激励反映在倒谱上是同样周期的冲激，借此，可从倒谱波形中估计出基音周期。对倒谱进行低时窗选，通过语音倒谱分析的最后一级，进行DFT后的输出即为平滑后的对数模函数，这个平滑的对数谱显示了特定输入语音段的谐振结构，即谱的峰值基本上对应于共振峰频率，对于平滑过的对数谱中的峰值进行定位，即可估计共振峰。对于倒谱计算与分析的设计实验可作如下训练： 1、复倒谱的几种计算方法： 2、最小相位信号法和递归法； 3、基音检测； 4、共振峰检测。 二、实验仪器或设备：windowsXP下的Matlab编程环境 三、总体设计（设计原理、设计方案及流程等） 1．复倒谱的几种计算方法： 在复倒谱分析中，z变换后得到的是复数，所以取对数时要进行复对数运算。这时存在相位的多值性问题，称为“相位卷绕”。 设信号为 则其傅里叶变换为 对上式取复对数为 则其幅度和相位分别为： 上式中，虽然 ， 的范围均在 内，但 的值可能超过 范围。计算机处理时总相位值只能用其主值 表示，然后把这个相位 主值“展开”，得到连续相位。所以存在下面的情况： （K为整数） 此时即产生了相位卷绕。下面介绍几种避免相位卷绕求复倒谱的方法。 最小相位信号法 这是解决相位卷绕的一种较好的方法。但它有一个限制条件：被处理的信号想x(n)必须是最小相位信号。实际上许多信号就是最小相位信号，或可以看作是最小相位信号。语音信号的模型就是极点都在z平面单位圆内的全极点模型，或者极零点都在z平面单位圆内的极零点模型。 设信号x(n)的z变换为X(z)=N (z)/ D(z) ，则有 根据z变换的微分特性有 若x(n)是最小相位信号，则 必然是稳定的因果序列。 由Hilbert变换的性质可知，任一因果复倒谱序列都可分解为偶对称分量和奇对称分量之和: 其中 这两个分量的傅里叶变换分别为 的傅里叶变换的实部和虚部。 从而可得 此即复倒谱的性质3，也就是说一个因果序列可由其偶对称分量来恢复。如果引入一个辅助因子g(n)，上式可写作 其中 最小相位信号法求复倒谱原理框图如下 递归法 这种方法仅限于是最小相位信号的情况。根据z变换的微分特性得 对上式求逆z变换，根据z变换的微分特性，有 所以 设x(n)是最小相位序列，而最小相位信号序列一定为因果序列 ，所以有 由于 及 可得递推公式 递归运算后由复倒谱定义 可知 如果x(n)是最大相位序列，则变为 其中 2、基音检测； 语音的倒谱是将语音的短时谱取对数后再进行IDFT得到的，所以浊音信号的周期性激励反映在倒谱上是同样周期的冲激。借此，可从倒谱波形中估计出基音周期。一般把倒谱波形中第二个冲激，认为是对应激励源的基频。下面给出一种倒谱法求基音周期的框图及流程图如下 3.共振峰检测 倒谱将基音谐波和声道的频谱包络分离开来。对倒谱进行低时窗选，通过语音倒谱分析系统的最后一级，进行DFT后的输出即为平滑后的对数模函数，这个平滑的对数谱显示了特定输入语音段的谐振结构，即谱的峰值基本上对应于共振蜂频率，对平滑过的对数谱中的峰值进行定位，即可估计共振峰。原理框图及流程图如下。 四、实验步骤（包括主要步骤、代码分析等） 1．倒谱MATLAB实现代码段 clear all; %倒谱 [s,fs,nbit]=wavread(beijing.wav); %读入一段语音 b=s; %将s转置 x=b(5000