具有矩阵变量的数字多项式

Question

p = np.polynomial.Polynomial((0, 0, 1))  # p(x) = x^2
x = np.array(
    [[2, 2],
     [2, 2]]
)
print(p(x))
    [[4. 4.]
     [4. 4.]]
print(x @ x)
    [[8 8]
     [8 8]]

尽管第一个print语句的结果，在numpy中有一种方法可以使用多项式，在矩阵点求值，并获得由原始矩阵乘法产生的结果，如第二个print所示，>

p = np.polynomial.Polynomial（（0，0，1））＃p（x）= x ^ 2 x = np.array（[[2，2]，[2，2]]）print（p （x））[[4。 4.] [4。 4.]] print（x @ x）[[8 8] [8 8]]在numpy中有没有办法...

Answer 1

使用库函数为矩阵x ^ n和整数x计算n：