1. 项目概述从玩具到数学与工程的桥梁万花尺这个很多人童年记忆中的绘图玩具其背后隐藏的数学原理远比我们想象的要深邃。它不仅仅是一个能画出繁复美丽图案的尺规更是“圆内旋轮线”这一经典几何曲线最直观、最生动的物理演示器。当我们用MATLAB来仿真万花尺时我们实际上是在搭建一座桥梁一端连接着趣味性的艺术创作和直观的几何理解另一端则通向严肃的工程应用比如异形齿轮的设计、机构运动学分析以及轨迹规划。这个项目的核心就是利用MATLAB强大的数值计算和图形可视化能力去动态地再现并探索万花尺的绘图过程。我们不只是想画出一张静态的、漂亮的旋轮线图而是要模拟笔尖动点相对于尺框大圆和齿轮小圆的精确运动关系观察参数变化如何实时影响最终图案的形态。这对于理解非圆齿轮、偏心齿轮传动甚至是某些机械臂末端的复杂轨迹生成都有着直接的启发意义。通过这个仿真无论是数学爱好者、机械专业的学生还是对运动仿真感兴趣的工程师都能获得一个低成本、高效率的探索平台亲手“调参”来感受几何与运动的魅力。2. 核心数学模型与原理拆解万花尺的数学本质是一个经典的“两个旋转运动合成”问题。我们可以将其抽象为一个纯粹的几何模型一个大圆定圆内部有一个小圆动圆紧贴其内壁做纯滚动。在小圆上或延伸出的连杆上固定一个点笔尖该点在运动过程中于大圆所在的平面内留下的轨迹就是圆内旋轮线也称为内摆线。2.1 坐标系与参数定义为了用MATLAB精确描述这个过程我们首先需要建立清晰的数学模型和坐标系。大圆定圆圆心位于坐标系原点 O(0, 0)半径为 R。它是我们的参考系。小圆动圆初始时刻小圆的圆心 C 位于大圆的正右方 (R - r, 0) 处半径为 r。小圆在大圆内部做无滑动的纯滚动。笔尖点 P位于小圆上或与小圆固连的连杆上。设其在小圆自身坐标系中的位置为 (a, b)。当 a r 且 b 0 时笔尖就在小圆周上当 a 和 b 不为零时笔尖就相当于在小圆外“伸”出了一根杆。关键的物理参数是两圆的半径比k R / r。这个比值决定了旋轮线是否闭合以及其“花瓣”的数量。2.2 运动学方程的推导推导笔尖 P 在世界坐标系 (O-XY) 中的坐标 (x, y) 随时间 t或转角变化的参数方程是整个仿真的基石。设小圆圆心 C 绕原点 O 的公转角速度为 ω_c通常我们关心的是相对运动可以设大圆静止。由于是纯滚动小圆上的接触点线速度相等。经过推导过程涉及弧长相等关系我们可以得到小圆绕自身圆心 C 的自转角速度 ω_s 与 ω_c 的关系。一个更直观的推导方法是追踪两个旋转运动公转小圆圆心 C 绕原点 O 的旋转。设公转角度为 θ。那么 C 点坐标为C ( (R-r)*cos(θ), (R-r)*sin(θ) )。自转由于纯滚动当小圆公转 θ 角时其滚过的弧长等于大圆上对应的弧长。这导致小圆自身的旋转角度 φ 与 θ 满足r * φ R * θ考虑方向后。因此φ (R/r) * θ k * θ。这里的负号是关键因为小圆在内侧滚动其自转方向与公转方向通常相反。所以更准确的关系是φ -k * θ。现在笔尖 P 在小圆坐标系圆心为C随小圆公转但不自转的坐标系中的坐标是 (a, b)。为了得到世界坐标我们需要先将这个向量进行“自转”变换旋转 -φ 角度因为小圆自身转了 φ相当于笔尖在小圆坐标系中转了 -φ然后再叠加“公转”的平移。最终笔尖 P 的世界坐标参数方程为x(θ) (R - r)*cos(θ) a*cos(φ) - b*sin(φ) y(θ) (R - r)*sin(θ) a*sin(φ) b*cos(φ) 其中φ - (R / r) * θ -k * θ参数 θ 从 0 变化到2*pi*NN 是使得轨迹闭合所需的最小公转圈数通常与半径比 k 是否为有理数有关。注意这个方程是仿真的核心。参数 a, b 给了我们极大的灵活性。当 (a, b) (r, 0) 时就是标准的内摆线当 a r 或 b ≠ 0 时画出的就是延长或偏移的内摆线这也是万花尺玩具能产生千变万化图案的原因。2.3 图案闭合性与“花瓣”数轨迹是否闭合取决于小圆公转多少圈后笔尖能回到起始点。这要求存在整数 m, n 使得m * θ 2π * n且同时满足相位条件。分析表明当半径比k R/r是一个有理数即可以表示为两个整数之比 p/q且 p, q 互质时轨迹是闭合的。并且图案的“花瓣”数或尖点数量与 p, q 的值直接相关。若k p/q为最简分数则小圆公转 q 圈后轨迹闭合。图案的“瓣数”通常为 p当笔尖在小圆上时或与 p、q 的奇偶性相关的其他值。若 k 为无理数则轨迹永不闭合最终会稠密地填满大圆内的一个环域。3. MATLAB仿真实现与核心代码解析有了数学模型我们就可以用MATLAB将其转化为生动的动画和图像。我们的目标是实现一个交互性强、可视化的仿真程序。3.1 基础静态绘图首先我们实现一个最基本的静态绘图函数用于快速验证公式和观察图案。function drawSpirograph(R, r, a, b, nRotations) % R: 大圆半径 % r: 小圆半径 % a, b: 笔尖在小圆坐标系中的坐标 % nRotations: 小圆公转的圈数 k R / r; % 半径比 theta linspace(0, nRotations*2*pi, 5000); % 公转角度5000个点保证曲线光滑 % 核心参数方程 phi -k * theta; % 小圆自转角度 x (R - r) * cos(theta) a * cos(phi) - b * sin(phi); y (R - r) * sin(theta) a * sin(phi) b * cos(phi); % 绘图 figure(Position, [100, 100, 800, 800]); plot(x, y, b-, LineWidth, 1.5); axis equal; grid on; title(sprintf(圆内旋轮线 (R%.1f, r%.1f, a%.1f, b%.1f), R, r, a, b)); xlabel(X); ylabel(Y); end调用示例drawSpirograph(100, 30, 35, 0, 10)会绘制一个经典的、花瓣密集的旋轮线。3.2 动态动画仿真静态图虽美但动态过程更能揭示运动本质。下面我们实现一个动画实时展示小圆滚动和笔尖画线的过程。function animateSpirograph(R, r, a, b, nRotations, speed) % speed: 动画速度控制越小越快但不应小于0.01 k R / r; thetaVec linspace(0, nRotations*2*pi, 1000); % 减少点数以提高动画流畅度 totalFrames length(thetaVec); % 预计算所有帧的坐标以提高效率 x zeros(1, totalFrames); y zeros(1, totalFrames); Cx (R - r) * cos(thetaVec); % 小圆圆心轨迹 Cy (R - r) * sin(thetaVec); for i 1:totalFrames theta thetaVec(i); phi -k * theta; x(i) Cx(i) a * cos(phi) - b * sin(phi); y(i) Cy(i) a * sin(phi) b * cos(phi); end % 初始化图形窗口 fig figure(Position, [100, 100, 900, 900]); ax axes(Parent, fig); axis(ax, equal); grid(ax, on); hold(ax, on); title(ax, sprintf(万花尺动态仿真 (R%.1f, r%.1f, a%.1f, b%.1f), R, r, a, b)); % 绘制静态的大圆 rectangle(Position, [-R, -R, 2*R, 2*R], Curvature, [1,1], ... EdgeColor, k, LineWidth, 1, LineStyle, --); % 初始化动态图形对象 hTrail plot(ax, nan, nan, b-, LineWidth, 1.5); % 轨迹线 hCircle plot(ax, nan, nan, r-, LineWidth, 1.5); % 小圆 hCenter plot(ax, nan, nan, ro, MarkerSize, 8, MarkerFaceColor, r); % 小圆圆心 hPen plot(ax, nan, nan, go, MarkerSize, 10, MarkerFaceColor, g); % 笔尖 hLine plot(ax, nan, nan, k-, LineWidth, 0.5); % 圆心到笔尖的连线 xlim(ax, [-R*1.2, R*1.2]); ylim(ax, [-R*1.2, R*1.2]); % 动画循环 for i 1:totalFrames theta thetaVec(i); phi -k * theta; % 当前小圆的位置 currentCx Cx(i); currentCy Cy(i); % 计算当前小圆轮廓上的点用于绘制红色小圆 alpha linspace(0, 2*pi, 100); circleX currentCx r * cos(alpha); circleY currentCy r * sin(alpha); % 当前笔尖位置 currentPenX x(i); currentPenY y(i); % 更新图形对象数据 set(hTrail, XData, x(1:i), YData, y(1:i)); set(hCircle, XData, circleX, YData, circleY); set(hCenter, XData, currentCx, YData, currentCy); set(hPen, XData, currentPenX, YData, currentPenY); set(hLine, XData, [currentCx, currentPenX], YData, [currentCy, currentPenY]); drawnow; pause(speed); % 控制帧率 end hold(ax, off); end调用示例animateSpirograph(50, 15, 20, 5, 5, 0.02)。你会看到红色小圆沿大圆内壁滚动绿色笔尖画出蓝色轨迹黑色线连接圆心和笔尖。实操心得在编写动画时预计算所有帧的数据如笔尖轨迹x,y至关重要。如果放在循环内实时计算尤其是总帧数多时会导致动画严重卡顿。drawnow命令强制刷新图形pause(t)控制播放速度。将图形对象hTrail,hCircle等的句柄保存起来在循环中只更新其‘XData’和‘YData’属性是制作流畅MATLAB动画的标准高效做法。3.3 创建交互式图形用户界面GUI为了让参数探索更方便我们可以创建一个简单的GUI使用滑块实时调整参数并观察图案变化。这里使用MATLAB传统的GUIDE或面向对象的uifigure方式皆可。以下是一个基于uifigure和回调函数的简洁示例框架。function spirographGUI() % 创建主窗口 fig uifigure(Name, 万花尺仿真实验室, Position, [100 100 1200 700]); % 创建坐标区域 ax uiaxes(fig, Position, [300 50 650 600]); title(ax, 万花尺图案); axis(ax, equal); grid(ax, on); hold(ax, on); % 创建控制面板 panel uipanel(fig, Title, 参数控制, Position, [20 50 250 600]); % 大圆半径 R 滑块 uilabel(panel, Text, 大圆半径 R:, Position, [20 520 100 22]); sliderR uislider(panel, Position, [20 500 200 3], ... Limits, [10 200], Value, 100, ... ValueChangedFcn, (src,event) updatePlot()); % 小圆半径 r 滑块 uilabel(panel, Text, 小圆半径 r:, Position, [20 460 100 22]); sliderr uislider(panel, Position, [20 440 200 3], ... Limits, [5 100], Value, 40, ... ValueChangedFcn, (src,event) updatePlot()); % 笔尖坐标 a 滑块 uilabel(panel, Text, 笔尖坐标 a:, Position, [20 400 100 22]); sliderA uislider(panel, Position, [20 380 200 3], ... Limits, [0 150], Value, 80, ... ValueChangedFcn, (src,event) updatePlot()); % 笔尖坐标 b 滑块 uilabel(panel, Text, 笔尖坐标 b:, Position, [20 340 100 22]); sliderB uislider(panel, Position, [20 320 200 3], ... Limits, [-50 50], Value, 0, ... ValueChangedFcn, (src,event) updatePlot()); % 公转圈数输入框 uilabel(panel, Text, 公转圈数:, Position, [20 280 100 22]); editRot uieditfield(panel, numeric, Position, [120 280 100 22], ... Value, 10, Limits, [1 100], ... ValueChangedFcn, (src,event) updatePlot()); % 绘图按钮 btnPlot uibutton(panel, push, Text, 重新绘图, ... Position, [80 200 100 30], ... ButtonPushedFcn, (src,event) updatePlot(true)); % 清空按钮 btnClear uibutton(panel, push, Text, 清空画布, ... Position, [80 150 100 30], ... ButtonPushedFcn, (src,event) clearPlot()); % 存储图形对象句柄 hPlot plot(ax, nan, nan, b-, LineWidth, 1.5); hBigCircle plot(ax, nan, nan, k--, LineWidth, 1); % 更新绘图函数嵌套函数可以访问GUI控件 function updatePlot(forceRecalc) if nargin 1 forceRecalc false; end persistent lastParams currentData if isempty(lastParams) lastParams zeros(1,5); end % 获取当前参数 R sliderR.Value; r sliderr.Value; a sliderA.Value; b sliderB.Value; nRot editRot.Value; currentParams [R, r, a, b, nRot]; % 如果参数未变且不是强制重算则跳过优化性能 if isequal(currentParams, lastParams) ~forceRecalc return; end lastParams currentParams; % 计算轨迹 k R / r; theta linspace(0, nRot*2*pi, 3000); phi -k * theta; x (R - r) * cos(theta) a * cos(phi) - b * sin(phi); y (R - r) * sin(theta) a * sin(phi) b * cos(phi); % 更新轨迹线 set(hPlot, XData, x, YData, y); % 更新大圆参考线 th linspace(0, 2*pi, 100); set(hBigCircle, XData, R*cos(th), YData, R*sin(th)); % 更新坐标轴范围 axisRange max(abs([x, y])) * 1.1; if axisRange 0 xlim(ax, [-axisRange, axisRange]); ylim(ax, [-axisRange, axisRange]); else xlim(ax, [-R, R]); ylim(ax, [-R, R]); end title(ax, sprintf(R%.1f, r%.1f, a%.1f, b%.1f, k%.3f, R, r, a, b, k)); drawnow; end function clearPlot() set(hPlot, XData, nan, YData, nan); drawnow; end % 初始化绘图 updatePlot(true); end运行spirographGUI()即可启动交互界面。拖动滑块图案会实时变化这极大地便利了参数探索。4. 从仿真到应用异形齿轮与工程启示万花尺仿真绝不仅仅是数学游戏其生成的复杂、精确的闭合曲线在机械工程领域尤其是非圆齿轮和异形齿轮的设计中有着直接的应用。4.1 旋轮线与齿轮节曲线在齿轮传动中两个齿轮的节圆或节曲线在啮合过程中做纯滚动。这与我们万花尺模型中大圆和小圆的关系完全一致。万花尺的小圆圆心轨迹即(R-r)*cos(θ), (R-r)*sin(θ)可以看作一个齿轮的节圆而笔尖P的轨迹则可以启发我们设计另一个齿轮的齿廓形状使得两个齿轮在传动时其中一个齿轮上的某个点相当于笔尖能走出预设的复杂轨迹。具体来说如果我们把万花尺模型中的笔尖P的轨迹固定下来将其视为一个“工件”需要被加工出的形状或者一个执行器需要走过的路径。那么小圆的圆周就可以被视作一个“刀具”或“生成齿轮”的节圆。当这个生成齿轮小圆沿着固定的内齿轮大圆做纯滚动时其上的“刀尖”笔尖P就会包络出我们想要的轨迹。这个过程在机械原理中称为“包络法”生成共轭齿廓。4.2 仿真用于齿轮设计验证MATLAB仿真在这里可以成为一个强大的前期验证工具参数化设计将齿轮的节曲线参数R, r、刀具参数a, b作为输入快速生成理论上的齿廓曲线。干涉检查通过高密度采样仿真啮合过程可以检查在传动过程中两个齿轮的齿廓是否会发生干涉碰撞。这可以通过计算两个齿廓曲线之间的最小距离来实现。运动规律分析对于非圆齿轮传动比是变化的。我们可以通过仿真轻松计算出主动轮转角与从动轮转角的关系即φ与θ的关系进而分析速度、加速度的变化是否满足设计要求。可视化啮合过程类似我们的动画可以制作两个异形齿轮啮合传动的动画直观展示传动过程和平顺性。4.3 一个简单的异形齿轮齿廓生成示例假设我们要设计一个内啮合齿轮副内齿轮大轮的节曲线是半径为R的圆外齿轮小轮的节曲线是半径为r的圆。我们希望小轮上的一个点能走出一个特定的内摆线轨迹。那么该轨迹就是小轮齿廓的等距线。我们可以利用仿真先画出理论接触点轨迹即旋轮线然后通过数学计算求法向等距线生成实际齿廓。MATLAB代码可以如下function generateGearProfile(R, r, a, b, toothWidth) % 生成基于旋轮线的异形齿轮齿廓简化示例 % toothWidth: 齿宽法向距离 % 1. 生成理论接触点轨迹旋轮线 theta linspace(0, 2*pi, 1000); k R / r; phi -k * theta; x_t (R - r) * cos(theta) a * cos(phi) - b * sin(phi); y_t (R - r) * sin(theta) a * sin(phi) b * cos(phi); % 2. 计算轨迹上各点的法向量简化计算使用数值微分 dx gradient(x_t, theta); dy gradient(y_t, theta); % 法向量为 (-dy, dx)并单位化 normals sqrt(dx.^2 dy.^2); nx -dy ./ normals; ny dx ./ normals; % 3. 生成内外两侧的等距线齿廓 x_inner x_t toothWidth/2 * nx; y_inner y_t toothWidth/2 * ny; x_outer x_t - toothWidth/2 * nx; y_outer y_t - toothWidth/2 * ny; % 4. 绘图 figure; plot(x_t, y_t, k-, LineWidth, 1, DisplayName, 理论接触线); hold on; plot(x_inner, y_inner, r-, LineWidth, 2, DisplayName, 齿廓一侧); plot(x_outer, y_outer, b-, LineWidth, 2, DisplayName, 齿廓另一侧); axis equal; grid on; legend; title(基于旋轮线的异形齿轮齿廓生成); end注意事项这是一个高度简化的示例。真实的齿轮设计需要考虑根切、顶隙、压力角、强度等一系列因素并且需要保证共轭齿廓的精确数学关系。此处的价值在于展示了如何将万花尺的几何模型与齿轮设计的基本思路包络法联系起来并用MATLAB进行快速原型验证。5. 常见问题、技巧与扩展探索在实际仿真和探索过程中你可能会遇到以下问题这里提供一些解决思路和进阶技巧。5.1 仿真精度与性能平衡问题轨迹有棱角不光滑。原因参数theta的采样点太少。linspace(0, 2*pi*nRotations, N)中的N太小。解决增加N例如从 1000 增加到 5000 或 10000。但要注意这会增加计算量和动画的存储负担。问题动画卡顿严重。原因每帧计算量太大或pause时间太短导致图形刷新跟不上。解决预计算如前面代码所示在循环前计算好所有数据。减少帧数动画的thetaVec点数可以比最终绘图点数少。使用drawnow limitrate在动画循环中用drawnow limitrate替代drawnow它可以阻止MATLAB刷新图形直到系统准备好能提升流畅度。升级图形对象更新方式对于极复杂的动画可以考虑使用animatedline对象逐步添加点效率更高。5.2 参数选择的艺术与规律如何得到对称、漂亮的图案闭合性确保R/r为有理数。尝试让R和r为整数比如 5:2, 7:3, 12:5。“花瓣”数量当b0时图案的尖点或环数大致与分子pR/r的最简分数形式相关。p为奇数时常得到p个尖点p为偶数时常得到p/2个环或p个瓣。笔尖位置的影响a控制图案的大小和复杂度a越接近r曲线越靠近中心a大于r会形成带有环的复杂图形。b会破坏对称性产生螺旋或偏移效果带来不对称的美感。一个实用技巧先固定R然后用一个循环去尝试不同的r和a批量生成图片再挑选喜欢的图案。MATLAB可以很方便地实现R 100; figure; for i 1:9 r 10 5*i; % 变化的r a r * (0.5 0.2*rand()); % 在r附近随机取a subplot(3,3,i); drawSpirograph(R, r, a, 0, lcm(round(R), round(r))/round(r)); % 确保闭合 title(sprintf(r%d, a%.1f, r, a)); end5.3 扩展方向超越经典万花尺多级万花尺让小圆内部再有一个更小的圆滚动形成更复杂的复合旋轮线。数学模型需要叠加更多层的旋转运动。椭圆旋轮线将大圆或小圆改为椭圆。运动方程将基于椭圆的弧长参数计算会复杂很多但能产生更独特的图案。3D 万花尺让笔尖在三维空间中运动。例如小圆在一个球体内或一个圆锥面上滚动。这需要引入三维旋转矩阵。与物理引擎结合使用 Simulink 或 Simscape Multibody 建立真正的多体动力学模型考虑质量、惯性、摩擦研究实际万花尺玩具的动力学特性。艺术化应用将生成的轨迹点坐标导出为 SVG 或 DXF 格式供激光切割、CNC雕刻或绘图机器人使用制作实体艺术品。5.4 调试与排查轨迹看起来不对首先检查参数方程中的正负号。自转角度φ前面的负号最容易出错它决定了小圆是顺滚还是逆滚。如果不确定可以先将b设为0a设为r此时笔尖应画出一个标准的多尖点内摆线。如果画出的是一条简单的弧线或圆大概率是φ的符号错了。动画中图形闪烁在更新图形对象属性前设置ax.NextPlot add即hold on并在创建图形对象后不要频繁创建和删除对象只更新其数据。GUI 响应慢确保回调函数updatePlot中的计算是高效的。避免在每次滑块移动时都重新计算大量数据。可以加入防抖逻辑比如使用pause和drawnow结合或者检查参数是否真的改变了再重算。这个项目就像一把钥匙打开了一扇连接数学、编程、艺术和工程的大门。从一行行代码中看到那些精妙的曲线逐渐浮现再到理解它背后严谨的几何约束和广阔的工程应用可能性这种体验远比单纯使用一个玩具或学习一个数学公式要深刻得多。我个人最享受的时刻就是在调整参数时心中根据有理数的规律预测图案的大致形状然后运行程序验证预测与视觉结果是否吻合——那一刻仿佛直接触摸到了数学的纹理。你不妨也从一组简单的参数开始运行起来然后大胆地拖动滑块去发现属于你自己的、独一无二的旋轮线世界。