Appel de variable dans une fonction déjà défini dans un prog principal

(Re-)Bonjour à tout le monde,

Voilà j'ai un autre problème, enfin disons que j'aimerai tester dans une première étape la qualité de l'optimisation avec ou sans contrainte sous Matlab quand je connais les vraies valeurs.

1er problème :

Je dispose d'un vecteur de taille n appelé u.
J'ai une fonction K dépendant du temps t et de cette valeur u qui elle est fixée dans une boucle quand je cherche le minimum de K.

J'ai crée ma fonction 'fun' K(t) mais elle dépend de u aussi (en fait de u(i) pour être précis).
Le problème c'est que je n'arrive pas à mettre tout ça en place pour avoir le minimum de ma fonction K et surtout la valeur de t pour laquelle j'attends ce minimum.

Code:

---

1 2 3

function K=mygen(t)   K=((u-m1)./S11)-t.*(S12./S11).*m1+t.*m2+(1./2)*((((u-m1)./S11)-t.*(S12./S11).^2).*S11+2.*((u-m1)./S11)-t.*(S12./S11).*t.*S12+(t.^2).*S22);

---

que j'appelle là dedans:

Code:

---

1 2 3 4 5

for i=1:n t0=0; u=u(i); [t(i),fval(i)] = fminunc(@mygen,t0); end

---

2ème problème:

Dans la réalité je ne peux que connaitre K(s,t) <- matrice et je vais chercher pas que le minimum, toujours pour u donné, de K(s,t).
Je cherche s et t tel que la dérivée première de K(s,t) en s soit égale à u(i) et la dérivée seconde de K(s,t) en t soit égale à 0.

Je n'arrive pas à formaliser le problème sous Matlab ...

:(:(:(

Bonjour tout le monde,

Alors j'ai remarqué que le problème venait de ma fonction 'fun' :

Code:

---

1 2 3

function K=mygen(t)   K=((10-m1)./S11)-t.*(S12./S11).*m1+t.*m2+(1./2).*((((10-m1)./S11)-t.*(S12./S11).^2).*S11+2.*((10-m1)./S11)-t.*(S12./S11).*t.*S12+(t.^2).*S22);

---

Car dedans j'utilise des variables que j'ai définis dans mon code principal mais je ne sais pas comment dire à l'ordinateur que m1 = l'ancien m1 sans avoir à retaper la valeur à chaque fois.

Voici le code erreur :

Code:

---

1 2 3 4 5 6 7

??? Undefined function or variable 'm1'.   Error in ==> mygen at 3 K=((10-m1)./S11)-t.*(S12./S11).*m1+t.*m2+(1./2).*((((10-m1)./S11)-t.*(S12./S11).^2).*S11+2.*((10-m1)./S11)-t.*(S12./S11).*t.*S12+(t.^2).*S22);   Error in ==> fminunc at 184 f = feval(funfcn{3},x,varargin{:});

---

Du coup si j'arrive à comprendre ça j'aurai résolu mon problème qui n'est pas lié à fmincon ou fminunc à prirori.

Hollow

salut,

une solution serait de tout mettre dans un seul fichier:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

%-- Declaration/affectation des variables m1 = 1; % exemple % ... le reste des variables   %-- Ta fonction MYGEN : elle verra les variables m1, ... mygen = @(t)( ((10-m1)./S11)-t.*(S12./S11).*m1+t.*m2+ ... (1./2).*((((10-m1)./S11)-t.*(S12./S11).^2).*S11+ ... 2.*((10-m1)./S11)-t.*(S12./S11).*t.*S12+(t.^2).*S22) );     %-- La boucle for i=1:n t0=0; u=u(i); [t(i),fval(i)] = fminunc(@mygen,t0); end

---

++ bonne chance

Bonjour,

Je vais essayer.
Je savais pas que c'était possible :oops::oops:

Merci encore je vous tiens au courant!

J'ai donc essayer, ça me parait la bonne technique mais j'ai encore un problème de code apparemment :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19

function f=myfunction(u,m1,m2,S11,S12,S21,S22)   ...   Tf=zeros(n,1); fval=zeros(n,1);   %-- Function MYGEN : elle verra les variables m1, ... mygen = @(t)(((u-m1)./S11)-t.*(S12./S11).*m1+t.*m2+(1./2).*((((u-m1)./S11)-t.*(S12./S11).^2).*S11+2.*((u-m1)./S11)-t.*(S12./S11).*t.*S12+(t.^2).*S22));     %-- La boucle for i=1:n t0=0; u=u(i); [Tf(i),fval(i)] = fminunc(@mygen,t0); end

---

Malheureusement j'ai ce message d'erreur lorsque je lance la fonction :

Code:

---

1 2 3 4

??? Error: File: myfunction.m Line: 118 Column: 32 "mygen" was previously used as a variable, conflicting with its use here as the name of a function or command.

---

Et de plus, je penses qu'il serait préférable pour moi de mettre mygen=@t dans la boucle car cela dépend de u(i) aussi.

Hollow

L'operateur '@' fourni un identifiant (handler) de la fonction. Dans notre cas, on dispose déja de cela, c'est mygen :

Code:

---

1 2 3 4

% ... [Tf(i),fval(i)] = fminunc(mygen,t0); % Sans @ % ...

---

Sinon, je ne vois pas l'intérêt de ta boucle FOR ... Tu minimise la même fonction en boucle.

Edit :

Si tu veux que u change dans ta fonction (donc dans ta boucle) :

Code:

---

1 2 3 4 5 6 7 8

%-- La boucle for i=1:n t0=0; u=u(i); %-- Sinon, la fonction MYGEN ne mettra pas à jour sa definition ;) mygen = @(t)(((u-m1)./S11)-t.*(S12./S11).*m1+t.*m2+(1./2).*((((u-m1)./S11)-t.*(S12./S11).^2).*S11+2.*((u-m1)./S11)-t.*(S12./S11).*t.*S12+(t.^2).*S22)); [Tf(i),fval(i)] = fminunc(mygen,t0); end

---

++

Salut,

Alors j'ai fait ça avant que tu m'envoies une solution plus "propre" :

Code:

---

1 2 3 4 5 6 7 8 9 10

Tf=zeros(n,1); fval=zeros(n,1);     %-- La boucle for i=1:n t0=0; [Tf(i),fval(i)] = fminunc(@(t)(((u(i)-m1)./S11)-t.*(S12./S11).*m1+t.*m2+(1./2).*((((u(i)-m1)./S11)-t.*(S12./S11).^2).*S11+2.*((u(i)-m1)./S11)-t.*(S12./S11).*t.*S12+(t.^2).*S22)),t0); end

---

Le problème c'est que même si u varie, et donc mygen varie aussi, il trouve toujours le même minimum ... or ce n'est pas le cas quand je fais le calcul explicite sur papier.

Alors soit il y a un problème parce que je pars de t0=0 et il s'arrête à un minimum local (ce qui m'étonnerait car ça voudrait dire que fminunc est pas très bien pensé lol) soit j'ai fait une boulette ce qui ne m'étonnerait pas du tout :mrgreen:

Je vais essayer de voir ce qui ne va pas là ...

PS: quand je cherche le minimum d'une fonction c'est bien équivalent à chercher quand la dérivée est nulle (avec les bonnes hypothèses de décroissance, croissance, etc) ?

EDIT

Même en changeant le t0 ça reste contant donc le problème vient demoi ce qui me rassure :D

Pour ton probléme de solution inchangée. Je soupçonne que u à un effet additif sur la fonction MYGEN (le minimum ne change pas dans ce cas). Fait un PLOT de MYGEN pour différents valeurs deu ...


Remarques, on sort de Matlab maintenant ;) .


++

J'ai fait des plot pour différentes valeurs de u et effectivement ça ne change rien du tout.

Or ça devrait!

Peut-être ai-je mal tapé ma fonction.

En tout cas tout ce qui concerne Matlab est résolu :king::king::king:
Merci encore à vous!

J'espère m'améliorer quand même ...

Hollow