计算lpcc的matlab程序

该matlab程序用于计算语音信号的lpcc参数

function[y]=lpcc（input）
p=12;
NJ=256;
for m=0:1:p
    r（m+1）=0;
    for i=1:1:（NJ-m）
        r（m+1）=r（m+1）+input（i）*input（i+m）;
    end
end
R=zeros（1p）;
R（1（1:p））=r（1（2:p+1））;
R=R‘;
E=zeros（1p）;
k=E;
a=zeros（pp）;
e=r（1）;          %Durbin算法初始化
k（1）=r（2）/r（1）;
a（11）=k（1）;
E（1）=（1-k（1）^2）*e;
l=2;
for i=2:1:12;    %算法开始
    m_t=sum（a（（i-1:-1:1）i-1）.*R（（1:1:i-1）1））;
    k（i）=（R（i）-m_t）/E（i-1）;
    a（ii）=k（i）;
    a（（1:1:i-1）l）=a（（1:1:i-1）i-1）-k（i）*a（（i-1:-1:1）i-1）;
    E（i）=（1-k（i）^2）*E（i-1）;
    l=l+1;
end
y=a（:12）‘;            %预测系数
A=zeros（112）;         %计算LPCC参数，使用12阶
A=y（1end:-1:1）;
c=zeros（112）;
c（1）=A（1）;
for n=2:1:12
    c（n）=A（n）+sum（（1:1:n-1）/12.*c（1:1:n-1）.*A（n-（1:1:n-1）））;
end      
y=c;            

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        788  2007-07-04 19:32  lpcc.m

----------- ---------  ---------- -----  ----

                  788                    1