如何访问numpy多项式基函数

dldeef67 于 2023-10-19 发布在其他

关注(0)|答案(2)|浏览(121)

我不知道如何用numpy.polynomial.legendre访问勒让德基函数L_i（x）。我只能得到它们的加权和，和c_i，L_i，其中c是常系数。如何访问底层的基函数？我一定是漏掉了什么。我在下面举一个简单的例子。我手工编写了正确的多项式，并将它们与我使用Numpy访问的内容进行比较。显然，加权和不是期望的单个基函数。如何使用numpy的polynomial模块使顶部图与底部图匹配？

import numpy as np
import matplotlib.pyplot as plt
f, ax = plt.subplots(2)
x = np.linspace(-1, 1, 1000)
def correct_legs(x, c):
   if c == 0: return(np.ones(len(x)))
   if c == 1: return(x)
   if c == 2: return(0.5 * (3 * x ** 2 - 1))
   if c == 3: return(0.5 * (5 * x ** 3 - 3 * x))
   if c == 4: return((1. / 8.) * (35 * x ** 4 - 30 * x ** 2 + 3))
   if c > 4 or c < 0 or type(c)!=int: return(np.nan)
for order in range(5):
   # Attempt to generate legendre polynomials with numpy
   ax[0].plot(x, np.polynomial.legendre.legval(x, np.ones(order + 1)), label=order)

   # plot propper legendre polynomials
   ax[1].plot(x, correct_legs(x, order))
ax[0].legend()
plt.savefig('simpleExample')

numpy

来源：https://stackoverflow.com/questions/57829124/how-to-access-numpy-legendre-polynomial-basis-functions

2条答案

按热度按时间

fiei3ece1#

当然，所需的解决方案取决于您希望如何使用多项式。为了“绘图”的目的，你可以看看Legendre系列类。它可以通过Legendre.basis方法生成所需域中的基函数。
对于绘图，它工作正常

import numpy as np
import matplotlib.pyplot as plt

for order in range(5):
   # Attempt to generate legendre polynomials with numpy
   x, y = np.polynomial.legendre.Legendre.basis(order, [-1, 1]).linspace(100)
   plt.plot(x, y, label=order)

plt.legend()
plt.plot()

赞(0）回复(0）举报 2023-10-19

vshtjzan2#

这是一种方法。

from numpy import cos, arange, pi
from scipy.special import legendre
import matplotlib.pyplot as plt

def P(l,chi ) :
  polynom = legendre(l)
  return polynom(chi)

CHI = arange(-1,1, 0.01)

plt.plot(CHI, P(1,CHI), lw=2, ls = '-', color = 'k')
plt.plot(CHI, P(2,CHI), lw=2, ls = '--', color = 'g')
plt.plot(CHI, P(3,CHI), lw=2, ls = '-.', color = 'b')
plt.plot(CHI, P(4,CHI), lw=2, ls = ':', color = 'r')

plt.legend(['$P_1(\chi)$','$P_2(\chi)$','$P_3(\chi)$','$P_4(\chi)$'])
plt.title(' Polynômes de Legendre ')
plt.xlabel(' $\chi = \cos\phi$ ')
plt.ylabel(' $P_l(\chi)$ ')

这是输出Output

赞(0）回复(0）举报 2023-10-19

我来回答

如何访问numpy多项式基函数

2条答案

相关问题

热门标签

最新问答