4

sympyの方程式の指数を単純化する方法

from sympy import symbols 
a,b,c,d,e,f=symbols('abcdef')
j=(a**b**5)**(b**10)
print j 
(a**(b**5))**(b**10) #ans even after using expand simplify 
# desired output 
a**(b**15)

sympyでそれが不可能な場合、どのモジュールをPythonにインポートする必要がありますか？

'b'を実数として定義した場合でも編集し、他のすべての記号も編集します

b = symbol（'b'、real = True）単純化された指数を取得しない指数が定数の場合にのみ単純化

a=symbols('a',real=True)
b=symbols('b',real=True)
(a**5)**10
 a**50  #simplifies only if exp are numbers
(a**b**5)**b**10


(a**(b**5))**b**10  #no simplification
4

3 に答える 3

7

（x m）n = x mnは、m、nが実数の場合にのみ真になります。

>>> import math
>>> x = math.e
>>> m = 2j*math.pi
>>> (x**m)**m      # (e^(2πi))^(2πi) = 1^(2πi) = 1
(1.0000000000000016+0j)
>>> x**(m*m)       # e^(2πi×2πi) = e^(-4π²) ≠ 1
(7.157165835186074e-18-0j)

AFAIK、sympyは複素数をサポートしているので、本物であることが証明できない限り、この単純化を行うべきではないと思いますb。

編集：xが正でない場合もfalseです。

>>> x = -2
>>> m = 2
>>> n = 0.5
>>> (x**m)**n
2.0
>>> x**(m*n)
-2.0

編集（gnibblerによる）：これは、ケニーの制限が適用された元の例です

>>> from sympy import symbols 
>>> a,b=symbols('ab', real=True, positive=True)
>>> j=(a**b**5)**(b**10)
>>> print j
a**(b**15)
于 2010-07-18T06:24:06.950 に答える
2 
a,b,c=symbols('abc',real=True,positive=True)
(a**b**5)**b**10
a**(b**15)#ans
于 2010-07-18T07:36:48.880 に答える
0

これは、このバグに関連している可能性があります。

于 2010-07-18T06:10:35.890 に答える