Dans MATLAB, sont des variables VRAIMENT en double précision par défaut?

Ils sont doubles. vpa() est tout simplement choisir d'afficher des chiffres non significatifs au-delà de la virgule flottante précision relative, où printf() et disp() tronquer ou zéro dehors.

Vous êtes seulement obtenir l'original de votre quatre chiffres de retour parce que le littéral que vous avez choisi d'initialiser num avec seulement se trouve être l'exact décimal à l'expansion d'un binaire double de la valeur, car il a été copier-collé à partir de la sortie de l'extension d'une réelle valeur double de l'autre question. Il ne fonctionnera pas pour la proximité d'autres valeurs, comme vous le montrez dans votre "BONUS" additif.

Plus précisément, tous les littéraux numériques dans Matlab produire des valeurs de type double. Ils se converti à l'binaire double de la valeur qui est la plus proche à la valeur décimale de la valeur qu'ils représentent. En effet, les chiffres dans un littéral au-delà de la limite de précision du type double sont ignorés en mode silencieux. Lorsque vous copiez et collez la sortie de vpa pour créer une nouvelle variable, comme les questions des autres affiches ont fait avec le e = ... déclaration, vous êtes à l'initialisation d'une valeur à partir d'un littéral, au lieu de traiter directement avec le résultat d'une expression précédente.

Les différences ici sont juste en sortie de mise en forme. Je pense que ce qui se passe est que vpa() est de prendre de la double précision binaire double et la traiter comme une valeur exacte. Pour une donnée binaire mantisse-valeur d'exposant, vous pouvez calculer l'équivalent décimal de façon arbitraire, nombre de décimales. Si vous avez une précision limitée ("largeur") dans la valeur binaire, comme vous le faites avec n'importe quel de taille fixe, type de données, que beaucoup de ces chiffres décimaux sont importantes. printf() et Matlab par défaut de l'affichage de gérer cela par la troncation de la sortie ou de l'affichage des chiffres non significatifs à 0. vpa() est ignorant les limites de la précision et de continuer à calculer autant de décimales que vous demandez.

Ces chiffres sont faux, dans le sens que si elles ont été remplacées par d'autres valeurs pour produire de la proximité d'une valeur décimale, ils seraient tous get "arrondi" à la même binaire valeur double.

Voici une manière de le montrer. Ces valeurs de x sont tous les mêmes lorsqu'ils sont stockés dans des chambres doubles, et seront tous représentés de la même façon par vpa().

x = [
    2.7182818284590455348848081484902650117874145507812500
    2.7182818284590455348848081484902650117874145507819999
    2.7182818284590455348848
    2.71828182845904553488485555555555555555555555555555
    exp(1)
    ]
unique(x)

Voici une autre façon de démontrer qu'on. Voici deux doubles qui sont très proches les uns des autres.

x0 = exp(1)
x1 = x0 + eps(x0)

vpa(x0) et vpa(x1) devrait produire des résultats qui diffèrent beaucoup passé le 16e chiffres de. Cependant, vous ne devriez pas être en mesure de créer un double de la valeur x tels que vpa(x) produit une représentation décimale qui tombe entre vpa(x0) et vpa(x1).

(Mise à JOUR: Amro points que vous pouvez utiliser fprintf('%bx\n', x) pour afficher une représentation exacte du sous-jacent valeur binaire au format hex. Vous pouvez l'utiliser pour confirmer les littéraux de la carte à la même double.)

Je soupçonne vpa() se comporte de cette façon, parce qu'il traite ses entrées comme des valeurs exactes, et polymorphically prend en charge d'autres Matlab types de la Symbolique boîte à outils qui ont plus de précision que doublé. Ces valeurs doivent être initialisés par des moyens autres que les littéraux numériques, c'est pourquoi sym() prend une chaîne de caractères en entrée et vpa(exp(1)) diffère de vpa(sym('exp(1)')).

Sens? Désolé pour le long-windedness.

(Remarque je n'ai pas la Symbolique boîte à outils, donc je ne peux pas tester vpa() moi-même.)