切线的魔法：用 SymPy 和 Manim 轻松搞定导数动画

合集 - manim动画(43)

1.【manim动画教程】-- 安装2023-03-282.【manim动画教程】-- 基本图形2023-03-293.【manim动画教程】-- 坐标系2023-04-104.【manim动画教程】-- 文本样式2023-04-075.【manim动画教程】-- 文字和公式2023-04-046.【manim动画教程】-- 图形样式2023-03-317.【manim动画教程】--相机2023-04-198.【manim动画教程】--高级动画效果2023-04-149.【manim动画教程】--常用动画效果2023-04-1210.【manim】之滚动字幕2022-12-0611.【manim】之圆规动画2023-01-3112.【manim】之目录动画2023-02-2813.Manim：动画制作背后的魔法2025-05-1814.manim 动画效果总结2025-05-2015.manim变换效果总结2025-05-2116.Manim动画渲染：从代码到屏幕的幕后故事2025-05-2217.轻松掌握Manim的.animate语法：让动画编程更简单2025-05-2518.Manim中三种函数图像类的比较2025-06-0519.掌握ChangeSpeed类：让数学动画速度随心而动！2025-06-1520.把数学对象画出来：Manim Mobject类库速查手册2025-08-1621.「亲手画出动态数学」：让数学可视化触手可及的Manim入门课2025-11-0722.让你的动画“活”过来：Manim 节奏控制指南 (Rate Functions)2025-11-2323.Manim进阶：用背景图片让你的数学视频脱颖而出2025-11-2624.Manim v0.19.1 发布啦！三大新特性让动画制作更丝滑2025-12-0225.🚀 Manim CE v0.20.0 发布：动画构建更丝滑，随机性终于“可控”了！02-2126.ManimCE v0.20.1 发布：LaTeX 渲染修复与动画稳定性提升03-0227.探索视觉的边界：用 Manim 重现有趣的知觉错觉03-0728.用Manim实现动态交点计算--从一个动点问题说起04-1629.Manim如何在数学公式中完美显示中文？04-1930.如何灵活设置公式中各个部分的颜色？04-2231.如何让多个动画“齐步走”？04-2632.告别手动计算，SymPy 初识与 Manim 联动05-1633.让数学公式自动推导05-17

34.切线的魔法：用 SymPy 和 Manim 轻松搞定导数动画05-19

35.填充与积累：积分与面积的可视化05-2136.轨迹的蓝图：方程求解与交点计算05-2337.Manim物理模拟：别自己写欧拉了！05-2538.用 SymPy 解决 Manim 曲线绘制速度不均的问题06-0139.一次函数图像工厂：用 SymPy 自动生成 y=kx+b 对比动画06-0340.用SymPy自动计算抛物线求根、判别式与顶点06-0541.用SymPy自动求解三角形构造与全等条件验证06-0842.用SymPy自动求解追及问题的方程06-1043.用SymPy自动因式分解：从面积拼图到代数恒等式06-13

收起

大家好，你有没有试过在 Manim 里做导数定义的动画？

就是那个经典的场景：画一条曲线，再画一条割线，然后让割线上的一个点无限逼近另一个点，最后变成切线。

这个过程的核心是计算割线的斜率 (f(x+h) - f(x)) / h，并观察当 h 趋近于 0 时，这个斜率是如何变化的。

听起来很简单，但实际操作起来，手动去推导极限、计算每一帧的坐标，不仅繁琐，还特别容易出错。

相信不少朋友都为此头疼过。

想象一下，我们要为函数 $f\left(x\right)=x^3-2x+1$ 做一个在 $x=1$ 处的切线动画。

1. 定义割线：我们需要两个点，$P(1,f(1\left)\right)$ 和 $Q(1+h,f(1+h\left)\right)$。
2. 计算斜率：$slope=\left(f\right(1+h)-f(1\left)\right)/h$。
3. 求极限：为了让动画平滑过渡到切线，我们需要知道当 $h\to 0$ 时，$slope$ 的精确值，也就是 $f^{\prime }\left(1\right)$。
4. 动态更新：在动画中，$h$ 是一个不断变小的值（比如从 1 变到 0.01），我们需要为每一个 h 实时计算 Q 点的坐标和割线的斜率。

如果手动来做，第2、3步就需要展开 $(1+h{)}^3-2(1+h)+1$，再减去 $f\left(1\right)$，化简，最后求极限。

对于复杂的函数，这简直是灾难！（比如函数$f\left(x\right)=sin\left(x^2\right)$）

而且，在代码里硬编码这些公式，一旦函数变了，所有计算都得重来。

这就是我们的痛点：动态、精准、自动化地处理符号计算。

SymPy 解决方案：让计算机做数学

SymPy 正是解决这个问题的完美工具，它可以把 x, h 当作真正的数学符号来处理，而不是具体的数字。

针对我们的需求，只需要两个核心函数：

• diff(f, x): 自动求导。直接告诉我们 f(x) 的导函数是什么。
• limit(expr, h, 0): 计算极限。可以验证我们的割线斜率在 h->0 时的确等于导数值。

下面看一段核心的 SymPy 代码，感受一下它的威力：

from sympy import symbols, diff, limit

# 定义符号变量
x = symbols('x')

# 定义我们的函数 f(x)
f = x**3 - 2*x + 1

# --- 核心操作 ---
# 自动求导，得到 f'(x)
f_prime = diff(f, x)
print(f"导函数 f'(x) = {f_prime}") 
# 输出: 导函数 f'(x) = 3*x**2 - 2

# 在 x=1 处的导数值
slope_at_1 = f_prime.subs(x, 1)
print(f"x=1 处的瞬时变化率 (斜率) = {slope_at_1}")
# 输出: x=1 处的瞬时变化率 (斜率) = 1

# 用极限来验证割线斜率
# 割线斜率表达式
secant_slope_expr = (f.subs(x, 1+h) - f.subs(x, 1)) / h
# 计算 h->0 时的极限
limit_slope = limit(secant_slope_expr, h, 0)
print(f"通过极限计算得到的斜率 = {limit_slope}")
# 输出: 通过极限计算得到的斜率 = 1

看！我们完全不用关心中间复杂的代数运算，SymPy 几行代码就帮我们完成了求导和极限验证，并且结果精确无误。

这为我们接下来的 Manim 动画提供了坚实的数学基础。

Manim 联动实战：让切线“动”起来

现在，我们将 SymPy 的计算能力嵌入到 Manim 动画中。

我们将使用 ValueTracker 来控制 h 的值，让它从一个较大的数（如1）逐渐减小到接近0。

在每一帧，Manim 都会调用 SymPy 重新计算 Q 点的位置和割线，从而实现动态效果。

下面是核心的代码：

from manim import *
from sympy import symbols, lambdify, diff

class DerivativeAnimation(Scene):
    def construct(self):
        # ========== SymPy 符号计算部分 ==========
        x_sym = symbols("x")
        f_sym = x_sym**3 - 2*x_sym + 1           # 原函数：f(x) = x³ - 2x + 1
        f = lambdify(x_sym, f_sym, "numpy")       # 转为 NumPy 函数供绘图

        f_prime_sym = diff(f_sym, x_sym)          # SymPy 自动求导：f'(x) = 3x² - 2
        x_p = 1                                    # 切点横坐标
        exact_k = float(f_prime_sym.subs(x_sym, x_p))  # 精确导数 f'(1) = 1

        # ========== Manim 坐标系与曲线 ==========
        ax = Axes(x_range=[-2, 3], y_range=[-3, 5])
        graph = ax.plot(f, color=YELLOW)          # 原函数曲线

        p_point = Dot(ax.c2p(x_p, f(x_p)), color=RED)  # 切点 P

        # ========== ValueTracker 驱动割线动态逼近 ==========
        h_tracker = ValueTracker(1)               # h 从 1 逐渐减小到 0.001

        # 割线：随 h 变化而重新绘制
        def get_secant_line():
            h_val = h_tracker.get_value()
            x_q = x_p + h_val
            k = (f(x_q) - f(x_p)) / h_val         # 割线斜率 Δy/Δx
            return ax.plot(
                lambda x: k * (x - x_p) + f(x_p), # 点斜式
                color=GREEN, x_range=[x_p - 1, x_q + 1]
            )

        secant_line = always_redraw(get_secant_line)

        # 切线：使用 SymPy 算出的精确导数
        tangent_line = ax.plot(
            lambda x: exact_k * (x - x_p) + f(x_p),
            color=PURPLE, x_range=[-0.5, 2.5]
        )

        # ========== 动画流程 ==========
        self.play(Create(ax), Create(graph), Create(p_point))
        self.play(Create(secant_line))

        # 核心：h → 0，割线动态逼近切线
        self.play(
            h_tracker.animate.set_value(0.001),
            run_time=5,
            rate_func=rate_functions.ease_in_out_quad,
        )

        # 对比展示精确切线
        self.play(Create(tangent_line))
        self.wait(1)

代码关键点解析

• lambdify: 连接 SymPy 和 Manim 的桥梁。它把 SymPy 的符号表达式 f_sym 转换成一个普通的 Python 函数 f，这个函数可以接受 NumPy 数组作为输入，正好符合 Manim ax.plot() 的要求。
• ValueTracker: Manim 中创建动态效果的核心。h_tracker 存储了 h 的当前值。
• always_redraw: 这个装饰器告诉 Manim，被它修饰的对象（如 q_point 和 secant_line）需要在每一帧都重新计算和绘制。它们内部的函数 get_q_point 和 get_secant_line 会读取 h_tracker 的最新值，并调用 f 函数来获取最新的坐标。
• 动态割线: 在 get_secant_line 中，我们虽然可以直接用两点式画线，但这里展示了如何利用 SymPy 的思想——通过计算斜率和截距来定义直线，逻辑更清晰。

效果展示说明

运行这段代码，你会看到以下动画效果：

1. 坐标系与函数登场：黄色的三次函数$f\left(x\right)=x^3-2x+1$被绘制出来。
2. 固定点 P：在$x=1$处，一个红色的点 P 被标记出来。
3. 动态点 Q 与割线：一个蓝色的点 Q 出现在 P 的右侧（因为初始 h=1），一条绿色的割线连接 P 和 Q。
4. 魔法时刻：动画开始，Q 点开始平滑地向 P 点移动（h 值从 1.5 逐渐减小到 0.01）。与此同时，绿色的割线也随之旋转。
5. 切线显现：当 Q 无限接近 P 时，割线几乎不再变化。此时，一条紫色的精确切线被绘制出来，你会发现它和最终的割线几乎完全重合！

整个过程直观地展示了导数作为瞬时变化率的几何意义，而这一切的精准性都由 SymPy 在幕后保证。

小结

我们已经成功地将 SymPy 的符号计算能力与 Manim 的动画渲染能力结合起来，解决了制作导数定义动画时的手动计算痛点。

通过 diff 和 limit，我们获得了精确的数学结果；

通过 ValueTracker 和 always_redraw，我们让这些结果在屏幕上“活”了起来。

这种 “SymPy 负责思考，Manim 负责表现” 的模式非常强大，可以应用到各种复杂的数学可视化场景中。