function [P_w,Q_w,H]=lsp(s,p,num)
%此函數(shù)的功能是求語音信號的線譜對參數(shù)
              
a=lpc_coefficients(s,p);              %求線性預(yù)測系數(shù)
a=a';                                 %將a轉(zhuǎn)化為行向量
p=length(a); 

P_c=[1,-a-a(p:-1:1),1];               %多項式P(Z)的系數(shù)
Q_c=[1,-a+a(p:-1:1),-1];              %多項式Q(Z)的系數(shù)

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);             %為線譜對系數(shù)

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