割圆术几何推导与圆周率计算
时间: 2025-10-07 22:35:58
好的，我们来一步步推导割圆术的几何原理。
割圆术（刘徽方法）的核心思想是用圆内接正多边形来逼近圆，通过倍增边数来计算圆周率 π。
1. 基本设定
设圆的半径为 $ R $，圆内接正 $ n $ 边形的边长为 $ a_n $，面积（或周长）可以用来近似圆的面积（或周长）。
刘徽从正六边形开始，然后逐步倍增边数：
$ n = 6, 12, 24, 48, 96, \dots $
2. 正 $ n $ 边形边长到正 $ 2n $ 边形边长
已知正 $ n $ 边形边长 $ a_n $，求正 $ 2n $ 边形边长 $ a_{2n} $。
几何关系
- 圆心 $ O $，正 $ n $ 边形一边 $ AB = a_n $。
- $ M $ 为 $ AB $ 的中点，则 $ OM \perp AB $。
- $ MB = \frac{a_n}{2} $，$ OB = R $。
- 由勾股定理：
$OM = \sqrt{R^2 - \left( \frac{a_n}{2} \right)^2 }$
- 线段 $ MC = R - OM = R - \sqrt{R^2 - \left( \frac{a_n}{2} \right)^2 } $。
求 $ a_{2n} $
$ a_{2n} $ 是正 $ 2n $ 边形的边长，即 $ AM $ 或 $ MB $ 对应的弦长？
不对，正 $ 2n $ 边形顶点包括 $ A, M, B, \dots $ 吗？
实际上，正 $ 2n $ 边形在 $ A $ 与 $ B $ 之间插入中点 $ M $ 作为新顶点，所以 $ AM $ 是正 $ 2n $ 边形的边长。
在三角形 $ OAM $ 中：
- $ OA = R $
- $ OM $ 已知
- $ AM = \sqrt{ (OA)^2 - (OM)^2 } $？ 不对，这样不对，因为 $ AM $ 是直角边。
更直接的方法：在直角三角形 $ OMB $ 中：
$MB = \frac{a_n}{2}, \quad OM = \sqrt{R^2 - \left( \frac{a_n}{2} \right)^2 }$
在直角三角形 $ AMB' $ 中？ 我们直接看三角形 $ AMC $ 吗？ 不，我们看三角形 $ AMO $ 不是直角三角形。
正确推导（用弦长公式）：
在圆上取 $ A, B $ 为 $ n $ 边形的相邻顶点，$ M $ 是弧 $ AB $ 的中点，则 $ AM $ 是正 $ 2n $ 边形的边长。
在等腰三角形 $ OAB $ 中：
$AB = a_n, \quad \angle AOB = \frac{2\pi}{n}.$
由余弦定理：
$a_n^2 = R^2 + R^2 - 2R^2 \cos\frac{2\pi}{n} = 2R^2 \left( 1 - \cos\frac{2\pi}{n} \right).$
利用半角公式：
$1 - \cos\frac{2\pi}{n} = 2\sin^2\frac{\pi}{n}.$
所以：
$a_n = 2R \sin\frac{\pi}{n}.$
同理，正 $ 2n $ 边形边长：
$a_{2n} = 2R \sin\frac{\pi}{2n}.$
用 $ a_n $ 表示 $ a_{2n} $
已知：
$a_n = 2R \sin\frac{\pi}{n}, \quad a_{2n} = 2R \sin\frac{\pi}{2n}.$
利用公式：
$\sin\frac{\pi}{2n} = \sqrt{ \frac{1 - \cos(\pi/n)}{2} }.$
而：
$\cos\frac{\pi}{n} = \sqrt{1 - \sin^2\frac{\pi}{n}} = \sqrt{ 1 - \frac{a_n^2}{4R^2} }.$
所以：
$a_{2n}^2 = (2R)^2 \cdot \frac{1 - \cos(\pi/n)}{2} = 2R^2 \left( 1 - \sqrt{ 1 - \frac{a_n^2}{4R^2} } \right).$
即：
$a_{2n} = \sqrt{ 2R^2 - 2R^2 \sqrt{ 1 - \frac{a_n^2}{4R^2} } }.$
这是精确的递推公式。
3. 刘徽的几何方法（简化推导）
刘徽的方法更直观，不直接用三角学，而是用勾股定理两次。
设 $ OA = R $，$ AB = a_n $，$ M $ 为 $ AB $ 中点，则 $ AM = a_{2n} $。
在直角三角形 $ OMB $ 中：
$OM = \sqrt{ OB^2 - MB^2 } = \sqrt{ R^2 - (a_n/2)^2 }.$
在直角三角形 $ AMC $ 中，其中 $ C $ 是 $ OA $ 与过 $ M $ 的垂线交点？ 不，直接看直角三角形 $ AMN $ 吗？
更直接：在直角三角形 $ AOM $ 中？ 不对，$ AOM $ 不是直角。
正确几何推导（勾股两次）：
设 $ OA = R $，$ AB = a_n $，$ M $ 是弧 $ AB $ 中点，则 $ AM = a_{2n} $。
1. 在直角三角形 $ OMB $ 中：
$OM = \sqrt{R^2 - \left( \frac{a_n}{2} \right)^2 }.$
2. 在直角三角形 $ AMC $ 中，点 $ C $ 是 $ OM $ 延长线与圆的交点？ 不，我们直接算 $ AM $：
考虑直角三角形 $ AMD $ 其中 $ D $ 是 $ M $ 在 $ OA $ 上的垂足？ 更简单的方法：
$ MC = R - OM $。
在直角三角形 $ AMC $ 中（$ C $ 是圆上点？ 不，我们定义 $ C $ 为 $ OA $ 与过 $ M $ 的垂线交点？ 这不对）。
其实经典割圆术推导（《九章算术》刘徽注）：
已知 $ a_n $，求 $ a_{2n} $ 公式为：
$a_{2n} = \sqrt{ \left( R - \sqrt{ R^2 - (a_n/2)^2 } \right)^2 + (a_n/2)^2 }.$
解释：
- 从 $ A $ 到 $ M $ 的水平距离是 $ a_n/2 $，
- 竖直距离是 $ R - OM = R - \sqrt{R^2 - (a_n/2)^2} $，
- 所以 $ AM = \sqrt{ (R - OM)^2 + (a_n/2)^2 } $。
化简：
$a_{2n}^2 = (R - OM)^2 + \left( \frac{a_n}{2} \right)^2$
$= R^2 - 2R\cdot OM + OM^2 + \frac{a_n^2}{4}.$
但 $ OM^2 + (a_n/2)^2 = R^2 $，所以：
$a_{2n}^2 = R^2 - 2R\cdot OM + R^2 = 2R^2 - 2R\cdot OM.$
即：
$a_{2n} = \sqrt{ 2R^2 - 2R \sqrt{ R^2 - (a_n/2)^2 } }.$
这与前面三角法导出的等价。
4. 从正六边形开始
取 $ R = 1 $（单位圆），正六边形边长 $ a_6 = 1 $（因为正六边形边长等于半径）。
- $ a_6 = 1 $
- $ a_{12} = \sqrt{ 2 - 2\sqrt{ 1 - (1/2)^2 } } = \sqrt{ 2 - 2\sqrt{3/4} } = \sqrt{ 2 - \sqrt{3} } $
- 继续迭代可得 $ a_{24}, a_{48}, a_{96} $。
5. 求圆周率
圆内接正 $ n $ 边形的周长 $ P_n = n \times a_n $。
当 $ n $ 很大时，$ P_n \approx 2\pi R $，所以：
$\pi \approx \frac{P_n}{2} \quad (\text{当 } R=1)。$
刘徽用面积来算更精确：
正 $ n $ 边形面积 $ S_n = \frac{n}{2} R^2 \sin\frac{2\pi}{n} $，但他用勾股关系逐次算。
面积递推（刘徽方法）：
从正 $ n $ 边形到正 $ 2n $ 边形，增加 $ n $ 个等腰三角形，每个底边 $ a_n $，高 $ R - OM $：
$S_{2n} = S_n + \frac{n}{2} \cdot a_n \cdot (R - OM).$
这样逐步提高精度。
6. 最终结果
刘徽计算到正 96 边形，得到：
$\pi \approx 3.141024 \quad (\text{96边})$
再算到 192 边形得 $ 3.14159 $（准确到小数点后 4 位）。
总结公式（单位圆 $ R=1 $）：
$a_{2n} = \sqrt{ 2 - 2\sqrt{1 - (a_n/2)^2} }$
$\text{周长 } P_{2n} = 2n \cdot a_{2n}$
$\pi \approx P_{2n} / 2$
这就是割圆术的几何推导核心。
