formant.m

matlab code for EMT field theory project stating mashgrid quiver

function [P, F, I] = formant(y)
% calculate the PSD using Yule-Walker's method
order = 12;
P = pyulear(y,order,[]);
P = 10*log10(P); % convert to DB
F = 0:1/128:1; 
[Pm,I] = pickmax(P);
I = I/128; % normalize the index