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

您需要阅读this。

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

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

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

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

虽然循环和双重表现奇怪

1 个答案: