Matlab: Medición de frecuencia

A partier de la inquietud planteada por el Lucho de medir la frecuencia de sus mediciones del osciloscopio salió el siguiente código:

function f = findfreq(t,u,fmin,fmax,varargin)
if nargin>4
if strcmpi(varargin{1},'fpoints')
frecs = linspace(fmin,fmax,varargin{2})*2*pi; % valores de frecuencia a probar
elseif strcmpi(varargin{1},'fstep')
frecs = [fmin:varargin{2}:fmax]*2*pi; % valores de frecuencia a probar
end
else
frecs = [fmin:fmax]*2*pi; % valores de frecuencia a probar
end
m = zeros(size(frecs,1),2);
for j = 1:length(frecs)
m(j,1) = sum((u.*sin(frecs(j)*t)).^2);
m(j,2) = sum((u.*cos(frecs(j)*t)).^2);
end
[a,b] = max(sum(m.^2,2),[],1);
f = frecs(b)/(2*pi);
plot(frecs,m)

ojalá le sirva a alguien, si necesitan alguna explicación del código, ¡¡¡pregunten poh!!!

Funcionó bastante bien con respecto a la medición de frecuencia del osciloscopio y las mediciones al ojímetro, asi que bien...

Lo que sí me entró la duda si funciona bien cuando las ondas no estan en fase con el seno o coseno (eso se trata de superar comparando con seno y coseno), si alguien encuentra que falla, que avise ;)

Bueno, eso sería por esta jornada Matlabiana