寻找最速降线.ppt

《寻找最速降线.ppt》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、寻找最速降线数学给我们一个用之不竭,充满真理的宝库,这些真理不是孤立的,而是以相互密切的关系并立着,而且随着科学的每一成功进展,我们会不断发现这些真理之间的新的接触点.──C.F.Gauss数学既不严峻,也不遥远,它和几乎所有的人类活动有关,又对每个真心对它感兴趣的人有益.──R.C.Buck介绍一类最优问题的求解新框架-变分方法连续,多元函数极值,积分等内容提要回顾微积分有关知识复习微分方程的求解的解析与数值方法最速降线求解的仿真方法1696年JohnBernoulli向他的兄长和其他数学家挑战性地提出了最速降线（捷线）问题:一质量为m的质点，在重

2、力作用下从定点A沿曲线下滑到定点B，AB试确定一条曲线，使得质点由A到B下滑时间最短.假定B比A低，不计摩擦力和其他阻力等因素.此问题导致数学新分支的产生.背景故事思考这是一个求最值的问题与求函数的极值一样吗？与求线性规划问题中的极值一样吗？它的数学形式怎样？历史1697年5月号“教师学报”接收了5篇解答报告贝努利约翰Bernoulli,Johann欧洲著名科学家族涉猎微积分、微分方程、解析几何、概率论以及变分法谁发现L’Hospital法则欧拉的指导者和老师更贡献于物理、化学和天文学瑞士的骄傲问题数学形式ABxyc设曲线为满足y(0)=

3、0,y(c)=H我们要求的是怎样的函数y(x)下滑的时间}{质点沿y=y(x)若使得T(y)取得最小值minT(y)近似方法如图建立坐标系,设A为原点,B为(c,H),将带状区域直线y=yk=kH/n把这区域ABxycyk-1xk-1ykxk分成n个带状小区域.在带状域yk-1

4、令解出故求解极值数值方法令为下列方程的解再将代入(*)式中，将用曲线连接即得拟合最速降线，再求出时间.functionm6_1(G,H,n)h=H/n;g=9.8;f=1.0;a=0;b=2/(sqrt(2*g*(n-1)/n*H));c=(a+b)/2;i=1;whileabs(f)>1e-10s=0;forj=1:nv=sqrt(2*g*j*h);s=s+v/sqrt(1.0-c^2*v^2);endf=c-G/(h*s);iff>0b=c;elsea=c;endc=(a+b)/2;i=i+1;endx(1)=sqrt(g*h/2)*c*h/sqrt(1

5、.0-c*c*2*g*h);T=sqrt((x(1)-a)^2+h^2)/sqrt(2*g*h)fork=2:nv=sqrt(2*g*k*h);x(k)=x(k-1)+c*v*h/sqrt(1.0-c*c*v*v);T=T+sqrt((x(k)-x(k-1))^2+h^2)/v;endplot(x,-(0.1:h:H),'*r')利用数学软件求近似最速降线和最短时间利用数学软件求解得到的曲线再作分析质点要走最快的路线(曲线)，应该如何变化？依然用从质点速度变化的角度考虑设质点从A1经直线l到达A2，质点速度在l的上侧为v1，下侧为v2，则质点如何运动才最省

6、时？A1A212ClOD如图,若A1,A2到l的垂足分为a,b,OD=c,质点经过l于C别为O,D,A1,A2到l的距离分别OC=x那么质点由A1到A2需时间惟一驻点满足也即这就是光学中的Snell折射定律A1A212ClODxabc－x建立数学模型分析：如图建坐标系，AB分割成小段,考虑在第kABxyck+1k层与k+1层质点在曲线上的下滑，依能量守恒律，可近似认为质点在每层内的速度不变，于是依辅助结论知注意上式对任何k成立，若用与x轴平行的直线将ABxyc令平行线的间距趋于零，我们就得到在曲线上任何一点其中为该点切线与铅垂线的夹角故导出导

7、出微分方程ABxyc又因于是得到一个引理设集合E0={g(x)C1│g(a)=g(b)=0}如果在[a,b]连续函数f(x)满足那么f(x)0对g(x)E0，总有另一种方法－变分法ABxyc设曲线为满足y(0)=0,y(c)=H在曲线上P(x,y）处质点速度为又设从A到P的弧长为s,则从而质点沿曲线由A到B需时间那么我们的问题成为求某个使得引进集合显然若是最速曲线函数,则于是函数在取得最小值故得设集合那么对依复合函数求导法注意第二项为了计算,记上式乘以这里于是导出注意从降线定义可知故1）可求解析解解法2）也可以用数值方法，例如欧拉法求解得到方程

8、为由于在原点y=0,可改写方程求解析解提示：解析解