如何简化方程式中的指数
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 =符号( 'B',真=真) 没有得到简化的指数 只有当指数是常数时才会简化
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
答案 0 :(得分:7)
(x m ) n = x mn 为真only if m, n are real。
>>> 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 supports complex numbers,所以我认为除非你能证明b
是真实的,否则不应该进行这种简化。
编辑:如果x不是正数,也是假的。
>>> x = -2
>>> m = 2
>>> n = 0.5
>>> (x**m)**n
2.0
>>> x**(m*n)
-2.0
编辑(通过gnibbler):这是应用Kenny限制的原始示例
>>> from sympy import symbols
>>> a,b=symbols('ab', real=True, positive=True)
>>> j=(a**b**5)**(b**10)
>>> print j
a**(b**15)
答案 1 :(得分:2)
a,b,c=symbols('abc',real=True,positive=True)
(a**b**5)**b**10
a**(b**15)#ans
答案 2 :(得分:0)
这可能与this bug有关。