Discussion :

[Stepsbysteps]

Bonsoir à tous,

J'ai deux questions pour vous :
1) Mon if ne marche pas alors que les deux variables sont égales
2) Je voudrai trouver une valeur approchée et non une valeur juste

La valeur 0.0679 étant la première valeur du vecteur y

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
function [beta0]=tableF(m1,m2)
 
b_cal=m1./(sqrt(m2));
 
x=0.1:0.01:10;
y=F(x);
 
val_rechercher=0.0679;
 
% Recherche dans F(x) la position valeur la plus proche de b_cal
    for i=1:length(y)
        val=y(i);
        if val == val_rechercher
            indice=i;
            break;
        end
 
    end
 
% Recherche de la valeur de cette position dans le vecteur des beta
beta0=x(indice);
 
 
 
end

Merci à ceux qui pourront m'aider

[duf42]

Bonjour,

Plutôt que de faire un test d'égalité avec ==, fait plutôt un test sur la valeur absolue de la différence

Voir la FAQ : Pourquoi 0.3-0.2-0.1 est-il différent de 0 ?

Dans ton cas:

Code :

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

1
2
 
if abs(val - val_rechercher) < tol

Avec tol, l'écart maximal autorisé pour considérer que les 2 valeurs sont égales.

Une petite remarque, dans MATLAB tu peux éviter d'utiliser la boucle FOR, en remplaçant les lignes 10 à 18 par:

Code :

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

1
2
 
indice = find(abs(y - val_rechercher) < tol, 1, 'first');

Et tu pourrais même sans doute améliorer encore un peu en utilisant l'indexation logique (à vérifier que ca fonctionne bien sur ton code):

Code :

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

1
2
 
beta0 = x(abs(y - val_rechercher) < tol);

Duf

[Stepsbysteps]

super merci je regarderai ca.

J'étais passé sur ca

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
function [beta0]=tableF(m1,m2)
 
b_cal=m1./(sqrt(m2));
 
x=0.1:0.01:10;
y=F(x);
 
val_rechercher=b_cal;
eps=0.010;
 
% Recherche dans F(x) la position valeur la plus proche de b_cal
    for i=1:length(y)
        val=y(i);
        binf=val_rechercher-eps;
        bsup=val_rechercher+eps;
        if (val >= binf) && (val <= bsup)
            indice=i;
            goto(24);
        end
 
    end
 
% Recherche de la valeur de cette position dans le vecteur des beta
beta0=x(indice);
 
 
 
end

[FLB]

Fait attention, eps est une variable native de matlab :

Code :

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

1
2
3
>> eps
ans =
     2.22044604925031e-016

Il vaut mieux en général ne pas y toucher...

[le fab]

Une petite remarque, dans MATLAB tu peux éviter d'utiliser la boucle FOR, en remplaçant les lignes 10 à 18 par:

Code :

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

1
2
 
indice = find(abs(y - val_rechercher) < tol, 1, 'first');

Et tu pourrais même sans doute améliorer encore un peu en utilisant l'indexation logique (à vérifier que ca fonctionne bien sur ton code):

Code :

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

1
2
 
beta0 = x(abs(y - val_rechercher) < tol);

Duf

autre solution :

Code :

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

[~,indice] =  min(abs(y - val_rechercher));

[duf42]

Code :

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

[~,indice] =  min(abs(y - val_rechercher));

C'est sans doute la meilleure solution en effet

[Stepsbysteps]

Merci beaucoup à vous.
Vous êtes au top.

Un petit rappel sur l'utilisation du goto? pour sortir d'une boucle par exemple

[Jean Dumoncel]

Bonjour,

Un petit rappel sur l'utilisation du goto? pour sortir d'une boucle par exemple

C'est très simple, le goto n'existe pas dans matlab.

Si tu veux sortir de la boucle avant la fin du compteur, tu peux utiliser l'instruction break.

[shaiHulud]

Pour sortir d'une boucle, il n'y a que break, return, ou throw une error

Attention, ceci ne marche pas non plus

Code :

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

1
2
 
for k=1:10;disp(k);k=1000;end % pas de sortie boucle apres k=1 !!!