自行车轮饰物的运动轨迹问答题

1、+理工大学暑期数学建模强化训练专题四自行车轮饰物的运动轨迹问题学员：曹阳许佳利 倪迪杭学院：通信工程学院时间：2010.08.19自行车轮饰物的运动轨迹问题本文就自行车轮饰物的运动轨迹问题，采用解析几何的方法建立数学模型， 求出了自行车在各种不同形状的道路上行驶时饰物和椭圆板中心的运动轨迹方 程，并且利用Matlab软件模拟仿真出了两者的运动轨迹。对于问题1和问题2,先运用解析几何方法求出自行车轮轴心的轨迹方程， 而后利用饰物始终绕车轮轴心作圆周运动建立参数方程，求出饰物的轨迹方程。 求出的曲线轨迹分别见图2、图4和图5。对于问题4,将“圆板”换为“椭圆板”，通过设定参数，结合坐标转换的 知识

2、，将转动过程中椭圆板中心的坐标用该参数表示，求出了其运动轨迹的参数方程。其轨迹图像见图 &关键词：运动轨迹，解析几何，抛物线，椭圆，坐标转换、问题的提出为了改变平淡的自行车外表，给自行车添加一分美妙的动感，同时，也为了 增加骑车人的“安全系数”，一些骑车人及自行车厂家在自行车的辐条上安装一 块亮丽的饰物。当有这种饰物的自行车在马路上驶过时，这种饰物就如游龙一样， 对街边的行人闪过一道波浪形的轨迹。这一波一闪的光亮游龙，也默默地维护着 骑车人的安全。建立数学模型解决以下问题：1、这轨迹是什么曲线？试画出它的图形。2、当这自行车又在一个抛物线形的拱桥上通过时，或是在一拱一拱的正弦 曲线（例

3、如山地摩托车赛场）上通过时，这饰物又画出一条曲中有曲的轨迹，这 轨迹是什么曲线？试画出它的图形。3、这种滚动中圆盘中心的运动轨迹是什么？4、将问题中“圆板”换为“凸形板”（例如椭圆板）时，其滚动轨迹会有什么结 果？二、问题的分析对于问题1、2,装有饰物的自行车在马路上行驶过程中，饰物会形成一道 道曲线轨迹。而随着路况的不同，如平坦的公路、抛物线形的拱桥、正弦曲线形 的山地摩托车赛场等，饰物会形成不同的曲线轨迹。不管形状的道路怎么样，不 管曲线轨迹有多复杂，都可以先运用解析几何方法求出自行车轮轴心的轨迹方 程，而后利用饰物始终绕车轮轴心作圆周运动建立参数方程，求出饰物的轨迹方程。饰物绕车轮轴心的

4、参数方程是较易得到的， 故问题的关键在于求出不同形状 的道路上自行车轮轴心的轨迹方程。求出饰物的轨迹方程后，设定合理的参数， 利用Matlab软件就可以模拟仿真出饰物的曲线轨迹。对于问题4,研究凸形板（椭圆板）中心的运动轨迹时，可以通过设定参数, 将转动过程中椭圆板中心的坐标通过坐标转换用参数表示，求出其运动轨迹方 程，设定合理的参数，利用 Matlab软件就可以模拟仿真出椭圆板中心的曲线轨 迹。三、模型假设1、自行车在行驶过程中车轮不打滑；2、自行车在行驶过程中速度保持不变；3、自行车在行驶过程中车轮始终与地面接触。1、符号说明Ra,b四、模型的建立及求解自行车轮的半径饰物距车轮轴心的距离自

5、行车车轮上的饰物 自行车轮的转动速度 定常量2、模型建立2.1模型1 （问题1）2.1.1模型的分析自行车在马路上行驶时，可以认为马路是平坦的，因此自行车轮轴心的运动 轨迹就是一条直线。求出这条直线方程后联立饰物绕车轮轴心的参数方程，就可以容易的得到饰物的运动轨迹。在检验模型的正确性时，可以利用图像求出几个 特殊点的坐标，与利用轨迹方程求出的坐标一一比较， 若两者求出的坐标都是相 同的，则验证了模型的正确性。2.1.2模型的建立与求解设自行车轮的半径为R，某一点为坐标原点建立坐标系，1所示，则饰物B的初始坐标车轮上的饰物B距离车轮轴心为r。以自行车轮上 假设饰物B的初始位置在y轴上，其示意图如

6、图Bo 为 0,R r。v设自行车轮的转动速度为，则自行车的速度为物B的绕自行车轮轴心旋转需满足参数方程：R。经过时间t后，饰xrsin t，yr cos tox Rt， y Ro所以，饰物B的坐标应满足参数方程：r sin t r cos ty根据一般自行车的规格，取R 软件画出其图像如图2所示：0.5m, rRt，R。0.4m,20rad / s ,使用 Matlab而自行车轮轴心的运动轨迹需满足参数方程:图2饰物在平坦道路上的轨迹2.1.3模型的验证用饰物B在特殊点的坐标来验证模型。假设自行车沿正x轴方向行驶，贝哨行车轮顺时针转动。饰物B转动一、 和- 时的坐标分别为B1、B2和Ba。2

7、 2饰物B转动一时的图像如图3所示。2由图像我们可以得到Bj的纵坐标为R。饰物B转动一时，车轮中心经过的2距离为一R，故Bi的横坐标为一R r。因此饰物B转动一时的坐标为2 2 2B R, R r。而把t T一 带入方程式求得的B的坐标Bi'R-R r，这2422与Bi的坐标是相同的。3用同样的方法，我们可以把饰物B转动 和3时的图像作出，然后利用图 像求得饰物B转动 和3时的坐标，这与由方程式计算出饰物 B转动2时的坐标都是相同的。从而验证了我们的模型的正确性。2.2模型2 （问题2）个类 抛物线。在求饰物的轨迹时，先利用拱桥的抛物线方程求出自行车轮轴心的轨迹， 然后联立饰物的坐标的

8、参数方程，就可以得到饰物的轨迹。当这自行车在一拱一拱的正弦曲线上通过时， 使用类似的方法，先求出自行 车轮轴心的轨迹，然后确定饰物的轨迹。2.2.2模型的建立与求解（1）当这自行车在一个抛物线形的拱桥上通过时，设拱桥的抛物线方程为：2.2.1模型的分析当这自行车在一个抛物线形的拱桥上通过时，自行车轮轴心的轨迹是Xi X0,2yiaX0bx0(a 0, b 0)。抛物线上任意一点的切线斜率ki为:人yi |x2ax0 b ,则该点的法向量的斜率k2为:1k2。ki假设自行车轮的轴心坐标为（X2,y2），则自行车轮的轴心的轨迹方程为:2 2 2(y2 yi) (X2 xi)R ,k2oX2为y2

9、yi解之得:X2Xiy2yiRk2F?Ro厂这就是自行车轮轴心的轨迹方程。抛物线上从原点到任意一点的距离为 s，则有:xs 0Ji (yi')2dx-1 ki Ji ki2 ln4abji b2 In 4aki Ji ki2b Ji b2又根据行驶路程与行驶时间之间的关系，有s vt Rt，令u t，则:su 一。R而经过时间t后，饰物B的绕自行车轮轴心旋转需满足参数方程:x r si nu,y r cosu。而自行车轮轴心的运动轨迹需满足参数方程:xX2,yy2。所以，饰物B的坐标应满足参数方程:x rsinu x2,y r cosu y2o这就是当这自行车在一个抛物线形的拱桥上通过

10、时饰物B的运动轨迹方程。取a 0.05, b 2，R 0.5m，r 0.4m，使用Matlab软件画出其图像如 图4所示：图4饰物在抛物线拱桥上的轨迹(2)当这自行车在一拱一拱的正弦曲线上通过时，设正弦曲线的方程为:X Xo,y1 asi nbXoO 0,b 0)。正弦曲线上任意一点的切线斜率k1为:kiyi lx Xab cosxo,则该点的法向量的斜率k2为:k2丄ki假设自行车轮的轴心坐标为(X2,y2)，贝哨行车轮的轴心的轨迹方程为:(y2 yi)2 (X2 xi)2 R2,k2oy2 yX2 X,解之得:X2Xiy2yiRR_ 。丄这就是自行车轮轴心的轨迹方程。设自行车在正弦曲线上行

11、驶的距离为 s，则有：X0 Ps x'j1 (y1)2dx o这个定积分较为复杂，故在求解s的时候，运用微积分的思想，将图形等步 长微分成n份，对每一部分求出面积后累加即可得到 s的解。n越大，求出的s就 却精确。故有：s 冷(y1)dxx0 'nJ1 (abcosi)2。i 1 n Vn同样令u t，根据s vtRt，得到饰物B的坐标应满足参数方程:x rsinuX2,图5饰物在正弦曲线上的轨迹1, b 0.6, n 100 , R 0.5m, r 0.4m，使用 Matlab 软件画出其 5所示：y r cosu y2o这就是当这自行车在一拱一拱的正弦曲线上通过时饰物 B的

12、运动轨迹方程。 取a图像如图模型3 (问题4)2.32.3.1模型的建立与求解将问题中的“圆板”换为“椭圆板”之后，考虑在水平道路上运动。 设初始状态时椭圆的方程为：2 2y a x221 a b 0 ，a b其参数形式为:x bsin , y a acos o则椭圆中心的坐标为0 0, a，坐标原点为01 0,0，如图6所示。图6设椭圆转动任意角度 后，椭圆中心为0' X2,y2 A Xo,yo。尝试通过坐标转换，找出椭圆转动角度 达式，即：，椭圆与X轴的切点为 后x2和y2关于的参数表X2X2y2y2则过切点A x0,y0的切线方程为:yo a y2 aX0gX 1。b2所以，该切

13、线的斜率k2为:k2b2 yoa2Xo atan 。从而可以推得:yo a2a帯Xo。b k2点A Xo,yo在椭圆上，故又满足:yo a2 a联立（1）、（2）可以得到：2Xo1。b2Xob2k2yoa这样就得到了 Xo和yo关于 的参数表达式。 由于点A Xo，yo在椭圆上，亦满足椭圆的参数方程:Xoyobsin o,0°a a cos故得到:令直线OA与椭圆长轴的夹角为tanarc cos 1，则有:Xoayoyo 。a故,arc tanXoa yo令坐标原点Oi和切点A之间的弧长OiA为s，则有:o722,s va Sin b cos d 。o运用微积分的思想求解s，将图形等

14、步长微分成n份，对每一部分求出面积 后累加即可得到s的解。n越大，求出的s就却精确。故有：7a2 sin2b2 cos2 d2.i o Sin n2i ocos。n切点A X0,y0到椭圆中心O 0,a的距离d2yo a观察图7,容易得到椭圆中心o'在坐标轴中满足:%o1a'1 1OAy2O'a' cossin所以，椭圆中心O在椭圆转动过程当中的轨迹方程为:x2s d sin取a 4.5 , b 2.5，使用Matlab软件画出其图像如图8所示:图8椭圆中心在水平道路上运动轨迹五、模型的评价本文通过对自行车运动过程中饰物以及圆板（凸形板）的运动进行了详细的 分析

15、，并通过一定的物理和数学思想及方法，借助于MATLAB件模拟了问题中所需的各点的轨迹。在解决问题二时，我们巧妙地运用了数值积分的方法， 把自 行车在抛物线轨道上通过时饰物的轨迹方程求解出来，从而模拟出了其轨迹。对于问题四，我们设定了椭圆参数方程，通过对其运动过程物理状态的分析，完成 了对其中心轨迹的求解。本文的缺点在于由于时间有限，对于问题 4的求解，只考虑椭圆板在水平直 线上的运动。七、模型的推广可以将问题4拓展为较为简单的凸形板在常见平滑曲线上的运动， 通过设定 参数，结合坐标转换的知识，将转动过程中简单的凸形板中心的坐标用该参数表 示，从而求出其运动轨迹的参数方程。参考文献1.同济大学数

16、学系高等数学第六版上册高等教育出版社 2007 年4月出版附录一一图2的程序:t=0:0.01:3;r=0.4;R=0.5;a=3;v=a*t w=v./R x =-r*si n(w.*t)+w.*R.*t ; y =R-r*cos(w.*t); plot(x,y)axis(0 7 0 1.2)附录二图4的程序：x0=-1:0.01:40;a=-0.05;b=2;R=0.5;r=0.4;x1=x0;y仁a*x0.A2+b*x0k1=2*a*x0+b;plot(x1,y1)hold onk2=-1./k1;%法线的斜率；x2=R./(k2.*sqrt(1+(1./k2).2)+x1;y2=R./

17、sqrt(1+(1./k2).A2)+y1%p Iot(x2,y2)hold onu=1/R*(1/(4*a)*(k1.*sqrt(1+k1.2)+Iog(k1+sqrt(k1.2+1)-1/(4*a)*(b*sqrt(1+b2)+Iog(b+sqrt(1 +b2);x=-r*s in( u)+x2;y=y2-r*cos(u);Plot(x,y);axis(-2 45 0 30)图5程序：x0=-2:0.01:18; a=1； b=0.6;R=0.5;r=0.4;x1=x0;y1=a*si n(b*x0);k1=a*b*cos(b*x0);plot(x1,y1)hold onk2=-1./k1

18、;%法线的斜率；x2=R./(k2.*sqrt(1+(1./k2).2)+x1;y2=R./sqrt(1+(1./k2).A2)+y1;%p Iot(x2,y2)hold onfor x3=-2:0.01:18s=0;for i=1:100s=s+sqrt(1+(a*b*cos(b*i*(x3+2)/100-2)A2)*(x3+2)/100end u=s/Ry1=a*si n(b*x3);k仁a*b*cos(b*x3);k2=-1./k1;x2=R./(k2.*sqrt(1+(1./k2).2)+x3;y2=R./sqrt(1+(1./k2)A2)+y1;x=-r*s in( u)+x2;y=y2-r*c