因此,我正在编写一个评估泰勒级数的脚本。但是,我希望它能够评估所有类型的功能。因此,例如,我尝试使用功能acot(x)

x = sy.Symbol('x')
f = acot(x)
...
func = taylor(f,0,3)
taylor_lambda = sy.lambdify(x, func, 'numpy')

上面的代码毫无异常(exception)地运行(例如,除非我使用acsch,但它不会运行)。

但是,当它到达这一行时:
plt.plot(x1,taylor_lambda(x1),label='taylor approximation')

我得到:
NameError: name 'acot' is not defined

我试图在lambdify调用中将numpy替换为sympy,但这似乎具有象征意义。某些情况(更罕见的功能)正在发生这种情况,而其他情况则没有。谢谢!

我的进口如下:
import sympy as sy
import numpy as np
from sympy.functions import *
from sympy import pi, E,acot
import matplotlib.pyplot as plt
import math

最佳答案

我之前确实有过类似的问题,并且已经设法解决了。

plt.plot(x1,taylor_lambda(x1),label='taylor approximation')

看起来还可以。我给我一个较旧的代码,可以正常工作,您可以进行比较。
from sympy.abc import x
from sympy import sin, series
from sympy.utilities.lambdify import lambdify

import numpy as np
import matplotlib.pyplot as plt


func = sin(x)/x
taylor = series(func, n=6).removeO()

evalfunc = lambdify(x, func, modules=['numpy'])
evaltaylor = lambdify(x, taylor, modules=['numpy'])

t = np.linspace(-7.5, 7.5 , 100)
plt.plot(t, evalfunc(t), 'b', label='sin(x)/x')
plt.plot(t, evaltaylor(t), 'r', label='Taylor')
plt.legend(loc='best')
plt.show()

关于python - (某些功能)未使用SymPy Lambdify定义,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/45273827/

10-12 21:22