Discussion :

[sylvain62]

Bonjour à tous,

Je poste ici, un problème ou plutôt une question.

Je demande à matlab (version 7.6.0.324 mais pareil sur Matlab version 5.2.0.3084), de me créer une matrice d'une seule colonne et de la remplir de valeurs. Je lui indique (il me semble) les valeurs que je veux.

Pourtant, certaines des valeurs que je retrouve dans cette matrice colonne, diffèrent de celles que j'avais demandées.

Concrètement voilà mon code :

Code :

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

a=(-110:0.2:70)';

Je créer un vecteur que je veux remplir (ici de longitudes) allant de -110 à 70 avec un incrément de 0.2 degrés. Matlab me le crée

Code :

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

1
2
3
4
>> whos
  Name        Size            Bytes  Class     Attributes
 
  a         901x1              7208  double

A mon sens, les valeurs inclues dans "a" seraient par exemple:
-110
-109.8
-109.6
...

Ce qu'elles sont en fait pour une grosse partie:

Code :

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

1
2
3
4
5
6
7
8
9
a(1:3)
 
ans =
 
  1.0e+002 *
 
  -1.100000000000000
  -1.098000000000000
  -1.096000000000000

Mais pas pour toute.Si je tape:

Code :

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

>> a(180)

j'obtiens:

Code :

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

1
2
3
ans =
 
 -74.199999999999989

Au lieu de -74.2, et ceci que ce soit le format long ou format short. Disons qu'en format short j'ai un chiffre arrondi, donc je ne m'en étais pas rendu compte tout de suite

Code :

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

1
2
3
4
5
6
>> format short
>> a(180)
 
ans =
 
  -74.2000

Mais après un test:

Code :

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

1
2
3
4
5
a(180)==-74.2
 
ans =
 
     0

Il semblerait que ma valeur, que je voulais égale à -74.2 ne l'est pas...

Quelqu'un a une idée de pourquoi Matlab fait ça?
J'ai cherché dans le forum, du côté des "arrondis" mais je n'ai pas vraiment trouvé d'éléments de réponse.

[Caro-Line]

Matlab 5.2

Réponse dans la
Pourquoi 0.3-0.2-0.1 est-il différent de 0 ?

Même si le titre parait loin de ta question

Et MATLAB n'y est pour rien le pauvre.

[sylvain62]

Cool Merci Caro-Line

Dans la FAQ... Pourtant j'y passe régulièrement, mais j'avoue "virgule flottante", n'avait pas attiré mon attention.

C'est pile les soucis quand on lit dans la FAQ:
"Il est nécessaire d'être particulièrement vigilant lorsque des tests logiques d'égalité sont effectués sur ces nombres.

A ma décharge, hier soir (je me suis pris la tête avec des test logiques pendant une partie de la soirée), je n'avais pas internet Donc pas de FAQ (hum hum je m'en sortirai comme ça...)

Donc si je peux redonner un conseil (déjà donner par Caro-Line, Dut et tant d'autres), aller toujours toujours voire la FAQ


et oui je travaille toujours de temps à autre avec Matlab 5...
Encore Merci Caro

[phryte]

Bonjour.
Une autre solution ?

Code :

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

1
2
3
4
5
6
7
clear
a=(-110:2e-1:70)';
a1=(floor((-110:2e-1:70)*100)/100)';
a(180)
a1(180)
a(180)+74.2
a1(180)+74.2

[kmaniche]

Bonsoir,

Pour enrichir d'avantage les solutions proposées.

a = -110 à 70 avec un pas de 0.2 => 70-(-110) = 180

le nombre d'éléments du vecteur a est de : 180/0.2 = 900

Ma solution :

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
 
>> a=linspace(-110, 70, 900) ;
>> numel(a)
 
ans =
 
   900
 
>> max(a)
 
ans =
 
    70
 
>> min(a)
 
ans =
 
  -110
 
>>

[Jerome Briot]

Code :

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

1
2
 
>> a=linspace(-110, 70, 900) ;

Euh... plutôt :

Code :

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

a=linspace(-110, 70, 901) ;

Mais ce qui ne change (presque) rien au problème :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
>> a=linspace(-110, 70, 901) ;
>> a(180)
 
ans =
 
  -74.2000
 
>> sprintf('%.15f',a(180))
 
ans =
 
-74.200000000000003

[kmaniche]

On peux régler avec Format

[Jerome Briot]

Ce n'est pas uniquement une question d'affichage...

L'outil informatique a des limites qu'il est bon d'avoir en tête lorsqu'on programme

Maintenant, ce n'est pas non plus une raison pour conclure que tous les calculs sont faux... ils sont simplement justes à une précision donnée.

[kmaniche]

L'outil informatique a des limites qu'il est bon d'avoir en tête lorsqu'on programme

Exacte.

[sylvain62]

Pour revenir à la solution de Phryte

Bonjour.
Une autre solution ?

Code :

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

1
2
3
4
5
6
7
8
 
clear
a=(-110:2e-1:70)';
a1=(floor((-110:2e-1:70)*100)/100)';
a(180)
a1(180)
a(180)+74.2
a1(180)+74.2

Effectivement ça fonctionne pour 74.2 mais pas pour tous, par exemple pour 74.4:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
a(179)-74.4
 
ans =
 
   -1.488000000000000e+002
 
a1(179)-74.4
 
ans =
 
   -1.488100000000000e+002

Donc pour mes tests logiques ça marche pas, je passe par des "valeurs proches" comme expliquées dans la FAQ. Et garder en tête cette limitation.
Comme disait Dut

justes à une précision donnée

En tous cas merci à tous pour vos réponses