Compréhension fonction optimset

Version imprimable

Bonjour,

Je dois travailler sur un fichier qui comprend cette fonction :

Code:

options=optimset('Display','iter','Diagnostics','on','MaxFunEvals',100000000,'MaxIter',1000,'TolFun',1e-7,'TolX',1e-10,'NonlEqnAlgorithm','lm','LineSearchType','cubicpoly','LargeScale','on','Hessian','on');

Lorsque j'exécute le fichier, le message d'erreur suivant s'affiche :
Code:

1 2 3 4 5 Error using optimset (line 199) The NonlEqnAlgorithm option is no longer valid. Set the Algorithm option instead. Error in modele_aeraulique_36_zones (line 138) options=optimset('Display','iter','Diagnostics','on','MaxFunEvals',100000000,'MaxIter',1000,'TolFun',1e-7,'TolX',1e-10,'NonlEqnAlgorithm','lm','LineSearchType','cubicpoly','LargeScale','on','Hessian
Pourriez-vous m'expliquer à quoi sert cette fonction et comment résoudre ce problème ?

Merci

18/07/2013, 16h15
le fab

salut

tu dois utiliser un programme d'optimisation qui a été développé sous une version antérieure de Matlab
la fonction optimset sert à définir les options pour l'optimisation (fonctions comme fmincon ou lsqnonlin par exemple)

là visiblement tu utilise un argument obsolète pour définir un algo de levenberg-marquardt (lm)
il faut donc enlever la paire d'arguments 'NonlEqnAlgorithm','lm' et la remplacer par ce qui va bien, soit probablement 'Algorithm','levenberg-marquardt'

tu aura probablement un soucis avec la paire d'arguments 'LineSearchType','cubicpoly'

peux tu en dire plus sur ton programme et surtout sur la fonction d'optimisation utilisant la variable options ainsi définie ?

Fabien

En effet, le message suivant s'affiche :

Code:

1
2
3
4
5
Warning: The LineSearchType option is no longer valid. It was
only used by the Gauss-Newton algorithm, which is no longer
used in Optimization toolbox solvers. 
> In optimset at 223
  In modele_aeraulique_36_zones at 138

Mon programme me permet d'étudier les débits d'air entre les différentes zones d'un bâtiment. Au début du programme, je définis les températures de chaque zone, les zones et leurs caractéristiques (surface des murs, hauteurs, emplacement des ventilations).

La fonction qui utilise la variable options est la suivante :

Code:

[x,fval,exitflag,output]=zonal(P,options);

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
function [x,fval,exitflag,output]=zonal(xstart,options)
fun = @residus;
%options=optimset('Display','iter','MaxFunEvals',100000000,'MaxIter',1000,'TolFun',1e-10,'TolX',1e-6,'NonlEqnAlgorithm','dogleg','LineSearchType','cubicpoly','LargeScale','on');   % Option to display output
[x,fval,exitflag,output] = fsolve(fun,xstart,options)  % Call optimizer
 
function F=residus(P)
 
global H
global M
global A
global T
global F
global debvent
 
 
nz=38;
g=9.81;
Cd=0.83;
R=287;
b=0.5;
 
%calcul des débits
[debit]=FLOW(P,T,b,R,nz,M,Cd,A,g,H);
 
%résidus
F=[-3*debvent -12*debvent -3*debvent 0 0 0 0 0 0 0 0 0 2*debvent 8*debvent 2*debvent 0 0 0 0 0 0 0 0 0 debvent 4*debvent debvent 0 0 0 0 0 0 0 0 0 0 0];
 
for i=1:38
    for j=1:38
        F(i)=F(i)-debit(i,j)+debit(j,i);
 
    end
end

Et la fonction FLOW est la suivante :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
function [debit]=FLOW(P,T,b,R,nz,M,Cd,A,g,h)
 
 
 
rho0=1.21;
T0=293.15;
 
for i=1:nz
    rho(i)=rho0*T0/T(i);
end
 
Cd=0.83;
 
for i= 1:nz
    for j=1:nz
        if abs(M(i,j))==1
            debit(i,j)=(Cd*((2*rho(i))^0.5)*A(i,j)*(abs(P(i)-P(j)-0.5*(rho(i)*g*h(i)+rho(j)*g*h(j))))^b)*(P(i)-P(j)-0.5*(rho(i)*g*h(i)+rho(j)*g*h(j)))/abs((P(i)-P(j)-0.5*(rho(i)*g*h(i)+rho(j)*g*h(j))));
        else
            debit(i,j)=0;
        end
    end
end

18/07/2013, 17h26
le fab

bref, si tu remplace le premier couple d'argument comme je te l'ai conseillé ey que tu vire le 2ème, ca devrais marcher

en effet le deuxième argument n'était pas compatible avec levenberg-marquardt

et tu as aussi probablement l'argument large-scale en trop
18/07/2013, 17h45
Rhius

J'ai fait les modifications et maintenant mon programme fonctionne.

Merci pour ton aide.
19/07/2013, 08h54
le fab

juste un remarque générale en optimisation

ton programme fonctionne, c'est bien
mais fonctionne t'il comme prévu et surtout comme avant (sur une ancienne version de matlab) depuis que tu as été obligé de changer des paramètres ?
c'est toujours bien quand on fait ce genre d'exercice de comparer la solution trouvée avec la solution qu'il proposait avant (si tu en a gardé une)

Fabien

Au passage, tu peux améliorer ton code comme ceci :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
function debit = FLOW(P,T,b,R,nz,M,Cd,A,g,h)
 
rho0 = 1.21;
T0 = 293.15;
 
rho = rho0*T0./T;
 
Cd = 0.83;
 
debit = zeros(nz);
 
for i = 1:nz
    for j = 1:nz
        if abs(M(i,j))==1
            debit(i,j) = (Cd*((2*rho(i))^0.5)*A(i,j)*(abs(P(i)-P(j)-0.5*(rho(i)*g*h(i)+rho(j)*g*h(j))))^b)*(P(i)-P(j)-0.5*(rho(i)*g*h(i)+rho(j)*g*h(j)))/abs((P(i)-P(j)-0.5*(rho(i)*g*h(i)+rho(j)*g*h(j))));       
        end
    end
end