第六篇 章 语音信号线性预测分析 数字语音处理及MATLAB仿真 .pptVIP

LSP特性实验仿真结果 1.多项式根分布图：16阶LPC系数构成的17阶对称和反对称多项式的根在单位圆上的分布图如下 LSP特性实验仿真结果 2. LSP轨迹图：连续20帧16阶LPC系数对应的LSP轨迹图如下 LSP特性实验仿真结果 3. LSF轨迹图：连续20帧16阶LPC系数对应的LSF轨迹图如下 LSP特性实验仿真结果 4. 一帧语音信号16阶LPC谱包络和相应LSF轨迹图： 6.6.2 LPC参数到LSP参数的转换及MATLAB实现 将P(z)和Q(z)中与LSP系数无关的两个实根取 掉，得到如下两个新的多项式 求解上两式等于零时 值，即得LSP系数 第一种方法：利用代数方程式求解。 求解LSP系数的几种方法 由 令 根据 得 第二种方法：离散傅立叶变换（DFT）方法 对 和 系数求离散傅立叶变换，得到 各点的值，搜索最小值的位置，即是零点所在。由于除了0和?之外，总共有p个零点，而且 和 的根是相互交替出现的，因此只要很少的计算量即可解得，其中N的取值取64～128就可以。 第三种方法：利用切比雪夫(Chebyshev)多项式求解 则 可以写作： 其中 是m阶的Chebyshev多项式 f(i)是由递推关系计算得到的 和 的每个系数。 用下面的递推关系 其中 。 多项式C(x)在x=cos?时的递推关系是： 其中初始值 第四种方法： 将0～?之间均分为60个点，以这60个点的频率值代入(6-41)、(6-42)式，检查它们的符号变化，在符号变化的两点之间均分为4份，再将这三个点频率值代入方程(6-41)、(6-42)，符号变化的点即为所求的解。这种方法误差略大，计算量较大，但程序实现容易。 教材给出从LPC参数到LSP参数转换的MATLAB程序，其中a_lsf_conversion.m为求解LSF的函数，a_lsf_main.m为主程序。由于MATLAB程序本身有求多项式根的函数，因此在求解 和 零点时直接调用即可，这极大简化了求解过程。 6.6.3 LSP参数到LPC参数的转换及MATLAB实现 已知量化和内插的LSP系数 ，可用下式计算 和 的系数 和 ： 以下的递推关系可利用qi，i=0、1…、p-1，来计算 其中 初始值 ， 把上面递推关系中的q2i-1替换为q2i，就可以得到 ， 根据 、 和 和 利用 得到LPC系数为 6.7 导抗谱对ISP分析 导抗谱对ISP（Immittance Spectral Pair）也是与LPC系数等价的一种参数，可提高LPC系数鲁棒性，由Yuval Bistritz 和Shlomo Peller在1993年提出。目前已经用于自适应多速率宽带（AMR-WB）语音编码算法中。 6.7.1 ISP的定义和特点 设线性预测逆滤波器为 用LPC系数构造的对称和反对称多项式如下 反映声门激励的导抗函数如下 ISP包括Ip(z)的极点和零点，以及一个反射系数 由于导抗函数的所有系数都是实数，因此其分子和 分母多项式的根将以共轭复数的形式出现，且分子 和分母多项式所有的根均位于单位圆上而且彼此轮 流出现。此处针对p为偶数进行讨论。这样 和 在单位圆上分别有p/2和p/2-1个共轭复根， 因此可得到如下的多项式 其中的cosωi，i=1，2，? p -1 ，是ISP前p-1个系数 在余弦域的表示，式中的kp是ISP的最后一个系数， 也称为常数增益。 ISP的p个参数如下： LPC滤波器稳定性得到保证的条件： 1.与前p-1个ISP系数相对应的频率按升序排列，即 2.常数增益满足 ，在AMR-WB中，取kp=ap 导抗谱频率ISF：与p个ISP系数相对应的频率，即 前p-1个ISP表现出与p个LSP相似的一些特性： 1．都在单位圆上且满足降序排列的特性。 2．与ISP对应的前p-1个ISF都满足升序排列特性，且ISP的第p个系数小于1，这使得与之对应的LPC滤波器的稳定性可以得到保证。因此ISP分析就是用p-1个离散频率和离散