Sympy符号行列式或evalf函数存在bug？不同计算顺序结果差异巨大

这是数值计算中典型的灾难性抵消问题，既不是操作错误，也不是Sympy本身的功能bug，是不同计算流程的数值稳定性差异导致的。

结果差异原因解释 #

• det1的计算逻辑是先把所有数值代入矩阵，得到纯数值矩阵后再调用行列式计算接口，内部使用LU分解这类数值稳定性极强的算法完成行列式计算，因此得到的结果误差在可接受范围内，符合预期。
• det3的计算逻辑是先得到符号行列式的完整表达式，再代入所有数值做精确符号运算，没有引入浮点误差，因此可以得到精确为0的正确结果。
• det2的问题出在计算顺序：先对符号矩阵做行列式展开，得到的符号表达式包含大量绝对值极高的项，这些项在理论上可以完全抵消得到0，但直接调用evalf时会用有限精度浮点数逐个计算每一项的值再累加，大项的浮点计算误差会完全覆盖掉抵消后的真实小值，最终得到完全错误的结果。哪怕调高evalf的精度，只要精度没有达到抵消需要的阈值，结果都不会正确。

解决方案 #

日常使用时优先选择det1的计算流程：先代入数值得到数值矩阵再计算行列式，这是数值行列式计算的最优实践，稳定性最高。如果确实需要先展开符号行列式再求值，可以先对det返回的符号结果调用simplify等符号化简接口，消去可抵消的项之后再调用evalf，就能得到正确的结果，示例如下：

# 先做符号化简再做数值求值
det_expr = sympy.simplify(sympy.det(matrix))
det2_fixed = det_expr.evalf(subs={z1:-1,z2:-1,z3:1,mv1:50E-6,mv2:90E-6,mv3:150E-6,R:8.314,T:298,D12:1E-9,D23:9E-9,D13:3E-9,D14:0.165171466666667,D24:0.0917619259259259,D34:0.0550571555555556,x1:0.3,x2:0.2,x3:0.5})

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