Gradient Search

dx6.m

achromat_setup;
load achromat;
show_contour3;
v = [0 3]';
p(1,:)=v';
for k=2:6
    [A fz] = calculate_derivatives(@achromat2,v);
    g = A'*fz/norm(A'*A);
    v = v -g;
    p(k,:)=v';
end
hold on;
plot(p(:,1),p(:,2),'ko-','LineWidth',1.5);
hold off;

Output

Search in direction opposite to gradient vector.

Maintained by John Loomis, last updated 6 Feb 2003