最速降线问题仿真的方法Matlab程序上海交通大学数学系.ppt

数学实验 寻找最速降线 一王海交通大学数学系三  数学给我们一个用之不竭,充满真理的宝庳,这 些真理不是孤立的,而是以相互密切的关系并立着, 而且随着科学的每一成功进展,我们会不断发现这 些真理之间的新的接触点 C.F. Gauss 数学既不严峻,也不遥远,它和几乎所有的人 类活动有关,又对每个真心对它感兴趣的人有益 R C. Buck  内容提要 回顾微积分有关知识 连续,多元函数极值,积分等 〉复习微分方程的求解的解析与数值方法 介绍一类最优问题的求解新框架一变分方法 最速降线求解的仿真方法  背景故事 1696年 John bernoulli向他的兄长和其他 数学家挑战性地提出了最速降线(捷线)问题: 一质量为m的质点,在重力作用下从定点A 沿曲线下滑到定点B 试确定一条曲线,使得质 点由A到B下滑时间最短. 假定B比A低,不计摩擦力 和其他阻力等因素. 此问题导致数学新分支的产生  思考 这是一个求最值的问题 与求函数的极值一样吗? 与求线性规划问题中的极值一样吗? 它的数学形式怎样? 历史 1697年5月号“教师学报”接收了5篇解答报告  贝努利·约翰 Bernoulli, Johann 欧洲著名科学家族 冫涉猎微积分、微分方程、解析几 何、概率论以及变分法 更贡献于物理、化学和天文学 谁发现 LHospital法贴 欧拉的指导者和老师 a 瑞士的骄傲  问题数学形式 A 设曲线为y=y(x),(x∈,c) 满足y(0)=0,y(c)=H 若T=(y)={质点沿y=y(x) B 下滑的时间 我们要求的是怎样的函数y(x) 使得T(y)取得最小值 minTo)  近似方法 如图建立坐标系,设A为原A1 点,B为(C,H),将带状区域 0yH用平行于x轴的 直线y=yk=kH/m把这区域 分成n个带状小区域 在带状域yyy,可近似认为vk=√28y 而曲线段近似认为是直线段,其长度 x2-x1)2+(△y)2  于是质点从A到B所需时间近似为 (x-x1)2+(Ay) (n-1元函数!) (4y是已知的!) 〉求这个函数的极小值,就得到问题的近似解 (为简单计,可取g=1000cm/s2) 可以使用数学软件来求极值,但所得曲线为 离散形式,无解析表达式  求解极值数值方法 OT 0 0xv[(x1-x1)2+h2 -x i)+h X -x x:;-X 可令 (x-x)2+h2 (x1-x1)2+h2 解出 CV h (i=1,2,…,n)(*) 故 xtch