Discussion :

[namklaw]

Bonjour à tous,

Je voudrais déclarer et utiliser une fonction à 2 variables en passant un seul paramètre, donc un vecteur de dimension 2 dont j'utiliserai les composantes pour définir ma fonction. Je voudrais également retourner la valeur du gradient et du Hessian de ma fonction j'ai écrit le code suivant mais cela ne donne pas de résultats concluants ^^'.

Code :

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

1
2
3
4
5
function [f,g,H]=Rosenbrock(x)
 
f=(x[1]-1)^2+100*((x[1]^2)-x[2])^2;
g=gradient(f,x);
H=hessian(f,x);

Pouvez vous m'éclairez ou me conseillez je vous en serai grandement reconnaissant

[Jerome Briot]

L'indexage des tableaux se fait à l'aide de parenthèses dans MATLAB :

Code :

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

f=(x(1)-1)^2+100*((x(1)^2)-x(2))^2;

[namklaw]

D'accord merci pour cette précision , mon programme compile et s’exécute mais ne m'affiche que le résultat de la fonction f en x il ne m'affiche ni mon gradient ni mon hessian, je n'ai pas utilisé la bonne syntaxe pour ces deux fonctions de matlab ?

[Gooby]

Bonjour,

Peux tu nous faire voir comment tu utilises ta fonction?

[namklaw]

Bonjour,

Et bien en fait je voudrais utiliser la fonction de rosenbrock pour tester le code d'une méthode quasi newton (BFGS). Mais pour le moment j'ai juste taper la commande suivante :

Code :

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

1
2
>> x=[1,2]
>> Rosenbrock(x)

Et ça me sort la bien la valeur de f en (1,2) mais il me dit que la fonction hessian n'est pas définie pour des arguments de type double et pour le gradient il me renvoie juste 0 qui correspond à la somme de mes 2 dérivées partielles mais pas les composantes de mon vecteur gradient: (-400 , 400).

[Gooby]

Pour la fonction Hessian, possèdes-tu bien la MATLAB Symbolic Math Toolbox, nécessaire à l'usage de la fonction?

Pour la fonction gradient, je ne comprends pas bien l'utilisation que tu en fais. Le gradient sur matlab s'applique à des vecteurs, là tu l'appliques à un double.

[namklaw]

Comment sait on si on possede MATLAB Symbolic Math Toolbox, j'utilise la version R2012a de matlab. Concernant mon gradient en fait je voudrais que g me renvoie le vecteur du gradient avec les composantes évaluées en x1 et x2 soit 1 et 2 dans mon exemple, la syntaxe de mon code est donc incorrect ? Si je comprend bien mon code applique le gradient à f évaluée en x1 x2 et n'effectue pas les dérivées partielles puis l'évaluation ?
Le problème est alors similaire pour le hessian et l'erreur n'est peut etre pas du a la toolbox non ?

[Gooby]

Pour connaître les toolbox installées voir la FAQ de développez.com - Comment connaitre la liste des Toolbox installées sur un ordinateur ?

Pour définir une fonction sur MATLAB: Création de fonctions anonymes

[namklaw]

Pour être un peu plus précis, le résultat auquel je veux aboutir est le suivant :

J'ai une fonction :

Code :

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

function unc_opt(fname,x0)

qui en théorie doit me donner le minimum de la fonction fname en utilisant une methode quasi newton avec comme point de depart x0.

Et je voudrais appliquer cette fonction/ méthode à ma fonction de Rosenbrock qui doit être déclaré de cette façon :

Code :

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

function [f,g,H]=fname(x)

Car dans l’idéal je dois pouvoir appliquer ma méthode quasi newton à d'autre fonction. Or ma fonction de rosenbrock à 2 composantes x1 et x2 c'est donc pour cela que je devais passer un vecteur en argument mais ensuite j'ai du mal me servir des composantes de ce dernier pour coder ma fonction et ensuite obtenir le gradient et le laplacien ? Pour être tout à fait honnête je suis un peu perdu je ne sais pas si l'utilisation d'un vecteur en argument pour x était vraiment la bonne démarche mais vu que x a 2 composantes je ne vois pas vraiment d'autre solution...

edit : Merci gooby, j'ai donc vérifié et j'ai bien la toolbox symbolic math