[matlab] la dichotomie

[invite373ba4d2]

bonjour


j'ai un problém avec ce code là ,




quand je fait par exemple

f = @(x) 5*x^4 - 2.7*x^2 - 2*x + .5;
l = .1;
p = 0.5;
k=0.00000000000000000001
s=3
x = d(f, l, p,k,s);

le programme affiche :

Undefined function or variable "n_iter".

Error in d (line 29)
while( abs(d-c)>epsilon && drapeau ==0 && n_iter<n_max)


!!!
-----

[Antoane]

Bonsoir,

Tu n'as pas défini la variables n_iter, qui permet de compter combien de fois on a parcouru la boucle while(), et d'assurer que l'on entre pas dans une boucle infinie.

[invite373ba4d2]

et comment je définie variable n_iter ??!! comment savoir combient itérations on peux faire?

[Antoane]

et comment je définie variable n_iter ??!! comment savoir combient itérations on peux faire?

As-tu compris comment elle fonctionne, à quoi elle sert ?
Peux-tu fournir le reste de ton code ?

[A voir en vidéo sur Futura]

[invite373ba4d2]

voilla le reste:

 n_iter=n_iter+1;
    end
end

[Antoane]

Pour une raison ou pour une autre, il est possible que la condition :
abs(d-c)<epsilon
ne soit jamais remplie.
Pour éviter de se retrouver alors à tourner dans la boucle while, on ajoute un compteur (n_iter) qui va permettre de compter combien de fois on parcours la boucle. Si ce nombre est trop grand (c'est à dire >n_max), on comprendra qu'il y a un problème et que l'algorithme n'a pas réussi à converger. Dans ce cas, on sortira de la boucle (même si la solution n'a pas été trouvée) et on affichera un message d'erreur.

[invite373ba4d2]

mais ya un autre problem , comment afficher c et d ? ( les valeurs des bords de l'intervalle [c;d] qui contient la solution)

j'ai trouvé d'autre code

function [x e] = mybisect(f,a,b,n)

format long

c = f(a); d = f(b); 
if c*d > 0.0

end
disp('           x                  y')
for i = 1:n

    x = (a + b)/2; 
    y = f(x);
    disp([    x     y])
    if y == 0.0    
        e = 0;
        break      
    end

    if c*y < 0
        b=x;
    else
        a=x;
    end
end

x = (a + b)/2;
e = (b-a)/2;

a ce que je peux inclure la partie for i=1:n , dans le code précédent ou non?

[Antoane]

Bonsoir,

Tant que tu n'auras pas compris le principe de la méthode, il sera inutile de chercher à faire un code fonctionnel.
Une base de départ : https://fr.wikipedia.org/wiki/Dichotomie. Google's thy friend.

Récupérer des morceaux de code et essayer de les combiner sans les comprendre ne te fera pas avancer, et a peut de chance de donner un résultat bon et fonctionnel.

mais ya un autre problem , comment afficher c et d ? ( les valeurs des bords de l'intervalle [c;d] qui contient la solution)

Que ce passe-t-il dans matlab lorsque tu entres

c, d

dans l'invite de commande ?
Et

help display

Et

help sprintf