C#与JavaScript湿球温度计算结果差异的原因排查

背景 #

在气象程序CumulusMX中，湿球温度的C#计算实现如下：

public static double CalculateWetBulbC(double tempC, double dewPointC, double pressureMb)
{
    double svpDP = SaturationVapourPressure1980(dewPointC);

    return (((0.00066 * pressureMb) * tempC) + ((4098 * svpDP) / (Sqr(dewPointC + 237.7)) * dewPointC)) / ((0.00066 * pressureMb) + (4098 * svpDP) / (Sqr(dewPointC + 237.7)));
}

public static double SaturationVapourPressure1980(double tempC)
{
    return 6.112 * Math.Exp(17.67 * tempC / (tempC + 243.5));
}

注：Sqr(a)表示计算参数的平方（即a*a）

对应的JavaScript计算实现为：

(0.00066 * Pressure[i][1] * Temperature[i][1] + 4098 * 6.112 * Math.exp(17.67 * Temperature[i][1] / (243.5 + Temperature[i][1])) 
       / Math.pow((Dewpoint[i][1] + 237.7), 2) * Dewpoint[i][1]) / (0.00066 * Pressure[i][1] + 4098 * 6.112 
       * Math.exp(17.67 * Temperature[i][1] / (243.5 + Temperature[i][1])) / Math.pow((Dewpoint[i][1] + 237.7), 2))

注：JavaScript将计算整合为单一表达式，[i][1]指向JSON数组中过去三天内某一分钟的数据，实现逐分钟计算；C#同样按该频率计算

目前观察到两者计算结果存在差异：

• wetbulb_1minwetbulb：CumulusMX的C#计算结果，每分钟一次
• Natte Bol Temperatuur：JavaScript计算结果，每分钟一次

经确认，C#与JavaScript的计算公式逻辑一致，输入数据完全相同。现需分析结果差异的产生机制及可能原因。

可能的原因 #

1. 浮点数精度差异 #

C#的double和JavaScript的Number（基于IEEE 754双精度浮点数）虽遵循同一标准，但运算细节可能存在细微偏差：

• 指数运算Math.Exp/Math.exp的底层实现不同，导致中间结果的精度差异
• 除法、乘法的运算顺序解析时，可能产生隐式的精度损失差别
• 平方运算：C#用a*a，JS用Math.pow(a,2)，极端数值下的精度表现可能不同

2. 变量类型与数据截断 #

• 输入数据的读取方式有差异：比如C#读取时可能做了类型转换截断，JS直接读取为浮点数，导致输入值看似一致但实际存在微小差别
• 数据存储/传输时的精度丢失，导致两边实际使用的输入值存在细微偏差

3. 运算顺序的隐式差异 #

虽然公式逻辑一致，但表达式的运算优先级在不同语言中可能存在易被忽略的差异：

• 乘法、除法的结合顺序不同，导致中间计算步骤的数值结果偏差
• JS表达式中重复计算了Math.exp(17.67 * Temperature[i][1] / (243.5 + Temperature[i][1]))，而C#通过svpDP只计算一次，重复计算时的精度偏差会传递到最终结果

4. 边界值处理差异 #

当输入值处于极端范围（如极低/极高温度、露点）时，两种语言的数学函数可能有不同的边界处理逻辑，导致计算结果偏差

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