Discussion :

[tenderstoune]

Bonjour

Compte tenu de la FAQ : Pourquoi 0.3-0.2-0.1 est-il différent de 0 ? j ai essayé de faire un simple test d'égalité mais j'arrive toujours pas à le faire

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
exemple1:
 
a=rand(4,3)
 
a =
 
    0.7711    0.5628    0.1346
    0.3355    0.1243    0.5207
    0.9937    0.5756    0.5977
    0.1816    0.6876    0.1447
 
>> find(abs(a-05628)<eps)
 
ans =
 
   Empty matrix: 0-by-1

exemple2

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
 
%fonction
function ID = get_ID( a,CF )
for i=1:size(a,1)
    B=a(i,3);
 
if (abs(B-CF)<eps)
    ID= a(i,4);
end
 
end
return;
%execution
 
a=topo(4,1,1,0)
 
a =
 
    0.8868    0.2778    0.2178    1.0000
    0.1109    0.7176    0.5436    2.0000
    0.7046    0.1638    0.9978    3.0000
    0.1720    0.6947    0.2719    4.0000
get_ID(a,0.2178)
%pas de réponse il entre pas dans la condition vérifié en debug

Merci d'avance

[Jean Dumoncel]

Bonjour,

L'affichage de a montre uniquement les 4 premières décimales. (pour changer cela, voir la fonction format).

Tu utilises un test qui a une précision de environ 10^(-16) (c'est le résultat par défaut de la fonction eps).

Donc soit tu connais les 16 premières décimales de a, soit tu utilises une précision plus grande :

Code :

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

1
2
3
4
epsilon = 10^-4;
find(abs(a-0.7711)<epsilon)
epsilon = 10^-5;
find(abs(a-0.7711)<epsilon)

[tenderstoune]

merci pour ta réponse et ton idée

merci encore une fois