Question

我正在将一些代码从Matlab移植到C ++。

在Matlab中

format long
D = 0.689655172413793 (this is 1.0 / 1.45)
E = 2600 / D
// I get E = 3.770000000000e+03

在C ++中

double D = 0.68965517241379315; //(this is 1.0 / 1.45)
double E = 2600 / D;
//I get E = 3769.9999999999995

对我来说这是一个问题，因为在这两种情况下我都必须向下舍入到0（Matlab的修复），在第一种情况下（Matlab）变为3770，而在第二种情况下（C ++）它变为3769。

我意识到这是因为C ++案例中有两个额外的最低有效数字“15”。鉴于Matlab似乎只能在双精度中存储多达15位有效精度数字（如上图所示 - 0.689655172413793），我怎样才能有效地告诉C ++忽略后面的“15”？

所有计算均以双精度完成。

Answer 1

你对C ++和MATLAB打印双值的不同方式感到困惑。 MATLAB的format long仅打印15 significant digits，而C ++打印17 significant digits。在内部都使用相同的数字：IEEE 754 64位浮点数。为了在MATLAB中重现C ++行为，我定义了一个anonymous function disp17，它打印出17位有效数字的数字：

>> disp17=@(x)(disp(num2str(x,17)))

disp17 = 

    @(x)(disp(num2str(x,17)))

>> 1.0 / 1.45

ans =

   0.689655172413793

>> disp17(1.0 / 1.45)
0.68965517241379315

你看到MATLAB和C ++中的结果是一样的，它们只是打印不同数量的数字。如果现在使用相同的常量继续使用两种编程语言，则会得到相同的结果。

>> D = 0.68965517241379315 %17 digits, enough to represent a double.

D =

   0.689655172413793

>> ans = 2600 / D %Result looks wrong

ans =

     3.770000000000000e+03

>> disp17(2600 / D) %But displaying 17 digits it is the same.
3769.9999999999995

打印17或15位数的背景：

A double requires 17 significant digits to be stored without precision loss，这是C打印的内容。
对于最多15位数字，任何数字都可以从字符串转换为double到字符串，并返回原始数字，这就是MATLAB所做的。

Answer 2

如果我理解您要实现的目标，使用ceil功能可能有所帮助：

ans = ceil(ans); /* smallest integral value that is not less than ans. */
// now ans in C++ is also be 3970.

这是usage reference。

Matlab与C ++双精度

2 个答案: