Lecture_01.ppt

;“诸位在校，有两个问题应该自己问问， 第一，到浙大来做什么？ 第二，将来毕业后做什么样的人？” ------- 竺可桢;这一讲的主要内容;数值计算（分析）是做什么的？;今有上禾三秉，中禾二秉，下禾一秉，实三十九斗； 上禾二秉，中禾三秉，下禾一秉，实三十四斗； 上禾一秉，中禾二秉，下禾三秉，实二十六斗。 问上、中、下禾实一秉各几何？ 答曰：上禾一秉九斗四分斗之一。中禾一秉四斗四分斗之一。下禾一秉二斗四分斗之三。-------《九章算术》;线性方程组的数值方法！;2、天体力学中的Kepler方程;全球定位系统：在地球的任何一个位置，至少可以同时收到4颗以上卫星发射的信号 ; 表示地球上一个接收点R的当前位置，卫星Si的位置为 ，则得到下列非线性方程组; 非线性方程组的数值方法!;4、已经测得在某处海洋不同深度处的水温如下： 深度（M） 466 741 950 1422 1634 水温（oC）7.04 4.28 3.40 2.54 2.13 根据这些数据，希望合理地估计出其它深度（如500米，600米，1000米…）处的水温 ;5、用比较简单的函数代替复杂的函数;6、人口预测 ;曲线拟合!;7、铝制波纹瓦的长度问题; 这个问题就是要求由函数f(x)=sin x给定的曲线从x=0到x=48英寸间的弧长L. 由微积分学我们知道,所求的弧长可表示为: ;8、生物化学反应的例子 ;通过建模可以得到如下方程组 ; ; 矩阵特征值问题的数值方法!;数学 模型 ; ;科学计算的重要性;本课程数值分析讲课范围 ;主要教材 颜庆津，数值分析。北航出版社，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;本门课程的特点;学习要求;其它要注意的几点;数值分析中的基本概念：误差; 误差的来源 ;误差的分类;误差的分类（续）;误差的分类（续）;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.