当前位置：首页 > news >正文

轨迹的蓝图：方程求解与交点计算

news 2026/6/29 19:30:30

做数学动画时，我经常遇到这样一个问题：

辛辛苦苦画好了两条曲线，y = sin(x)和y = x/2，想让它们交点的位置亮起一个发光点，标注出坐标。结果发现：我根本不知道交点到底在哪。

你可能会怎么做？

方案 A：用NumPy生成一堆点，然后暴力遍历找最接近的位置。
方案 B：打开Desmos或GeoGebra，估计交点坐标，再手动写进代码里。

这两种方法我都试过，结果总是不尽人意。要么精度不够，动画一放大就露馅；

要么过程繁琐，每次改函数都得重新估算。更别提当函数变得更复杂时，比如加入指数或对数项，肉眼基本就失效了。

有没有一种方法，能让 Python自动算出精确的解析解，或者至少给出一个高精度的数值解，

然后直接喂给Manim？还真有，就是SymPy。

1. 痛点场景还原

想象你要做一个动画：展示正弦曲线y = sin(x)和直线y = x/2的交点，并在交点处放置一个脉动的光点。

靠肉眼看，交点大概在x ≈ 1.895附近。但这显然不是一个精确值。即便你用numpy的数值方法去逼近，也需要手动写寻根算法，步骤繁琐且不够优雅。

# 以前你可能这样干（甚至更原始） import numpy as np # 手动定义寻根函数，还要选初始猜测值 def f(x): return np.sin(x) - x/2 # 需要手动调用 root 求解器，还有可能不收敛 from scipy.optimize import root result = root(f, 1.8) # 初始猜测还得靠蒙 print(result.x) # 输出 [1.89549427]

这种方式的问题在于：

流程割裂：需要在一个环境里求解，再手动复制数字到Manim脚本里
不精确：数值解有精度损失，且解析解（如果存在）被完全忽略
不可维护：换个函数就要重新跑一遍流程

SymPy 可以优雅地解决这一切，并且能和 Manim 无缝衔接。

2. SymPy 解决方案介绍

SymPy可以像在纸上推公式一样处理数学表达式。

对于交点问题，它提供了三大利器：

方法	适用场景	特点
`solve(Eq, x)`	有解析解的方程	返回精确表达式
`solveset()`	复杂方程，尤其是周期函数	返回解集，支持无穷多解
`nsolve()`	没有解析解的情况	数值求解，需要初始猜测

2.1 符号求解：solve 和 solveset

对于sin(x) = x/2这种超越方程，没有初等函数形式的解析解。但我们可以用 SymPy 验证这一点，并获取符号形式的解集表示：

import sympy as sp # 定义符号变量 x = sp.Symbol('x') # 定义方程 sin(x) = x/2 equation = sp.Eq(sp.sin(x), x/2) # 尝试求解（对于超越方程，可能返回未求值的形式） symbolic_solutions = sp.solve(equation, x) print(symbolic_solutions) # 因为没有解析解，返回空列表或用数值表示 # 使用 solveset 可以更好地处理这类方程 solution_set = sp.solveset(equation, x, domain=sp.S.Reals) print(solution_set) # 返回一个包含数值解的条件集合

对于有解析解的方程，solve就非常强大了：

# 有解析解的例子：x^2 - 5x + 6 = 0 equation_quad = sp.Eq(x**2 - 5*x + 6, 0) exact_roots = sp.solve(equation_quad, x) print(exact_roots) # 输出 [2, 3] —— 精确解！

2.2 数值求解：nsolve

对于没有解析解的超越方程，nsolve是真正的主力。它基于牛顿法等高精度算法，只需一个初始猜测就能给出高精度数值解：

# 数值求解 sin(x) = x/2 # 初始猜测 1.8，因为从图像看交点大概在那个位置 numerical_solution = sp.nsolve(equation, 1.8) print(numerical_solution) # 输出 1.89549426703398 print(type(numerical_solution)) # <class 'mpmath.ctx_mp_python.mpf'>

关键点：nsolve返回的是mpmath的高精度浮点数，可以轻松转换成 Python 的float供 Manim 使用。

# 转换成普通浮点数 x_intersection = float(numerical_solution) y_intersection = float(sp.sin(numerical_solution)) print(f"交点坐标：({x_intersection:.6f}, {y_intersection:.6f})") # 输出：交点坐标：(1.895494, 0.947747)

3. Manim 联动实战

现在我们把 SymPy 的计算结果直接传递给 Manim，实现“计算即所得”的动画效果。

下面是核心代码部分：

from manim import * import sympy as sp class IntersectionDemo(Scene): def construct(self): # ========== SymPy 数值求解交点 ========== x = sp.Symbol('x') f_expr = sp.sin(x) # 曲线1：y = sin(x) g_expr = x / 2 # 曲线2：y = x/2 equation = sp.Eq(f_expr, g_expr) # 方程：sin(x) = x/2 # nsolve：牛顿法数值求解（超越方程无解析解） x_sol = float(sp.nsolve(equation, 1.8)) # 初始猜测 1.8 y_sol = float(f_expr.subs(x, x_sol)) # 代入求 y # 符号表达式 → Python 函数（供绘图） f = lambda t: float(f_expr.subs(x, t)) g = lambda t: float(g_expr.subs(x, t)) # ========== Manim 可视化 ========== axes = Axes( x_range=[-5, 5, 1], y_range=[-2, 2, 0.5], x_length=8, y_length=5, axis_config={"color": BLUE} ) sin_graph = axes.plot(f, color=YELLOW, stroke_width=3) # sin 曲线 line_graph = axes.plot(g, color=GREEN, stroke_width=3) # 直线 # 交点（使用 SymPy 算出的精确坐标） dot = Dot(axes.c2p(x_sol, y_sol), color=RED, radius=0.1) dot.set_z_index(10) label = MathTex(f"({x_sol:.2f}, {y_sol:.2f})", font_size=30, color=RED).next_to(dot, UR) # ========== 动画播放 ========== self.play(Create(axes)) self.play(Create(sin_graph), Create(line_graph)) self.play(GrowFromCenter(dot), Write(label)) self.wait(1)