function h_c=lpc_cepstrum(s,p)
%本函數的功能是求lpc倒譜h_c,它為聲道沖激響應h(n)的復倒譜
%p為選定模型的階數
n=length(s);
a=lpc_coefficients(s,p);                                                    %獲得線性預測的系數

h_c(1)=a(1);
%下面為i<=p的情況
for i=2:p                                                        
    k=1:i-1;
    h_c(i)=a(i)+sum(a(1:i-1)'.*(ones(1,i-1)-k/i).*h_c(i-1:-1:1)); 
end
%下面為i>=p的情況
for i=p+1:n
    k=1:p;
    h_c(i)=sum(a(1:p)'.*(ones(1,p)-k/i).*h_c(i-1:-1:i-p)); 
end