BigDecimal实现浮点数精确运算

BigDecimal实现浮点数精确运算

1.精度缺失

1.1 加减乘除

​ Java在对浮点数进行四则运算时，会出现明显的精度缺失问题。

1.2 具体原因

​ 我们发现，计算出来的值和我们预期结果不一致。原因在于我们的计算机是二进制的。浮点数没有办法使用二进制进行精确表示。

​ 计算机的 CPU 表示浮点数由两个部分组成：指数和尾数，这样的表示方法一般都会失去一定的精确度，有些浮点数运算也会产生一定的误差。

​ 浮点运算很少是精确的，只要是超过精度能表示的范围就会产生误差。往往产生误差不是因为数的大小，而是因为数的精度。因此，产生的结果接近但不等于想要的结果。

2.BigDecimal

2.1 API简介

​ Java在java.math包中提供的API类BigDecimal，用来对超过16位有效位的数进行精确的运算。双精度浮点型变量double可以处理16位有效数。在实际应用中，需要对更大或者更小的数进行运算和处理。float和double只能用来做科学计算或者是工程计算，在商业计算中要用java.math.BigDecimal。BigDecimal所创建的是对象，我们不能使用传统的+、-、*、/等算术运算符直接对其对象进行数学运算，而必须调用其相对应的方法。方法中的参数也必须是BigDecimal的对象。构造器是类的特殊方法，专门用来创建对象，特别是带有参数的对象。

2.2 DFormat

​ 以下就是我们依据BigDecimal实现的能精确浮点数的四则运算development工具类

package Unit;

import java.math.BigDecimal;
import java.math.RoundingMode;

public class DFormat {
    //相减
    public static Double sub(Double v1, Double v2) {
        BigDecimal b1 = new BigDecimal(v1.toString());
        BigDecimal b2 = new BigDecimal(v2.toString());
        return b1.subtract(b2).doubleValue();
    }

    //相加
    public static Double add(Double v1, Double v2) {
        BigDecimal b1 = new BigDecimal(v1.toString());
        BigDecimal b2 = new BigDecimal(v2.toString());
        return b1.add(b2).doubleValue();
    }

    //相乘
    public static Double mul(Double v1, Double v2) {
        BigDecimal b1 = new BigDecimal(v1.toString());
        BigDecimal b2 = new BigDecimal(v2.toString());
        return b1.multiply(b2).setScale(6,RoundingMode.HALF_UP).doubleValue();
    }
    //相除
    public static Double div(Double v1, Double v2) {
        BigDecimal b1 = new BigDecimal(v1.toString());
        BigDecimal b2 = new BigDecimal(v2.toString());
        return b1.divide(b2,6,RoundingMode.HALF_UP).doubleValue();
    }
}

3.取舍模式

3.1 模式介绍

RoundingMode.CEILING：取右边最近的整数

RoundingMode.DOWN：去掉小数部分取整，也就是正数取左边，负数取右边，相当于向原点靠近的方向取整

RoundingMode.FLOOR：取左边最近的正数

RoundingMode.HALF_DOWN:五舍六入，负数先取绝对值再五舍六入再负数

RoundingMode.HALF_UP:四舍五入，负数原理同上

RoundingMode.HALF_EVEN:这个比较绕，整数位若是奇数则四舍五入，若是偶数则五舍六入

3.2 模式选择

​ 经过多次的计算与对比测试，我们选择了RoundingMode.HALF_UP取舍模式，并保留了小数点后六位的数字。这种计算模式是与计算器计算结果最近似的情况，并且计算的误差控制在10^-3~10^-6之内可以近似忽略。在下面的第四部分会具体体现这一点。

4.计算测试

4.1 计算结果

4.1.1 运算参数

g/kg, g/kg

,

, ，，

，

4.1.2 运算公式

公式1.1

样例（1步运算）：

计算器：

程序：

公式1.2

样例（5步运算）：

计算器：

程序：

公式2

样例（3步运算）：

计算器：

程序：

公式3

样例（5步运算）：

计算器：

程序：