数值绪论辩析.ppt

任课教师： 王 永 丽 E-mail: wangyongli@ Mobile: Office: 信息楼 J13-108 应用问题举例 2、天体力学中的Kepler方程 全球定位系统：在地球的任何一个位置，至少可以同时收到4颗以上卫星发射的信号 3、全球定位系统（Global Positioning System, GPS) 4、已经测得在某处海洋不同深度处的水温如下： 深度（M） 466 741 950 1422 1634 水温（oC）7.04 4.28 3.40 2.54 2.13 根据这些数据，希望合理地估计出其它深度（如500米，600米，1000米…）处的水温. 5、人口预测 下面给出的是中国1900 年到2000年的人口数， 我们的目标是预测未来 的人口数（数据量较大时） 6、铝制波纹瓦的长度问题 这个问题就是要求由函数f(x)=sin x 给定的曲线从x=0到x=48英寸间的弧长L. 由微积分学我们知道,所求的弧长可表示为: 威尔金森（James Hardy .Wilkinson，1919-1986） Wilkinson是数值分析和数值计算的开拓者和奠基人。1940 年，开始研究弹道的数学模型与数值计算。 1946 年成为Turing 的助手，协助设计 Pilot ACE 计算机。1969年他当选为英国皇家学会院士；1970年工业和应用数学会(s1am)授予他冯·诺伊曼奖；1987年他获得美国数学会的chauvenet奖。著名的美国阿尔贡国家实验室曾聘威尔金森为荣誉高级研究员并两次向他授奖。 Wilkinson在数值分析研究领域作出了杰出贡献，是数值计算的早期开拓者，其工作加速了数字计算机 ( 在科学计算中 ) 的使用。他研究的主要问题是线性代数方程组和矩阵特征值问题的数值解法，特别是他的向后误差分析法 (backward error analysis)的创造性工作奠定了数值分析和数值计算早期的理论基础。 1975 年 J. H. Wilkinson成为第五位图灵奖获得者。 教材 现代科学与工程计算 孟大志 刘伟（高等教育出版社） 参考书目 ? 数值分析 孙志忠 袁慰平等（东南大学出版社,第二版） ? 应用数值方法 使用MATLAB和C语言 Robert J.Schilling Sandra L.Harris （机械工业出版社） ? 数值分析（第4版） 李庆杨 王能超 易大义 华中科技大学出版社 ? 数值分析与科学计算 Jeffery J.Leader 著，张威，刘志军，李艳红等译，（清华大学出版社) ?数值分析精品课程网站：/na/html/index_common.htm §2 算 法 一、算法的概念 描述算法可以有不同的方式。例如，可以用日常语言 和数学语言加以叙述，也可以借助形式语言（算法语言） 给出精确的说明，也可以用框图直观地显示算法的全貌。 定义：由基本运算及运算顺序的规定所构成的完整的 解题步骤，称为算法。 例：求解二元一次联立方程组 用行列式解法：首先判别 （1）如果 ，则令计算机计算 输出计算的结果x1，x2。 （2）如果D= 0，则或是无解，或有无穷多组解。 是否为零，存在两种可能： 令 通过求解过程，可以总结出算法步骤如下： S2 计算 S3 如果 则输出原方程无解或有无穷多组解的信息; 否则 S1 输入 S4 输出计算的结果 输入 D=a11a22-a12a21 D=0 开始 输出 x1, x2 结 束 No 输出无解信息 Yes 二、算法优劣的判别 ? 计算量的大小 ? 存贮量 ? 逻辑结构 例：用行列式解法求解线性方程组: n阶方程组，要计算n + 1个n阶行列式的值， 总共需要做n! (n - 1) (n + 1) 次乘法运算。 n=20 需要运算多少次？ n=100? 一、误差的来源与分类 从实际问题中抽象出数学模型 —— 模型误差 例:质量为m的物体，在重力作用下,自由下落， 其下落距离s 与时间t 的关系是： 其中 g 为重力加速度。 §3 误 差 通过测量得到模型中参数的值 —— 观测误差 求近似解 —— 方法误差 (截断误差） 例如，当函数 用Taylor多项式 近似代替时，数值方法的截断误差是 与0之间。 在 机器字