MATLAB设置vpa100位精度后fprintf与double仅输出16位问题咨询

你给出的测试代码如下：

clear;
clc;
syms x;
f=log(x)*sin(x^2);
a=vpa(subs(f,x,2),100)
fprintf('a=%.100f\n',a);
doublea=double(a);
fprintf('a=%.100f\n',doublea);

以下是对应问题的解答：

问题1：已通过vpa设置100位精度，为何使用fprintf输出时仅能保留16位有效小数？ #

MATLAB原生的fprintf函数不支持直接输出vpa/sym类型的高精度数值，当你把vpa类型变量传入fprintf时，MATLAB会先隐式将其转换为double类型，再执行格式化输出。double类型的精度上限只有15~17位十进制有效数字，因此即使你指定了%.100f的格式，16位之后的数值都是无效的补位0，且前16位之后的精度已经丢失。

问题2：将vpa生成的变量a转换为double类型后，也仅保留16位有效数字的原因是什么？ #

double类型是遵循IEEE 754标准的双精度浮点数，固定占用64位存储空间，其中用于存储有效数字的尾数部分仅占52位，换算为十进制有效数字的上限就是15~17位。超过这个范围的数值无法被double类型存储，因此无论vpa变量原始精度是多少，转换为double后都会被截断，最多保留17位有效数字。

问题3：若需使用超过16位精度的计算，该情况是否会导致计算误差？应当如何解决？ #

会产生计算误差。如果在高精度计算流程中中途将vpa类型转换为double，或者使用了仅支持double类型的运算函数，精度会被强制截断到double的上限，和目标高精度结果之间会出现偏差，多步运算后误差还会累计放大。
解决方法如下：

• 全程使用sym或vpa类型完成所有计算流程，不要中途转换为double类型
• 输出高精度结果时不要使用fprintf，直接用disp输出vpa变量，或先通过char()函数将vpa变量转换为字符串后再自定义格式输出，示例代码：

  % 输出完整100位精度结果
  disp(['a=' char(a)])
  

• 计算时优先使用符号数学工具箱中支持vpa类型的运算函数，避免调用仅支持double类型的原生数值函数。

内容的提问来源于stack exchange，提问作者Ongky Denny Wijaya