这份文档是国际文凭（IB）课程中**数学分析与方法（Mathematics: Analysis and Approaches）高级别（HL）**的内部评估（Exploration）。以下是其主要内容的概述：
主要内容概述引言与背景- 主题：在曲线 y=1x y = \frac{1}{x} y=x1 下拟合一系列圆形，探索其几何与数值规律。
- 个人背景：
- 作者受Desmos Marbleslides游戏（需调整函数收集星星）和Numberphile视频（介绍Kissing Number Problem）启发，对图形复杂性和优化问题产生兴趣。
- Kissing Number Problem涉及优化N维空间中与中心球相邻的球数，激发了作者将圆形拟合与曲线结合的灵感。
- 选择理由：
- 选择 y=1x y = \frac{1}{x} y=x1 因其对称性（沿 y=x y = x y=x 和 y=−x y = -x y=−x）和简单性，聚焦第一象限正值部分。
目的- 主要目标：
- 推导一个通用的圆方程，使圆形序列满足以下条件：
- 相邻圆相交一次。
- 圆与曲线 y=1x y = \frac{1}{x} y=x1 相切。
- 圆与 x x x-轴或 y y y-轴相切。
- 通过数学方法（几何与数值）确定圆的中心和半径序列。
- 研究范围：仅分析第一象限，利用对称性简化问题。
研究过程与分析- 欧几里得几何方法（Euclidean Geometric Approach）：
- 切点与切线：
- 设圆与曲线交点为 (s,1s) \left(s, \frac{1}{s}\right) (s,s1)，计算切线斜率 −1s2 -\frac{1}{s^2} −s21，得切线方程 y=−1s2x+2s y = -\frac{1}{s^2}x + \frac{2}{s} y=−s21x+s2。
- 法线方程（过圆心）为 y=s2x−s3+1s y = s^2 x - s^3 + \frac{1}{s} y=s2x−s3+s1。
- 圆的参数：
- 圆心 (h(s),f(s)) (h(s), f(s)) (h(s),f(s))，半径 f(s) f(s) f(s)。
- h(s)=2s−s2+1s2 h(s) = 2s - \sqrt{s^2 + \frac{1}{s^2}} h(s)=2s−s2+s21，f(s)=s2+1s2⋅tan(arctan(1s2)2) f(s) = \sqrt{s^2 + \frac{1}{s^2}} \cdot \tan\left(\frac{\arctan\left(\frac{1}{s^2}\right)}{2}\right) f(s)=s2+s21⋅tan(2arctan(s21))。
- 圆方程：(x−h(s))2+(y−f(s))2=f(s)2 (x - h(s))^2 + (y - f(s))^2 = f(s)^2 (x−h(s))2+(y−f(s))2=f(s)2（Equation 1）。
- 验证：以 s=1.21 s = 1.21 s=1.21 测试，确认方程有效。
- 坐标几何方法（Coordinal Geometric Approach）：
- 首圆分析：
- 首圆与 x x x-轴、y y y-轴及 y=1x y = \frac{1}{x} y=x1 相切，交点为 (1,1) (1, 1) (1,1)。
- 使用半径为1的圆，内接正方形边长 2 \sqrt{2} 2，交点为 (2+22,2+22) \left(\frac{2+\sqrt{2}}{2}, \frac{2+\sqrt{2}}{2}\right) (22+2,22+2)。
- 缩放得首圆半径 r1=2−2 r_1 = 2 - \sqrt{2} r1=2−2，圆心 (2−2,2−2) (2 - \sqrt{2}, 2 - \sqrt{2}) (2−2,2−2)。
- 第二方法：
- 应用勾股定理直接计算，确认 r1=2−2 r_1 = 2 - \sqrt{2} r1=2−2。
- 验证：在Desmos中绘制首圆，满足条件。
- 数值方法（Numerical Approach）：
- 圆心与半径关系：
- 观察到半径等于圆心 y y y-坐标，即 rn=yn r_n = y_n rn=yn。
- 相邻圆关系：(xn+1−xn)2=4yn+1yn (x_{n+1} - x_n)^2 = 4 y_{n+1} y_n (xn+1−xn)2=4yn+1yn（Equation 2）。
- 圆心到曲线的距离：
- 距离 p(s)=(xn−s)2+(yn−1s)2 p(s) = \sqrt{(x_n - s)^2 + \left(y_n - \frac{1}{s}\right)^2} p(s)=(xn−s)2+(yn−s1)2，求极值得 yn=s2(xn−s)+1s y_n = s^2 (x_n - s) + \frac{1}{s} yn=s2(xn−s)+s1（Equation 3）。
- 结合圆方程 (xn−s)2+(yn−1s)2=yn2 (x_n - s)^2 + \left(y_n - \frac{1}{s}\right)^2 = y_n^2 (xn−s)2+(yn−s1)2=yn2（Equation 4），解得：
- xn=2s−s2+1s x_n = 2s - \frac{\sqrt{s^2 + 1}}{s} xn=2s−ss2+1（Equation 5）。
- yn=s3−ss4+1+1s y_n = s^3 - s \sqrt{s^4 + 1} + \frac{1}{s} yn=s3−ss4+1+s1（Equation 6）。
- 递推关系：
- 推导九次多项式 f(s) f(s) f(s)（见PAGE13），根 sn+1 s_{n+1} sn+1 由 (xn,yn) (x_n, y_n) (xn,yn) 确定。
- 使用牛顿法（Newton's Method）求解：
- snk+1=snk−f(snk)f′(snk) s_n^{k+1} = s_n^k - \frac{f(s_n^k)}{f'(s_n^k)} snk+1=snk−f′(snk)f(snk)（Equation 7）。
- 以 s1=1 s_1 = 1 s1=1 开始，计算 s2=1.58313 s_2 = 1.58313 s2=1.58313、s3=2.08115 s_3 = 2.08115 s3=2.08115 等（Python代码验证）。
- 通解探索（General Equation）：
- 迭代过程：
- 计算前九个圆的 sn s_n sn 和 rn r_n rn（见表，PAGE18）。
- 绘制 rn r_n rn 对 n n n，R2=0.9999 R^2 = 0.9999 R2=0.9999，但无法准确预测第10个圆。
- 发现：
- 从第三个圆起，圆心近似沿 y=12x y = \frac{1}{2x} y=2x1 分布，半径约为曲线的“半径”。
结论与扩展- 成果：
- 成功拟合圆形序列，满足三条件。
- 首圆半径 2−2 2 - \sqrt{2} 2−2，后续通过牛顿法迭代计算。
- 反思：
- 九次多项式无闭合解，需数值方法。
- y=12x y = \frac{1}{2x} y=2x1 的拟合提示几何规律。
- 扩展：
- 探索六边形填充法求精确解。
- 计算圆总面积，判断收敛性及与曲线面积的比率。
附录与参考文献- 附录：Python代码解释。
- 参考文献：包括Kissing Number Problem、圆方程和牛顿法来源。
科目确认这份文档明确标注为数学分析与方法（Mathematics: Analysis and Approaches）HL的Exploration。依据： - 标题：注明“IB Diploma Mathematics AA HL Exploration”。
- 内容深度：涉及微分、九次多项式、牛顿法等，超出SL范围，符合HL的复杂性要求。
- 方法多样性：结合几何、坐标和数值分析，体现HL的探索性。
总结这份Exploration通过几何与数值方法，成功在 y=1x y = \frac{1}{x} y=x1 下拟合一列圆形，展示了作者对数学建模和优化问题的深入理解。核心内容包括圆方程推导、首圆计算及迭代求解，体现了HL水平的高阶分析能力。
文档链接地址：
| 
发表于 