【信号处理】基于matlab实现语音信号变速

1 简介

一种语音变速的方法包含下列步骤:接收原始语音信号;计算出该原始语音信号的基本周期;依据该基本周期定义搜寻范围;在该原始语音信号的每个该搜寻范围内寻找最大值;依据该最大值将该原始语音信号划分为多个语音区段;依据变速指令对该原始语音信号的每个该语音区段进行变速运算,而取得变速语音信号;及输出该变速语音信号.上述方法的好处是,将语音信号以最细小的基本周期单位进行后续的变速处理.

2 部分代码

%LPC转LSP函数function lsf=lpctolsf(a)a=a(:);                       % 将a转换为列向量% 如果a不是实数，输出错误信息：LSF不适用于复多项式的求解if ~isreal(a) ,      error('Line spectral frequencies are not defined for complex polynomials.');end% 如果a(1)不为1,将矩阵a的每个元素除以a(1)再赋给矩阵if a(1) ~= 1.0,    a=a./a(1); end% 如果a的根不在单位圆内,显示错误信息并返回if (max(abs(roots(a))) >= 1.0),    error ('The polynomial must have all roots inside of the unit circle. ');    return;end%求对称和反对称多项式的系数p=length(a) - 1;              % 求对称和反对称多项式的阶次a1=[a;0];                     % 给行矩阵a再增加一个元素0的行 a2=a1(end:-1:1);              % a2的第一行为a1的最后一行，最后一行为al的第一行P1=a1+a2;                     % 按式(3-82)求对称多项式的系数Q1=a1-a2;                     % 按式(3-83)求反对称多项式的系数%如果阶次p为偶数次，按式(3-89)和式(3-90)从P1去掉实数根z=-1，从Q1去掉实数根z=1%如果阶次为奇数次，从Q1去掉实数根z=1及z=-1if rem(p,2),                  % 求解P除以2的余数，如果P为奇数次。余数为1，否则为0   Q=deconv(Q1,[1 0 -1]);     % 奇数阶次，从Q1去掉实数根z=1及z=-1   P=P1;else                          % P为偶数阶次执行下面的操作   Q=deconv(Q1,[1 -1]);       % 从Q1去掉实数根z=1   P=deconv(P1, [1 1] );      % 从P1去掉实数根z=-1endrP=roots(P);                  % 求去掉实根后的多项式P的根rQ=roots(Q) ;                 % 求去掉实根后的多项式Q的根aP=angle(rP(1: 2:end) ) ;     % 将多项式P的根转换为角度(为归一化角频率)赋给apaQ=angle(rQ(1: 2:end));       % 将多项式Q的根转换为角度(为归一化角频率)赋给aQ lsf= sort([aP; aQ]);          % 将P、Q的根(归一化角频率)按从小到大的顺序排序后即为lsf​​

3 仿真结果

4 参考文献

[1]张昊. 语音数字信号增强处理及其Matlab实现[J]. 通信技术, 2009(1):3.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。

​