4

SymPy で、次の方法で一般的なガウス分布を統合すると:

x = Symbol("x", real=True)
y = Symbol("y", real=True)
s = Symbol("s", real=True)

gaussian = exp(-((x-y)**2)/(2*(s**2)))
nfactor = simplify(integrate(gaussian, (x,-oo,oo)))

私はnfactorと にPiecewise依存しperiodic_argumentていpolar_liftます。以下によって作成されたものと同じオブジェクト:

Abs(periodic_argument(polar_lift(s)**(-2), oo))

pi/2私が望むソリューションよりも小さくなければなりません。この条件だけを持つ良い方法はありますか - Mathematica ではAssumingandRefineまたはSimplifyを使用するかもしれませんが、私は sympy が初めてで、ここで何をすべきかわかりません。

4

2 に答える 2

4

Symbol答えには、バインドされた の平方根が含まれsます。

s適切に単純化するために、肯定を宣言します。に：

s = Symbol('sigma', real=True, positive=True)
于 2013-05-17T00:52:00.437 に答える