view dsp/formants.m @ 36:9e7be347b3a0

Renamed sequence classes to avoid clashes with seq methods;
Fixed default slicing dimension while retaining behaviour of window.m;
Updated use of sequences in dsp/synth.

% formants - Compute formants from LPC coefficients
function F=formants(a,M,fs)

const=fs/(2*pi);
if all(isfinite(a))
	rts=roots(a);

	fx = const*angle(rts);
	bw = -2*const*log(abs(rts));

	F = log(sort(fx(fx>90 & fx<4000 & bw<700)));
	N = length(F);
	if N<M, F=[F;nan(M-N,1)];
	else F=F(1:M); end
else
	F=nan(M,1);
end
% only first 3?