如何將正交多項式分配給函數？

Question

我正在嘗試 plot Laguerre 的 class 分配的正交多項式。 我想創建 10 個函數，每個 function 通過索引 i 分配給一個特定的多項式。

for (i in 1:10){
  opl[i] <- function(x) {opl[3]}
}

然后使用curve()到 plot 它。 但它不起作用。 laguerre.polynomials() function 將多項式作為列表提供給您，我認為問題在於我的循環無法通過索引從列表中提取項目並將其分配給 function。

關於如何做的任何想法？

Answer 1

您可以使用as.function將多項式轉換為函數，例如通過

library(orthopolynom)
library(ggplot2)

opl <- laguerre.polynomials(10)

opl_functions <- lapply(opl, as.function)

# x interval
x <- seq(-1, 1, 0.05)

# plot the first two polynomials
ggplot(data.frame(x), aes(x = x, y = y)) +      # basic graphical object
    geom_line(aes(y = opl_functions[[1]](x)), colour = "red") +  # first layer
    geom_line(aes(y = opl_functions[[2]](x)), colour = "blue") # second layer

# and so on ...

opl_functions的第 i 個元素是第 i 個多項式，具體取決於x 。 然后可以使用它來 plot 多項式。