如何在 R 中制作这个函数？

Question

我目前正在尝试编写一个函数，通过 E = e21.44 * (e(8.74901^10−20 / (1.3806488*10−23 *T)。等式中的 T 等于温度在开尔文。我知道如果我要插入 20ºC 的温度，该函数应该给我大约 0.829 的输出。但是，当我运行该函数（见下文）时，我得到 2.628622e-121。我相信这是方法我在说明如何将温度从摄氏温度转换为开尔文温度是不正确的，我目前坚持使用它。

TempKelvin <- 273.15
Energy = function(TempKelvin = 273.15 {
exp(24.11) * exp(-8.74901*10^20 / (1.3806488*10^-23 * (T + TempKelvin)))
}
Energy(T = 28)

Answer 1

有几个问题：

• 对Energy函数的调用使用T参数，但该函数（如问题中所定义）没有T参数
• T在 R 中也表示TRUE 。使用不同的名称。 除非被覆盖，否则在数字上下文中使用时T将被解释为 1。 例如，在 R 中，默认情况下T+8等于9 。
• 问题中的公式使用 8.74901 但代码使用该值的负数。 哪个是正确的？
• 对 R 常数使用适当的指数表示法，例如 1.38064852e−23
• 玻尔兹曼常数的一些数字（乘以温度的数字）是不正确的。
• 尽管添加 273.15 是正确的摄氏到开尔文转换公式，但我注意到如果我们使用 273 而不是更准确的 273.15，那么我们确实得到了假设的答案。

请注意，因为答案仅取决于celsius和 273.15 的总和，所以我们可以获得使用 273 而不是 273.15 的效果，如下所示：

celsius + 273 
= (celsius - 0.15) + 273.15

也就是说，传递函数 20-0.15 而不是 20 有效地使用 20 和 273 而不是使用 20 和 273.15 。

Energy <- function(celsius) {
  exp(21.44) * exp(-8.74901e−20 / (1.38064852e−23 * (celsius + 273.15)))
}

# match assumed answer (despite assumed answer being slightly wrong)
Energy(20 - 0.15)  
## [1] 0.8289531

# correct answer
Energy(20)
## [1] 0.8381777

Answer 2

欢迎来到 SO！ 如果给你摄氏度，那么先转换成开尔文，然后再计算公式如何，然后你就可以从那里缩小问题的范围。 我猜-8.7..与8.7..存在括号问题或符号问题。

Energy <- function(temp_c) {
  temp_k <- temp_c + 273.15
  energy <- exp(21.44) * exp(-8.74901*10^20 / (1.3806488*10^-2 * temp_k))
  print(energy)
  }
Energy(28)
#[1] 0