python - 使用 Numpy 创建勒让德公式

因此，我尝试使用 numpy.polynomial.legendre 命令生成 P2 到 Pn 多项式公式。
我想输入 2 它给我 p2 = 1/2 *(-1 +3x**2) 或者如果输入是 3 它给我 P3 公式。

这样我就可以给出 x 值来计算每个 Pn 并使用我的一些类方法来计算进一步的误差来找到根。

我设法使用以下方法制作情节:

 numpy.polynomial.legendre.legval (x, np.identity(10))

最佳答案

我认为您正在寻找函数 scipy.special.legendre 。

#Build the polynomial
>>> import scipy.special as sp
>>> sp.legendre(2)
poly1d([ 1.5,  0. , -0.5])

#Compute on an interval from -1 to 1
>>> sp.legendre(2)(np.linspace(-1,1,10))
array([ 1.        ,  0.40740741, -0.03703704, -0.33333333, -0.48148148,
       -0.48148148, -0.33333333, -0.03703704,  0.40740741,  1.        ])

关于python - 使用 Numpy 创建勒让德公式，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/23419402/