Matlab qui "change mes valeurs", problème de format?

Version imprimable

19/01/2010, 10h25
sylvain62

Matlab qui "change mes valeurs", problème de format?
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. 8O

Concrètement voilà mon code :

Code:

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:

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:

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:

>> a(180)

j'obtiens:
Code:

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:

1 2 3 4 5 6 >> format short >> a(180) ans = -74.2000
Mais après un test:
Code:

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.
19/01/2010, 10h39
Caro-Line

Matlab 5.2 8O

Réponse dans la :faq:
:arrow: 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.
19/01/2010, 10h49
sylvain62

Cool Merci Caro-Line

Dans la FAQ...:oops: 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:lol: 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 :)

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

Bonjour.
Une autre solution ?

Code:

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

19/01/2010, 19h15
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:

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 >>
19/01/2010, 19h22
Jerome Briot
Citation:
Envoyé par kmaniche

Code:

1 2 >> a=linspace(-110, 70, 900) ;
Euh... plutôt :

Code:

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

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

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
;)
19/01/2010, 21h04
kmaniche

On peux régler avec Format
19/01/2010, 22h53
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.
20/01/2010, 14h44
kmaniche

Citation:

Envoyé par Dut

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

Exacte.
29/01/2010, 11h55
sylvain62
Pour revenir à la solution de Phryte
Citation:
Bonjour.
Une autre solution ?

Code:

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:

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

Citation:

justes à une précision donnée

En tous cas merci à tous pour vos réponses