将湿球温度数学公式转C#代码时遇运算符错误，求助

问题背景 #

需要转换的湿球温度计算公式：

Tw = T × arctan[0.151977 × (RH% + 8.313659)^(1/2)] + arctan(T + RH%) - arctan(RH% - 1.676331) + 0.00391838 ×(RH%)^(3/2) × arctan(0.023101 × RH%) - 4.686035

写出的C#代码：

double temp = 25;
double rh = 40;

var Tw = (temp * Math.Atan(0.151977 * (rh + 8.313659)) + Math.Atan(temp + rh) - Math.Atan(rh - 1.676331) + 0.00391838 * rh) ^ 3 / 2 * Math.Atan(0.023101 * rh) - 4.686035;

遇到错误：

^ cannot be applied to double and int.

错误原因及修正 #

1. 运算符误用 #

C#中^是按位异或运算符，并非数学中的幂运算符号。计算幂必须使用Math.Pow(base, exponent)方法。

2. 运算优先级错误 #

原代码的括号完全不符合公式逻辑：

• 公式中(RH% + 8.313659)^(1/2)是对rh+8.313659开平方，而非对前面的大段表达式做幂运算
• (RH%)^(3/2)是对rh做3/2次方运算，原代码错误地将前面所有项的和拿来做幂运算

修正后的代码 #

double temp = 25;
double rh = 40;

// 拆分公式为子项，确保运算逻辑与公式完全对齐
double term1 = temp * Math.Atan(0.151977 * Math.Pow(rh + 8.313659, 0.5));
double term2 = Math.Atan(temp + rh);
double term3 = Math.Atan(rh - 1.676331);
double term4 = 0.00391838 * Math.Pow(rh, 1.5) * Math.Atan(0.023101 * rh);
double Tw = term1 + term2 - term3 + term4 - 4.686035;

说明 #

• Math.Pow(rh + 8.313659, 0.5)等价于公式中的(RH% + 8.313659)^(1/2)（平方根）
• Math.Pow(rh, 1.5)等价于(RH%)^(3/2)（3/2次方）
• 拆分子项既保证了运算优先级正确，也提升了代码可读性，便于后续维护和调试

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