home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 4.ddi / OPTIM.DI$ / MINIMAX.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  2.5 KB  |  86 lines
  1. function [x,OPTIONS]=minimax(FUN,x,OPTIONS,VLB,VUB,GRADFUN,P1,P2,P3,P4,P5,P6,P7,P8,P9,P10)
  2. %MINIMAX Finds a minimax solution of a function of several variables.
  3. %    MINIMAX solves the following problem:
  4. %    min (max {FUN(X} )  where FUN and X can be vectors or matrices.
  5. %     X 
  6. % 
  7. %    X=MINIMAX('FUN',X0) starts at X0 and finds a minimax solution to 
  8. %    the function which is described in FUN (usually an M-file: FUN.M).
  9. %    The function 'FUN' should return functions to be minimized: F=FUN(X). 
  10. %
  11. %    X=MINIMAX('FUN',X,OPTIONS) allows a vector of optional parameters to 
  12. %    be defined. For more information type HELP FOPTIONS.
  13. %    
  14. %    X=MINIMAX('FUN',X,OPTIONS,VLB,VUB) defines a set of lower and upper
  15. %    bounds on the design variables, X, so that the solution is always in 
  16. %    the range VLB < X < VUB. 
  17. %
  18. %    X=MINIMAX('FUN',X,OPTIONS,VLB,VUB,'GRADFUN') allows a function 
  19. %    'GRADFUN' to be entered which returns the partial derivatives of the 
  20. %    functions at X stored column-wise:  GRADS = GRADFUN(X). 
  21. %
  22. %    X=MINIMAX('FUN',X,OPTIONS,VLB,VUB,GRADFUN,P1,P2,..) allows
  23. %    coefficients, P1, P2, ... to be passed directly to FUN:
  24. %    [F,G]=FUN(X,P1,P2,...). Empty arguments ([]) are ignored. 
  25.  
  26. %    Copyright (c) 1990 by the MathWorks, Inc.
  27. %    Andy Grace 7-9-90.
  28.                      
  29. %
  30.  
  31. if nargin < 3, OPTIONS=[]; end
  32. if nargin < 4, VLB=[]; end
  33. if nargin < 5, VUB=[]; end
  34. if nargin < 6, GRADFUN=[]; end
  35. lenopt = length(OPTIONS); 
  36. if ~lenopt, OPTIONS=0; end
  37.  
  38. xnew=x(:);
  39. OPTIONS=foptions(OPTIONS);
  40. OPTIONS(7) = ~OPTIONS(7);
  41. neqcstr=OPTIONS(15); 
  42.  
  43. if ~any(FUN<48)
  44.     evalstr1 = ['[f,g]=',FUN,];
  45.         evalstr1=[evalstr1, '(x'];
  46.         for i=1:nargin - 6
  47.                 evalstr1 = [evalstr1,',P',int2str(i)];
  48.         end
  49.         evalstr1 = [evalstr1, ');'];
  50. else
  51.         evalstr1=[FUN,';'];
  52. end
  53.  
  54. eval(evalstr1)
  55. ncstr=length(f(:));
  56. WEIGHT = ones(ncstr,1);
  57. GOAL = zeros(ncstr,1);
  58.  
  59. evalstr2='';
  60.  
  61. if length(GRADFUN)
  62.     if ~any(GRADFUN<48) % Check alphanumeric
  63.         evalstr2 = ['[gf,gg]=',GRADFUN,'(x'];
  64.         for i=1:nargin - 6
  65.             evalstr2 = [evalstr2,',P',int2str(i)];
  66.         end
  67.         evalstr2 = [evalstr2, ');'];
  68.         gfun  = 'goalgra';
  69.     else
  70.         evalstr2=[GRADFUN,';'];
  71.     end
  72. else
  73.     gfun = [];
  74. end
  75.  
  76. evalstr = ['[xnew, OPTIONS] = constr(''goalfun'',[xnew;0],OPTIONS,VLB,VUB,gfun,neqcstr,evalstr1,evalstr2,WEIGHT,GOAL,x'];
  77.  
  78. for i=1:nargin - 6
  79.     evalstr = [evalstr,',P',int2str(i)];
  80. end
  81. evalstr = [evalstr, ');'];
  82.  
  83. eval(evalstr)
  84. OPTIONS(7) = ~OPTIONS(7);
  85. x(:) = xnew(1:length(xnew)-1);
  86.