数值计算方法与算法.pptVIP

数值计算方法与算法 教学目标 掌握常用的数值计算方法 掌握计算方法的数学原理 学会选择恰当的计算方法 学会使用流行的计算软件 教学计划 教材及参考 数值计算方法与算法，张韵华、奚梅成、陈效群编，科学出版社，2006。 科学计算导论，Michael T. Heath著，张威、贺华、冷爱萍译，清华大学出版社， 2005。 数值计算方法解题指导，张韵华编，科学出版社，2003。 网络教学资源 联系方式 第0章 绪论 误差的类型 绝对误差＝真实值－近似值 相对误差＝绝对误差／真实值 误差的来源 原始误差、截断误差、舍入误差 一些例子： 计算地球的体积 计算 计算 如何减小计算误差？ 选择好的算法、提高计算精度 范数的定义 满足非负性,齐次性,三角不等式的实函数 常用的向量范数 常用的矩阵范数 矩阵的谱半径 例：计算矩阵 的范数和谱半径。 例：范数在误差估计中的应用 作业一、145页习题6第1,2题. 作业二、利用公式 编写两个计算ex的C程序(一个用单精度, 另一个用双精度).令x=±1,±5,±10,±15,±20, 比较它们和库函数exp(x)之间的运行时间和计算误差. 思考如何减小误差？ 第1章 插值 作业一、习题1第2,4,6,8,10,12,14,16题。 作业二、在半圆 上随机选取10个点, 构造插值多项式, 画出函数图像, 并比较3种插值方法的计算误差。 作业三、思考3种插值系数之间的关系。比较3种插值方法的优缺点和应用范围。 第2章 数值微分和数值积分 数值微分 差商法 向前差商 向后差商 中心差商 插值法在 x 附近取点(xi,f(xi))构造插值多项式φ。 样条法在 x 附近取点(xi,f(xi))构造样条函数φ。 f’(x)≈φ’(x) 例：用中心差商公式计算f’(xi)。 例：用向后差商公式计算f’’(0.2), f’’(0.4)。 例：设xi=x0+i*h, i=1,...,n。计算φ’(xk)。 解： 误差估计 前后差商中心差商插值微分 数值积分 插值法 若积分公式对任意m次多项式都取等号，则称积分公式具有至少m阶的代数精度。 插值型积分公式的代数精度≥n。 当积分节点 x0,...,xn 给定时，代数精度≥n的积分公式唯一。 例：设xi=a+i*h, i=0,...,n, h=(b-a)/n。 计算Newton-Cotes积分 解： 特别，当n=1,2时，积分公式分别称为 梯形公式 Simpson公式 误差估计 特别，梯形公式和Simpson公式的误差为 代数精度＝1 代数精度＝3 复化数值积分 梯形公式 Simpson公式 Richardson外推法 我们要计算 假设 则 有比　　 和　　 更高的精度。 误差估计 Romberg积分公式 等分的梯形公式， 瑕积分 重积分 Gauss-Legendre积分 定理：假设 满足 则插值积分公式具有2n+1阶的代数精度。 证明： 课本21页性质1.3：若f(x)为m次多项式， 则f[x0,...,xn,x]为m-n-1次多项式。 求多项式空间在内积下的标准正交基。 解法1：对任意基作Gram-Schmidt正交化。 解法2：对任意度量方阵作相合对角化。 解法3：求解正交关系的线性方程组。 解法4：Legendre多项式 作业一、习题2第1~11题。 作业二、使用各种方法计算 的值，其中 0x1，要求误差1e-9。 第3章 曲线拟合的最小二乘法 第4章 非线性方程求根 第5章 解线性方程组的直接法 第6章 解线性方程组的迭代法 作业一、145页习题3、6、7。 作业二、用迭代法计算10阶Hilbert矩阵H的逆矩阵A以及AH-I的元素平方和。 第7章 计算矩阵的特征值和特征向量 作业一、167页习题1(3)、2(2)、3(4)。 作业二、计算10阶Hilbert矩阵的正交相似标准形以及过渡矩阵。 第8章 常微分方程数值解 作业、191-192页习题1--6。 插值 拟合 微积分 非线性方程 线性方程组 特征值问题 常微分方程 结 束 问题：求解n元线性方程组AX=B。 　　　方法？速度？精度 ？存储？ 下三角方程组 上三角方程组 n(n-1)/2次加减法，n(n+1)/2次乘除法。 Gauss消元法解一般方程组 (2n3+3n2-5n)/6次加减法，(n3+3n2-n)/3次乘除法。 追赶法解三对角方程组 3n-3次加减法，5n-4次乘除法。 线性方程组解的精度 矩