1. 从玩具到数学万花尺的奇妙世界如果你小时候玩过那种带齿轮的塑料尺在纸上转出各种繁复而规律的曲线那你一定对万花尺不陌生。它看起来只是个简单的绘图玩具几块钱就能买到一套由一个大尺环和几个带孔的小齿轮组成。但就是这个不起眼的小玩意儿其背后隐藏的数学原理和生成的轨迹之美足以让数学爱好者和图形学研究者着迷。今天我们不聊怎么玩而是深入探讨一下“万花尺轨迹性质研究”这个听起来有点学术但实际上充满趣味和挑战的话题。简单说我们要搞清楚为什么两个简单的圆尺环内壁和齿轮外沿相对滚动能画出如此复杂又对称的图案这些图案遵循什么数学规律我们能否预测和控制最终画出的图形这篇文章适合所有对几何、数学可视化、编程绘图或者单纯对“美”的生成规律感兴趣的朋友。我们将从最基础的原理拆解开始一步步推导出核心参数方程并用代码实现可视化最后探讨一些高级的轨迹性质和应用联想。2. 核心原理拆解两个圆的“舞蹈”要理解万花尺首先得抛开“尺”和“齿轮”的物理概念将其抽象为一个纯粹的几何问题。本质上万花尺描绘的是一个动圆小齿轮在一个定圆尺环内部或外部沿内壁纯滚动时动圆上某一点笔尖所在的孔所留下的轨迹。这是一个经典的“摆线”问题家族中的一员。2.1 关键参数定义为了进行定量研究我们必须先明确几个核心参数。假设我们有一个固定的定圆尺环半径为R。一个动圆齿轮半径为r在定圆内部紧贴内壁做纯滚动。在动圆上距离其圆心为d的地方有一个点这就是笔尖的位置d可以小于、等于或大于r。当d r时笔尖就在动圆的圆周上当d r时笔尖在动圆内部当d r时笔尖在动圆外部像一根延长臂。这三个参数(R, r, d)就完全决定了一个万花尺系统的几何构型。而轨迹的形态则取决于动圆滚动的过程。2.2 从物理滚动到参数方程动圆在定圆内纯滚动意味着没有滑动。我们可以想象两圆接触点的线速度必须时刻相等。这导出了一个关键关系动圆圆心绕定圆圆心的公转角速度与动圆绕自身圆心自转的角速度之比是固定的。设动圆圆心O绕定圆圆心O转动的角度为θ公转角。由于是纯滚动动圆上与定圆的接触点走过的弧长必须相等。定圆上走过的弧长为R * θ而动圆上由于滚动转过的弧长对应其自转的角度φ乘以半径r即r * φ。因此有R * θ r * φ由此可得自转角φ (R / r) * θ。注意这里有一个至关重要的符号当动圆在定圆内部滚动时其自转方向与公转方向相同。所以动圆上某一点相对于其自身圆心的位置会随着自转角φ而旋转。现在考虑笔尖点P。它的位置可以看作两个旋转的叠加动圆圆心O绕定点O旋转了θ角。点P绕动圆圆心O旋转了φ角方向与公转相同。在复平面或二维坐标系中我们可以轻松写出点P的坐标(x, y)关于参数θ的参数方程设定圆圆心O为原点(0, 0)。 动圆圆心O的坐标为( (R - r) * cosθ, (R - r) * sinθ )。因为两圆内切圆心距为R - r。 笔尖P相对于动圆圆心O的向量为( d * cos(φ), d * sin(φ) )其中φ (R / r) * θ。 因此点P的绝对坐标方程为x(θ) (R - r) * cosθ d * cos( (R/r) * θ ) y(θ) (R - r) * sinθ d * sin( (R/r) * θ )这就是内摆线Hypotrochoid的参数方程。当d r时笔尖在动圆圆周上轨迹称为内摆线Hypocycloid当d ! r时就是更一般的内旋轮线内次摆线。万花尺在尺环内画图对应的就是内摆线族。注意如果动圆在定圆外部滚动方程会变为外摆线Epitrochoid其圆心距为(R r)且自转角φ与公转角θ的比值关系为φ (R / r) * θ但方向通常相反或通过符号体现。万花尺套装中如果使用尺环外沿画出的就是外摆线。3. 轨迹形态的定性分析参数比R/r的决定性作用拿到参数方程后我们最关心的是R,r,d这三个数如何影响最终图案的样式、对称性和复杂度其中半径比R/r是最核心的控制器。3.1 有理数与图案的闭合性R/r这个比值决定了轨迹是否闭合以及闭合所需的“圈数”。如果R/r是一个有理数即可以表示为两个互质整数的比R/r p/qp, q为整数那么轨迹将是闭合的。因为公转θ每增加2π * q自转φ将增加(R/r) * 2π * q 2π * p。这意味着笔尖点P相对于动圆圆心的位置经历了整数p个完整的周期同时动圆圆心也绕行了整数q圈。此时点P会精确地回到起始位置形成一个闭合图形。例如若R/r 5/2则p5, q2。动圆公转2圈笔尖相对于动圆自转了5整圈图案闭合。通常我们称其有q2个“花瓣”或“循环节”但更准确地说其旋转对称阶数为q如果p, q互质。若R/r 4即4/1则p4, q1。动圆公转1圈笔尖自转4圈图案闭合且具有4阶旋转对称性。如果R/r是一个无理数那么无论公转多少圈笔尖相对于动圆的位置永远不会精确重复轨迹将永不闭合最终会以某种方式稠密地填满一个环形区域。这在物理万花尺上几乎看不到因为齿轮齿数都是整数比值为有理数但在数学模拟中可以实现产生类似“毛线团”或“环形密铺”的效果。3.2p与q对图形样式的具体影响假设R/r p/q已约分那么q值通常决定了图案的“主干”或“突出尖角”的数量。例如q3常常产生三角形或三瓣花结构q4产生方形或四瓣结构q5产生五角星或五瓣结构。q可以理解为图案的“基础对称阶数”。p值影响了图案内部的复杂度和“缠绕”次数。p越大在q个主干之间产生的环线、内旋结构就越复杂。p和q需要结合看。例如R/r5/2(p5, q2) 会画出一个类似五角星但在两个主轴方向来回缠绕的图案而R/r5/3则会画出更复杂的三瓣图案。3.3 距离d的作用从星形到玫瑰形参数d笔尖到动圆圆心的距离控制着轨迹的“胖瘦”和风格。d r笔尖在动圆内部。画出的曲线相对“圆润”尖角不明显更偏向于在环形区域内平滑变化的曲线有时像一朵饱满的花。d r笔尖在动圆周上。此时画出的是标准的内摆线Hypocycloid轨迹会经过动圆与定圆的切点图形通常有非常尖锐的尖角比如著名的三尖瓣线R/r3就是一个三角形但顶点是尖的。d r笔尖在动圆外部。这相当于给动圆加了一个“延长臂”。画出的曲线会有交叉的回环图形更加张扬、复杂经常出现类似“恒星”放射状的光芒或 loops环视觉效果非常华丽。d可以看作一个“振幅”调制器它放大了自转部分的贡献。当d接近0时轨迹几乎就是动圆圆心的轨迹——一个半径为(R-r)的圆。随着d增大轨迹开始偏离这个圆呈现出丰富的内部结构。4. 从方程到图形编程实现与可视化探索理论分析之后最好的理解方式就是亲手把它画出来。我们可以用任何具备绘图功能的编程语言来实现这里以 Python 的 Matplotlib 库为例因为它简单直观。4.1 基础绘图代码实现import numpy as np import matplotlib.pyplot as plt def draw_hypotrochoid(R, r, d, theta_max20*np.pi, num_points5000): 绘制内摆线万花尺轨迹 参数 R: 定圆半径 r: 动圆半径 d: 笔尖距离动圆圆心的距离 theta_max: 参数θ的最大值决定画多少圈 num_points: 采样点数量 # 生成参数θ数组 theta np.linspace(0, theta_max, num_points) # 计算参数方程 x (R - r) * np.cos(theta) d * np.cos((R - r) / r * theta) y (R - r) * np.sin(theta) d * np.sin((R - r) / r * theta) # 绘图 plt.figure(figsize(8, 8)) plt.plot(x, y, linewidth0.5) plt.axis(equal) # 保证x,y轴比例相同图形不变形 plt.title(fHypotrochoid: R{R}, r{r}, d{d}) plt.grid(True, alpha0.3) plt.show() # 示例1一个经典的5瓣花图案 (R/r 5, d r) draw_hypotrochoid(R5, r1, d0.7, theta_max10*np.pi) # 示例2一个尖锐的星形 (R/r 7/2, d r) # 注意这里R7, r2, 比值7/23.5 draw_hypotrochoid(R7, r2, d2, theta_max4*np.pi*2) # 公转q2圈闭合这段代码定义了核心的绘图函数。关键点在于theta_max的设置。为了确保画出闭合图形我们需要让θ至少经历2π * qq是R/r分数形式的分母。通常设置theta_max为2 * np.pi * q的若干倍以确保图形完整。4.2 交互式探索与参数影响观察静态代码只能看一个结果。要真正体会参数的影响最好能进行交互式调整。我们可以使用ipywidgets库在 Jupyter Notebook 中创建滑动条或者简单写一个循环来生成参数网格图。# 参数研究固定R5观察不同r和d的影响 R_fixed 5 d_values [0.3, 1.0, 2.0] # d r, r, r r_values [1, 2, 3, 4] # r 必须小于 R fig, axes plt.subplots(len(d_values), len(r_values), figsize(15, 10)) fig.suptitle(万花尺轨迹参数研究 (R5), fontsize16) for i, d in enumerate(d_values): for j, r in enumerate(r_values): if r R_fixed: # 避免非法值 axes[i, j].axis(off) continue theta np.linspace(0, 20*np.pi, 3000) # 注意使用 (R-r)/r 作为比值这是从几何推导中来的。 # 更通用的形式是 k R/r那么公式中的系数是 (k-1) (R-r)/r ratio (R_fixed - r) / r x (R_fixed - r) * np.cos(theta) d * np.cos(ratio * theta) y (R_fixed - r) * np.sin(theta) d * np.sin(ratio * theta) ax axes[i, j] ax.plot(x, y, linewidth0.7) ax.set_title(fr{r}, d{d}\nR/r≈{R_fixed/r:.2f}) ax.axis(equal) ax.axis(off) # 关闭坐标轴让图形更清晰 plt.tight_layout() plt.show()运行这段代码你会得到一个4x3的图形矩阵。横向看r变化R/r比值在变图形的对称性和瓣数在剧烈变化。纵向看d变化图形的“风格”从圆润到尖锐再到带有回环的复杂星形一目了然。这种可视化是理解参数作用最直接的方式。4.3 实操中的注意事项与技巧theta_max的设置艺术如果R/r是有理数p/q理论上theta_max 2π * q就能画完一个完整周期。但在实际数值计算中由于浮点数精度问题q可能很大例如R100, r33q33或者你根本不知道q是多少。一个稳妥的做法是设置一个较大的theta_max如20π或50π并观察图形何时开始重复覆盖。也可以通过计算R/r的分数近似来估算q。采样点数量num_points对于非常复杂的图形p和q都较大需要更多的采样点来保证曲线光滑。否则图形会出现明显的折线段。一个经验法则是num_points至少为theta_max / (2π) * 100即每圈公转采样100个点。性能与精度权衡当theta_max很大用于绘制无理数比近似图形时计算量会增大。如果不需要极高精度可以适当减少num_points。另外使用numpy的向量化运算如上例远比用for循环逐个计算点要快得多。闭合性判断在程序中如何自动判断图形是否闭合一个简单的方法是计算轨迹起点和终点的距离。但更可靠的方法是检查参数方程中的角度部分(R/r) * theta。如果(R/r)是有理数则cos和sin函数会在theta为2π * q的整数倍时周期重复。你可以计算(R/r) * 2π除以2π的小数部分是否接近0考虑浮点误差。5. 超越基础高级性质与扩展思考当我们能够熟练生成各种内摆线后可以进一步研究其更深层的数学性质和扩展应用。5.1 轨迹的对称性分析万花尺轨迹的对称性非常丰富。除了之前提到的旋转对称阶数为q当R/r p/q且p, q互质时许多图形还具有反射对称。例如当p和q均为整数且q 1时图形通常关于q条射线对称这些射线从中心点出发夹角为2π/q。此外如果p是偶数图形可能还具有关于中心的点对称即旋转180度重合。研究对称性不仅有趣在图形学和艺术设计中也很有用。你可以通过对称性来优化绘图算法只计算一个对称单元然后通过旋转复制得到完整图形大大提高效率。5.2 奇异点与尖点当d r标准内摆线时轨迹上会出现尖点Cusp。在尖点处曲线的方向发生突变导数切线斜率不连续。尖点的位置和数量有明确的规律对于R/r p/q约分后标准内摆线恰好有p个尖点。例如R/r3的三尖瓣线有3个尖点R/r4的四尖星形线有4个尖点。当d r时曲线是光滑的没有尖点。当d r时曲线会出现自交点Loop但通常不是尖点。研究这些特殊点的坐标和形成条件是微分几何的内容能帮助我们更深刻地理解曲线的局部行为。5.3 从二维到三维空间万花尺一个自然的扩展是如果让动圆在一个球面或其他曲面上滚动笔尖的轨迹会怎样或者我们增加一个维度的运动。例如让笔尖的距离d随时间周期性变化或者让定圆的半径R缓慢变化。这可以产生出极其复杂和美丽的三维螺线或Lissajous-like图形。用三维绘图库如 Matplotlib 的 3D 轴或 Mayavi可以将其可视化。其参数方程会变为x(θ) (R(θ) - r) * cosθ d(θ) * cos(k * θ) y(θ) (R(θ) - r) * sinθ d(θ) * sin(k * θ) z(θ) A * sin(ω * θ) # 例如增加一个垂直方向的简谐运动这已经进入了参数化艺术设计的领域。5.4 与傅里叶级数的深刻联系这是一个非常美妙的数学事实任何闭合、平滑的平面周期曲线都可以用有限或无限的傅里叶级数来表示。而内摆线/外摆线族恰好对应着一种特殊形式的傅里叶级数——只有两项谐波且频率比为有理数。观察我们的参数方程x(θ) A * cos(θ) B * cos(kθ)y(θ) A * sin(θ) B * sin(kθ)其中A R - r,B d,k R/r。这正是两个做匀速圆周运动的点的合成。第一个点以角速度1、半径A旋转第二个点以角速度k、半径B旋转。它们的合运动轨迹就是万花尺曲线。因此万花尺是理解二维傅里叶合成Lissajous 图形的一种特例的绝佳物理模型。当k为有理数时合成运动是周期的闭合图形当k为无理数时是非周期的但轨迹会稠密地覆盖一个环面区域。这个视角将玩具、古典几何和现代信号处理联系了起来。在数字图像处理中用少量几个傅里叶描述子来近似一个封闭轮廓其思想与此异曲同工。6. 实际应用联想与创作启发虽然万花尺本身是一个数学玩具但其原理和生成的图形在多个领域都有启发或直接应用。艺术与设计这是最直接的应用。万花尺图案具有强烈的数学美和韵律感被广泛应用于装饰艺术、纺织品设计、Logo设计、建筑纹样如玫瑰窗等领域。通过编程控制参数可以批量生成一系列相关但又各不相同的图案作为设计素材库。数学教育它是展示参数方程、周期运动合成、有理数与几何闭合性、对称性等概念的完美教具。学生可以通过调整参数即时看到图形变化建立直观理解。机械工程与机构学内摆线和外摆线是某些特殊齿轮如行星齿轮齿廓的理论基础。历史上外摆线齿轮曾用于钟表制造因为它能提供平稳的传动。研究点在这些曲线上的运动规律对分析一些连杆机构或凸轮机构的运动学有帮助。计算机图形学与屏幕保护程序万花尺算法简单视觉效果却非常丰富是编写演示程序Demo或屏幕保护程序的经典题材。通过加入颜色渐变根据θ或曲率变化、线宽变化、多重轨迹叠加用不同的d同时画多个点可以创造出令人惊叹的视觉效果。物理模拟在微观世界一些电子或粒子在电磁场中的运动轨迹有时也能用摆线族来近似描述。虽然不完全相同但研究思路有相通之处。对我个人而言研究万花尺最大的乐趣在于“可控的复杂性”。用寥寥几个参数R, r, d就能驾驭从简单圆环到复杂星云般的无穷图案谱系。它像一把钥匙打开了一扇门门后连接着古典几何、数论和现代计算机图形学。下次当你再看到万花尺画出的图案时希望你能认出其中p和q的舞蹈以及d赋予它的独特性格。不妨自己写几行代码试试改变几个数字你就能成为这个奇妙几何世界的创造者。