Problèmes numériques dans Matlab: même entrée, même code -> sortie différente?

Question

Je rencontre des problèmes lorsque je compare les résultats de différentes exécutions de mon logiciel Matlab avec la même entrée. Pour réduire le problème, je l'ai fait ce qui suit:

• enregistrer toutes les variables pertinentes en utilisant la méthode
• appel d'une méthode save() Matlab qui calcule quelque chose
• enregistrer toutes les variables de sortie pertinentes à l'aide de nouveau save()

Sans changer la méthode appelée, je l'ai fait une autre course avec

• charge les variables enregistrées Abov e et comparer avec les variables d'entrée courantes en utilisant isequal()
• appeler à nouveau ma méthode avec les variables d'entrée en cours
• charger les variables out enregistrées ci-dessus et comparer.

Je ne peux pas croire que la comparaison dans la dernière "ligne" détecte de légères différences. Les calculs comprennent des nombres simples et doubles de précision, l'erreur est de l'ordre de 1e-10 (la sortie est un nombre double). La seule explication possible que je puisse imaginer est que Matlab perde de la précision lors de l'enregistrement des variables (ce que je considère très improbable, j'utilise le format binaire Matlab par défaut) ou qu'il y a des calculs comme a = b + c + d, qui peut être calculé comme = (b + c) + d ou a = b + (c + d), ce qui pourrait conduire à des différences numériques. Savez-vous quelle pourrait être la raison des observations décrites ci-dessus?

Merci beaucoup!

Answer 1

il semble vraiment être causé par le mélange simple/double dans les calculs. Depuis que je suis passé en double précision, le problème ne s'est plus produit. Merci à tous pour vos pensées.

Answer 2

ces erreurs pourraient être arrondies. vous pouvez trouver la précision en virgule flottante de votre système comme ceci:

>> eps('single') 

ans = 

    1.1921e-07 

Sur mon système ce rapports 10^-7 qui expliquerait les écarts de votre commande

Answer 3

Pour obtenir des résultats reproductibles, en particulier si vous utilisez toutes les fonctions de génération aléatoire (directement ou indirectement), vous devez restaurer le même état au début de chaque course:

%# save state (do this once) 
defaultStream = RandStream.getDefaultStream; 
savedState = defaultStream.State; 
save rndStream.mat savedState 

%# load state (do this at before each run) 
load rndStream.mat savedState 
defaultStream = RandStream.getDefaultStream(); 
defaultStream.State = savedState;