Lpcc.rar_LPCC中matlab实现

clear all; close all; [x,fs,bit]=wavread('f:\voice\a.wav'); framelength=512; H=hamming(framelength); P=14; m=20; for n=1:framelength; y(n)=x((m-1)*framelength+n); end y=y.*H'; R=autocorrelate(y,framelength); E(1)=R(1); K(1)=R(2)/R(1); a(1,1)=K(1); E(2)=(1-K(1)^2)*E(1); m=2; while(m<=P) K(m)=parcor(E,R,a,m); a(m,m)=K(m); for n=1:m-1; a(m,n)=a(m-1,n)-K(m)*a(m-1,m-n); end E(m+1)=(1-K(m)^2)*E(m); m=m+1; end for m=1:P; A(m)=a(P,m); end G=1; for m=1:framelength; sum=0; for n=1:P; sum=sum+A(n)*exp(-j*2*pi*n*m/framelength); end H1(1,m)=G/(1-sum); end plot(f,20*log10(abs(H1(1:framelength/2))),'r');