数值计算与方法总复习 .ppt

复化Simpson求积公式 将 [a, b] 分成 2n 个小区间，在每两个相邻小区间 上用Simpson求积公式，再将这些小区间上的数值积分累加起来，就得到区间[a, b]上的数值积分。这种方法称为复化Simpson积分。 根据Simpson求积公式，有 将 2n 个小区间上的积分累加起来，得 复化Simpson求积公式 截断误差为： 收敛性： 复合Simpson公式分解 系数首尾为1，奇数点为4，偶数点为2 例 计算 使计算结果有5位有效数字。 分别用复化梯形和复化Simpson求积公式，求分点数。 解 由复化梯形求积公式得： 由复化Simpson求积公式得： 数值微分 x-h x x+h B C A T f (x) 差商型求导公式 自然，而又简单的方法就是，取极限的近似值，即差商。 导数与差商的关系 x-h x x+h B C A T f (x) 从截断误差的角度看，步长越小，计算结果越准确；从舍入误差的角度来看，步长不宜太小。下面将介绍具有较高精度的数值微分公式。 差商型求导公式的余项 插值型求导公式 由于采取的是 n 次Lagrange插值多项式，而高次插值会产生 Runge现象，因此实际应用中多采用低次插值型求导公式。 （1）两点数值微分公式 (2) 三点数值微分公式 同样，针对m也可扩展。 第七章 常微分方程数值解初步 欧拉折现公式及其变形 梯形公式 龙格-库塔格式 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 一般最小二乘拟合多项式 对于离散数据: (xk, yk), k=1,2,…,m, 用 n (nm) 次多项式来拟合曲线。设多项式 的系数是下述极小值问题的解： 一阶必要条件： 直接计算易得 故 或 称为正规方程组。可表示为 例5.1.4 用二次多项式函数拟合如下数据: xi -3 -2 -1 0 1 2 3 yi 4 2 3 0 -1 -2 -5 解 设 p(x) = a0 + a1x + a2 x2, 形成正规方程组： m =7. 约定 直接计算有： 方法一 一般地， 为定义在X上的广义多项式，记为 定义残差的平方和： 最小二乘问题为：求解极小值问题 正规方程组便可化为： 将其表示成矩阵形式 例5.1.6 用给定数据，求经验公式 f (x) = a + bx3. x = -3 -2 -1 2 4 y = 14.3 8.3 4.7 8.3 22.7 解 约定 直接计算得 法方程 于是法方程组为： 所求经验公式为：f (x) =10.675 + 0.137 x3. 函数的最佳平方逼近 上述问题等价于求多元函数 的最小值。 由多元函数取极值的必要条件 得 于是有 上述方程组称为正规方程组。也可以写为 例5.2.1 解 取基函数为 建立正规方程组： 根据内积公式，可得 正规方程组为： 所求拟合函数为： 第六章 数值微积分 基本概念 代数精度 算法及其收敛性 Newton-Cotes 求积公式 复化梯形公式和复化辛普森公式 两点微分公式和三点微分公式 定义 设有近似式 则称该近似式具有m次的代数精度。 代数精度也称 代数精确度 代数精度 容易验证：梯形求积公式的代数精度为1， Simpson求积公式的代数精度为3. 例 试确定下面积分公式中的参数使其代数精确度尽量高. 解 因此 所以该积分公式具有3次代数精确度。 例 求以下微分公式的代数精度 解 故代数精度为2. 对于一般的n阶插值型求积公式 由误差公式 可知，其代数精度至少为n. 对于Newton－Cotes求积公式，我们有： 定理 定理 Newton-Cotes 求积公式 在微积分中，定积分是 Riemann 和的极限，即 数值积分就是取定积分极限中的有限项的和，即 其中，xk 称为积分节点,Ak 称为积分系数。 我们的任务就是确定积分系数Ak ,使 I [ f ]≈In[ f ]. 最常用的方法就是用插值多项式近似代替被积函数 f (x) 来确定Ak. 得数值积分公式： 截断误差为： （1