Example 1.  Use the Nelder-Mead method to find the minimum of   [Graphics:Images/NelderMeadMod_gr_94.gif].  

Solution 1.

Enter the function  [Graphics:../Images/NelderMeadMod_gr_95.gif] and graph the surface [Graphics:../Images/NelderMeadMod_gr_96.gif].

[Graphics:../Images/NelderMeadMod_gr_97.gif]


[Graphics:../Images/NelderMeadMod_gr_98.gif]

[Graphics:../Images/NelderMeadMod_gr_99.gif]



[Graphics:../Images/NelderMeadMod_gr_100.gif]


[Graphics:../Images/NelderMeadMod_gr_101.gif]

[Graphics:../Images/NelderMeadMod_gr_102.gif]



[Graphics:../Images/NelderMeadMod_gr_103.gif]


[Graphics:../Images/NelderMeadMod_gr_104.gif]

[Graphics:../Images/NelderMeadMod_gr_105.gif]

Start with the three vertices  

[Graphics:../Images/NelderMeadMod_gr_106.gif]

The function [Graphics:../Images/NelderMeadMod_gr_107.gif] takes on the values  

[Graphics:../Images/NelderMeadMod_gr_108.gif]


[Graphics:../Images/NelderMeadMod_gr_109.gif]

The function values must be compared to determine [Graphics:../Images/NelderMeadMod_gr_110.gif], [Graphics:../Images/NelderMeadMod_gr_111.gif] and [Graphics:../Images/NelderMeadMod_gr_112.gif];  

 

[Graphics:../Images/NelderMeadMod_gr_113.gif]

The vertex  [Graphics:../Images/NelderMeadMod_gr_114.gif]  will be replaced.  The points [Graphics:../Images/NelderMeadMod_gr_115.gif] and [Graphics:../Images/NelderMeadMod_gr_116.gif] are constructed.  

[Graphics:../Images/NelderMeadMod_gr_117.gif]



[Graphics:../Images/NelderMeadMod_gr_118.gif]

The function value  [Graphics:../Images/NelderMeadMod_gr_119.gif]  is less than  [Graphics:../Images/NelderMeadMod_gr_120.gif],  so the situation is Case (i).  
Since  
[Graphics:../Images/NelderMeadMod_gr_121.gif],  we have moved in the right direction, and the vertex [Graphics:../Images/NelderMeadMod_gr_122.gif] must be constructed:

[Graphics:../Images/NelderMeadMod_gr_123.gif]



[Graphics:../Images/NelderMeadMod_gr_124.gif]

The function value  [Graphics:../Images/NelderMeadMod_gr_125.gif]  is less than  [Graphics:../Images/NelderMeadMod_gr_126.gif],  and we replace [Graphics:../Images/NelderMeadMod_gr_127.gif] with [Graphics:../Images/NelderMeadMod_gr_128.gif] and the new triangle has vertices

 

[Graphics:../Images/NelderMeadMod_gr_129.gif]

The process continues and generates a sequence of triangles that converges down on the solution point  [Graphics:../Images/NelderMeadMod_gr_130.gif].  The computer implementation of the Nelder-Mead algorithm is given below

[Graphics:../Images/NelderMeadMod_gr_131.gif]


[Graphics:../Images/NelderMeadMod_gr_132.gif]


[Graphics:../Images/NelderMeadMod_gr_133.gif]


[Graphics:../Images/NelderMeadMod_gr_134.gif]
[Graphics:../Images/NelderMeadMod_gr_135.gif]

Let us compare this answer with Mathematica's built in procedure FindMinimum.

[Graphics:../Images/NelderMeadMod_gr_136.gif]


[Graphics:../Images/NelderMeadMod_gr_137.gif]
[Graphics:../Images/NelderMeadMod_gr_138.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004