English
- 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
