语音合成与识别

开发平台：

Matlab

lsp.m：源码内容

							function [P_w,Q_w,H]=lsp(s,p,num)
%此函数的功能是求语音信号的线谱对参数
              
a=lpc_coefficients(s,p);              %求线性预测系数
a=a';                                 %将a转化为行向量
p=length(a); 
P_c=[1,-a-a(p:-1:1),1];               %多项式P(Z)的系数
Q_c=[1,-a+a(p:-1:1),-1];              %多项式Q(Z)的系数
P_root=roots(P_c);                    %多项式P(Z)=0的根
Q_root=roots(Q_c);                    %多项式Q(Z)=0的根
P_r_angle=sort(angle(P_root));        %多项式P(Z)=0的根的辐角,并按升序排列
Q_r_angle=sort(angle(Q_root));        %多项式Q(Z)=0的根的辐角,并按升序排列
P_w=P_r_angle(p/2+1:p);               
Q_w=Q_r_angle(p/2+2:p+1);             %为线谱对系数
for k=1:num
    H(k)=2^((1-p)/2)/(sin(pi/num*k).^2*prod((cos(2*pi/num*k*ones(p/2,1))-cos(Q_w)).^2)+cos(pi/num*k).^2*prod((cos(2*pi/num*k*ones(p/2,1))-cos(P_w)).^2));
end   

						