数值分析—第一讲误差.pptVIP

;“诸位在校，有两个问题应该自己问问， 第一，到浙大来做什么？ 第二，将来毕业后做什么样的人？” ------- 竺可桢;这一讲的主要内容;数值计算（分析）是做什么的？;1、天体力学中的Kepler方程;全球定位系统：在地球的任何一个位置，至少可以同时收到4颗以上卫星发射的信号 ; 表示地球上一个接收点R的当前位置，卫星Si的位置为 ，则得到下列非线性方程组; 非线性方程组的数值方法!;3、已经测得在某处海洋不同深度处的水温如下： 深度（M） 466 741 950 1422 1634 水温（oC）7.04 4.28 3.40 2.54 2.13 根据这些数据，希望合理地估计出其它深度（如500米，600米，1000米…）处的水温 ;4、铝制波纹瓦的长度问题; 这个问题就是要求由函数f(x)=sin x给定的曲线从x=0到x=48英寸间的弧长L. 由微积分学我们知道,所求的弧长可表示为: ;数学 模型 ; ;科学计算的重要性;本课程数值分析讲课范围 ;主要教材 颜庆津，数值分析。北航出版社，2006。 数值逼近参考书 李庆扬，王能超，易大义：数值分析。清华大学出版社，2001。 李岳生，黄友谦：数值逼进。北京人民教育出版社，1979。 王德人、杨忠华，数值逼近引论，高等教育出版社，1990 王仁宏，数值逼近，北京：高等教育出版社，1999 徐萃薇：计算方法引论。北京高等教育出版社，1985。 数值代数参考书 曹志浩，数值线性代数，上海：复旦大学出版社，1996. 徐树方，矩阵计算的理论与方法，北京大学出版社，1995. 微分方程数值解参考书 李立康、於崇华、朱政华，微分方程数值解法，复旦大学出版社，1999 陆金甫、关? 治，微分方程数值解法（第二版），北京：清华大学出版社，2004 综合类（数值分析与科学计算、习题、实验等）参考书 蔡大用，数值分析与实验学习指导，北京：清华大学出版社，2001 Numerical Recipes(数值方法库） in C/Matlab/Fortran/C++, 其他 /wiki/Numerical_analysis Software：IMSL，NAG，MATLAB;本门课程的特点;学习要求;其它要注意的几点;数值分析中的基本概念：误差; 误差的来源 ;误差的分类（1/4）;误差的分类（2/4）;误差的分类（3/4）;S4 ;§1.2.2 误差与有效数字 /* Error Significant Digits */;;有效数字 定义：如果近似值a的误差限是某一位数的半个单位， 则称a准确到小数点后n位，并从第一个非零的数字到这一位的所有数字均为有效数字。 例：π＝3.1415926535, 3.1416有五位有效数字,误差限为0.00005。 例： a=0.003400±0.5E-5近似值准确到小数点后 五位，有三位有效数字。;§1.2.3 函数求值的的误差估计 /*Error Estimation for Functions*/; f 的条件数在某一点是小\大，则称 f 在该点是好条件的 /* well-conditioned */ \坏条件的 /* ill-conditioned */。;特别情况[u=f（x）]：;误差的四则运算;例 求积分;两种算法与真实值的比较;说明; 在我们今后的讨论中，误差将不可回避， 算法的稳定性会是一个非常重要的话题。 ;4. 数值运算中误差分析的方法与原则;§5 Remarks;避免相近的两数相减 (会耗失许多有效数字,可以用数学公式化简后再做). 例 各有五位有效数字的23.034与22.993相减. 23.034-22.993=0.041 0.041只有两位有效数字,有效数字的耗失,说明准确度减小,因此,在计算时需要加工计算公式,以免这种情况发生. 例 当 x 较大时,计算 ;防止”大数”吃”小数” 当两个绝对值相差很大的数进行加法或减法运算时,绝对值小的数有可能被绝对值大的数吃掉从而引起计算结果很不可靠. 例:求 两者结果不同,因为计算机计算时做加减法要 “对阶”,“对阶”的结果使大数吃掉了小数.产生了误差.为了避免由于上述原因引起的计算结果严