Sympy计算复合函数导数后，如何去除doit()残留的Subs对象？

要解决这个问题，我们可以通过递归遍历表达式，将所有Subs对象替换为直接作用在目标函数上的导数形式，下面是通用的解决方案：

1. 重现问题 #

先确认你的场景：

from sympy import *
t = symbols('t')
p = Function('p')
x = Function('x')

v = diff(x(p(t)), t)
a = diff(v, t)
print(a.doit())
# 输出包含Subs的结果：
# Derivative(p(t), t)**2*Subs(Derivative(x(_xi_3), _xi_3, _xi_3), (_xi_3,), (p(t),)) + Derivative(x(p(t)), p(t))*Derivative(p(t), t, t)

2. 通用替换方法 #

我们可以写一个递归函数来处理所有嵌套的Subs对象，确保表达式里不再有占位符变量：

def remove_subs(expr):
    # 处理Subs对象：替换占位变量为实际参数并求值
    if isinstance(expr, Subs):
        deriv_expr = expr.expr
        placeholder_var = expr.variables[0]
        target_arg = expr.point[0]
        return deriv_expr.subs(placeholder_var, target_arg).doit()
    # 递归处理复合表达式（加法、乘法、幂次等）
    elif expr.is_Add or expr.is_Mul or expr.is_Pow:
        return expr.func(*[remove_subs(arg) for arg in expr.args])
    # 处理导数对象，递归处理其内部参数
    elif isinstance(expr, Derivative):
        cleaned_args = [remove_subs(arg) for arg in expr.args]
        return Derivative(*cleaned_args)
    # 其他类型直接返回
    else:
        return expr

3. 测试效果 #

调用这个函数处理你的结果：

clean_a = remove_subs(a.doit())
print(clean_a)
# 输出符合要求的结果：
# Derivative(p(t), t)**2*Derivative(x(p(t)), p(t), p(t)) + Derivative(x(p(t)), p(t))*Derivative(p(t), t, t)

4. 简洁版实现 #

如果不需要处理特别复杂的嵌套场景，也可以用SymPy的replace方法结合匿名函数快速实现：

clean_a = a.doit().replace(
    lambda e: isinstance(e, Subs),
    lambda e: e.expr.subs(e.variables[0], e.point[0]).doit()
)

原理说明 #

SymPy在计算复合函数高阶导数时，会用占位符变量（比如_xi_1）来临时表示中间变量，最终生成Subs对象。我们的方法就是把这些占位符替换回实际的p(t)，并让SymPy将导数转换为直接对x(p(t))求导的形式，这个逻辑适用于任意复杂的复合函数导数表达式。

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