Discussion :

[Audinet]

Bonjour,

Dans le cadre d'un projet pour l'université, je dois écrire un programme pour calculer le minimum d'une fonction donnée par l'utilisateur.

L'algorithme donné est le suivant:


Voici mon code:

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
20
21
22
23
24
25
26
27
28
29
30
31
32
function r=f(x)
    r=x*x;
endfunction
 
function res= minimum(f,l,T,bet,Niter,N)
 
    x0=rand(1,1);
    x=[];
    for i=1:N
        x(1)=x0;
        A=x(i)-l
        B=x(i)+l
        y=grand(1,1,'bet',A,B)
        if f(y)<=f(x(i))then 
            x(i+1)=y;
        elseif f(y)>=f(x(i)) then 
            r=rand(1,1);
            p=exp(-(f(y)-f(x(i)))/T);
            if (r<p)
                then
                x(i+1)=y;
            else 
                x(i+1)=x(i);
            end
        end
        if N==Niter then 
            T=bet*T;
            l=bet*l;
        end
    end
end
endfunction

A la ligne 14

if f(y)<=f(x(i))then

Scilab me renvoie:

Opérateur, virgule ou point-virgule manquant.
at line 10 of function minimum called by :
endfunction

Voilà, je ne sais pas trop ce qui coince ici? J'ai bien le point virgule à la fin de l'expression après le premier "then". La fonction f(x)->x^2 est définie et marche correctement. Il y a surement autre chose que je ne comprends pas..

Merci en avance pour toute aide possible

[Audinet]

if f(y)<=f(x(i))then

A la ligne 14 Scilab me renvoie:

Opérateur, virgule ou point-virgule manquant.
at line 10 of function minimum called by :
endfunction

Voilà, je ne sais pas trop ce qui coince ici? J'ai bien le point virgule à la fin de l'expression après le premier "then". La fonction f(x)->x^2 est définie et marche correctement. Il y a surement autre chose que je ne comprends pas..

Merci en avance pour toute aide possible

Rebonjour,

Je viens de m'aperçoir qu'il y a une erreur dans mon code. La ligne 14 est

Code :

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

if (f(y)<=f(x(i)))then

Donc pas d'oubli de parenthèse.. J'obtiens la même erreur de Scilab.

Quel est l'élément qui manque ou l'erreur dans le programme qui fait que Scilab pense qu'il manque qqch? Je ne vois pas.

Merci encore

[Jerome Briot]

Attention à la syntaxe du code Scilab…

Le mot-clé then doit être précédé d'un espace et doit se trouver sur la même ligne que la condition
Il y a aussi un end en trop dans le code

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
20
21
22
23
24
25
26
27
28
29
30
function r=f(x)
    r=x*x;
endfunction
 
function res= minimum(f,l,T,bet,Niter,N)
 
    x0=rand(1,1);
    x=[];
    for i=1:N
        x(1)=x0;
        A=x(i)-l
        B=x(i)+l
        y=grand(1,1,'bet',A,B)
        if f(y)<=f(x(i)) then
            x(i+1)=y;
        elseif f(y)>=f(x(i)) then
            r=rand(1,1);
            p=exp(-(f(y)-f(x(i)))/T);
            if (r<p) then
                x(i+1)=y;
            else 
                x(i+1)=x(i);
            end
        end
        if N==Niter then 
            T=bet*T;
            l=bet*l;
        end
    end
endfunction

[Audinet]

Bonsoir et merci,

Effectivement avec les espaces et le END en mois ça fonctionne.

Merci

Bonne soirée