Discussion :

[kitts]

Bonjour,

Codant en Matlab un algorithme de type gradient projeté, j'ai à coder deux traitements pour le calcul de deux réels.

f est non-linéaire de Rn dans R, passée en paramètre d'entrée du programme, x et y sont connus avant les calculs de max / min.

1) amax = max{a tel que f(x+a*y) <= 0}
2) b = min{norme(x+a*y) tel que 0<=a<=amax}

J'ai essayé des traitements avec
1)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
 
function [c,ceq] = mycon(f,a,x,y)
         c = @(a)(feval('f',x-a*y));
         ceq = [];
end
 
! dans le programme principal
amax = fmincon(@(a) feval('f',x-a*k), 1, [],[], [],[],0,[], @(a)mycon(f,a,x,y));

2)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
! dans le programme principal
[b,fb] = fmincon(@(a) norm(x+a*y,2), amax/2 ,[],[],[],[], 0, amax);

Le traitement 1) me met une erreur : je ne dois pas avoir tout compris pour prendre en compte la contrainte f(...)<=0
Le traitement 2) me met "Warning: Large-scale (trust region) method does not currently solve this type of problem, switching to medium-scale (line search)."

Si quelqu'un peut m'aider, peut-être que je n'utilise pas la fonction Matlab la plus appropriée pour faire ce que j'ai à coder...
Merci.

[christophe_halgand]

Quelle est ton erreur ?

[kitts]

Sur le traitement 1), cela me fait à l'éxecution :

Warning: Large-scale (trust region) method does not currently solve this type of problem, switching to medium-scale (line search).
> In fmincon at 260
In methode_gradient at 102
Warning: To support parenthesis notation for invocation, "indexing" a scalar function handle by ":" will continue to work in R14, but may be illegal or may work differently in R15. To prepare for R15, and to avoid this warning, use cell arrays of function handles instead of arrays. For more information, type 'help function_handle' and see the section at the end entitled Note on Backward Compatibility.
> In fmincon at 406
In methode_gradient at 102
Warning: non-scalar function handles will not be supported past R14
> In optim\private\nlconst at 175
In fmincon at 498
In methode_gradient at 102
Warning: non-scalar function handles will not be supported past R14
> In optim\private\nlconst at 176
In fmincon at 498
In methode_gradient at 102
??? Error using ==> vertcat
The following error occurred converting from double to function_handle:
Attempt to execute SCRIPT function_handle as a function.

Error in ==> optim\private\nlconst at 176
c = [ Aeq*XOUT-Beq; nceq; A*XOUT-B; ncineq];

Error in ==> fmincon at 498
[X,FVAL,lambda,EXITFLAG,OUTPUT,GRAD,HESSIAN]=...

Error in ==> methode_gradient at 102
amax = fmincon(@(a) feval('f',x-a*k), 1, [],[], [],[], 0,[], @(a)mycon(f,a,x,y));