No BSD License  

Highlights from
simps

5.0

5.0 | 3 ratings Rate this file 15 Downloads (last 30 days) File Size: 9.37 KB File ID: #102

simps

by

 

The method is based on iterations of full-dimensional simplex calls

| Watch this File

File Information
Description

SIMPS (StrategySimplex)-Constrained minimizer

The method is based on iterations of full-dimensional simplex calls (Nelder-Mead direct search method), each time followed by a series of two-dimensional simplex calls (local improvements by subspaces).

The inner full- and two-dimensional Nelder-Mead minimizers are realized through the calls to the internal AMOEBA function, based essentially on the Matlab's FMINS non-linear simplex implementation. The constraints are implemented by penalizing the target function.

The method provides clear advantage over the plain nonlinear simplex, and it has been proved to be specially useful for target functions with plenty of narrow local minima - standard traps for analitically based minimizers.

The method is not limited to continuos functions and does not require derivatives.

Additional help and information is available by calling HELP SIMPS, as well as looking into the documented code of M-files as included. A compete example GO.M is included, along with an exemplar target function FUN.M.

Authors:
Zeljko Bajzer (bajzer@mayo.edu) and Ivo Penzar (penzar@mayo.edu)
Mayo Clinic and Foundation, Rochester, Minnesota, USA

MATLAB release MATLAB 5.2 (R10)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (4)
08 Dec 2013 Richi

More specifically, the approach is different for the else cases for fxr<fv(n+1). In your implementation, there is either a small contraction or small reflection or a multi-contration (shrink). In the original paper, shrinking can be combined with small reflection and small contration. However, this is only a minor difference from the original paper and will of course work.

07 Dec 2013 Richi

Really good work. What is the basis for the approach? I saw, that the implementation is different from, e.g. Numerical Recepies. Is there a specific reason for this?

13 Jan 2012 LIAO xj  
21 May 2009 Chris Men

Excellent work.

Contact us