%求解倒譜程序
function lpcc=lpc2lpcc(lpc1)
%  lpc1 實際為預(yù)測系數(shù)Ai;
n_lpc=8;n_lpcc=12;
lpcc=zeros(1,n_lpcc);
lpcc(1)=lpc1(1);   %  h(1)=a(1);  
%%%%%%%%%%
for n=2:n_lpc      % h(n),1<=n<=p; h(n)=a(n)+ sigma(.........)
lpcc(n)=lpc1(n);   
for l=n-1          % 這里的循環(huán)是分步計算n＝2、3、4......的時候，計算h（n）的值，
lpcc(n)=lpcc(n)+lpc1(l)*lpcc(n-l)*(n-l)/n;  %需要循環(huán)完成后才能得到一個完整的lpcc(n)； 
end
end

%%%%%%%%%%%%%%%%%%%%%%%%%
for n=n_lpc+1:n_lpcc  % n>p的時候，h(n)＝sigma(...)的部分，也就是（<n）時的后面部分的表達(dá)式。
lpcc(n)=0;
for l=1:n_lpc         % 見課本p100頁。
lpcc(n)=lpcc(n)+lpc1(l)*lpcc(n-l)*(n-l)/n;
end
end
lpcc=-lpcc;