虽然循环和双重表现奇怪

Question

试图以数​​字方式(第一次)求解微分方程.我的程序要求一堆参数,然后计算将一瓶葡萄酒冷却到给定温度所需的时间.这是我的主要内容:

import java.util.Scanner;

public class Step2_lab11 {
    public static int TAO = 50;
    public static double DELTA_MINUTES = 0.1;

    public static void main(String[] args) {

        System.out.println("VINETS AVKYLNINGSTID \n");
        System.out.println("Ange vinets temperatur:");

        Scanner userIn = new Scanner(System.in);
        double wineTemp = userIn.nextDouble();

        System.out.println("Vinets önskade temperatur:");
        double preferredTemp = userIn.nextDouble();

        System.out.println("Kylens/frysens temperatur:");
        double chillTemp = userIn.nextDouble();

        WineChiller wineChiller = new WineChiller();

        double elapsedTime = 0.0;

        while(wineTemp > preferredTemp) {

            elapsedTime = elapsedTime + DELTA_MINUTES;
            double dT = wineChiller.getChillingTime(TAO, DELTA_MINUTES, chillTemp, preferredTemp, wineTemp);
            wineTemp = wineTemp - dT;
            System.out.println(elapsedTime);


        }
    }

}

这是WineChiller.java文件:

public class WineChiller {
    public WineChiller() {


    }

    public double getChillingTime(int TAO, double DELTA_MINUTES, double chillTemp, double preferredTemp, double wineTemp) {

        double dT = (wineTemp - chillTemp) * DELTA_MINUTES  / TAO;
        return dT;
    }

}

while循环的Syso部分产生这个(wineTemp = 25,preferredTemp = 16,chillTemp = 5)

0.1
0.2
0.30000000000000004
....
29.300000000000146
29.400000000000148
29.50000000000015
29.60000000000015
29.700000000000152
29.800000000000153
29.900000000000155

不知道为什么它会添加随机小数.我也认为(但我不是100%)正确的答案应该是30分钟,而不是29.9分.我在这里错过了一些明显的逻辑错误吗？

Answer 1

你需要阅读这个.

这就是二进制数和IEEE浮点表示的工作原理.

你不能再用二进制表示精确的0.1,而不是十进制的1/3.

这就是为什么在使用double或float时不应该比较值的原因; 你需要对差异的绝对值进行容忍.

看起来像一个简单的Euler积分的一阶ODE,用于集总质量的瞬态热传递.确保您了解时间步长选择如何影响稳定性和准确性.