使用Manim绘制secp256k1椭圆曲线时，点偏离曲线的问题求助

核心问题定位 #

你遇到的点不在曲线上的问题，本质是混淆了「实数域椭圆曲线」和「有限域secp256k1椭圆曲线」的差异，或者代码中存在坐标计算/绘制范围的错误，具体分两种场景：

场景1：你想绘制实数域的x³+7=y²曲线 #

排查步骤

1. 检查点的坐标计算：x=1时，y²=1³+7=8，因此y必须是±√8≈±2.828，而非硬编码的整数或其他错误值。
2. 确认曲线绘制范围：如果Manim的y_range设置过小（比如[-2,2]），会导致曲线未显示到±2.828的位置，看起来点不在曲线上。
3. 修正后的示例代码：

from manim import *
import math

class RealEllipticCurve(Scene):
    def construct(self):
        # 定义实数域椭圆曲线方程
        curve_eq = lambda x, y: x**3 + 7 - y**2
        # 设置包含目标点的坐标范围
        curve = ImplicitFunction(
            curve_eq,
            x_range=[-5, 5],
            y_range=[-4, 4],
            color=BLUE
        )
        self.add(curve)
        
        # 计算x=1对应的正确y值
        x_val = 1
        y_squared = x_val**3 +7
        y1 = math.sqrt(y_squared)
        y2 = -math.sqrt(y_squared)
        
        # 添加曲线上的点
        point1 = Dot(point=(x_val, y1, 0), color=RED, radius=0.08)
        point2 = Dot(point=(x_val, y2, 0), color=RED, radius=0.08)
        self.add(point1, point2)
        
        # 添加坐标轴与标签
        axes = Axes(
            x_range=[-5,5,1],
            y_range=[-4,4,1],
            axis_config={"include_numbers": True}
        )
        self.add(axes)
        self.add(Text(f"({x_val}, {y1:.2f})", font_size=20).next_to(point1, UP))
        self.add(Text(f"({x_val}, {y2:.2f})", font_size=20).next_to(point2, DOWN))

场景2：你想绘制有限域GF(p)上的secp256k1曲线 #

secp256k1的定义域是有限域GF(p)（p=2^256-2^32-977），所有运算都是模p的，和实数域曲线完全无关：

1. 有限域中点的计算：x=1时，需要找到满足y² ≡ 8 mod p的y值，这是一个超大整数，直接放到实数坐标系里肯定不在实数曲线上。
2. 可视化方式：有限域的点是离散的，不能用ImplicitFunction绘制连续曲线，需要先计算点坐标，再通过映射缩小到可视范围。
3. 示例代码（简化版）：

from manim import *

# secp256k1的有限域参数
p = 2**256 - 2**32 - 977

def mod_sqrt(a, p):
    # 针对p≡3 mod4的Tonelli-Shanks简化算法（secp256k1满足该条件）
    assert pow(a, (p-1)//2, p) == 1, "该x值在GF(p)中无对应y值"
    return pow(a, (p+1)//4, p)

class Secp256k1FiniteField(Scene):
    def construct(self):
        # 计算GF(p)中x=1对应的y值
        x = 1
        y_sq = (x**3 +7) % p
        y1 = mod_sqrt(y_sq, p)
        y2 = (-y1) % p
        
        # 将超大整数映射到[-5,5]的可视范围
        def map_coord(val):
            scaled = (val % 2**32) / (2**32) * 10 -5
            return scaled
        
        mapped_x = map_coord(x)
        mapped_y1 = map_coord(y1)
        mapped_y2 = map_coord(y2)
        
        # 添加映射后的点
        point1 = Dot(point=(mapped_x, mapped_y1, 0), color=RED)
        point2 = Dot(point=(mapped_x, mapped_y2, 0), color=RED)
        self.add(point1, point2)
        
        # 添加坐标轴
        axes = Axes(x_range=[-5,5], y_range=[-5,5], axis_config={"include_numbers": True})
        self.add(axes)

总结排查流程 #

1. 先明确你要绘制的是实数域曲线还是有限域secp256k1曲线，两者逻辑完全不同。
2. 实数域场景：重点检查点的坐标计算是否正确、曲线绘制范围是否覆盖目标点。
3. 有限域场景：必须用模运算计算点坐标，再通过映射缩小到可视范围，不能用连续曲线展示。

内容的提问来源于stack exchange，提问作者Radim Bača