Question

在R中如何表示多项式表达式并使用数字向量对象进行多项式数学运算？例如：

x1 <- c(2,1)  # 2 + x
x2 <- c(-1,3)  # -1 + 3*x

想要：

x1 * x2 # to return -2 + 5*x + 3*x^2

注意：我今天早上回答了一个问题，然后海报明显将其删除了（让我想知道这是否是家庭作业。）所以我从记忆中重新发布了这个问题。

Answer 1

可以使用outer直接乘以系数，然后汇总结果

x1 <- c(2,1)  # 2 + x
x2 <- c(-1,3)  # -1 + 3*x
tmp <- outer(x1, x2)
tapply(tmp, row(tmp) + col(tmp) - 1, sum)
# 1  2  3 
#-2  5  3

x1 <- c(2, 1) # 2 + x
x2 <- c(-1, 3, 2) # -1 + 3*x + 2*x^2
tmp <- outer(x1, x2)
tapply(tmp, row(tmp) + col(tmp) - 1, sum) # should give -2 + 5*x + 7*x^2 + 2*x^3
# 1  2  3  4 
#-2  5  7  2

如评论中所述，代码中的“-1”不是必需的。在提出帮助我的解决方案时，因为它允许我将outer输出中的每个位置映射到最终向量中的最终位置。如果我们改为'-2'，那么它将映射到结果多项式中x的指数。但我们真的不需要它，所以像下面这样的东西也会起作用：

tmp <- outer(x1, x2)
tapply(tmp, row(tmp) + col(tmp), sum)

Answer 2

使用多项包：

 require(polynom)
# Loading required package: polynom
# From the example for as.polynomial
 p <- as.polynomial(c(1,0,3,0))
 p
# 1 + 3*x^2 

 x1 <- c(2,1)
 x2 <- c(-1,3)
 px1 <- as.polynomial(x1)
 px2 <- as.polynomial(x2)

 px1*px2
# -2 + 5*x + 3*x^2 
 prod.p <- .Last.value
 str(prod.p)
# Class 'polynomial'  num [1:3] -2 5 3
 unclass(prod.p)
# [1] -2  5  3

Answer 3

多项式乘法是系数的卷积

convolve(c(2,1),rev(c(-1,3)),type="open")
#[1] -2  5  3

如何用R中的数字向量表示多项式

3 个答案: