欧拉极限与公理跳跃：计算实验报告
本实验报告旨在通过四个定量化的数值实验，将您的三篇文章（从圆周率到曲率、形式系统，再到可计算物理的递进思想）转化为可视化的实验语言。实验锚定欧拉公式与恒等式，向您展示无理数与超越数如何在一项项可数无穷的计算与定义跳跃（Axiom Jump）中逐步展开。
实验一：割圆曲率与“代数税”的 2:1 渐近法则
1. 实验目的与背景（关联第一弹）
第一篇文章指出，当我们用离散的折线（多边形）去逼近连续的圆时，多边形的曲率并不是均匀的，而是在顶点处呈狄拉克 $\delta$ 函数奇异性凝聚，在边上为 0。这一拓扑深渊导致了“代数税”——即逼近误差的产生。本实验通过计算内接与外切正 $N$ 边形的面积误差，验证其二次收敛律与渐近误差比。
2. 数学模型
设单位圆半径为 $R=1$，真实面积为 $\pi$。
- 内接正 $N$ 边形面积为：$S_{in}(N) = \frac{N}{2} \sin\left(\frac{2\pi}{N}\right)$，面积误差 $E_{in}(N) = \pi - S_{in}(N)$。
- 外切正 $N$ 边形面积为：$S_{out}(N) = N \tan\left(\frac{\pi}{N}\right)$，面积误差 $E_{out}(N) = S_{out}(N) - \pi$。
- 渐近分析预测，当 $N \to \infty$ 时：
$$E_{in}(N) \approx \frac{2\pi^3}{3N^2}$$
$$E_{out}(N) \approx \frac{\pi^3}{3N^2}$$
两者之比 $\frac{E_{in}(N)}{E_{out}(N)} \to 2$。
3. 实验数据与结果
运行多边形边数 $N$ 从 4 到 1024。对误差进行对数拟合（$\log E = \log a + b \log N$）：
- 内接误差拟合斜率：$-1.9866$（理论值 $-2$）
- 外切误差拟合斜率：$-2.0297$（理论值 $-2$）
- 极限误差比（当 $N=1024$ 时）：$1.99999$（完美收敛于 $2$）
4. 图像分析
- 左图展示了离散多边形顶点处的红色“曲率刺针” ($\delta$ 奇异性)。即使 $N \to \infty$，这些刺针也无法平滑地融入圆的均匀曲率 $1$，这在物理上代表了维度的不连续，表明了“圆”必须通过一次公理跳跃（定义极限）来创造。
- 右图展示了在对数双轴下误差完美呈现斜率为 $-2$ 的直线，证明了曲率在局部空间呈二次展宽的代数必然性（代数税）。
实验二：欧拉公式在复平面的“泰勒级数行走”与圆的融化
1. 实验目的与背景（关联第一、三弹）
欧拉恒等式 $e^{i\pi} + 1 = 0$ 是将刚性的勾股定理在连续、旋转的复平面中进行大一统的杰作。泰勒级数拥有无限多项，这正是因为三角函数在有曲率变化的连续环境中运行。本实验通过可视化 $e^{i\pi}$ 泰勒展开路径在复平面上的“向量行走”，展示离散的有理数折线段如何通过可数无穷项的求和，一步步“融化”为平滑的圆，并收敛到超越数坐标 $-1$。
2. 数学模型
对 $e^{i\pi}$ 进行泰勒级数求和：
$$z_M = \sum_{n=0}^{M} \frac{(i\pi)^n}{n!}$$
我们跟踪前 $M=20$ 项的复数部分 $z_M$ 及其模长与目标值 $1$（即勾股定理 $x^2 + y^2 = 1$ 的守恒性）的偏差。
3. 实验结果
- 当项数 $M=20$ 时，计算得到的复数坐标与理论目标 $-1 + 0i$ 的绝对误差仅为 $3.568 \times 10^{-9}$。
- 在有限的 $M$ 步中，每一次折线移动的端点坐标都是代数数（有理数的代数运算结果），且其对应的三角形模长均不是完美的 $1$。
4. 图像分析
- 左图清晰展示了向量在复平面上的“直角转弯折线”（每一次转弯都是 $i$ 乘以实数，呈 $90^\circ$ 刚性直角）。随着项数 $M$ 增加，折线形成一个逐渐收缩的螺旋，最终在 $M \ge 10$ 时几乎融化为完美的圆弧，精准锁定在 $-1$ 点。
- 右图显示了级数项数与绝对误差的半对数图，误差呈超指数（阶乘分母级）的断崖式下跌，生动展示了“可数无穷”计算如何填平有理与无理的深渊。
实验三：欧拉常数 $\gamma$ 的“测度夹缝”几何叠加
1. 实验目的与背景（关联第一弹）
欧拉-马歇罗尼常数 $\gamma$ 诞生于离散求和（调和级数，阶梯）与连续积分（对数曲线，斜坡）的级数偏差中。它是一个由于离散与连续碰撞而残留的几何测度。本实验展示调和级数与对数在 $N \to \infty$ 下的偏差收敛，并将这无穷多个“局部弯角余项”平移叠放到单位正方形内，可视化常数 $\gamma$ 的实体形态。
2. 数学模型
$$\gamma = \lim_{N \to \infty} \left( \sum_{k=1}^N \frac{1}{k} - \ln N \right)$$
在单位区间 $x \in [0, 1]$ 内，第 $k$ 个偏差余项可以表示为函数 $Y_k(x)$ 与 $Y_{k-1}(x)$ 之间的差：
$$Y_k(x) = Y_{k-1}(x) + \left( \frac{1}{k} - \frac{1}{x+k} \right)$$
这些层层叠加的阴影区域总面积即为 $\gamma$。
3. 实验结果
- 计算至 $N=1000$ 时，$\gamma$ 的数值逼近值为 $0.577716$，与理论真实值 $0.577216$ 的绝对误差为 $0.0005$。
4. 图像分析
- 左图展示了随着调和级数项数增加，其差值快速抚平，最终紧紧贴着水平虚线 $\gamma$ 运行。
- 右图是极具独创性的几何叠加图。它将调和级数矩形在曲线 $1/x$ 上方的“离散突起弯角”依次搬移到单位区间 $[0, 1]$ 中进行垂直叠加（第 1 层到第 8 层）。这些参差不齐的有理数弯角经过无穷叠加，填满了虚线下的面积，其总面积收敛为 $\gamma$。这不仅是个极限过程，更是离散折线向连续曲线过渡时流下的“代数眼泪”。
实验四：丢番图逼近谱系与无理数测度的动力学映射
1. 实验目的与背景（关联第二弹）
第二篇文章探讨了无理数的“无理程度”并非一概而论。在丢番图指数的光谱中，黄金比例和代数无理数最难被有理数逼近，圆周率 $\pi$ 难度稍低，而李维尔数极易被有理数逼近。这一现象在双曲模块面的负曲率几何中有着完美的对应：轨道在尖点处的逃逸深度决定了无理数测度的大小。本实验计算这三类数 Continued Fraction（连分数）逼近的误差，可视化其数值光谱。
2. 数学模型
对于数 $x$，其连分数渐近分数为 $p_k / q_k$。我们根据定义检验丢番图误差特征：
$$\left| x - \frac{p_k}{q_k} \right| \sim \frac{1}{q_k^{\mu(x)}}$$
- 代数无理数（如 $\sqrt{2}$）：由罗斯定理可知，其无理数测度 $\mu = 2$。
- 圆周率 $\pi$：超越数，目前已知的无理数测度上限为 $7.103$。
- 李维尔数 $L$：$\mu = \infty$。
3. 实验数据与结果
计算这三个数前若干项渐近分数的误差 $\log_{10} |x - p_k/q_k|$ 与分母对数 $\log_{10} q_k$：
- $\sqrt{2}$ 逼近分母达到 $195,025$ 时，误差为 $9.29 \times 10^{-12}$。
- $\pi$ 逼近分母达到 $364,913$ 时，误差为 $1.61 \times 10^{-12}$。
- 李维尔数 在分母仅为 $1,000,000$ 左右时，由于阶乘级跳跃，其逼近误差在计算机的双精度浮点极限下瞬间归 $0$（达到不可区分的底线）。
4. 图像分析
- 蓝线 ($\sqrt{2}$) 紧紧贴着黑色虚线（$1/q^2$，Dirichlet 逼近底线）下行，表现出极强的反抗性，说明在双曲空间中，它的测地线紧紧缠绕在紧致区域内，绝不落入代表有理数的“尖点黑洞”。
- 绿线 ($\pi$) 下行趋势略微陡峭，表现出比 $\sqrt{2}$ 更弱的无理度。
- 红线（李维尔数）表现为断崖式的急速下坠。对应在双曲面上，它的测地线几乎一瞬间便深插进尖点无限深处，展现了无限逼近有理数网格的倾向，这也是其柯尔莫哥洛夫定义复杂度低、极度有序的几何表现。
结论：公理跳跃的认识论沉淀
这四个实验生动地呈现了您三篇文章中的核心哲学观点：
1. 连续不是现成的物理密度：当我们用勾股定理或割圆多边形来界定世界时，我们在曲率的边缘必然产生“代数税”的漏洞。
2. 无穷级数是计算的缝合线：泰勒展开与调和求和中，可数无穷的计算像一根根缝合线，将破碎、离散的有理数坐标拼接成了看似光滑连续的复平面圆弧。
3. 公理跳跃的符号封装：当计算面临斯贝克序列、蔡廷常数 $\Omega$ 的不平凡边界或超越极限时，系统并不在无穷的泥潭中死锁，而是通过做出一次次伟大的“公理跳跃”，强行定义 $\pi$、$\gamma$、$e$ 以及虚数单位 $i$。
4. 欧拉公式的大一统：欧拉恒等式将这些通过跳跃所诞生的符号算子在复平面的旋转中进行大一统。无理数与超越数在这套机制下，不仅仅是数轴上的点，更是离散观察者的心智与连续几何曲率交互时所雕刻出来的物理和算术映像。