Comment corriger ce code?

Version imprimable

J'ai besoin de faire une fonction qui permet d'appliquer la méthode de dichotomie afin de rechercher
la racine approchée d'une fonction non linéaire.J'ai écrit le code suivant mais il y a des erreurs.Qu'est ce que
je dois rectifier?

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
function [r, fr, compt]=Dichotomie(x, f, e, Nmax,a,b)
compt=1;
pas=0.01;
a(1)=a;
b(1)=b;
x0(1)=(a+b)/2;
dist=2*e;
while ((compt<Nmax) && ( dist>e ))
    u=f((a(compt)-a)/pas+1);
    v=f((x0(compt)-a)/pas+1);
    if(u.*v<0)
        b(compt+1)=x0(compt); a(compt+1)= a(compt);
    end
    if(u.*v>0)
        a(compt+1)=x0(compt); b(compt+1)=b(compt);
    end
    if(u.*v==0)
        dist=0;
    end
    compt=compt+1;
    x0(compt)=(a(compt)+b(compt))/2;
        dist=abs(x0(compt)-x0(compt-1)); 
end
r=x0(compt) ;
fr=dist;
fprintf('La racine de l''équation est: %f obtenue aprés %f itérations\n ',r,compt);

27/04/2008, 16h41
Caro-Line

Comme on est dimanche et que je n'ai pas le courage d'essayer de faire tourner ton code, peux-tu nous dire quel est le message d'erreur et à quelle ligne il correspond ?
Et aussi de quel type et quelle taille sont les entrées de ta fonction ?
Un exemple de ces entrées serait aussi le bienvenu.

Voici les erreurs qui apparaissent
Code:

1 2 3 4 ??? Input argument "a" is undefined. Error in ==> Dichotomie at 4 a(1)=a;

27/04/2008, 17h29
Dam2227

Je crois que le message d'erreur est assez clair. Matlab ne reconnait pas la variable a. Es-tu sûr de l'avoir bien définie auparavant et de la passer correctement en argument dans ta fonction?

Salut,

Essayes ceci :

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
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
 
function [r, fr, compt]=Dichotomie(f, a, b, e, Nmax)
% [r, fr, compt]=Dichotomie(f, a, b, e, Nmax)
% Dichotomie function
 
if nargin < 5
    Nmax = 100;
end
 
if nargin < 4 
    e = 1.0e-8; % standard 
end
 
if nargin < 3 
   error('Nombre d''arguments insuffisant!') 
end
 
if f(a)*f(b) > 0 
   error('Condition : f(a)*f(b) < 0 insatisfaite!!') 
end
 
if a >= b 
   error('Condition : a < b insatisfaite!!') 
end
 
compt=1;
fr = 2*e;
 
% Boucle principale
while compt < Nmax && fr > e
 
    if  f(a) == 0
        r = a;
        fr = b-a;
        break
    elseif f(b) == 0
        r = b;
        fr = b-a;
        break
    elseif f(a)*f((a+b)/2) <= 0
        b = (a+b)/2;
    else
        a = (a+b)/2;
    end
    r = (a+b)/2;
    fr = b-a;
    compt=compt+1;
end
 
fprintf('La racine de l''équation est: %f obtenue aprés %f itérations\n ',r,compt);

problème de dichotomie

J'ai besoin de rechercher la racine d'une équation par la méthode de dichotomie.Cette équation est f(x)=x-e.sin(x) =0, x dans [1,10].
Pour créer cette fonction et la traiter j'ai utilisé ce code du fichier test.m:

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
clear all;
close all;
clc;
pas=0.01;
% On définit la précision de la longueur de l'intervalle dans lequel on désire avoir la solution 
e=10^(-3);
% On définit un nombre d'itérations maximal aprés lequel le traitment sera interrompu quelque soit le cas.  
Nmax=100;
% On définit les bornes de l'intervalle 
a=1;
b=10;
% On définit la variable de la fonction
 x1=1:pas:10;
 % On définit la fonction
 L=length(x);
for i=1:L
    f1(i)=x(i)-exp(1)*sin(x(i));
end
 
% On visualise la fonction
figure(1),plot(x,f)
% On définit r: la solution estimée,fr son image par la fonction f1,compt: le nbre d'itérations
[r, fr, compt]=Dichotomie(x, f, e, Nmax,a,b);
 
% Fin du code

Voici le code de la fonction de dichotomie du fichier Dichotomie.m:

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
function [r, fr, compt]=Dichotomie( f, e, Nmax,a,b)
compt=1;
pas=0.01;
a(1)=a;
b(1)=b;
x0(1)=(a+b)/2;
% On initialise un detecteur de précision >e
dist=2*e;
 
while ((compt<Nmax) && ( dist>e ))
    u=f((a(compt)-a)/pas+1);
    v=f((x0(compt)-a)/pas+1);
    if(u.*v<0)
        b(compt+1)=x0(compt); a(compt+1)= a(compt);
    end
    if(u.*v>0)
        a(compt+1)=x0(compt); b(compt+1)=b(compt);
    end
    if(u.*v==0)
        dist=0;
    end
    compt=compt+1;
    x0(compt)=(a(compt)+b(compt))/2;
        dist=abs(x0(compt)-x0(compt-1)); 
end
r=x0(compt) ;
fr=x0(compt)-exp(1)*sin(x0(compt));

Voici la reponse de MATLAB:

Citation:

??? Subscript indices must either be real positive integers or logicals.

Error in ==> Dichotomie at 12
v=f((x0(compt)-a)/pas+1);

Error in ==> test at 23
[r, fr, compt]=Dichotomie(f, e, Nmax,a,b);

29/04/2008, 11h24
LordPeterPan2
Code:

1 2 3 4 ??? Subscript indices must either be real positive integers or logicals. Error in ==> Dichotomie at 12 v=f((x0(compt)-a)/pas+1);
Le message te dis texto qu'un indice doit être strictement positif et entier ou logique.

Dans l'expression incriminé l'indice est : (x0(compt)-a)/pas+1

Donc je pencherai plutôt pour dire que (x0(compt)-a)/pas ne doit pas être entier.
29/04/2008, 13h28
Caro-Line

Cette erreur est dans la :faq:
Subscript indices must either be real positive integers or logicals.