Sympy部分分式分解残差异常：与Octave/Wolfram结果不符问询

核心问题本质是部分分式分解的表示存在等价变形空间，极点一致但残差不同，通常是不同工具的运算策略、精度处理逻辑差异导致，具体原因如下：

• 浮点数与符号运算的处理差异
  Sympy以符号运算为核心，会优先将浮点数系数转换为高精度有理数或近似符号形式再执行分解；而Octave的residue()、Wolfram Alpha多基于数值算法直接计算残差。这种处理逻辑的差异会导致残差的数值呈现形式不同，但分解式整体与原函数是代数等价的——将不同工具的分解式通分后，结果应与原函数一致（考虑浮点数精度误差）。

• 重极点的分解形式差异
  若原函数存在重极点，部分分式分解存在多种等价写法。例如二阶重极点，Sympy可能将残差表示为导数形式的简洁项，而其他工具可能展开为多个线性项的组合，虽然残差数值看似不同，但整体分解式的数学意义完全一致。

• apart方法的参数影响
  你使用的apart(Z, full=True).doit()中，full=True会强制展开所有可能的分式，但Sympy的分解逻辑优先选择符号上更简洁的形式，而非匹配数值工具的输出格式。尝试去掉full=True，直接调用apart(Z)，可能得到更接近其他工具的残差结果。

• 等价性验证方法
  要确认残差差异是否为等价变形，可通过以下操作验证：

  1. 将Sympy得到的分解式通分，与原函数Z做差，结果应趋近于0（允许浮点数精度范围内的误差）。
  2. 代入具体的s值（如s=1），计算不同分解式的结果，应与原函数的计算值一致。

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