数值计算办法及算法.pptxVIP

数值计算方法与算法;第0章 绪论; 数学建模 数值计算;误差的类型 绝对误差＝真实值－近似值 相对误差＝绝对误差／真实值 误差的来源 原始误差、截断误差、舍入误差;一些例子： 计算地球的体积 计算 计算 如何减小计算误差？ 选择好的算法、提高计算精度 范数的定义 满足非负性,齐次性,三角不等式的实函数;常用的向量范数 常用的矩阵范数 矩阵的谱半径 例：计算矩阵 的范数和谱半径。 例：范数在误差估计中的应用;第1章 插值;函数逼近 用未知函数f(x)的值构造近似函数φ(x)。要求误差小、形式简单、容易计算。 常用的函数逼近方法 插值：φ(xi)=yi, i=0,1,…,n. 拟合：||φ(x)-f(x)||尽可能小 通常取 φ(x) = a0φ0(x) + … + anφn(x)，其中 {φi(x)}为一组基函数。;多项式插值 给定平面上n+1个插值点(xi,yi), 构造n次多项式φ(x), 满足φ(xi)=yi, i=0,1,…,n.;???项式插值;Lagrange插值;Newton插值;差商表;差商的性质 以x0,…,xn为节点的n次插值多项式φ(x)的首项系数等于f[x0,…,xn]。 证明：分别以x0,…,xn-1和x1,…,xn为节点构造n-1次插值多项式φ1(x)和 φ2(x)，则有对n用归纳法。 f[x0,…,xn]与x0,…,xn的顺序无关。;误差估计： 证明：设 ，则 　 有n+2个零点。 根据中值定理，存在 于是 。;Hermite插值 给定平面上n+1个插值点(xi,yi,mi), 构造2n+1次多项式φ(x), 满足φ(xi)=yi, φ’(xi)=mi, i=0,1,…,n.;单项式基函数;误差估计： 证明：设 ，则 　 有2n+3个零点。根据中值定理，存在 于是 。;Runge现象：并非插值点取得越多越好。 解决办法：分段插值;三次样条插值 给定平面上n+1个插值点(xi,yi), 构造分段三次多项式φ(x), 满足φ(xi)=yi, φ’(x)可微，φ”(x)连续。;第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多项式;第3章 曲线拟合的最小二乘法;曲线拟合 对区间 I 上的连续函数 f ，构造特定类型的函数 φ 使 φ≈f 。 对离散数据序列 (xi, yi), i=1,2,…,m，构造特定类型的函数 φ 使 φ(xi)≈yi 。 最小二乘法 求 φ 使 最小。 求 φ 使 最小。;多项式拟合 其中 是标准正交基， 。 求 使 最小。; 奇异值分解 Moore-Penrose 广义逆 矛盾方程组 　　　　的解;其他类型的离散数据拟合 ;第4章 非线性方程求根;问题 求f(x)