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.