数值分析重要公式最终修改版03.docVIP

1、Doolittle分解法 设A的各阶主子式非奇异，那么 A = LDU = L(DU) = LU1 (Doolittle分解, L为单位下三角阵) = (LD)U= L1U (Crout分解, U为单位上三角阵)。 例1 用Doolittle分解法解方程组: 解 记 ，设 由矩阵相等得 ，，从而 由 ， 解得：；由 ， 解得： (1) Jacobi迭代格式 ， ，，为初始向量（矩阵形式） A=D+L+U例 设 其中 ，试求 ，， 解：， ，，， 求 2 Gauss–Seidel迭代法 （1）迭代格式 ， ，，为初始向量（矩阵形式） 例题0 构造收敛的Gauss–Seidel迭代格式说明收敛理由 行变换系数矩阵主对角线元素按行严格占优迭代格式收敛， 任选；k=0,1,2…. 例题1：确定 并求代数精度 令 求 例题2 三次样条函数 【0,1】且求在【1,2】表达式 设 求A B C D 例题3 梯形法求初值 公式 例题4合曲线 ， A= b= a= 2散点拟合例1 已知 0 1 2 3 1 2 5 7 9 4 2 1 试求它的最小二乘拟合曲线（。 解：（1）取直角坐标系，描点，由图可知，这些点位于一条双曲线附近。取 ，即，； ，=1.842857, =1.310408，=16， =11.542857 (3) 解方程组 ， 例5设，，求它在上的最佳平方逼近多项式。 解： ，， ，，，，，， 解方程组 ，， 例题6Householder做正交相似变换化为三对角矩阵 取 因为 = 方法二1），，，， （2） 得 HAH 例题7判断收敛及收敛速度 因0所以有m阶收敛速度1．Newton迭代格式 ， 为初值， 令 所以在只有一个实根 至少平方收敛 例题8angrange拉格朗日插值多项式注： n=1时，线性插值公式， , 或 n=2时，二次插值（或抛物线插值）公式： Newton插值多项式插值公式余项（即截断误差）= 例题9 复化梯形公式计算 将[a, b] n等分，，，，在每个子区间上应用梯形公式 得复化梯形公式： ， 复化梯形公式的截断误差： ， （3） 复化抛物线公式： ， 复化抛物线公式的截断误差： ， 例1 若用复化梯形公式和复化抛物线公式计算积分，问积分区间要等分多少才能保证有5位有效数字? 解 由余项公式 因为 f (x) = f ((x) = … = e x, b – a = 1 故 或 由 | R | ( 10 – 4/2, 得 n ( (e(10 4/6)1/2 =68 (梯形公式), 或 n ( (e(104/1440)1/4 =2.1(抛物线公式, 因区间是2n, 实际上还要乘以2, 故6等分即可). 例题10逐次分半复化梯形公式截断误差不超过步长 为绝对误差限满足时就停止计算m=0时 m=1, m=2, m=3, m=4, 例题11 Euler法求初值 原始Euler法=局部截断误差 向后Euler法截断误差 求最佳逼近元素的步骤： 所以 1 1/2 1/3 2/3 2/5（2）求解方程组（4）得： （3）将代入（1）式即得 Simpson公式截断误差 3次代数精度 梯形公式截断误差 1次代数精度