复习_仔细做几遍-副本.pptVIP

以下是复习内容, 要非常熟练!带计算器 1. Langrange 例1 拉格朗日余项定理 例2 已知sin0.32=0.314567, sin0.34=0.333487, 用三角分解 更简单的Cholesky分解 不动点的收敛性 说明: 定理2 例2 讨论迭代公式 8. 局部收敛 例 收敛阶 定理4 简单迭代法的敛速 例3 简单迭代法求 9.超定方程组 定理 例题 对应超定方程组: 用LLT分解做: 10.Gauss-Seidel迭代法 Gauss-Seidel公式 例3 用G-S迭代法解 11. 迭代法的收敛性p190 例 对于Ax=b, 13.证明中点公式 所以 x是最小二乘解 ATAx=ATb 对于||b-Ax||=min, 假设为y 实际计算时 而是 ATA的逆, (1) 计算C=ATA; d=ATb; (2) 计算C的 Cholesky分解: C=LLT; 则原方程组的解等价于 x=(LLT)-1d =L-TL-1d 先解 Ly=d得y; 再解x= 并非计算 L-Ty. 记住 正则化方法 用二次多项式拟合数据 x -2 -1 0 1 2 y 0 1 2 1 0 解 设多项式为: 基函数: f2(x)=x, f3(x)=x2; y=c0+c1x+c2x2; 矩阵A为: f1(x)=1, A= 最小二乘问题为: ||b-Ax||2= min. 左乘矩阵AT 即: 效果图 A= C=ATA = =LLT 其中L= d=ATb = ATAx=ATb LLTx=d Ly=d得y= 再解LTx=y x= 故拟合曲线为: y= 1.657 -0.429x2 在A=M-N中, 取M为A的 即 下三角部分 D-L, 又称: 高斯-塞德尔迭代法 =-L+D-U 取M=D-L,N=U, 故G-S公式为: 迭代矩阵 k=0,1,2,... 以上是 GS法的 矩阵表示方法. 换成(k+1) 可从 即GS法的 雅可比公式 J修改而得: 计算格式. 解 x1={3+x2+x3}/5 x1 ={3+x2 +x3 }/5 x2 ={3+x1 +x3 }/5 x2 ={3+x1 +x3 }/5 x3 ={3+x1+x2 }/5 x3 ={3+x1 +x2 }/5 (k+1) (k) (k) (k+1) (k+1) (k) (k+1) (k+1) (k+1) 取初始向量 x(1)=( )T x(9)= x(2)= x(0)=(0 0 0)T, 迭代9次: …… 0.864000 0.600000 0.720000 x1 ={3+x2 +x3 }/5 x2 ={3+x1 +x3 }/5 x3 ={3+x1 +x2 }/5 (1) (0) (0) (1) (1) (0) (1) (1) (1) x1 ={3+x2 +x3 }/5 x2 ={3+x1 +x3 }/5 x3 ={3+x1 +x2 }/5 (2) (1) (1) (2) (2) (1) (2) (2) (2) (0.916800 0.956160 0.974592)T (1.00000 1.00000 1.000000)T x1 ={3+x2 +x3 }/5 x2 ={3+x1 +x3 }/5 x3 ={3+x1 +x2 }/5 (9) (8) (8) (9) (9) (8) (9) (9) (9) 本节讨论J和GS的收敛性. 将 x(k+1)=(I-D-1A)x(k)+D-1b x(k+1)=(D-L)-1Ux(k)+(D-L)-1b 统一 写成 x(k+1)= k=0,1,2,... 设Ax=b, B x(k)+ f x(k+1)=Bx(k)+f 收敛 (1) 充要条件: r(B)1; (2) 充分: ||B||1; (3) 充分: A对角占优; (4) 对称正定条件 则 Th7 Th7下面 Th9 Th10 解 由于 不好判断. 得 r( (D-