在R中编写一个函数来求解Plank方程

2023-05-17 admin

我正试着用R写出我的第一个函数，用普朗克函数计算不同温度下的发射度。我可以像下面这样手动完成，温度从200到310K。

pi <- 3.141593
h <- 6.626068963e-34
c <- 2.99792458e+8
lambda <- 4 * 1e-6
k <- 1.38e-23
t <- c (200:310)
a <- (2*pi*(c^2)*h)/(lambda^5)
b <- exp((h*c)/(lambda*k*t))
B <- a * (1/(b-1))

其中B是我想要的值的向量。

下面是用R编写函数的一种尝试：

P_function <- function(t, pi = 3.141593, h = 6.626068963e-34, c = 2.99792458e+8,
             lambda = 4 * 1e-6, k = 1.38e-2) {
         ((2*pi*(c^2)*h)/(lambda^5)) *((1/(exp((h*c)/(lambda*k*t))-1)))
}

现在对于不同的t值(200-300K)，我该如何实现这个函数？

转载请注明出处:http://www.nali5.com/article/20230517/1370247.html