gusucode.com > 信号处理工具箱 - signal源码程序 > signal\signal\signal\private\frefine.m
function rs = frefine(a,v,rs) % f = frefine(a,v,rs); % refine local minima and maxima of H using Newton's method % % H : H = cos(w*v)*a; % rs : initial values for the extrema of H % % Copyright (c) 1988-98 by The MathWorks, Inc. % $Revision: 1.1 $ a = a(:); v = v(:)'; w = rs(:); m = length(a)-1; for k = 1:5 H1 = -sin(w*v) * (v'.*a); H2 = -cos(w*v) * ((v.^2)'.*a); w = w - H1./H2; end rs(:) = w;