2

我发现自己想使用作为 mpmath 包的一部分提供的近似值，但对它们应该做什么感到困惑：

http://docs.sympy.org/dev/modules/mpmath/calculus/approximation.html

sympy 表达式和 sympy.mpmath 表达式之间到底有什么区别？

如果我想要一个符号表达式的泰勒近似值而不了解 mpmath 包在做什么，我可以执行以下操作：

#Imports
import sympy
import sympy.parsing
import sympy.parsing.sympy_parser
import Library_TaylorApproximation

#Create a sympy expression to approximate
ExampleStringExpression = 'sin(x)'
ExampleSympyExpression = sympy.parsing.sympy_parser.parse_expr(ExampleStringExpression)


#Create a taylor expantion sympy expression around the point x=0
SympyTaylorApproximation = sympy.series( 
    ExampleSympyExpression,
    sympy.Symbol('x'),
    1, 
    4,
    ).removeO()

#Cast the sympy expressions to python functions which can be evaluated:
VariableNames = [str(var) for var in SympyTaylorApproximation.free_symbols]
PythonFunctionOriginal =  sympy.lambdify(VariableNames, ExampleSympyExpression)
PythonFunctionApproximation = sympy.lambdify(VariableNames, SympyTaylorApproximation)

#Evaluate the approximation and the original at a point:
print PythonFunctionOriginal(2)
print PythonFunctionApproximation(2)

#>>> 0.909297426826
#>>> 0.870987413961

但是，如果我尝试根据文档对 mpmath 做同样的事情：

TaylorCoefficients = sympy.mpmath.taylor(ExampleSympyExpression, 1, 4 )
print 'TaylorCoefficients', TaylorCoefficients

#>>> TypeError: 'sin' object is not callable

我可以尝试在其中填充 python 函数（这是可调用的）：

TaylorCoefficients = sympy.mpmath.taylor(PythonFunctionOriginal, 1, 4 )
print 'TaylorCoefficients', TaylorCoefficients

#>>> TaylorCoefficients [mpf('0.8414709848078965'), mpf('0.0'), mpf('0.0'), mpf('0.0'), mpf('-8.3694689805155739e+57')]

但是以上没有任何意义，因为我知道不能对 python 函数进行导数。

我可以调用 mpmath 函数sin：

TaylorCoefficients = sympy.mpmath.taylor(sympy.mpmath.sin, 1, 4 )
print 'TaylorCoefficients', TaylorCoefficients
#>>> TaylorCoefficients [mpf('0.8414709848078965'), mpf('0.54030230586813977'), mpf('-0.42073549240394825'), mpf('-0.090050384311356632'), mpf('0.035061291033662352')]

但是我也不能按照我想要的方式对其进行操作 -> 就像我想要的那样

SinTimesCos = sympy.mpmath.sin*sympy.mpmath.cos
TaylorCoefficients = sympy.mpmath.taylor(SinTimesCos, 1, 4 )
print 'TaylorCoefficients', TaylorCoefficients
#>>> TypeError: unsupported operand type(s) for *: 'function' and 'function'

究竟什么是 mpmath 函数？

它不是一个 sympy 表达式，也不是一个 python 函数。如何对任意表达式进行操作？

看来我不能在文档中采用任意 sympy 表达式的近似值。 http://docs.sympy.org/dev/modules/mpmath/calculus/approximation.html

如何对任意 sympy 表达式进行任意近似（ Pade / Cheby Chev / Fourier ）？

编辑：

所以我正在寻找的一个例子是以下近似值：

#Start with a sympy expression of (a, b, x)
expressionString = 'cos(a*x)*sin(b*x)*(x**2)'
expressionSympy = sympy.parsing.sympy_parser.parse_expr(expressionString)

#Do not want to decide on value of `a or b` in advance.
#Do want approximation with respect to x:

wantedSympyExpression = SympyChebyChev( expressionSympy, sympy.Symbol('x') )

结果可以是作为a和的函数的系数表达式列表b：

wantedSympyExpressionCoefficients = [ Coef0Expression(a,b), Coef1Expression(a,b), ... , CoefNExpression(a,b)]

或者结果可能是整个 sympy 表达式本身（它本身是a,的函数b）：

wantedSympyExpression = Coef0Expression(a,b) + Coef1Expression(a,b) *(x**2) + ... + CoefNExpression(a,b) (x**N)

请注意a和b不是在执行近似之前选择的。

4

2 回答 2

1

mpmath 函数是普通的 Python 函数。他们只是在任意精度的算术中进行数学运算。

但是以上没有任何意义，因为我知道不能对 python 函数进行导数。

您不能以符号方式取导数，但您可以通过多次评估函数并使用数值微分技术来计算导数的近似值。这就是这样sympy.mpmath.taylor做的。引用文档：

使用高阶数值微分计算系数。该函数必须可以评估为任意精度。

于 2015-11-12T00:56:14.583 回答
1

如果您有一个 SymPy 表达式并希望将其评估为任意精度，请使用evalf，例如

sympy.sin(1).evalf(100)

您可以在评估之前使用sin(x).evalf(100, subs={x:1})替换为x。在后台使用 mpmath，因此这将为您提供与 mpmath 相同的结果，但不必直接使用 mpmath。1evalf

于 2015-11-12T17:24:56.903 回答