melp_lsf2lpc.m

语音编码中的MELP的经典估计

%LSF TO LPC
%INPUT: 
%        lsfs(LSF parameters)
%OUTPUT:
%        f(LPC parameters)
function f=melp_lsf2lpc(lsfs)
lsfs=pi*lsfs/4000;
for i=1:5
   for j=1:2
      w(j,i)=lsfs((i-1)*2+j);
   end
end
w=cos(w);
for i=1:2
   temp=[1 -2*w(i,5) 1];
   P=temp;
   for j=1:4
      temp=[1 -2*w(i,j) 1];
      w0=[P 0 0;0 P 0;0 0 P];
      P=temp*w0; 
   end
   Q(i,:)=P;
end
f=([Q(1,:) 0]+[0 Q(1,:)]+[Q(2,:) 0]-[0 Q(2,:)])/2;
f=f(2:11);