Matlab任意精度工具箱使用困惑及128位精度实现咨询

精度参数与存储位的对应关系 #

浮点数的十进制有效数字位数和二进制存储位存在固定转换关系：

• 十进制有效数字位数 ≈ 二进制尾数位数 × log₁₀(2) ≈ 二进制尾数位数 × 0.3010
• 64位双精度：尾数52位，对应约15-17位十进制有效数字
• 128位四精度：尾数112位，对应约33-36位十进制有效数字

工具箱的precision参数（十进制位数）本质是指定有效数字位数，设置为34左右即可匹配128位四精度的运算精度。

HPF工具箱快速示例 #

% 设置精度为34位十进制有效数字（对应128位四精度）
hpf('precision', 34);

% 创建128位精度复数
z_hpf = hpf('1.234567890123456789012345678901234') + hpf('9.876543210987654321098765432109876')*1i;

% 执行递归运算核心操作（exp、+、*）
result_hpf = exp(z_hpf) + z_hpf * hpf('2.718281828459045235360287471352662');
disp(result_hpf);

MPT工具箱快速示例 #

% 设置34位十进制有效数字，匹配128位四精度
mp.Digits(34);

% 创建128位精度复数
z_mpt = mp('1.234567890123456789012345678901234') + mp('9.876543210987654321098765432109876')*1i;

% 执行核心运算
result_mpt = exp(z_mpt) + z_mpt * mp('2.718281828459045235360287471352662');
disp(result_mpt);

关键说明 #

• 不要混淆“小数点后位数”和“有效数字位数”：工具箱的precision指的是从第一个非零数字开始的十进制有效数字总数，而非小数点后的位数。
• 128位四精度的核心是二进制存储的尾数长度（112位），通过设置34位十进制有效数字，工具箱会自动分配对应的二进制存储资源，实现128位精度的运算。

内容的提问来源于stack exchange，提问作者CitizenInsane