如何在不使用SymPy的情况下用Python实现符号导数求解？

嘿，已经搞定逆波兰式转换这一步真的超棒——这可是符号求导的核心基础之一！接下来咱们一步步把求导逻辑搭起来，其实本质就是把每个运算符的求导规则，对应到逆波兰式的栈操作上，咱们拆解细说：

逆波兰式天生适合栈操作，求导也可以沿用这个思路：遍历逆波兰式的每个元素，遇到操作数就直接计算它的导数压栈；遇到运算符就弹出栈顶的两个操作数（注意顺序：第一个弹出的是右操作数，第二个是左操作数），然后应用对应的求导规则，把组合后的导数表达式（逆波兰式形式）压回栈。最后栈里剩下的就是整个表达式的导数逆波兰式。

先搞定最基本的单个元素求导，这是所有复杂表达式的基础：

• 常数（比如2、5.7）：导数为0
• 自变量x：导数为1
• 幂项x^n（n为常数）：导数为n*x^(n-1)（如果n=1，结果简化为1；n=0的话原项是常数，导数为0）

针对你需要支持的+、-、*、/、^，每个都有对应的求导逻辑，咱们逐个说：

• 加法（+）：用和的导数等于导数的和，即 d/dx(A+B) = dA/dx + dB/dx
  操作逻辑：弹出右操作数B、左操作数A，分别求它们的导数dA、dB，把dA dB +（逆波兰式形式）压回栈
• 减法（-）：差的导数等于导数的差，即 d/dx(A-B) = dA/dx - dB/dx
  操作逻辑：弹出B、A，求dA、dB，把dA dB -压回栈
• 乘法（*）：用乘积法则，即 d/dx(A*B) = A*dB/dx + B*dA/dx
  操作逻辑：弹出B、A，求dA、dB，把A dB * B dA * +压回栈
• 除法（/）：用商法则，即 d/dx(A/B) = (dA*B - A*dB)/(B²)
  操作逻辑：弹出B、A，求dA、dB，把dA B * A dB * - B 2 ^ /压回栈
• 幂运算（^）：这里要分三种情况简化（避免不必要的复杂表达式）：
  1. 底数是x，指数是常数n：直接用n*x^(n-1)（对应逆波兰式n x n-1 ^ *）
  2. 底数是常数c，指数是x：导数为c^x * ln(c)（对应逆波兰式c x ^ c ln *）
  3. 一般情况（A、B都含x）：用对数求导法，即 d/dx(A^B) = A^B * (dB*lnA + B*dA/A)，对应逆波兰式A B ^ dB A ln * B dA A / * + *

1. 标记元素类型：先给逆波兰式的每个元素打上标签，区分是「常数」「变量x」「运算符」，这样处理时能快速判断逻辑
2. 递归求导函数：写一个函数，输入是逆波兰式子列表，输出是该子式的导数逆波兰式：
   • 初始化空栈
   • 遍历每个元素：
     • 如果是常数/变量：直接计算导数，把导数的逆波兰式元素压栈
     • 如果是运算符：弹出栈顶两个元素（B和A），递归求dA和dB，再根据运算符规则组合成新的逆波兰式压栈
3. 表达式简化：最后可以加一步简化逻辑，比如把0 + x改成x，x*1改成x，x^1改成x，让结果更简洁

比如你提到的表达式x^2 + 2，逆波兰式是x 2 ^ 2 +，求导过程：

1. 遍历x：导数是1，压栈 → [1]
2. 遍历2：导数是0，压栈 → [1, 0]
3. 遍历^：弹出0（B=2，常数）和1（A=x的导数），应用x^n规则，生成2 x 1 ^ *压栈 → [2 x 1 ^ *]
4. 遍历2：导数是0，压栈 → [2 x 1 ^ *, 0]
5. 遍历+：弹出0和2 x 1 ^ *，组合成2 x 1 ^ * 0 +，简化后就是2x，完全符合预期！

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