学习Euler,是认识数学的最好途径。高斯1654最全版.pptVIP

MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK 学习Euler,是认识数学的最好途径。 ― 高斯 纯粹数学与应用数学是理解世界及其发展的一把主要钥匙。 ― 里约热内卢宣言 * 数学实验 导弹跟踪问题 上海交通大学数学科学学院 * KUHKUK 实验目的 微分方程数学建模 回顾微分方程解析解法 介绍Euler方法 介绍改进Euler方法 介绍仿真方法 * 实际问题 某军的一导弹基地发现正北方向120 km处 海面上有敌艇一艘以90 km/h的速度向正东方向 行驶.该基地立即发射导弹跟踪追击敌艇,导弹 速度为450 km/h，自动导航系统使导弹在任一 时刻都能对准敌艇. 试问导弹在何时何处击中 敌艇？ （追踪问题） * 120 P(x, y) 当t =0时,导弹位于原点O, 敌艇位于(0,H)点, H =120 (km). 当时刻t ,导弹位于P(x(t),y(t)) 敌艇位于(90t, H)点，根据题设可建立如下数学模型. O 建立坐标系 x * 时刻t 导弹的位置为P(x(t),y(t)),其 Vw = 450 (km/h) Ve = 90 (km/h) 数学模型 速度可由水平分速度与垂直分速度合成 导弹方向指向敌艇,故有 * 导出一阶微分方程组 * 解析方法 技巧性较强: 消去 t ,化为二阶方程 连同初始条件 设 降阶得到一阶可分离变量方程 * 导弹击中艇时y=H, 得到此时其位置和时刻 进而解出 * Euler 方法十分简单, 就是用差商代替微商 数值方法 (Euler方法) 第k+1步的值之间的关系式 即采用 通常取Δt为常数τ,就得到由第 k 步的值到 * 伟大的欧拉（Leonhard?Euler，1707～1783 ） ? 著作辉煌：涉及几乎每一数学领域 ? 数学天才：约翰·伯努利指引 13岁进巴塞尔大学 ? 毅力惊人：失明17年，完成书和约400篇论文 ? 非凡记忆和心算力：前100个素数的前六次幂 886本书和论文 ，汇成100巨册 ? 风格高尚：支持年青的拉格朗日 ? 热爱生活 * ? 瑞士法郎上的欧拉 ? 邮票上的欧拉 * 解析法降阶后方程的Euler格式 和初始条件 由方程 ? Euler 迭代格式 h=H/n为步长 Matlab 文件m4_1 此例告诉我们高阶方程可化为一阶方程组 * 设 t = tk+1 时导弹位置为(xk+1 ，yk+1) 计算到 yk－1 H, yk≥H 停止，取L ≈ xk 原模型的 Euler 迭代格式 * 利用Matlab 参见教材中介绍的m文件并运行 则 L≈x6= 29.195 T ≈0.3244 例如 m4_3,当τ = 0.05 k tk xk yk 1 0.05 0.000 00 22.500 00 2 0.10 1.037 36 44.976 07 3 0.15 3.412 05 67.350 41 4 0.20 7.646 15 89.448 43 5 0.25 14.867 90 110.757 96 6 0.30 29.194 80 128.107 02 * 当τ = 0.001 L≈x278=25.049 T ≈0.27833 τ 0.1 0.05 0.005 0.0001 n 3 6 56 278 L 22.67495 29.19480 25.66731 25.04935 T 0.25194 0.32439 0.28519 0.27833 表1是对应不同的τ ，用Euler法所得相应的步长推进次数和计算结果 表1 * 关于微分方程数值方法 大多数微分方程－即使形式十分简单， 也可能无法用解析方法求解，因此数值方法 极其重要。 * 改进Euler法（预测－校正法） 一维例子 方程与初始条件 Euler 格式 改进Euler 格式 预测值 校正值 * 改进Euler法的几何解释 若 的解为x(t) , 代入方程从tk→tk+1 积分导出 O tk f t f =f (x(t),t) tk+1 Euler法 xk+1- xk=τ f (xk ,tk) 改进Euler法 xk+1- xk=τ[ f (xk ,tk)+ * MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK MKKLL JHKJK MKK