聚贤 發表於 2026-4-16 14:56:00

用Manim实现动态交点计算--从一个动点问题说起

<p>大家好,今天想和大家分享一个在制作<code>Manim</code>动画时非常实用的话题:如何<strong>动态计算两条直线的交点</strong>。</p>
<p>对于动点问题,比如初中数学中经典的<strong>“时钟模型”</strong>或<strong>“将军饮马”</strong>及其变种等等,硬编码坐标肯定不行,因为交点坐标是随动点变化的。</p>
<p>下面,我们结合 <code>Python</code> 的符号计算库 <code>SymPy</code> 和 <code>Manim</code> 的更新器(<code>Updater</code>),来实现真正的动态交点计算。</p>
<h1 id="1-从一个初中题目说起">1. 从一个初中题目说起</h1>
<p>先来看一道经典的初中动点问题:</p>

<p><img src="https://img2024.cnblogs.com/blog/83005/202604/83005-20260416145553478-1968406202.png" alt="" loading="lazy"></p>
<p>这个题目中,当点$ E <span class="math inline">\(和\)</span> F <span class="math inline">\(在线段\)</span> AD <span class="math inline">\(上移动时,点\)</span> H <span class="math inline">\(和\)</span> G $都随之变换。</p>
<p>题目如何解答我们不用管,我们的重点是如何实现点$ E <span class="math inline">\(和\)</span> F <span class="math inline">\(移动时,实时的更新点\)</span> H <span class="math inline">\(和\)</span> G $。</p>
<h1 id="2-解决方案用sympy解方程组">2. 解决方案:用Sympy解方程组</h1>
<p>这里简单说明下,<code>Sympy</code>是一个<code>Python</code>的符号计算库,可以像数学课本那样进行代数运算。</p>
<p>它能<strong>解方程</strong>、<strong>求导积分</strong>、<strong>化简表达式</strong>,最重要的是能<strong>精确求解方程组</strong>,返回的是精确的数学解而不是近似值。</p>
<p>我们实现这个动画效果时,是根据两个直线的方程来求解它的交点的,如果没有<code>Sympy</code>,我们要手动推导直线方程、联立求解,代码会非常冗长且容易出错。</p>
<p>从后面的代码你可以看出,<code>Sympy</code>让我们只需"翻译"数学表达式,就能得到精确结果,大大简化了代码逻辑。</p>
<p>我们的思路很简单:</p>
<ul>
<li>已知<strong>两个点</strong>的坐标,可以求出<strong>直线方程</strong></li>
<li>已知两条<strong>直线方程</strong>,联立求解得到<strong>交点坐标</strong></li>
</ul>
<p>在<code>Python</code>中,用<code>Sympy</code>这个符号计算库来实现再合适不过了。</p>
<h2 id="21-第一步定义求直线方程的函数">2.1. 第一步:定义求直线方程的函数</h2>
<pre><code class="language-python">from sympy import Symbol, solve

def get_line(p1, p2):
    """已知两点,求直线y = kx + b的k和b"""
    k = Symbol("k")
    b = Symbol("b")
   
    expr1 = p1 * k + b - p1
    expr2 = p2 * k + b - p2
   
    ret = solve((expr1, expr2), dict=True)
    return {"k": ret, "b": ret}
</code></pre>
<h2 id="22-第二步定义求交点的函数">2.2. 第二步:定义求交点的函数</h2>
<pre><code class="language-python">def cross_points(l1, l2):
    """已知两条直线方程,求交点坐标"""
    x = Symbol("x")
    y = Symbol("y")
   
    expr1 = l1["k"] * x + l1["b"] - y
    expr2 = l2["k"] * x + l2["b"] - y
    ret = solve((expr1, expr2), dict=True)
   
    return np.array((float(ret), float(ret), 0))
</code></pre>
<p>有了这两个函数,我们就可以在<code>Manim</code>中动态更新交点了。</p>
<h2 id="23-完整的manim动画实现">2.3. 完整的Manim动画实现</h2>
<pre><code class="language-python">from manim import *
import numpy as np
from sympy import Symbol, solve

class DynamicCrossPoint(Scene):
    def construct(self):
      # 定义矩形的顶点坐标
      points = {
            "A": np.array([-2.5, 2, 0]),
            "B": np.array([-2.5, -3, 0]),
            "C": np.array(),
            "D": np.array(),
      }
      
      # 初始动点的位置
      points["E"] = np.array([-0.52, 2, 0])   # E在AB上
      points["F"] = np.array()    # F在CD上,且AE=CF
      
      # 画矩形
      rectangle = Polygon(
            points["A"], points["B"],
            points["C"], points["D"],
            stroke_width=3, color=GREEN
      )
      self.play(Create(rectangle))
      
      # 创建初始的点和线
      d_e = Dot(points["E"], radius=0.05, color=BLUE)
      d_f = Dot(points["F"], radius=0.05, color=BLUE)
      d_h = Dot(points["A"], radius=0.05, color=YELLOW)# 初始随便放
      
      l_bf = Line(points["B"], points["E"], color=BLUE, stroke_width=2)
      l_ce = Line(points["C"], points["F"], color=BLUE, stroke_width=2)
      l_bd = Line(points["B"], points["D"], color=GREEN, stroke_width=2)
      
      self.play(Create(VGroup(l_bf, l_ce, l_bd, d_e, d_f, d_h)))
      
      # 核心部分:设置更新器
      # F点随着E点移动(保持AE=CF的关系)
      d_f.add_updater(
            lambda z: z.become(
                Dot(points["D"] - (d_e.get_center() - points["A"]),
                  radius=0.05, color=BLUE)
            )
      )
      
      # H点是BF和CE的交点
      d_h.add_updater(
            lambda z: z.become(
                Dot(
                  cross_points(
                        get_line(points["B"], d_e.get_center()),# BF
                        get_line(points["C"], d_f.get_center()),# CE
                  ),
                  radius=0.05, color=YELLOW
                )
            )
      )
      
      # 更新线段
      l_bf.add_updater(
            lambda z: z.become(
                Line(points["B"], d_e.get_center(), color=BLUE, stroke_width=2)
            )
      )
      
      l_ce.add_updater(
            lambda z: z.become(
                Line(points["C"], d_f.get_center(), color=BLUE, stroke_width=2)
            )
      )
      
      # 让E点动起来
      self.play(d_e.animate.shift(LEFT * 1.5), run_time=3)
      self.play(d_e.animate.shift(RIGHT * 3), run_time=6)
      self.play(d_e.animate.shift(LEFT * 1.5), run_time=3)
      
      # 清理更新器
      for mob in :
            mob.clear_updaters()
      
      self.wait()
</code></pre>

<p><img src="https://img2024.cnblogs.com/blog/83005/202604/83005-20260416145553527-1128468260.gif" alt="" loading="lazy"></p>
<h2 id="24-为什么不直接用manim的几何变换">2.4. 为什么不直接用Manim的几何变换?</h2>
<p>可能有朋友会问:<code>Manim</code>不是有<code>.move_to()</code>、<code>.shift()</code>这些方法吗?为什么非要用<code>Sympy</code>算交点?</p>
<p><strong>答案</strong>是:当动点之间没有简单的几何变换关系时,比如,在这个例子中,交点的位置是由两条动态直线决定的,没有办法通过简单的位移或旋转得到。</p>
<p>这时,数学计算就是最直接的方法。</p>
<h2 id="25-进阶思考">2.5. 进阶思考</h2>
<p>上面介绍的方法,通用性很强:</p>
<ul>
<li>可以用来求<strong>直线与圆</strong>的交点</li>
<li>可以用来求<strong>两条曲线</strong>的交点</li>
<li>甚至可以用来求<strong>动点轨迹</strong>的方程</li>
</ul>
<p>只要你能写出方程,<code>Sympy</code>就能帮你解出来。</p>
<h1 id="3-结语">3. 结语</h1>
<p>今天这个例子,我们学会了:</p>
<ol>
<li>用<code>Sympy</code>求解直线交点</li>
<li>在<code>Manim</code>中动态更新交点</li>
<li>处理多个相互关联的动点</li>
</ol>
<p>希望这篇文章能帮助你在制作<code>Manim</code>动画时,更自如地处理动态几何问题。</p><br><br>
来源:https://www.cnblogs.com/wang_yb/p/19877766
頁: [1]
查看完整版本: 用Manim实现动态交点计算--从一个动点问题说起