SymPy求解实变量tanh方程出现异常解及简化问题咨询

为什么会出现含负数对数的解？ #

SymPy的solve函数默认会先在复数域内求解方程的通解，再尝试代入变量的实约束，但不会自动完全过滤掉不符合实数域要求的解。
对于tanh(a)=b，实变量下b的有效定义域是(-1,1)，此时-(b+1)/(b-1) = (b+1)/(1-b)是正数，sqrt(-(b+1)/(b-1))是正实数，但第一个解里的-sqrt(...)是负数，对负数取自然对数在实数域无意义——这个解其实是复数域内的解，SymPy没有自动剔除它。

如何将解简化为a=atanh(b)的形式？ #

有两种常用方法：

1. 直接调用反双曲正切函数：因为atanh(b)本身就是tanh(a)=b的实解，直接使用sp.atanh(b)即可得到简洁形式。
2. 对solve结果化简或指定实域求解：
   • 对solve得到的对数形式解用sp.simplify()处理，示例代码：

     solutions = sp.solve(my_Eq, a)
     simplified_sol = [sp.simplify(sol) for sol in solutions]
     

     化简后符合实约束的解会转化为atanh(b)，无意义的复数解则会保留或被识别为非实解。
   • 调用solve时指定domain=sp.Reals，强制在实数域求解：

     sp.solve(my_Eq, a, domain=sp.Reals)
     

     此时会直接返回[atanh(b)]。

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