欢迎光临散文网 会员登陆 & 注册

美丽的摆线及其数学原理

2022-07-25 12:41 作者:现代微积分  | 我要投稿

最近看了有位昵称为“帆雨动画”的up主的科普视频,其中的许多直观的科普动画让观众赏心悦目。而从中挖掘出美妙的数学原理也是一件很有趣的事。下面主要来讲讲有关摆线的数学原理。

(1)摆线、次摆线

摆线:

我们可将运动过程分解为圆上一点匀速圆周,同时圆心向右做匀速直线运动

设圆方程为:x%5E2%2By%5E2%3Dr%5E2,圆上一点从(0,-1)开始顺时针转动,其坐标为:(rcos(-%5Calpha%20%2B%5Cfrac%7B3%5Cpi%20%7D%7B2%7D%20)%2Crsin(-%5Calpha%20%2B%5Cfrac%7B3%5Cpi%20%7D%7B2%7D%20))

(-rsin%5Calpha%20%2C-rcos%5Calpha%20)

转过的弧长为l%3D%5Calpha%20r

将圆心至(%5Calpha%20r%2Cr),此时动点平移至(%5Calpha%20r-rsin%5Calpha%20%2Cr-rcos%5Calpha%20)

摆线的参数方程为:%5Cleft%5C%7B%5Cbegin%7Bmatrix%7D%20%0A%20x%3Dr(%5Calpha%20-sin%5Calpha%20)%5C%5Cy%3Dr(1-cos%5Calpha%20)%0A%5Cend%7Bmatrix%7D%5Cright.%20

其中α为参数

如上图所示,取过该动点及动圆圆心的直径,可知该直径上的点在运动过程中的相对位置是不变的

取以圆心为起点,上述圆周上该动点为终点的向量(-rsin%5Calpha%20%2C-rcos%5Calpha%20)

将其进行伸缩得:(-%5Clambda%20rsin%5Calpha%20%2C-%5Clambda%20rcos%5Calpha%20)

其中λ∈[-1,1]

伸缩后向量终点坐标为:(%5Calpha%20r-%5Clambda%20rsin%5Calpha%20%2Cr-%5Clambda%20rcos%5Calpha%20)

次摆线的参数方程为:

%5Cleft%5C%7B%5Cbegin%7Bmatrix%7D%20%0A%20x%3Dr(%5Calpha%20-%5Clambda%20sin%5Calpha%20)%5C%5Cy%3Dr(1-%5Clambda%20cos%5Calpha%20)%0A%5Cend%7Bmatrix%7D%5Cright.%20

其中α为参数

下面是对圆上一直径上的几个点的轨迹示踪图:

美丽的摆线和次摆线

其中圆周上的红点和紫点的轨迹为摆线,其余的为次摆线

ps:上图用desmos制作,desmos可添加a=rgb(数字1,数字2,数字3)添加颜色(a换成其他字母也行),颜色的参数可参考网址:https://tool.oschina.net/commons?type=3

取一个周期2π研究运动轨迹长度:

%5Cbegin%7Barray%7D%0A%5C%5CL%3D%5Cint_%7B0%7D%5E%7B2%5Cpi%20%7D%5Csqrt%7B(%5Cfrac%7Bdx%7D%7Bd%5Calpha%20%7D%20)%5E2%2B(%5Cfrac%7Bdy%7D%7Bd%5Calpha%20%7D)%5E2%7D%20d%5Calpha%20%0A%5C%5C%3D%5Cint_%7B0%7D%5E%7B2%5Cpi%20%7D%20%5Csqrt%7Br%5E2(1-%5Clambda%20cos%5Calpha%20)%5E2%2Br%5E2(%5Clambda%20sin%5Calpha)%5E2%7D%20d%5Calpha%20%0A%5C%5C%3Dr%5Cint_%7B0%7D%5E%7B2%5Cpi%20%7D%20%5Csqrt%7B%7B%5Clambda%20%7D%5E2-2%5Clambda%20cos%5Calpha%20%2B1%7D%20d%5Calpha%20%0A%5Cend%7Barray%7D


特别地,当λ=0时,L=2πr,即圆心的运动路径长;

当λ=1时,

%20%5Cbegin%7Barray%7D%0A%5C%5CL%3Dr%5Cint_%7B0%7D%5E%7B2%5Cpi%20%7D%20%5Csqrt%7B2-2cos%5Calpha%20%7D%20d%5Calpha%20%0A%5C%5C%3Dr%5Cint_%7B0%7D%5E%7B2%5Cpi%20%7D%20%5Csqrt%7B4sin%5E2(%5Cfrac%7B%5Calpha%20%7D%7B2%7D%20)%7D%20d%5Calpha%20%0A%5C%5C%3D2r%5Cint_%7B0%7D%5E%7B2%5Cpi%20%7D%5Cleft%20%7C%20sin(%5Cfrac%7B%5Calpha%20%7D%7B2%7D%20)%20%5Cright%20%7C%20d%5Calpha%20%0A%5Cend%7Barray%7D

∵α∈[0,2π]

∴sin(α/2)≥0

%5Cbegin%7Barray%7D%0A%5C%5CL%3D2r%5Cint_%7B0%7D%5E%7B2%5Cpi%20%7Dsin(%5Cfrac%7B%5Calpha%20%7D%7B2%7D)d%5Calpha%20%0A%5C%5C%3D2r%5B-2cos(%5Cfrac%7Bx%7D%7B2%7D%20)%5D%5E%7B2%5Cpi%20%7D_%7B0%7D%0A%5C%5C%3D2r%5B2-(-2)%5D%0A%5C%5C%3D8r%0A%5Cend%7Barray%7D

即摆线路径长为8r;

λ=-1时同上;

若λ∈(-1,1),则轨迹为次摆线,积分无初等解(涉及到椭圆积分,跟椭圆周长的积分式类似),不能求出精确长度


(2)外摆线、内摆线

在此之前想先提及下所谓的“硬币悖论”:“有大小相同的两枚硬币,固定一枚硬币不动,将另一枚硬币贴着固定的硬币滚动一周回到原位,问运动的硬币转了几圈?”

首先有一点是确定的,也就是滚动过程中下图所示的两段弧相等

也就是说滚动一周动圆的圆周才能被完整地接触一次

而为什么我们看到的却是转了2圈呢?我认为用参考系来解释合理很多。

下面分享一位网友的思路图:

大致意思如下,转动过程中,相接触过的弧PM=弧PM',因此动圆相对于切点P转过的弧度为θ(即∠PC'M'),也就是若我们始终认为过P点的公切线是水平线,那么就是“一圈”。

而图中圆心相对于x轴的偏转角为2θ(即∠C'Ox+∠OC'M'),即圆心相对x轴的偏转角是相对过P点切线的偏转角的2倍,因此以定圆或x轴为参考系那么就是"两圈"。

所以不用为答案是“一圈”还是“两圈”纠结了,关键是参考系的选取

只要明白“滚过的弧长相等”的含义以及参考系的选取思路就清晰了,下文的求解需用到“相对于切点的转动”这一论述。


下面回到主题上,


外摆线:

其中滚过的弧长相等(标红的两段弧相等)

设定圆半径为r1,动圆半径为r2

定圆方程为:x%5E2%2By%5E2%3Dr_1%5E2

考虑到旋转用复数较好表示坐标,于是建立复平面进行研究

切点对应复数为r_1e%5E%7Bi%5Calpha%20%7D

由于两圆外切,则动圆圆心对应复数为(r_1%2Br_2)e%5E%7Bi%5Calpha%20%7D

由(r,0)开始滚过的弧长为r_1%5Calpha%20

则动圆相对于切点转过的弧度为%5Cfrac%7Bl%7D%7Br_2%7D%20%3D%5Cfrac%7Br_1%5Calpha%20%7D%7Br_2%7D%20

取由动圆圆心指向切点的向量,对应复数:

r_1e%5E%7Bi%5Calpha%20%7D-(r_1%2Br_2)e%5E%7Bi%5Calpha%20%7D%3D-r_2e%5E%7Bi%5Calpha%20%7D

将该向量逆时针旋转%5Cfrac%7Br_1%5Calpha%20%7D%7Br_2%7D%20得:

-r_2e%5E%7Bi%5Calpha%20%7D%5Ccdot%20e%5E%7B%5Cfrac%7Br_1%5Calpha%20%7D%7Br_2%7D%20i%7D%3D-r_2e%5E%7B(1%2B%5Cfrac%7Br_1%7D%7Br_2%7D%20)%5Calpha%20i%7D

故动点坐标对应的复数为:

(r_1%2Br_2)e%5E%7Bi%5Calpha%20%7D-r_2e%5E%7B(1%2B%5Cfrac%7Br_1%7D%7Br_2%7D%20)%5Calpha%20i%7D

用欧拉公式化为三角形式得:

%20%5Cbegin%7Barray%7D%0A%5C%5C%3D(r_1%2Br_2)(cos%5Calpha%20%2Bisin%5Calpha%20)-r_2(cos%5B(1%2B%5Cfrac%7Br_1%7D%7Br_2%7D)%5Calpha%5D%2Bisin%5B(1%2B%5Cfrac%7Br_1%7D%7Br_2%7D)%5Calpha%5D)%0A%5C%5C%3D(r_1%2Br_2)cos%5Calpha%20-r_2cos%5B(1%2B%5Cfrac%7Br_1%7D%7Br_2%7D)%5Calpha%5D%2Bi((r_1%2Br_2)sin%5Calpha%20-r_2sin%5B(1%2B%5Cfrac%7Br_1%7D%7Br_2%7D)%5Calpha%5D)%0A%5Cend%7Barray%7D

实部和虚部分别对应参数方程的x,y值

外摆线参数方程为:

%5Cleft%5C%7B%5Cbegin%7Bmatrix%7D%20%0A%20x%3D(r_1%2Br_2)cos%5Calpha%20-r_2cos%5B(1%2B%5Cfrac%7Br_1%7D%7Br_2%7D)%5Calpha%5D%0A%5C%5Cy%3D(r_1%2Br_2)sin%5Calpha%20-r_2sin%5B(1%2B%5Cfrac%7Br_1%7D%7Br_2%7D)%5Calpha%5D%0A%5Cend%7Bmatrix%7D%5Cright.%20

其中α为参数

作出轨迹示踪图:

取不同的r1和r2可绘制出各种美丽的图案

美丽的外摆线

内摆线:

其中滚过的弧长相等(标红的两段弧相等)

设定圆半径为r1,动圆半径为r2(r1>r2)

定圆方程为:x%5E2%2By%5E2%3Dr_1%5E2

切点对应复数为r_1e%5E%7Bi%5Calpha%20%7D

由于两圆内切,则动圆圆心对应复数为(r_1-r_2)e%5E%7Bi%5Calpha%20%7D

由(r,0)开始滚过的弧长为r_1%5Calpha%20

则动圆相对于切点转过的弧度为%5Cfrac%7Bl%7D%7Br_2%7D%20%3D%5Cfrac%7Br_1%5Calpha%20%7D%7Br_2%7D%20

取由动圆圆心指向切点的向量,对应复数:

r_1e%5E%7Bi%5Calpha%20%7D-(r_1-r_2)e%5E%7Bi%5Calpha%20%7D%3Dr_2e%5E%7Bi%5Calpha%20%7D

将该向量顺时针旋转%5Cfrac%7Br_1%5Calpha%20%7D%7Br_2%7D%20得:

r_2e%5E%7Bi%5Calpha%20%7D%5Ccdot%20e%5E%7B-%5Cfrac%7Br_1%5Calpha%20%7D%7Br_2%7D%20i%7D%3Dr_2e%5E%7B(1-%5Cfrac%7Br_1%7D%7Br_2%7D%20)%5Calpha%20i%7D

故动点坐标对应的复数为:

(r_1-r_2)e%5E%7Bi%5Calpha%20%7D%2Br_2e%5E%7B(1-%5Cfrac%7Br_1%7D%7Br_2%7D%20)%5Calpha%20i%7D

用欧拉公式化为三角形式得:

%20%5Cbegin%7Barray%7D%0A%5C%5C%3D(r_1-r_2)(cos%5Calpha%20%2Bisin%5Calpha%20)%2Br_2(cos%5B(1-%5Cfrac%7Br_1%7D%7Br_2%7D)%5Calpha%5D%2Bisin%5B(1-%5Cfrac%7Br_1%7D%7Br_2%7D)%5Calpha%5D)%0A%5C%5C%3D(r_1-r_2)cos%5Calpha%20%2Br_2cos%5B(1-%5Cfrac%7Br_1%7D%7Br_2%7D)%5Calpha%5D%2Bi((r_1-r_2)sin%5Calpha%20%2Br_2sin%5B(1-%5Cfrac%7Br_1%7D%7Br_2%7D)%5Calpha%5D)%0A%5Cend%7Barray%7D

实部和虚部分别对应参数方程的x,y值

内摆线参数方程为:

%5Cleft%5C%7B%5Cbegin%7Bmatrix%7D%20%0A%20x%3D(r_1-r_2)cos%5Calpha%20%2Br_2cos%5B(1-%5Cfrac%7Br_1%7D%7Br_2%7D)%5Calpha%5D%0A%5C%5Cy%3D(r_1-r_2)sin%5Calpha%20%2Br_2sin%5B(1-%5Cfrac%7Br_1%7D%7Br_2%7D)%5Calpha%5D%0A%5Cend%7Bmatrix%7D%5Cright.%20

其中α为参数

作出轨迹示踪图:

取不同的r1和r2可绘制出各种美丽的图案

美丽的内摆线

拓展一下,观察曲线参数方程的解析式可知,摆线和次摆线可由一“匀速圆周运动”和一“匀速直线运动”叠加绘制出;而外摆线和内摆线可由两“匀速圆周运动”叠加绘制出,多个匀速圆周运动叠加绘图也是傅里叶变换绘图中所用到的。还有许多美丽的几何有待去研究,还有许多数学知识有待去探索和发掘,让我们尽情在数学世界里遨游吧~

美丽的摆线及其数学原理的评论 (共 条)

分享到微博请遵守国家法律