寻找最速降线培训课件.pptVIP

一个引理 设集合E0={g(x)?C1 │g(a) =g(b)=0} 如果在[a,b]上连续函数 f (x)满足 那么 f (x) ? 0 对?g (x) ?E0 ，总有 另一种方法－变分法 A B x y b 设曲线为 满足 y(0)=0, y(b)=H 在曲线上P(x,y）处质点速度为 又设从A到P的弧长为s,则 从而质点沿曲线由A到B需时间 那么我们的问题成为 求某个 使得 引进集合 显然若 是最速曲线函数,则 于是函数 在 取得最小值 故得 设集合 那么对 依复合函数求导法 注意第二项 为了计算 ,记 上式乘以 这里 于是导出 注意从降线定义可知 故 1）可求解析解 解法 2）也可以用数值方法，例如欧拉法求解 得到方程为 上海交通大学数学科学学院 数学实验 寻找最速降线 数学给我们一个用之不竭,充满真理的宝库,这些真理不是孤立的,而是以相互密切的关系并立着,而且随着科学的每一成功进展,我们会不断发现这些真理之间的新的接触点. ── C. F. Gauss 数学既不严峻,也不遥远,它和几乎所有的人类活动有关,又对每个真心对它感兴趣的人有益. ── R. C. Buck ? 介绍一类最优问题的求解新框架-变分方法 连续,多元函数极值,积分等 内容提要 ? 回顾微积分有关知识 ? 复习微分方程求解的解析与数值方法 ? 最速降线求解的仿真方法 1696年John Bernoulli向他的兄长和其他 数学家挑战性地提出了最速降线（捷线）问题: 一质量为m的质点，在重力作用下从定点A 沿曲线下滑到定点B， A B 试确定一条曲线，使得质 点由A到B下滑时间最短. 假定B比A低，不计摩擦力 和其他阻力等因素. ? 此问题导致数学新分支的产生. 背景故事 思考 这是一个求最值的问题 ? 与求函数的极值一样吗？ ? 与求线性规划问题中的极值一样吗？ ? 它的数学形式怎样？ 历史 1697年5月号“教师学报”接收了5篇解答报告 贝努利? 约翰 Bernoulli,Johann ? 欧洲著名科学家族 ? 涉猎 微积分、微分方程、解析几 何、概率论以及变分法 ? 谁发现 L’Hospital 法则 ? 欧拉的指导者和老师 更贡献于物理、化学和天文学 ? 瑞士的骄傲 问题的数学形式 A B x y b 设曲线为 满足 y(0)=0, y(b)=H 我们要求的是怎样的函数y(x), 下滑的时间} {质点沿 y=y(x) 若 使得T(y) 取得最小值 求 minT(y) 近似方法 如图建立坐标系,设A为原 点, B为(b,H), 将带状区域 直线 y=yk=k H/n 把这区域 A B x y b yk-1 xk-1 yk xk 分成 n个带状小区域. 在带状域yk-1yyk ,可近似认为 而曲线段近似认为是直线段,其长度 0 y H用平行于 x 轴的 ? 求这个函数的极小值, 就得到问题的近似解 ( n -1元函数! ) 于是质点从A到B所需时间近似为 ? 可以使用数学软件来求极值，但所得曲线为 离散形式，无解析表达式 ( 是已知的！) (为简单计,可取g =1000cm/s2) 令 可令 解出 故 求解极值数值方法 令 为下列方程的解 再将 代入(*)式中，将 用曲线连接即得拟合最速降线，再求出时间 . function m6_1(b,H,n) h=H/n;g=9.8;f=1.0;a=0; d=2/(sqrt(2*g*(n-1)/n*H));c=(a+d)/2; i=1; while abs(f)1e-10 s=0; for j=1:n v=sqrt(2*g*j*h);s=s+v/sqrt(1.0-c^2*v^2); end f=c-b/(h*s); if f0 d=c;else a=c; end c=(a+d)/2;i=i+1;end x(1)=sqrt(g*h/2)*c*h/sqrt(1.0-c*c*2*g*h); T=sqrt((x(1)-a)^2+h^2)/sqrt(2*g*h) for k=2:n v=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; end plot(x,-(0.1:h:H),*r) 利用数学软件求近似最速降线和最短时间 利用数学软件求解得到的曲线 再作分析 质点要走最快的路线(曲线)，应该如何变化？ ? 依然用从质点速度变化的角度考虑 设质点从A1经直线 l 到达A