How to isolate the exponent in Sympy（如何在 Sympy 中“匹配”大公式）

Question

我正在尝试以 F*exp(E) 的形式提取类似于高斯 PDF 的公式中的指数。 SymPy 中的函数 match() 适用于简单的情况，但不适用于更复杂的公式，如下所示。

from __future__ import division
from sympy import *

δ, b, σ, m, s, μ, x = symbols('δ b σ m s μ x', real=True)

gauss = 1/σ/sqrt(2*pi)*exp( -(x-μ)**2/σ/σ/2 )
big = sqrt(2*σ**2*x + 2*s**2)* exp(
        b**2/(2*x*(σ**2*x + s**2)) - b**2/(2*x*s**2) - b*μ/(σ**2*x + s**2) + b*δ/s**2 + μ**2*x/(2*σ**2*x + 2*s**2) - μ**2/(2*σ**2) + μ*δ/σ**2 - x*δ**2/(2*s**2) - δ**2/(2*σ**2)
    ) / (2*sqrt(pi)*σ*s)

F, E = Wild('F', exclude=[exp,]), Wild('E', exclude=[exp,])
pprint( gauss)
print("match:")
pprint( gauss.match(F*exp(E)))
pprint( big )
print("match:")
pprint( big.match(F*exp(E)))  # Fails to match and returns None.

big 的公式显然是 F*exp(E) 的形式，但 match() 对它不起作用。如何从类似高斯的公式中提取指数和归一化因子？

如何在 exp() 的方程中分离出指数？

Answer 1

仅对Wild F 使用排除。然后E 将匹配表达式的exp 部分。这适用于高斯但不是很大——我不知道为什么。因此，或者，您仍然可以仅在 Wild F 上使用排除并使用 m = foo.match(E*F) - 指数将为 E.exp.

另一种无需匹配即可破解此 Mul 的方法是使用 as_independent

coeff, e = foo.as_independent(exp)
expo = e.exp

How to isolate the exponent in Sympy（如何在 Sympy 中“匹配”大公式）

How to isolate the exponent in Sympy (how to `match` large formulas in Sympy)

python

symbolic-math

sympy