寻找最速降线.

Slides:



Advertisements
Similar presentations
一、 一阶线性微分方程及其解法 二、 一阶线性微分方程的简单应用 三、 小结及作业 §6.2 一阶线性微分方程.
Advertisements

第五节 全微分方程 一、全微分方程及其求法 二、积分因子法 三、一阶微分方程小结. 例如 所以是全微分方程. 定义 : 则 若有全微分形式 一、全微分方程及其求法.
目录 上页 下页 返回 结束 习题课 一、导数和微分的概念及应用 二、导数和微分的求法 导数与微分 第二章.
第九章 常微分方程数值解法 §1 、引言. 微分方程的数值解:设方程问题的解 y(x) 的存在区间是 [a,b] ,令 a= x 0 < x 1
2.8 函数的微分 1 微分的定义 2 微分的几何意义 3 微分公式与微分运算法则 4 微分在近似计算中的应用.
第八章 第四节 机动 目录 上页 下页 返回 结束 一个方程所确定的隐函数 及其导数 隐函数的微分法.
2.6 隐函数微分法 第二章 第二章 二、高阶导数 一、隐式定义的函数 三、可微函数的有理幂. 一、隐函数的导数 若由方程 可确定 y 是 x 的函数, 由 表示的函数, 称为显函数. 例如, 可确定显函数 可确定 y 是 x 的函数, 但此隐函数不能显化. 函数为隐函数. 则称此 隐函数求导方法.
2.5 函数的微分 一、问题的提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、微分的求法 六、小结.
第二章 导数与微分 一. 内 容 要 点 二. 重 点 难 点 三. 主 要 内 容 四. 例 题与习题.
第二章 导数与微分. 二、 微分的几何意义 三、微分在近似计算中的应用 一、 微分的定义 2.3 微 分.
全微分 教学目的:全微分的有关概念和意义 教学重点:全微分的计算和应用 教学难点:全微分应用于近似计算.
第三节 微分 3.1 、微分的概念 3.2 、微分的计算 3.3 、微分的应用. 一、问题的提出 实例 : 正方形金属薄片受热后面积的改变量.
§3.4 空间直线的方程.
一、曲面及其方程 二、母线平行于坐标轴的柱面方程 三、以坐标轴为旋转轴的旋转曲面 四、小结
《解析几何》 -Chapter 3 §7 空间两直线的相关位置.
第八章 向量代数 空间解析几何 第五节 空间直线及其方程 一、空间直线的点向式方程 和参数方程 二、空间直线的一般方程 三、空间两直线的夹角.
3.4 空间直线的方程.
第七章 空间解析几何 §5 空间直线及其方程 一、空间直线的一般方程 二、空间直线的对称式方程与参数方程 三、两空间直线的夹角
圆复习.
例7-1 荡木用两条等长的钢索平行吊起,钢索的摆动规律为j= j 0sin(pt/4)。试求当t=0和t=2s时,荡木中点M的速度和加速度。
第三章 函数逼近 — 最佳平方逼近.
第二章 二次函数 第二节 结识抛物线
例题 教学目的: 微积分基本公式 教学重点: 牛顿----莱布尼兹公式 教学难点: 变上限积分的性质与应用.
恰当方程(全微分方程) 一、概念 二、全微分方程的解法.
高等数学电子教案 第五章 定积分 第三节 微积分基本定理.
第五节 微积分基本公式 、变速直线运动中位置函数与速度 函数的联系 二、积分上限函数及其导数 三、牛顿—莱布尼茨公式.
一、原函数与不定积分 二、不定积分的几何意义 三、基本积分公式及积分法则 四、牛顿—莱布尼兹公式 五、小结
第二节 微积分基本公式 1、问题的提出 2、积分上限函数及其导数 3、牛顿—莱布尼茨公式 4、小结.
第四章 函数的积分学 第六节 微积分的基本公式 一、变上限定积分 二、微积分的基本公式.
微积分基本定理 2017/9/9.
第四章 一元函数的积分 §4.1 不定积分的概念与性质 §4.2 换元积分法 §4.3 分部积分法 §4.4 有理函数的积分
第5章 定积分及其应用 基本要求 5.1 定积分的概念与性质 5.2 微积分基本公式 5.3 定积分的换元积分法与分部积分法
第四节 一阶线性微分方程 线性微分方程 伯努利方程 小结、作业 1/17.
第三节 格林公式及其应用(2) 一、曲线积分与路径无关的定义 二、曲线积分与路径无关的条件 三、二元函数的全微分的求积 四、小结.
多元函数微分学学习辅导 一、内容提要 二、典型例题 首页 上页 返回 下页 结束.
初中数学 九年级(下册) 5.3 用待定系数法确定二次函数表达式.
用函数观点看方程(组)与不等式 14.3 第 1 课时 一次函数与一元一次方程.
计算机数学基础 主讲老师: 邓辉文.
§2 求导法则 2.1 求导数的四则运算法则 下面分三部分加以证明, 并同时给出相应的推论和例题 .
双曲线的简单几何性质 杏坛中学 高二数学备课组.
§7.2 直线的方程(1) 1、经过两点P1(x1,y1),P2(x2,y2)的斜率公式: 2、什么是直线的方程?什么是方程的直线?
高等数学 西华大学应用数学系朱雯.
本节内容 平行线的性质 4.3.
2.1.2 空间中直线与直线 之间的位置关系.
第4章 非线性规划 4.5 约束最优化方法 2019/4/6 山东大学 软件学院.
第二十二章 曲面积分 §1 第一型曲面积分 §2 第二型曲面积分 §3 高斯公式与斯托克斯公式.
第三单元 第3课 实验 多元函数的积分 实验目的:掌握matlab计算二重积分与三重积分的方法,提高应用重积分解决有关应用问题的能力。
线段的有关计算.
第四章 一次函数 4. 一次函数的应用(第1课时).
⑴当∠MBN绕点B旋转到AE=CF时(如图1),比较AE+CF与EF的大小关系,并证明你的结论。
§1体积求法 一、旋转体的体积 二、平行截面面积为已知的立体的体积 三、小结.
第五节 对坐标的曲面积分 一、 对坐标的曲面积分的概念与性质 二、对坐标的曲面积分的计算法 三、两类曲面积分的联系.
复习: 若A(x1,y1,z1) , B(x2,y2,z2), 则 AB = OB - OA=(x2-x1 , y2-y1 , z2-z1)
第三单元 第2课 实验 一元函数的积分 实验目的:掌握matlab求解有关不定积分和定积分的问题,深入理解定积分的概念和几何意义。
3.1 变化率与导数   3.1.1 变化率问题 3.1.2 导数的概念.
抛物线的几何性质.
相似三角形存在性探究 嘉兴市秀洲区王江泾镇实验学校 杨国华
第一节 不定积分的概念与性质 一、原函数与不定积分的概念 二、不定积分的几何意义 三、基本积分表 四、不定积分的性质 五、小结 思考题.
直线和圆的位置关系 ·.
学习任务三 偏导数 结合一元函数的导数学习二元函数的偏导数是非常有用的. 要求了解二元函数的偏导数的定义, 掌握二元函数偏导数的计算.
高中数学选修 导数的计算.
轴对称在几何证明及计算中的应用(1) ———角平分线中的轴对称.
§2-2 点的投影 一、点在一个投影面上的投影 二、点在三投影面体系中的投影 三、空间二点的相对位置 四、重影点 五、例题 例1 例2 例3
选修1—1 导数的运算与几何意义 高碑店三中 张志华.
24.4弧长和扇形面积 圆锥的侧面积和全面积.
第四节 向量的乘积 一、两向量的数量积 二、两向量的向量积.
Volterra-Lotka方程 1925年, A. Lotka(美)和V. Volterra(意)给出了第一个两物种间的捕食模型。
5.1 相交线 (5.1.2 垂线).
第一模块 向量代数与空间解析几何 第六节 二次曲面与空间曲线 一、曲面方程的概念 二、常见的二次曲面及其方程 三、空间曲线的方程
§3.1.2 两条直线平行与垂直的判定 l1 // l2 l1 ⊥ l2 k1与k2 满足什么关系?
Presentation transcript:

寻找最速降线

数学给我们一个用之不竭,充满真理的宝库,这些真理不是孤立的,而是以相互密切的关系并立着,而且随着科学的每一成功进展,我们会不断发现这些真理之间的新的接触点. ── C.F.Gauss 数学既不严峻,也不遥远,它和几乎所有的人类活动有关,又对每个真心对它感兴趣的人有益. ── R.C.Buck

内容提要  回顾微积分有关知识 连续,多元函数极值,积分等  复习微分方程的求解的解析与数值方法  介绍一类最优问题的求解新框架-变分方法  最速降线求解的仿真方法

1696年John Bernoulli向他的兄长和其他 背景故事 1696年John Bernoulli向他的兄长和其他 数学家挑战性地提出了最速降线(捷线)问题: 一质量为m的质点,在重力作用下从定点A 沿曲线下滑到定点B, A B 试确定一条曲线,使得质 点由A到B下滑时间最短. 假定B比A低,不计摩擦力 和其他阻力等因素.  此问题导致数学新分支的产生.

思考 历史 这是一个求最值的问题  与求函数的极值一样吗?  与求线性规划问题中的极值一样吗?  它的数学形式怎样?  与求函数的极值一样吗?  与求线性规划问题中的极值一样吗?  它的数学形式怎样? 历史 1697年5月号“教师学报”接收了5篇解答报告

贝努利 约翰 Bernoulli,Johann  欧洲著名科学家族  涉猎 微积分、微分方程、解析几 何、 概率论以及变分法 更贡献于物理、化学和天文学  谁发现 L’Hospital 法则  欧拉的指导者和老师  瑞士的骄傲

问题数学形式 设曲线为 满足 y(0)=0, y(c)=H 若 {质点沿 y=y(x) 下滑的时间} 我们要求的是怎样的函数y(x) A B x y c 设曲线为 满足 y(0)=0, y(c)=H {质点沿 y=y(x) 若 下滑的时间} 我们要求的是怎样的函数y(x) 使得T(y) 取得最小值 minT(y)

近似方法 如图建立坐标系,设A为原 点, B为(c,H), 将带状区域 0< y <H用平行于 x 轴的 yk-1 xk-1 yk xk 如图建立坐标系,设A为原 点, B为(c,H), 将带状区域 0< y <H用平行于 x 轴的 直线 y=yk=kH/n 把这区域 分成 n个带状小区域. 在带状域yk-1<y<yk ,可近似认为 而曲线段近似认为是直线段,其长度

于是质点从A到B所需时间近似为 (n -1元函数!) ( 是已知的!)  求这个函数的极小值, 就得到问题的近似解 (为简单计,可取g =1000cm/s2)  可以使用数学软件来求极值,但所得曲线为 离散形式,无解析表达式

求解极值数值方法 令 可令 解出 故

令 为下列方程的解 再将 代入(*)式中,将 用曲线连接即得拟合最速降线,再求出时间 .

利用数学软件求近似最速降线和最短时间 function m6_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; 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-G/(h*s); if f>0 b=c;else a=c; c=(a+b)/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; plot(x,-(0.1:h:H),'*r')

利用数学软件求解得到的曲线

再作分析 质点要走最快的路线(曲线),应该如何变化?  依然用从质点速度变化的角度考虑 设质点从A1经直线 l 到达A2,质点速度在l 的  依然用从质点速度变化的角度考虑 设质点从A1经直线 l 到达A2,质点速度在l 的 上侧为v1,下侧为v2,则质点如何运动才最省时? 如图,若A1,A2到l 的垂足分 A1 A2 1 2 C l O D 别为O,D, A1,A2 到l的距离分别 为a, b, OD =c, 质点经过l于C OC =x 那么质点由A1到A2需时间

惟一驻点满足 A1 A2 1 2 C l O D x a b c-x 也即 这就是光学中的 Snell 折射定律

建立数学模型 若用与x 轴平行的直线将 分析:如图建坐标系, AB 分割成小段, 考虑在第k 层与k +1层质点在曲线上的下 y c k+1 k 层与k +1层质点在曲线上的下 滑,依能量守恒律,可近似 认为质点在每层内的速度不 变,于是依辅助结论知 注意上式对任何k成立,

故导出 令平行线的间距趋于零,我们就得到在曲线 上任何一点 A B x y c  其中 为该点切线与铅垂线 的夹角

导出微分方程 A B x y c   又因  于是得到

一个引理 设集合E0={g(x)C1 │g(a) =g(b)=0} 如果在[a,b]连续函数 f(x)满足 对g (x) E0 ,总有 那么 f (x)  0

另一种方法-变分法 设曲线为 满足 y(0)=0, y(c)=H 在曲线上P(x,y)处质点速度为 又设从A到P的弧长为s,则 B x y c 设曲线为 满足 y(0)=0, y(c)=H 在曲线上P(x,y)处质点速度为 又设从A到P的弧长为s,则 从而质点沿曲线由A到B需时间

设集合 那么我们的问题成为 求某个 使得 引进集合 显然若 是最速曲线函数,则 于是函数 在 取得最小值 故得

为了计算 ,记 那么对 依复合函数求导法 注意第二项

于是导出 上式乘以 这里

得到方程为 注意从降线定义可知 故 解法 1)可求解析解 2)也可以用数值方法,例如欧拉法求解

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

最速降线问题仿真方法Matlab程序 function cycloid(G,H,n) if nargin==2 %两个参数则默认n为100 end g=9.8; h=H/n; minc=0;maxc=1/sqrt(2*g*h*n); x=0;y=0; while abs(G-x)>1e-4 x=0; c=(minc+maxc)/2; %二分法求c值 for j=1:n y=j*h; v=sqrt(2*g*y); x=x+c*v*h/sqrt(1-c^2*v^2); gx(j)=x; gy(j)=y;

if x<G %判断最后一个点与所给点的位置情况 minc=c; else maxc=c; end T=0; for j=1:n v=sqrt(2*g*j*h); if j==1 s=sqrt(gx(1)^2+h^2); s=sqrt((gx(j)-gx(j-1))^2+h^2); T=T+s/v; plot(gx,-gy); T

取G=H=10,n=100

实验任务 1. 分别用数值方法和解析方法求出的最速降线 的曲线和下降时间,将两种结果比较 (设c=π/2, H=1) 2. 在一条直线 l 的上侧有两个点A,B,试找出一条从 A 到B的曲线,使得这曲线绕l 旋转所得的旋转面 的面积最小.设直线l与点A,B在xy 平面,l为x轴, A为(0,(e+e-1)/2), B为(3,(e2+e-2)/2)

3. 圆柱面方程为 用曲线连接面上A(0,0,1), B(1,3,0)两点,求使得 AB 弧长最短的曲线(短程线) 4. 在第3题中,将曲面改为 求在曲面上连接A(1,0,1),B(0,2,2)的最短弧线 (建议以数值和解析两种方法求解加以比较)

5. 若求出最速降线的曲线方程,试将质点从曲线 上任何一点无摩擦地滑到最低点,试求下滑所 需时间 注:若选择本次实验,必须完成任务1、2

谢谢各位!