lpc_cepstrum.m

关于双门限法的语音端点检测

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