语音信号的倒谱分析

语音信号的倒谱分析根据语音信号的产生模型，语音信号S(Z)是一个线性非移变因果稳定系统V(Z)受到信号E(Z)激励后所产生的输出。在时域中，语音信号s(n)是该系统的单位取样响应v(n)和激励信号e(n)的卷积。在语音信号数字处理所涉及的各个领域中，根据s(n)来求得v(n)和e(n)具有非常重要的意义。由卷积信号求得参与卷积的各个信号的过程称为解卷过程。语音信号的倒谱分析解卷算法可以分为两大类：第一类是首先为线性系统V(Z)建立一个模型，然后对模型参数按照某种最佳准则进行估计，这种方法称为参数解卷方法。采用的模型可以分为全极点模型（AR模型）和零极点模型（ARMA模型），如果采用最小均方误差准则对AR模型进行估计，就得到线性预测编码算法(LPC)。第二类算法称为非模型解卷。同态信号处理完成解卷任务就是其中最重要的一种。语音信号的倒谱分析对信号进行分析得出它的倒谱参数的过程称为同态处理。对语音信号的某一帧同样可以分析出它的短时倒谱参数，总的说来，无论对于语音通信、语音合成或语音识别，倒谱参数所含的信息比其他参数多，也就是语音质量好，识别正确率高。但其缺点是运算量比其他参数大，尽管如此，倒谱分析方法仍不失为一种有效的语音信号的分析方法。同态分析的基本原理有很多客观物理现象中的信号，其中各组成分量的组合，并不是按照加法组合原则组合起来的，如图像信号、地震信号、调制信号、语音信号等，它们都不是加性信号，而是乘积性或卷积性组合的信号。显然，这时不能用线性系统来处理，而必须用满足该组合规则的非线性系统来处理。但是非线性系统地分析非常困难。同态信号处理法就是设法将非线性问题转化为线性问题来处理的一种方法。按照被处理的信号来分类，大体上可以分为乘积同态信号处理和卷积同态信号处理。由于语音信号可以视为声门激励信号和声道响应信号的卷积结果。我们仅讨论卷积同态信号处理系统的问题。卷积同态信号处理系统同态系统可以分解为两个特征系统（即特征系统和逆特征系统）（指取决于信号的组合规则）和一个线性系统（仅取决于处理要求）卷积同态信号处理系统卷积同态信号处理系统由于加性信号的Z变换结果仍为加性信号，所以倒谱这种时域信号，是可以用线性系统来处理的，经线性处理之后，如欲在恢复出语音信号，则可以采用逆特征系统来实现，即特征系统的逆运算。即将线性系统输出的加性倒谱信号：卷积同态信号处理系统特征系统与逆特征系统的组成语音信号的倒谱利叶变换来代替。变换都可以利用正负福都包含单位圆，正反的收敛域问题中，绝大多数数字信号处理运算。而逆特征系统则为其逆谱域中加运算，的卷积运算转换为复倒特征系统将离散时域中域。的离散时域称为复倒谱的复倒谱。复倒谱所处是同样。的复倒谱，其英文为是ZZYZYZXZXnynystrumComplexCepnxnx)(ˆ),(),(ˆ),()()(ˆ)()(ˆdwjwnjwXjwXFnxjwXjwXjwnnxnxFjwXNNn)exp()(expˆ21)](expˆ[)(ˆ)](expln[)(expˆ)exp()()]([)(exp121征系统求得复倒谱的另一个特语音信号的倒谱。谱称为倒频谱，简称为倒求得倒谱的特征系统CepstrumncdwjwnjwCjwCFncjwXjwCjwnnxnxFjwXNNn)()exp()(exp21)](exp[)(])(expln[)(exp)exp()()]([)(exp121复倒谱经过正逆两个特征系统变换后，序列可以还原为本身。但是倒谱经过正逆两个特征系统变换后，序列不可以还原为本身。由序列的复倒谱求倒谱的方法由于偶对称序列的DTFT是实函数，奇对称序列的DTFT是虚函数。)(ˆ)(ˆ)(ˆnxnxnxoe的形式：之和和一个奇对称表示为一个偶对称序列首先将。它的倒谱求出，那么可以由的复倒谱如果已知一个实序列)(ˆ)(ˆ)(ˆ)()(ˆ)(ˆ)(nxnxnxncnxnxnxoe由序列的复倒谱求倒谱的方法)](expIm[)](expRe[)](exp[)(expln)](expln[)(expˆjwXjjwXjwXjArgjwXjwXjwX因此)](ˆ)(ˆ[21)()(nxnxnxnce相位倒谱的概念假设则称p(n)为相位倒谱。)]](exp[[)(1jwXArgFnp)](ˆ)(ˆ[21)(ˆ)(nxnxnxnpo。则包含两个方面的特征而的相位函数的特征，表现的是的模函数的特征，表现的是不难看出，)(ˆ)(exp)()(exp)(cnxjwXnpjwXn已知倒谱求复倒谱的方法要想由倒谱求复倒谱，首先复倒谱必须满足一定的条件，比如是因果序列则因此)()(ˆ)(ˆnunxnx0)(ˆ210)(ˆ0)(ˆ21)](ˆ)(ˆ[21)(ˆ)(nnxnnxnnxnxnxnxnce000)(0)(2)(ˆnnncnncnx已知倒谱求复倒谱的方法如果复倒谱是一个反因果序列：则可以推导出：只有当x(n)是一个因果最小相位序列是其复倒谱序列才是一个因果稳定序列。这要求x(n)应满足两个条件：1x(n)=x(n)u(n)；2X(Z)=Z[x(n)]的零极点都应该在单位圆内。)()(ˆ)(ˆnunxnx0)(20)(00)(ˆnncnncnnx语音信号倒谱和复倒谱的性质根据语音信号产生的模型，在z域中语音信号S(Z)等于激励信号E(Z)和声道传输函数V(Z)的乘积，即S(Z)=E(Z)V(Z)。经过同态系统后可以得到：先讨论声门激励信号。除了人们发清音时，声门激励是能量较小、频谱均匀分布的白噪声之外；发浊音时，声门激励是以基调周期为周期的周期脉冲序列)(ˆ)(ˆ)(ˆnvnens语音信号倒谱和复倒谱的性质语音信号倒谱和复倒谱的性质由上式可以得出以下结论：一个周期冲激的有限长度序列，其复倒谱也是一个同周期长度的周期冲激序列，只是其长度变为无限长度、振幅随着K值的增加而衰减，衰减速度比原来序列要快，显然，周期冲激序列的倒谱的这些性质对于语音信号的分析是很有用的，这意味着除了原点之外，可以用“高时窗”来从语音信号的倒谱中提取浊音激励信号的倒谱，从而使倒谱法提取音调成为现实。取其他值非零值nNNNnneppp,0,3,2,,)(ˆ语音信号倒谱和复倒谱的性质语音信号倒谱和复倒谱的性质语音信号倒谱和复倒谱的性质语音信号倒谱和复倒谱的性质语音信号倒谱和复倒谱的性质语音信号倒谱和复倒谱的性质在清音情况下，e(n)具有噪声特性，因而其复倒谱也没有明显的峰起点，且分布范围很宽，从低时域延伸到高时域。而v(n)的复倒谱仍然只分布在低时域中。求复倒谱的一种有效的递推算法前提：x(n)是最小相位序列。因为10)()(NnnZnxZX1011)]([)()()()('NnnnnxZZZnxnZdZZdXZX))(ln()(ˆZXZX)()('ˆ)('ZXZXZX求复倒谱的一种有效的递推算法)]([))(ˆ()]([nxZnxnZnnxZ)()(ˆ)()}(ˆ{))((lnxlxlnxnxnnxnl)0()1()1(ˆ1xxxn时，)1()0()1(ˆ21)0()2()2(ˆ2xxxxxxn时，)0()1()2(ˆ32)0()2()1(ˆ31)0()3()3(ˆ3xxxxxxxxxx时，)()0()(ˆ)()0()()(ˆ)0()(ˆ)()(ˆ)()(1010lnxxlxnlxnxnxxnxlnxlxnlnxnlnl可推导出：).0(ln)0(ˆ:)0(,ln)0(ˆxxAxAx所以，可以得出，对于因果序列而言，由前面的推导可知，语音信号的线性预测分析LinearPrediction1947年维纳提出；1967年板仓等人应用于语音分析与合成；语音信号处理与分析的核心技术提供了预测功能；提供了声道模型和声道模型的参数估计方法；基本思想：语音样本之间存在相关性，一个语音信号的样本可以用过去若干个样本的线性组合来逼近；Piiinsans1)()(ˆPiiPiiinsainsansnsnsne01)()()()(ˆ)()(线性预测原理线性预测是目前分析语音信号的最有效的方法之一，分析的结果是得到一组信号的全极点模型参数，所以又称为信号参数模型法。这个方法的基本思想是将被分析信号模型化，即用有限数目的模型参数来描述信号中的信息，具体来说，将被分析信号s(n)视为某系统（即模型）的输出，而系统的输入，在s(n)为确定性信号是采用单位取样序列。在s(n)为随机信号是采用白噪声序列。语音信号的模型常用来产生合成语音，所以称为合成滤波器PiiizaGzUzSzH11)()()()1}({Piai求解滤波器参数和G的过程就是线性预测的分析过程。线性预测原理在基于参数模型的谱估计方法和系统辨识中，常常假定系统的传递函数是有理函数，也就是变量Z的有理分式，这种有理分式有三种情况：只有零点没有极点的情况，称为滑动平均模型。即MA模型只有极点没有零点的情况，称为自回归模型。即AR模型既有零点又有极点的情况，称为自回归滑动平均模型。即ARMA模型线性预测原理全极点模型的参数估计十分简单，只需很小的几个极点就可以相当好的估计一种频谱或一种系统的频率响应，因此传递函数相当于一个递归数字滤波器。即IIR滤波器线性预测法正是基于全极点模型的假定，采用时域均方最小误差准则来估计模型参数的。参数解卷的通用模型)(/)()(ZAZGZV)(ne)(ZE)(ns)(ZS（未知）（未知）（已知）解卷的的问题。，从而全部解决便可以求得和根据便可求得。由于，那么未知的系统正确的估计出这些参数可能根据已知的。如果能有一种算法，都是实数，且和来表示。且域传输函数可以用其方程描述，那么其特性且可以用一个线性差分移变因果稳定系统，未知系统是一个线性非产生的。如果假设这个激励一个未知的系统是一个未知的序列假设一个已知序列E(Z))()(),()()(V(Z))(1)(,)(),(/)()()(Z)(e(n)s(n)000ZSZVZSZVZEnsaagZaZAZgZGZAZGZVZVnvijPiiiQjjj对模型的限制为了得到一种高效的求解方法。令G(Z)=1，模型中只含有极点不含有零点。这种模型称为全极点模型；对未知序列e(n)加以限制，表示成Ge(n)的形式，其中e(n)是一个周期脉冲序列或者高斯白噪声序列；系数G是非负实数，用来控制输出序列的幅度。斯白噪声序列。是周期脉冲序列还是高出并且判断和的基础上，求得其次，在估计得到和它的各个系数的阶数，也就是估计模型为首先正确估计以将参数解卷问题归结做了这些限制之后，可为周期。。。而对于其他值，时，其中，当，可以表示为：是一个周期脉冲序列时当)(),(GV(Z);~PA(Z))(/1)(0)(1)(0n)()(e(n)1neneaaZAZVNnnrNnnepprp由全极点模型构成的参数解卷模式)(/1)(ZAZV)(ne)(nsG1Z1Z1Z1Z1a2a1papa周期序列或高斯白噪声序列虽然采取了一些限制，但是这种全极点模型的适应面还是很宽的。特别适合于语音信号处理领域利用线性预测方法估计全极点模型的参数的命题现在讨论如何用线性预测的方法来解决全极点模型的参数估计问题，从而完成解卷的任务。参数解卷的命题是：如果得到语音序列s(n)，并且知道它是由一个周期脉冲序列或高斯白噪声序列激励一个全极点模型产生的。那么，采用什么样的算法，才能在最佳意义上估计出模型参数，激励信号e(n)和增益系数G。线性预测算法能按照最小均方误差的准则给出问题的解。利用线性预测的方法解决全极点模型的