frefine.m

来自「有关matlab的电子书籍有一定的帮助希望有用」· M 代码 · 共 20 行

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.5 $
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;