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