lpc_shorttime_filte.m

该代码实现了am

function lpc_shorttime_filte

for i=1:8
LARpp(1,i)=(LARc(1,i)-B(1,i))./A(1,i)
end
LARpppre=LARpp
LARppp(1,:)=0.75*LARpppre+0.25*LARpp(1,:)
LARppp(2,:)=0.5*LARpppre+0.5*LARpp(1,:)
LARppp(3,:)=0.25*LARpppre+0.75*LARpp(1,:)
LARppp(4,:)=LARpp(1,:)


for k=1:4
    for i=1:length(LARppp)
	if 	abs(LARppp(k,i))<0.675
		rp(k,i)=LARppp(k,i);
	elseif 	abs(LARppp(k,i))>=0.675 & abs(LARppp(k,i))<1.225
		rp(k,i)=sign(LARppp(k,i))*(0.500*abs(LARppp(k,i))+0.337500);
	elseif	abs(LARppp(k,i))>=1.225 & abs(LARppp(k,i))<=1.625
		rp(k,i)=sign(LARppp(k,i))*(0.125*abs(LARppp(k,i))+0.796875);
	end
end
end

[d(1:13),tt]=latcfilt([1,rp(1,:)],S(1:13))
[d(14:27),tt]=latcfilt([1,rp(2,:)],S(14:27))
[d(28:40),tt]=latcfilt([1,rp(3,:)],S(28:40))
[d(41:160),tt]=latcfilt([1,rp(4,:)],S(41:160))