惯性聚合 高效追踪和阅读你感兴趣的博客、新闻、科技资讯
阅读原文 在惯性聚合中打开

推荐订阅源

C
CXSECURITY Database RSS Feed - CXSecurity.com
Google Online Security Blog
Google Online Security Blog
The Last Watchdog
The Last Watchdog
S
Security @ Cisco Blogs
Help Net Security
Help Net Security
Security Archives - TechRepublic
Security Archives - TechRepublic
Exploit-DB.com RSS Feed
Exploit-DB.com RSS Feed
W
WeLiveSecurity
cs.CV updates on arXiv.org
cs.CV updates on arXiv.org
H
Hacker News: Front Page
人人都是产品经理
人人都是产品经理
aimingoo的专栏
aimingoo的专栏
Vercel News
Vercel News
Microsoft Azure Blog
Microsoft Azure Blog
小众软件
小众软件
Project Zero
Project Zero
T
Tailwind CSS Blog
V
Vulnerabilities – Threatpost
P
Privacy & Cybersecurity Law Blog
Know Your Adversary
Know Your Adversary
Last Week in AI
Last Week in AI
腾讯CDC
Schneier on Security
Schneier on Security
cs.AI updates on arXiv.org
cs.AI updates on arXiv.org
WordPress大学
WordPress大学
量子位
L
Lohrmann on Cybersecurity
J
Java Code Geeks
Cyberwarzone
Cyberwarzone
Recent Announcements
Recent Announcements
IT之家
IT之家
博客园_首页
罗磊的独立博客
博客园 - 聂微东
奇客Solidot–传递最新科技情报
奇客Solidot–传递最新科技情报
Application and Cybersecurity Blog
Application and Cybersecurity Blog
F
Fortinet All Blogs
博客园 - Franky
P
Palo Alto Networks Blog
V2EX - 技术
V2EX - 技术
让小产品的独立变现更简单 - ezindie.com
让小产品的独立变现更简单 - ezindie.com
Martin Fowler
Martin Fowler
N
News and Events Feed by Topic
C
Cybersecurity and Infrastructure Security Agency CISA
B
Blog RSS Feed
Cisco Talos Blog
Cisco Talos Blog
TaoSecurity Blog
TaoSecurity Blog
H
Heimdal Security Blog
G
GRAHAM CLULEY
Cloudbric
Cloudbric

博客园_首页

Plist 二进制格式 Milvus 和 PGVector,哪个更好? OpenClaw 已过时?在 VS Code 中运行 Hermes Agent! 第30篇文章:一个大三计科生的自白 Manim如何在数学公式中完美显示中文? Docker 部署 RocketMQ 5 并发编程核心概念辨析 C#事务处理最佳实践:别再让“主表存了、明细丢了”的破事发生 CLI 是什么?为什么大厂突然集体卷命令行? 【从0到1构建一个ClaudeAgent】协作-自主Agent UIImageView 设置图片不生效的原因排查 最小二乘问题详解20:无先验约束下的增量式SFM自由网平差 痞子衡嵌入式:大话双核i.MXRT1180之XIP应用里借助MU实现可靠Flash IAP的方法 AI Chat 封装, SemanticKerne.AiProvider.Unified 已发布 Windows下右键编辑js文件无法打开记事本——在注册表中使用环境变量 在后台服务中使用 Scoped 服务,为什么总是报错? H200 安装驱动并使用sglang启动模型 wireshark 抓包Trap上报告警内容 我用 AI 辅助开发了一系列小工具(2):图片压缩工具 [A Primer On MC and CC] 2.1 Memory Consistency 1 - 指令重排序和 SC 模型 Oracle数据库SCN推进技术详解与实践指南 玩转控件:封装个带图片的Label控件 Claude Code 4.7 真正该升级的不是模型,而是你的工作流 前端小白一句话,AI 帮我做了个颜值拉满的桌面媒体播放器。当代码不再是门槛,一句话编程就是现实。 5. WorkBuddy: 小龙虾的灵魂三件套,让你的小龙虾不只是工具 SQLite 分片方案实战:三种分片策略的深度对比 告别简陋 UI!一款基于 Fluent Design 和基于 WinUI 的开源免费、现代化的 Avalonia UI 控件库 关于二进制排列组合枚举的总结 AI开发-python-LangGraph框架(3-27-LangGraph从零实现大模型智能决策工作流) ElasticSearch主分片和副本分片概念详解 【002】HTTPS 粗解:证书、TLS 握手与对后端配置的影响 Hermes Agent 一周暴涨五万 Star,但我劝你别急着追 明明连接的是Redis的DB0,为什么能查到DB3的数据? 【从0到1构建一个ClaudeAgent】协作-Agent团队 熟悉电子元器件之后,电子小白下一步该怎么走? MAF快速入门(23)通过C#类定义Skills .NET 高级开发 | 手写一个对象映射框架 FastAPI数据库ORM怎么选?我肝了三个Demo后,终于不再纠结了 mysqldump 参数拾遗:在遗忘与铭记之间 C# .NET 周刊|2026年3月5期 Claude code入门 - 陈彦斌 一文学习入门 ThingsBoard 开源物联网平台 GitHub 热门项目 | 2026年04月16日 如何为GIT设置全局勾子,为每次提交追加信息 Number.isFinite和isFinite与isNaN()和Number.isNaN的区别 PortSwigger SQL注入LAB2 推荐一个测试人必备的Skills,从功能到性能全搞定(附详细实操和安装下载方式) 筑基期:掌握Odoo基础核心知识点02(Odoo XML 开发方式详解) GLM模型这么火,咱们用vllm也咧一个呗! 深入理解 AbortController:从底层原理到跨语言设计哲学 字符串学习笔记 多租户系统框架的基础模块设计和分析设计 Apache SeaTunnel Zeta 为什么能做到“又快又稳”? AI开发-python-LangGraph框架(3-26-LangGraph基本概念及第一个简单样例) Vue 3 组件通信,别只会用 Props 和 Emits 了,这几个狠活儿你得看看 ElasticSearch7.X版本配置密码 用Manim实现动态交点计算--从一个动点问题说起 团结引擎+Addressable+Instant Game打包抖音小游戏 function call 实战:让 LLM 自动判断 pod 异常、调用日志工具并完成故障分析 bubseek —— 让 Agent 的足迹,变成团队的洞察 通过 C# 读取并导出 PDF 书签 如何用 GitHub Actions 实现 Steam 自动化发布 【从0到1构建一个ClaudeAgent】并发-后台任务 .NET 高级开发 | 定制 ASP.NET Core 框架 电子小白:什么是运算放大器(运放) zero2Agent:面向大厂面试的 Agent 工程教程,从概念到生产的完整学习路线 堆上的ORW HC32F460 USB CDC通信异常:非对齐访问异常排查 20260413-Hyperbridge 攻击事件:发生在默克尔山上的验证绕过 那些喊着AI 要淘汰你的人,正在靠你的焦虑赚大钱! 深度学习进阶(八)Swin Transformer 最小二乘问题详解19:带先验约束的增量式SFM优化与实现 SnapTranslate 3.0 正式发布:全局划词翻译 + 完整英语学习闭环,一站式搞定查词、记词、复习 工作的意义、工作的困难认知再思考 .NET + AI 进阶实战:基于类的技能开发 - 打造可治理的 Agent 能力模块 【从0到1构建一个ClaudeAgent】规划与协调-技能 上周热点回顾(4.6-4.12) 电子小白的工具三件套:面包板、杜邦线、万能板 单表五亿数据的查询优化 | Mysql、StarRocks 2. WorkBuddy:从“我是谁”到“帮我干活” C# 如何减少代码运行时间:7 个实战技巧 基于HelixToolkit.SharpDX 渲染3D模型 - 笺上知微 从零开始的双臂具身VLA起源及现阶段发展综述 - SkyXZ 记对 xonsh shell 的使用, 脚本编写, 迁移及调优 - pluvium27 受够了Vibe Coding的失控?换个起点,让AI事半功倍 从开始配置漏洞环境到漏洞复现流程 - 難しい 关于10年工作经验的程序员对OpenClaw的实战经验分享以及看法 - 虚无境 Any metadata 的内存布局 C# .NET 周刊|2026年3月2期 - InCerry 我帮你测过了,测试圈排名第二的 Skill 依然很牛逼 Skill Discovery | 无监督技能发现的经典工作总结 - MoonOut 上下文工程是什么?过时了么?一文讲明白! - 一枫说码 开了 TUN 模式还是直连?90% 的人都踩过这个坑 AScript扩展多种脚本语言 - rockey627 AI 学习笔记:Agent 的记忆机制 你能被装进一个文件里吗?——7 万人把同事"蒸馏"成了 AI - 我没有三颗心脏 Claude Code 通关手册(七):给 AI 装上技能包——Skills 完全指南 - 暮色之狐 在浏览器中快速编辑代码:VSCode Web 集成实践 - Newbe36524 蒸馏自己 skill?基于 Deepseek 的蒸馏器,丐版蒸馏方式,简单便捷 - To_Carpe_Diem Spring AI Aliababa和AgentScope,哪个更好? - 苏三说技术
告别手动计算,SymPy 初识与 Manim 联动
wang_yb · 2026-05-16 · via 博客园_首页

下面是我正在做的一个抛物线演示动画。

需求很简单:展示一个二次函数 $ y = x^2 - 2x - 1 $ 的图像,并在上面标注几个关键点。

问题来了

  • 当我想调整函数参数时(比如把 $ -2x $ 改成 $ -3x $),所有点的坐标都要手动重算
  • 计算 $ x=1.5 $ 时的函数值?掏出计算器 → $ 1.5^2 - 2×1.5 - 1 = -1.75 $ → 再手动填回代码
  • 顶点坐标?求导 → 令导数等于0 → 解方程 → 计算 y 值 → 再填回代码
  • 对称轴和 x 轴的交点?求根公式 → 计算器 → 填代码

一个参数改动,我要重新计算七八个坐标值。这哪是在做动画,分明是在做数学作业!

直到我发现了 SymPy 这个神器。

SymPy 是什么?为什么 Manim 动画需要它?

简单说,SymPy 是一个 Python 的符号计算库

别被"符号计算"这个词吓到,用大白话讲就是:

让计算机帮你"列式子、解方程、求导数",而不是你自己手算。

数值计算 vs 符号计算

看一个更直观的对比,你会立刻明白符号计算的强大:

import math
import sympy as sp

# ========== 场景:计算 sin(π/3) 的精确值 ==========
# 数值计算 - 得到近似小数
result_num = math.sin(math.pi / 3)
print(f"数值计算: {result_num}")  
# 输出: 0.8660254037844386  ← 这是近似值,不知道它等于 √3/2

# 符号计算 - 得到精确表达式
x = sp.Symbol('x')
result_sym = sp.sin(sp.pi / 3)
print(f"符号计算: {result_sym}")  
# 输出: sqrt(3)/2  ← 精确的数学表达式!


# 场景1:求平方
print("\n=== 求 (sin(π/3))² ===")
# 数值计算 - 精度损失
square_num = result_num ** 2
print(f"数值: {square_num}")  
# 输出: 0.7499999999999999  ← 本应是 0.75,有浮点误差!

# 符号计算 - 精确化简
square_sym = result_sym ** 2
print(f"符号: {square_sym}")  
# 输出: 3/4  ← 精确值!

关键对比总结

特性 数值计算 (math) 符号计算 (sympy)
sin(π/3) 0.86602540378... √3/2
平方后 0.749999999999... 3/4
能否继续代数运算 ❌ 只能数值近似 ✅ 可代入方程、求导、化简
浮点精度问题 ⚠️ 存在误差累积 ✅ 完全精确

符号计算的灵活性体现在

  1. 保持数学形式√3/20.866... 更有数学意义
  2. 自动化简(√3/2)² 自动变成 3/4
  3. 代数兼容:可以继续解方程、求导、积分,保持精确形式

这对 Manim 动画尤为重要——你不仅需要坐标值,更需要数学关系的可视化,而符号计算保留了这种关系!

避免累积误差

符号计算在累积的计算中,能够有效的降低误差。

比如公式:$ I_n =1-n\times I_{n-1} $其中 $ I_0 = e-1 $。分别累积计算以后:

$ n $ 符号计算 数值计算 (模拟8位小数精度) 误差分析
0 $ e - 1 $ 0.71828183 初始误差: $ \approx 1.5 \times 10^{-9} $
1 $ 2 - e $ 0.28171817 误差微小
2 $ 2e - 5 $ 0.43656366 误差开始累积
3 $ 16 - 6e $ 0.30860902
4 $ 24e - 65 $ 0.23687292
5 $ 326 - 120e $ 0.18276460 误差开始显现
6 $ 1956e - 5315 $ 0.15054840
7 $ 13692 - 5040e $ 0.12145720
8 $ 109536e - 298325 $ 0.10364240
9 $ 985824 - 2691360e $ 0.08385840
10 $ 26913600e - 73309365 $ 0.07515840 偏差明显
11 $ 296049600 - 807408000e $ 0.09173440
12 $ 9688896000e - 26384952005 $ -0.08345440 灾难性错误:符号反转!
13 $ 342938611200 - 1258293216000e $ 2.10490720 完全失控

使用Sympy的话,可以在需要某一步结果的时候再代入$ e $去具体计算出来,不会累积误差。

SymPy 核心入门:把变量当作"符号"

在 SymPy 中,我们首先要定义符号变量

import sympy as sp

# 定义符号 - 告诉 SymPy "x 是一个数学变量,不是具体的数"
x = sp.Symbol('x')
y = sp.Symbol('y')

# 现在可以构建表达式了
expr = x**2 - 2*x - 1  # y = x² - 2x - 1
print(expr)  # 输出: x**2 - 2*x - 1

核心操作:代入求值 .subs()

有了表达式,我们可以轻松计算任意 x 对应的 y 值:

# 计算 x=1.5 时的函数值
result = expr.subs(x, 1.5)
print(result)  # 输出: -1.75000000000000
print(float(result))  # 转换为浮点数: -1.75

自动求导和解方程

这才是真正解放双手的功能:

# 求导数
derivative = sp.diff(expr, x)  # 对 x 求导
print(derivative)  # 输出: 2*x - 2

# 解方程:导数=0(找顶点)
vertex_x = sp.solve(derivative, x)[0]  # 解得 x=1
vertex_y = expr.subs(x, vertex_x)      # 代入求 y
print(f"顶点坐标: ({vertex_x}, {vertex_y})")  # 输出: (1, -2)

# 解方程:y=0(找与x轴交点)
roots = sp.solve(expr, x)
print(f"与x轴交点: {roots}")  # 输出: [1 - sqrt(2), 1 + sqrt(2)]

看到了吗? 原本需要手动计算的所有值,现在 SymPy 全自动搞定了!

SymPy 和 Manim 结合示例

现在我们把 SymPyManim 结合起来,做一个参数可调的抛物线动画。

核心代码示例

from manim import *
import sympy as sp

class AutoParabola(Scene):
    def construct(self):
        # ========== SymPy 自动计算部分 ==========
        x = sp.Symbol('x')
        a, b, c = 1, -2, -1                      # 抛物线参数:y = ax² + bx + c
        expr = a * x**2 + b * x + c              # SymPy 符号表达式

        # 自动求顶点:令导数为 0
        derivative = sp.diff(expr, x)            # 求导:2ax + b
        vertex_x = float(sp.solve(derivative, x)[0])
        vertex_y = float(expr.subs(x, vertex_x))

        # 自动求与 x 轴交点
        roots = sp.solve(expr, x)                # 解方程 ax² + bx + c = 0
        root_points = [(float(r), 0) for r in roots if r.is_real]
        
        # ========== Manim 可视化部分 ==========
        axes = Axes(x_range=[-2, 4, 1], y_range=[-3, 3, 1], axis_config={"color": BLUE})
        
        # 用 SymPy 表达式直接作为绘图函数
        parabola = axes.plot(
            lambda x_val: float(expr.subs(x, x_val)),  # SymPy 实时计算 y 值
            color=YELLOW, stroke_width=3,
        )
        
        # 顶点(使用 SymPy 算出的坐标)
        vertex_dot = Dot(axes.c2p(vertex_x, vertex_y), color=RED)
        vertex_label = MathTex(
            f"({vertex_x:.1f}, {vertex_y:.1f})", font_size=24, color=RED
        ).next_to(vertex_dot, UP)

        # x 轴交点
        root_dots = VGroup(*[
            Dot(axes.c2p(rx, ry), color=GREEN) for rx, ry in root_points
        ])

        # ========== 动画播放 ==========
        self.play(Create(axes))
        self.play(Create(parabola))
        self.play(Create(vertex_dot), Write(vertex_label))
        self.play(Create(root_dots))
        self.wait(2)

代码核心解析

关键点1:无缝衔接

# SymPy 计算出的值是符号类型,需要转为 float 给 Manim 使用
vertex_x = float(sp.solve(derivative, x)[0])

关键点2:动态函数映射

# 用 lambda 将 SymPy 表达式"翻译"成 Manim 能理解的数值函数
parabola = axes.plot(
    lambda x_val: float(expr.subs(x, x_val)),
    color=YELLOW,
)

关键点3:坐标系转换

# 数学坐标 → 屏幕坐标
vertex_dot = Dot(axes.c2p(vertex_x, vertex_y))

效果展示说明

运行这段代码后,你会看到:

  1. 坐标轴自动建立,范围根据函数特点自适应
  2. 抛物线精确绘制,形状由 SymPy 实时计算
  3. 红色顶点自动标注在正确位置,坐标值精确显示
  4. 绿色交点标记出抛物线与 x 轴的交点
  5. 白色虚线标出对称轴位置

最神奇的是:如果你想改成 $ y = 2x^2 + 3x - 1 $,只需要修改第 11 行的参数:

a, b, c = 2, 3, -1  # 其他代码完全不用动!

所有的点、线、标注都会自动更新到正确位置!修改后:

小结

今天我们解决了 Manim 动画制作中的一大痛点:手动计算坐标

通过 SymPy 的符号计算能力,我们实现了:

  • 表达式精确计算:告别计算器
  • 适合数学思维表达:将公式推导直接映射成代码
  • 自动求导找顶点:告别手算求导
  • 自动解方程找交点:告别求根公式
  • ✅ 等等... ...

核心代码模板

import sympy as sp
x = sp.Symbol('x')
expr = x**2 - 2*x - 1  # 你的表达式
y_value = float(expr.subs(x, x_value))  # 计算任意点的值