用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]