lpc2lpcc.m

来自「语音信号端点检测方面的源程序」· M 代码 · 共 22 行

%求解倒谱程序
function lpcc=lpc2lpcc(lpc1)
%  lpc1 实际为预测系数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          % 这里的循环是分步计算n＝2、3、4......的时候，计算h（n）的值，
lpcc(n)=lpcc(n)+lpc1(l)*lpcc(n-l)*(n-l)/n;  %需要循环完成后才能得到一个完整的lpcc(n)； 
end
end

%%%%%%%%%%%%%%%%%%%%%%%%%
for n=n_lpc+1:n_lpcc  % n>p的时候，h(n)＝sigma(...)的部分，也就是（<n）时的后面部分的表达式。
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;