ltermprediction.m

wav声音文件的通过多脉冲激励编码器

function [v,bvalue]=LTermPrediction(errorpre,error,fr,frmsize,pitch_marks)

v=zeros(1,frmsize);%long term prediction error
%caculate b* in long term prediction
bnumerator=0;
bdenominator=0;
for i=1:frmsize
    if (i-pitch_marks)<1
        bnumerator =bnumerator+errorpre(i)*error(frmsize+(i-pitch_marks));
        bdenominator=bdenominator+errorpre(frmsize+(i-pitch_marks))^2;
    else
        bnumerator =bnumerator+error(i)*error(i-pitch_marks);
        bdenominator=bdenominator+error(i-pitch_marks)^2;
    end
end
bvalue=bnumerator/bdenominator;
for i=1:frmsize
    if (i-pitch_marks)<1 
        v(i)=error(i)-bvalue*error(frmsize+(i-pitch_marks));
    else
        v(i)=error(i)-bvalue*error(i-pitch_marks);
    end
end