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

推荐订阅源

K
Kaspersky official blog
Engineering at Meta
Engineering at Meta
D
DataBreaches.Net
Stack Overflow Blog
Stack Overflow Blog
Microsoft Security Blog
Microsoft Security Blog
Y
Y Combinator Blog
B
Blog RSS Feed
GbyAI
GbyAI
P
Proofpoint News Feed
aimingoo的专栏
aimingoo的专栏
MyScale Blog
MyScale Blog
D
Docker
阮一峰的网络日志
阮一峰的网络日志
Cyber Security Advisories - MS-ISAC
Cyber Security Advisories - MS-ISAC
Recorded Future
Recorded Future
美团技术团队
The Register - Security
The Register - Security
V
Visual Studio Blog
H
Hackread – Cybersecurity News, Data Breaches, AI and More
T
Tailwind CSS Blog
爱范儿
爱范儿
freeCodeCamp Programming Tutorials: Python, JavaScript, Git & More
T
The Blog of Author Tim Ferriss
博客园 - 司徒正美
量子位
B
Blog
F
Fortinet All Blogs
Martin Fowler
Martin Fowler
博客园 - 【当耐特】
MongoDB | Blog
MongoDB | Blog
A
About on SuperTechFans
I
InfoQ
钛媒体:引领未来商业与生活新知
钛媒体:引领未来商业与生活新知
有赞技术团队
有赞技术团队
雷峰网
雷峰网
大猫的无限游戏
大猫的无限游戏
J
Java Code Geeks
L
LangChain Blog
Latest news
Latest news
S
SegmentFault 最新的问题
Exploit-DB.com RSS Feed
Exploit-DB.com RSS Feed
Blog — PlanetScale
Blog — PlanetScale
Threat Intelligence Blog | Flashpoint
Threat Intelligence Blog | Flashpoint
Cisco Talos Blog
Cisco Talos Blog
F
Full Disclosure
C
Cisco Blogs
D
Darknet – Hacking Tools, Hacker News & Cyber Security
W
WeLiveSecurity
T
Tenable Blog
T
Tor Project blog

博客园_首页

Linux实操--组管理、权限管理和定时任务 Java + EasyExcel 实现单个接口导出多个Excel Mem0 源码解析系列(二):提示词工程的深度剖析 Openclaw TaskFlow究竟是什么?和普通Skill技能有什么区别 博文阅读密码验证 - 博客园 嘉立创开源:应该是全网MicroPython教程最多的开发板 Hermes Agent 集成实践:从协议到生产 2026年AI编程工具横评:Cursor、Codex、Claude Code、Zed、Windsurf Java程序员必看的RAG入门教程 2026 AI效率神器:Superpowers + Claude Code 保姆级教程 本地大模型部署全攻略:从 0 到 1 玩转 Ollama 【从0到1构建一个ClaudeAgent】内存管理-上下文压缩 .NET 高级开发 | 设计、实现一个事件总线框架 电子小白入门之NE555 3. WorkBuddy:隐藏玩法,一键召唤专家,让 AI 以"专家身份"给你干活 和AI一起搞事情#3:Claude Teammate 游戏开发翻车实录 【OpenClaw】通过 Nanobot 源码学习架构---(7)Memory C# .NET 周刊|2026年3月3期 我在 Debian 11 上把 K8s 单机搭起来了,过程没你想的那么顺(/opt 目录版) 深度学习进阶(七)Data-efficient Image Transformer CLI+Skill搭建浏览器AI自动化框架,告别一切重复枯燥任务 告别Token账单无底洞:OpenClaw本地部署,重塑企业数据主权的唯一解 FastAPI+Vue:文件分片上传+秒传+断点续传,这坑我帮你踩平了! SBTI 爆火后,我做了个程序员版的 CBTI。。已开源 + 附开发过程 多模态检索开始进入工程期:用 Sentence Transformers 搭建可落地的 Multimodal RAG 100多行代码实现一个最简单的Agent(用ReAct) Claude Code 通关手册(八):推荐 5 个 Hooks,代码质量提升 3 倍 老板:“有人截图了!”。安全部门:“收到,马上查暗水印!” - why技术 技术之外,皆是人间 C#/.NET/.NET Core技术前沿周刊 | 第 69 期(2026年4.01-4.12) Snack JSONPath 项目架构分析 Claude Code Buddy 小析:一个非核心功能,如何体现产品的细节完成度 AI新时代下的图床管理方案-Cloudflare图床+MCP+Skills方案指南 化繁为简:顺丰速运App如何通过 HarmonyOS SDK实现专业级空间测量 从零实现富文本编辑器#13-React非编辑节点的内容渲染 AI开发-python-langchain框架(3-23-OpenAI Functions风格Tool Calling智能助手) .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 PbootCMS 网站内容数量多导致访问慢?这些实用优化方案帮你提速! - 家兴网络技术工作室 上下文工程是什么?过时了么?一文讲明白! - 一枫说码 网站漏洞怎么发现并修复?一篇实用指南(附完整流程) - 家兴网络技术工作室 开了 TUN 模式还是直连?90% 的人都踩过这个坑 Github日报|2026年04月12日 - AI一族 AScript扩展多种脚本语言 - rockey627 AI 学习笔记:Agent 的记忆机制 你能被装进一个文件里吗?——7 万人把同事"蒸馏"成了 AI - 我没有三颗心脏 Claude Code 通关手册(七):给 AI 装上技能包——Skills 完全指南 - 暮色之狐 在浏览器中快速编辑代码:VSCode Web 集成实践 - Newbe36524 蒸馏自己 skill?基于 Deepseek 的蒸馏器,丐版蒸馏方式,简单便捷 - To_Carpe_Diem Spring AI Aliababa和AgentScope,哪个更好? - 苏三说技术 Etsy 把 1000 个 MySQL 分片迁进 Vitess:425TB 数据背后的真正问题不是性能,而是运维规模 MicroPython LVGL基础知识和概念:底层渲染与性能优化 - FreakStudio 数据库草图算法 Python 潮流周刊#146:CPython 引入 Rust 的进展 - 豌豆花下猫 最小生成树 - mofei1116 红日靶场七:从外网入口、容器逃逸到 AD 接管的完整利用链复盘 - YouDiscovered1t 分享四款开源且实用的 Kafka 管理工具 - 追逐时光者 vLLM 权重加载机制全解析:从挑战到理想架构 LCT 学习笔记 - ACehomoxue Avalonia UI 12.0.0 正式发布:架构演进和性能飞跃 - 张善友 当 AI Agent 把调用链拉长,延迟开始成为一门生意 conhost.exe 无法显示 U+2717 - 145a 太秀了,我把自己蒸馏成了 Skill!已开源 - 程序员鱼皮 ASP.NET Core 内存缓存实战:一篇搞懂该怎么配、怎么避坑 基于 Ghostty 带有分割标签页和为 Claude 编程设计的通知终端 - BugShare AI 焊死入口:教育的“操作系统级”重塑 - 郝hai 初级Java开发工程师使用sql脚本编写代码的过程是简单而且不糊涂 - CoderOilStation Claude Code通关手册(六):MCP协议完全指南 - 暮色之狐 边框灯光环绕动画特效实现指南 - Newbe36524 开源:子木蒸馏版的 SEO 审计工具 seo-audit-skill v1.0 我所理解的Python元模型 【从0到1构建一个ClaudeAgent】规划与协调-TodoWrite - 程序员Seven Claude 和 Codex 在审计 Skill 上性能差异探究 - ACai_sec AScript如何实现中文脚本引擎 - rockey627 【渗透测试】HTB Season10 Garfield 全过程wp - dynasty_chenzi Android 开发者为什么必须掌握 AI 能力?端侧视角下的技术变革 树状数组正确性证明 - AC-wyr 你的 AI 焦虑,可能比 AI 本身更危险——ATM 机没有消灭银行柜员,但恐慌消灭了你的判断力 - 我没有三颗心脏 一个拉胯的分库分表方案有多绝望?整个部门都在救火! - 冰河团队 动态规划入门必学之走方格问题 - Ofnoname PostgREST 与 PostgreSQL 角色权限配置全解析(生产级实践) - SheepDog1998 使用 UEFI 图形输出协议 GOP 在屏幕上显示图像的方法 - 阿源- Claude Code通关手册(五):组建你的AI专家团队,子代理系统 - 暮色之狐 一个程序员到架构师的催婚路之感悟(整整10年后的催婚相亲感悟) - MisterLip 用 Agent Skill 自动生成工作周报 - 赵康
告别手动计算,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))  # 计算任意点的值