Discussion :

[d0n32]

Bonjour à toutes et à tous !

J'ai la fonction suivante, où Vg_ref est un vecteur 2x1 :

Code :

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

1
2
3
4
5
6
7
8
function y = Fobj2(Vg_ref)
 
Fvp = lambda_vp*n_vp*(alpha*sum(Vp_ref-Vp_mes)      -sum(Svp_vg*(Vg_ref-Vg_mes)))^2;
Fqg = lambda_qg*n_qg*(alpha*sum(Qg_max*q_ref-Qg_mes)-sum(Sqg_vg*(Vg_ref-Vg_mes)))^2;
Fvg = lambda_vg*n_vg*(alpha*sum(Vg_ref-Vg_mes)      -sum(Vg_ref-Vg_mes))^2;
 
y = Fvp + Fqg + Fvg;
end

Je souhaiterais simplement plotter cette fonction, disons sur l'intervalle [17;23]x[17;23]. J'ai essayé des tas de méthodes (plot3, mesh, surf...), mais je n'ai jamais réussi à obtenir autre chose que le message d'erreur suivant :

>> surf(X,Y,Fobj2);
??? Input argument "Vg_ref" is undefined.

Pourtant, c'est bien comme ça que je dois déclarer ma fonction dans une fonction de plot, non ?

Bref, si vous voyez où est l'erreur... Merci !

[soft001]

bonsoir,

J’ai un petit souci pour votre fonction.
Il faut que vous définissez les variables lambda_vp, lambda_qg,… ou bien les considérés comme des paramètres de votre fonction.
Il faut suivre ces étapes ;
définir vos variables :
lambda_vp=….
lambda_qg =….
….
Et après faire appel à votre fonction : surf(X,Y,Fobj2(les paramètre de votre fonction));
le résultat de votre fonction doit être une matrice.

[d0n32]

Bonjour,

Merci de votre réponse. Concernant les paramètres qui n'apparaissent pas dans la fonction, c'est par souci de simplification de la fonction aux yeux des membres du forum. J'ai fait le nécessaire pour les exprimer dans mon code, sans le faire apparaître dans mon précédent post. Désolé donc de ne pas avoir été assez clair !

Ensuite, je comprends bien ce que vous me conseillez de faire. Cependant, ma fonction doit accepter en argument un vecteur 2x1. Je ne peux simplement mettre, en suivant votre exemple :

surf(X,Y,Fobj2([X Y]'));

car le vecteur d'argument est alors le vecteur :

Code :

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

1
2
[17 18 19 20 21 22 23;
  17 18 19 20 21 22 23]

Et il manque manifestement des points.
Enfin, ma fonction ne renvoyant pas, comme vous l'avez justement fait remarqué, une matrice mais un scalaire, la méthode la plus adaptée doit être plot3 plutôt que surf ou mesh. Mais cela ne change rien au problème initial.

Merci !

[Jerome Briot]

L'idée est de faire ceci :

Code :

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

1
2
3
4
5
6
[X,Y] = meshgrid(17:23);
 
Z = Fobj2(X,Y);
 
figure
surf(X,Y,Z);

Maintenant, il ne reste plus qu'à modifier la fonction pour gérer directement des tableaux plutôt que des scalaires (voir les opérateurs .* et ./)

Une autre solution moins efficace consiste à calculer les termes de Z dans une boucle FOR-END

Code :

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

1
2
3
4
5
6
7
8
9
[X,Y] = meshgrid(17:23);
 
Z = zeros(size(X);
for i = 1:numel(Z)
   Z(i) = Fobj2([X(i) Y(i)]);
end
 
figure
surf(X,Y,Z);

[d0n32]

Merci beaucoup ! C'est bien l'idée que j'avais eu au départ (celle avec la boucle for/end), seulement, je n'avais pas eu l'idée d'utiliser meshgrid pour définir mes mailles. J'avais plutôt essayé de tout faire à la main..