10-11数值计算辅导.pptVIP

其中 U的第一行 L的第一列 U的第r行 L的第r列 对于线性方程组 系数矩阵非奇异,经过LU分解后 线性方程组可化为下面两个三角形方程组 2、矩阵的范数和条件数 (4) 根据算子范数的不同也有不同的条件数: 总复习提要及例题 一、误差 1、误差限的定义及求法 绝对误差限或误差限, 注 例 解: 一般地,凡是由精确值经过四舍五入得到的近似值, 其绝对误差限等于该近似值末位数字的半个单位。 二、插值法 （一）Lagrange插值 Lagrange线性插值基函数为 Lagrange线性插值多项式为 1、线形插值 其中 2、一般插值 例: 解: 且 作线性插值 解: Lagrange插值基函数为 Lagrange线性插值多项式为 （二）Newton插值 1、差商表 规定函数值为零阶差商 差商表 例 给定数据表f(x)=lnx数据表 xi 2.20 2.40 2.60 2.80 3.00 f(xi) 0.78846 0.87547 0.95551 1.02962 1.09861 1.构造差商表 2.用二次Newton差商插值多项式，近似计算f(2.65)的值3.写出四次Newton差商插值多项式N4(x) 解:差商表 N4(x)= 0.78846 +0.43505(x-2.20） - 0.087375(x-2.20)(x-2.40) +0.0225(x-2.20)(x-2.40)(x-2.60) -0.00755(x-2.20)(x-2.40)(x-2.60)(x-2.80) 只须解矩阵方程 由 因此求最小二乘解转化为 三、曲线拟合 一种简单常用的情况 基函数之间的内积为 已知一组实验数据如下，求它的拟合曲线. 例题 解 从图中看到各点在一条直线附近，故可选择线性函数作拟合曲线， 将所给数据在坐标纸上标出，见图3-4. 图3-4 令 这里 故 四、数值积分 （一）Newton-Cotes求积公式及余项 Newton-Cotes公式是指等距节点下使用Lagrange插值 多项式建立的数值求积公式 至少具有n次代数精度并且n为偶数时至少具有n+1次代数精度. 一般地 , Newton-Cotes公式的截断误差为 记为 余项为 1.梯形公式及其余项 具有1次代数精度 记为 Simpson公式的余项为 Simpson公式具有3次代数精度 2.Simpson公式及其余项 记为 Cotes公式的余项为 Cotes公式具有5次代数精度 3.Cotes公式及其余项 例 用梯形公式、Simpson公式和Cotes公式求积分 的近似值。 解 I?T=1/2(4+2)=3 I?S=1/6(4+12.8+2)=3.13333 I?C=1/90(28+…+14)=3.14212 （二）复化求积公式及余项 （一）复合求积公式 各节点为 复合梯形公式 则有误差估计 要使|I-Tn|? ,只要 如果记M2= 复合Simpson公式 如果记M4= 则有误差估计 要使|I-Sn|? ,只要 的数据表 已知函数 分别用复化梯形公式、复化Simpson公式计算积分 解 例 0.8414710 1 0.9361556 5/8 0.9896158 1/4 0.8771926 7/8 0.9588511 1/2 0.9973978 1/8 0.9088517 3/4 0.9767267 3/8 1 0 ?(xk) xk ?(xk) xk ?(xk) xk 的近似值 于是有 对复化梯形公式,若使|I-Tn|10-6, 只要 实际上,S3=0.9460838. 故应取n=167. 对复化Simpson公式,若使|I-Sn|10-6, 只要 故只需取n=3. 五、解线性方程组的直接法 1、LU分解 * *