Python中如何用sympy定义函数导函数？附代码及符号变量疑问

你代码的核心问题 #

你当前代码无法正常求导的核心原因是混用了数值计算库math和符号计算库SymPy的逻辑，两者无法直接兼容。

关于x = sp.Symbol('x')的作用 #

这是SymPy符号计算的基础声明：

• 普通Python里给变量赋值是存具体数值，比如x=3就是让x绑定数字3；而sp.Symbol('x')是生成一个代表数学未知变量x的特殊对象，告诉SymPy后续用这个x拼接的内容是数学公式，不是要立刻算出数值结果的算术式，只有声明了符号变量，SymPy才能完成求导、积分、公式化简这类符号运算。
• 你之前写的f(x)函数里全用了math.cos、math.pi这类math库的接口，这类接口只能处理浮点数输入，识别不了SymPy的符号对象，自然没法被sp.diff识别求导；另外你写的derv函数里调用sp.diff(f)时没有传入要求导的目标变量，也是语法错误。

可直接运行的正确实现 #

直接用SymPy内置的数学函数写符号表达式即可，不需要套普通Python函数的壳，求导得到的符号表达式还可以直接转成支持数值代入计算的可调用函数，参考代码如下：

import sympy as sp

# 声明符号变量x
x = sp.Symbol('x')

# 常量定义
W = 15 #kN/m
E = 70 # Gpa
I = 52.9*10**(-6) #m**4
L = 3 #m 
e = 0.01
xi = 1.8 
y = 9

# 用SymPy内置函数写原函数的符号表达式，不要用math库的函数/常量
f = ( ( y*3*(sp.pi**4)*E*I/(W*L) ) - ( 48*(L**3)*sp.cos(sp.pi*x/(2*L)) ) + ( 48*(L**3) ) + ( (sp.pi**3)*(x**3) ) )/(3*L*(sp.pi**3))**(sp.Rational(1,2))

# 对x求导，得到导函数的符号形式
f_deriv = sp.diff(f, x)

# 输出导函数公式
print("导函数表达式：", f_deriv)

# 如果需要代入具体数值计算导函数结果，可以转成数值计算函数
f_deriv_calc = sp.lambdify(x, f_deriv, 'numpy')
# 示例：计算x=1位置的导数值
print("x=1处的导数值：", f_deriv_calc(1))

关键注意点 #

• 做符号运算阶段，统一用sp.pi、sp.cos这类SymPy自带的数学常量、函数，不要混用math库的对应接口，math库仅支持浮点数计算，无法处理符号对象。
• 写分数次幂时优先用sp.Rational(分子,分母)代替直接写浮点数，比如把1/2写成sp.Rational(1,2)，避免Python提前计算浮点数带来的精度损失。
• 调用sp.diff时必须明确传入求导的目标符号变量，否则多变量场景下会出现计算错误。

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