Discussion :

[justin0203]

Bonsoir,

je suis entrain faire un programme sur matlab pour trouver le minimum d'un fonction

voici le programme que j'ai fais il m'affiche comme erreur

[x,y,k] = dicho(1e-50,1000,2,1)
??? Reference to a cleared variable a0.

Error in ==> dicho at 7
a(1)=a0;

fichier func.m

Code :

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

1
2
3
function [ f ] = func( x )
f = x-1-log(x);
end

fichier dicho.m

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
33
34
35
36
37
38
39
function [ a,b,k ] = dicho(a0,b0,tolerance,epsilon)
 
clear all;
% initialisation
 
k=1;
a(1)=a0;
b(1)=b0;
 
 
%Calcule de milieu du domiane x(k)
 
    while abs(b(k)-a(k))>tolerance
        x(k) = 0.5*(a(k) + b(k));
        c(k) = x(k) - (epsilon/2);
        d(k) = x(k) + (epsilon/2);
 
        if func(c(k)) <= func(d(k))
 
            % on prends la partie droite
 
            a(k+1) = a(k);
            b(k+1) = b(k);
        else
 
            % on prends la partie gauche
 
            a(k+1) = c(k);
            b(k+1) = b(k);
        end
        k=k+1;
    end
    k
    a(k)
    b(k)
 
 
 
end

Merci en avance

Cordialement,

Justin

[Jean Dumoncel]

Bonjour,

il ne faut pas mettre de clear all au début de ta fonction. clear efface toutes les variables, d'ou le message d'erreur.

[justin0203]

Bonsoir,

merci de votre réponse effectivement si j'efface toutes les variables ca ne fonctionnera pas! =)

Parcontre ca mets trop de temps pour compiler! et pour donner les resultas

[123djimy]

Bonjour,

Code :

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

1
2
3
4
5
6
7
8
9
 if func(c(k)) <= func(d(k))
            % on prends la partie droite        
            a(k+1) = c(k);
            b(k+1) = b(k);
   else
            % on prends la partie gauche
            a(k+1) = a(k);
            b(k+1) = d(k);
   end

[membreComplexe12]

clear all :
il ne doit jamais être mis dans les fonctions. Il est mis par contre dans le programme/fichier principal pour effacer ce que tu as pu faire auparavant.

a part cela :
ton code me semble bien mais je n'aurais pas fait des vecteurs pour les ak, bk ... j'aurais juste écrasé la précédente valeur. (mais après ça veut pas dire que ta solution ne marche pas (il y a tjs plusieurs solution en programmation) mais moi c'est ma préférence)

en algorithmie on donne souvent une relation de ce genre:

a(k+1)=a(k)

mais quand tu programmes ceci ça veut pas forcement dire que tu dois faire un vecteur "a" et que la colonne (k+1) sera la meme que la colonne (k) ça veut dire que la prochaine valeur de a sera égale à la précédente:

a=a

sinon à part ça la structure du programme m'a l'air très bien avec rien de superflu ou d'inutile