[理学]数值分析 第1章 绪论.ppt

应用问题举例 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、用比较简单的函数代替复杂的函数 误差为最小，即距离为最小 （在不同的度量意义下） 6、人口预测 下面给出的是中国1900 年到2000年的人口数， 我们的目标是预测未来 的人口数（数据量较大时） 7、铝制波纹瓦的长度问题 这个问题就是要求由函数f(x)=sin x 给定的曲线从x=0到x=48英寸间的弧长L. 由微积分学我们知道,所求的弧长可表示为: 教材 现代科学与工程计算 孟大志 刘伟（高等教育出版社） 参考书目 ?数值分析 孙志忠 袁慰平等（东南大学出版社,第二版） ? 应用数值方法 使用MATLAB和C语言 Robert J.Schilling Sandra L.Harris （机械工业出版社） ?数值分析基础教程 李庆扬 编 （高等教育出版社） ? 现代数值分析 李庆扬、易大义、王能超 编著 （高等教育出版社） ?数值分析与科学计算 Jeffery J.Leader 著，张威，刘志军，李艳红等译，（清华大学出版社) §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之间。 在 机器字长有限 —— 舍入误差 用计算机、计算器和笔算，都只能用有限位 ? = 3.1415926… 小数来代替无穷小数或用位数较少的小数来 代替位数较多的有限小数，如： 四舍五入后…… 在数值计算方法中，主要研究截断误差和舍入误差 （包括初始数据的误差）对计算结果的影响！ 二、 误差的概念 1、绝对误差与绝对误差限 例 :若用以厘米为最小刻度的尺去量桌子的长，大约为1.45米，求1.45米的绝对误差。 1.45米的 绝对误差=？ 不知道！ 是近似值 的绝对误差,简称为误差。 定义：设 是准确值，为 的一个近似值，称 但实际问题往往可以估计出 不超过某个正数 ， 即 ，则称 为绝对误差限，有了绝对误差限 就可以知道 的范围为 即 落在 内。 在应用上，常常采用下列写法来刻划 的精度。 2、相对误差与相对误差限 定义：设 是准确值， 是近似值，是近