Discussion :

[hollowdeadoss]

(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 :

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

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 :

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

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 ...

[hollowdeadoss]

Bonjour tout le monde,

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

Code :

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

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 :

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

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

[mr_samurai]

salut,

une solution serait de tout mettre dans un seul fichier:

Code :

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

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

[hollowdeadoss]

Bonjour,

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

Merci encore je vous tiens au courant!

[hollowdeadoss]

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

Code :

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

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 :

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

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

[mr_samurai]

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

Code :

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

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 :

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

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

++