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