scilab/modules/optimization/demos/neldermead/fminsearch_tolx.sce

   1 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab\r
   2 // Copyright (C) 2009 - DIGITEO - Michael Baudin\r
   3 // Copyright (C) 2010 - DIGITEO - Allan CORNET\r
   4 // Copyright (C) 2012 - Scilab Enterprises - Adeline CARNIS\r
   5 //\r
   6 // This file must be used under the terms of the CeCILL.\r
   7 // This source file is licensed as described in the file COPYING, which\r
   8 // you should have received as part of this distribution.  The terms\r
   9 // are also available at\r
  10 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt\r
  11 \r
  12 \r
  13 function demo_fmin_tolx()\r
  14 \r
  15     mprintf(_("Running optimization ...\n"));\r
  16 \r
  17     function y = banana (x)\r
  18         y = 100*(x(2)-x(1)^2)^2 + (1-x(1))^2;\r
  19     endfunction\r
  20 \r
  21     x0 = [-1.2 1]\r
  22     opt = optimset ( "TolX" , 1.e-2 );\r
  23     [x , fval , exitflag , output] = fminsearch ( banana , x0 , opt );\r
  24 \r
  25     //\r
  26     // Display results\r
  27     //\r
  28     mprintf("x = %s\n",strcat(string(x)," "));\r
  29     mprintf("fval = %e\n",fval);\r
  30     mprintf("exitflag = %d\n",exitflag);\r
  31     mprintf("output.message:\n");\r
  32 \r
  33     for i =1:3\r
  34         mprintf(output.message(i));\r
  35         mprintf("\n");\r
  36     end\r
  37 \r
  38     mprintf("output.algorithm:%s\n",output.algorithm);\r
  39     mprintf("output.funcCount:%d\n",output.funcCount);\r
  40     mprintf("output.iterations:%d\n",output.iterations);\r
  41 \r
  42     //\r
  43     // Load this script into the editor\r
  44     //\r
  45     m = messagebox(_("View Code?"), "Question", "question", [_("Yes") _("No")], "modal")\r
  46     if(m == 1)\r
  47         filename = 'fminsearch_tolx.sce';\r
  48         dname = get_absolute_file_path(filename);\r
  49         editor ( dname + filename, "readonly" );\r
  50     end\r
  51 endfunction\r
  52 \r
  53 demo_fmin_tolx();\r
  54 clear demo_fmin_tolx;\r
  55 \r
  56 \r
  57 \r
  58 \r
  59 \r
  60 \r
  61 \r\r