mgif_es.m

来自「这是一个用于语音信号处理的工具箱」· M 代码 · 共 28 行

%
% To estimate first five formants frequencies and bandwidths
% using lpc analysis
% This estimation is rough initial values for manual adjustment of
% formant frequency and bandwidth.
%
% Author : Minkyu Lee
% Date 10-Oct-1994
%Modified by KArthik on May 31 1997

ydata_h=hamming(length(ydata_m)).*ydata_m;
lpc_cof=lpc(ydata_h,12);
rt=roots(lpc_cof);
rt1=[];

for i=1:12,
	if angle(rt(i)) > 0
		rt1=[ rt1 ; rt(i) ];
	end
end

[frt fi]=sort(angle(rt1));

for i=1:5
	ff(i)=round( frt(i)/(2*pi*Ts) );
	fb(i)=round( -2*log(abs(rt1(fi(i))))/(Ts*2*pi) );
end