SymPy积分触发类型错误但Mathematica运行正常的差异与原因分析

你遇到的类型错误是SymPy符号积分接口对矩阵输入的适配缺失导致的：integrate接口默认对标量表达式设计，你直接传入矩阵形式的被积函数时，内部处理Piecewise分段函数的分支逻辑时，会出现「某一分支返回标量0、另一分支返回矩阵对象」的情况，SymPy尝试合并两个分支的返回值时无法完成矩阵和标量的加法操作，直接触发报错。
你之前查阅的矩阵替换标量相关问题与本次报错无关。

• Mathematica的NIntegrate是数值积分接口，首先会识别输入的向量/矩阵结构，自动对每个分量单独执行积分操作，全程保持数据结构类型一致，不会出现标量与矩阵混合的问题。
• SymPy的integrate默认走符号积分路径，没有内置矩阵输入的逐分量处理逻辑，在进入积分计算步骤前，会先完成分段表达式的分支合并校验，此时就先触发了类型冲突，根本不会走到后续的积分计算环节。

调整积分逻辑，手动对被积矩阵的每个分量单独执行积分，再合并为矩阵即可，同时建议将代码中用到的np.pi替换为sym.pi，避免numpy数值类型与SymPy符号类型的隐性冲突，修改后代码如下：

# 原k_u定义调整，把np.pi换成sym.pi
k_u = (2 * sym.pi)/lambda_u

# 调整后的积分函数
def rad(n, beta, betap):
    integrand = n.cross((n-beta).cross(betap))
    # 逐分量积分后重构矩阵
    res = Matrix([integrate(comp, (t, t_start, t_end)) for comp in integrand]).reshape(integrand.shape)
    # 若需要直接得到数值结果，取消下一行注释
    # res = sym.N(res)
    return res

执行上述代码后可得到与Mathematica一致的计算结果：Matrix([[0.00150421, 0, 0]])

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