Sympy未识别表达式为实数：共轭相乘未获预期简化结果

你遇到的问题是SymPy默认没有充分利用变量的正实数属性，导致sqrt(exp(gamma*t)-1)*conjugate(sqrt(exp(gamma*t)-1))无法自动简化为exp(gamma*t)-1。以下是几种可行的解决方式：

方法1：用refine明确根号内表达式的正性 #

通过refine函数结合Q.positive假设，让SymPy确认根号内的表达式是正实数，这样共轭运算后就能正确简化：

from sympy import symbols, sqrt, exp, conjugate, refine, Q

gamma, t = symbols('gamma t', real=True, positive=True)
expr = sqrt(exp(gamma*t) - 1)
# 给根号内的表达式添加正性假设
expr_refined = refine(expr, Q.positive(exp(gamma*t) - 1))
result = expr_refined * conjugate(expr_refined)
print(result.simplify())  # 输出 exp(gamma*t) - 1

方法2：使用real_root替代sqrt #

real_root专门处理实数域的根运算，对于正实数输入，它的共轭就是自身，相乘后直接得到预期结果：

from sympy import symbols, exp, conjugate, real_root

gamma, t = symbols('gamma t', real=True, positive=True)
expr = real_root(exp(gamma*t) - 1, 2)  # 等价于正实数域的平方根
result = expr * conjugate(expr)
print(result)  # 直接输出 exp(gamma*t) - 1

方法3：定义表达式时指定正实数属性 #

在创建sqrt表达式时，直接通过参数声明它是正实数，这样共轭运算会自动识别并简化：

from sympy import symbols, sqrt, exp, conjugate, simplify

gamma, t = symbols('gamma t', real=True, positive=True)
expr = sqrt(exp(gamma*t) - 1, real=True, positive=True)
result = expr * conjugate(expr)
print(simplify(result))  # 输出 exp(gamma*t) - 1

原因说明 #

SymPy的simplify默认不会强制应用所有变量假设，尤其是sqrt作为多值函数，默认会保留复数域的一般形式。只有明确告知上下文是正实数域，SymPy才会按照实数规则进行简化。

内容的提问来源于stack exchange，提问作者james_12890