2013数值分析课件.ppt

8.1 引言 物理、力学和工程技术中很多问题在数学上都归结为求矩阵的特征值问题。例如，振动问题(大型桥梁或建筑物的振动、机械的振动、电磁震荡等)，物理学中的某些临界值的确定。它们都归结为下述数学问题。 8.2 幂法及反幂法 一、幂法 幂法是一种求实矩阵A的按模最大的特征值λ1及其对应的特征向量x1的方法。特别适合于大型稀疏矩阵。 上述迭代法是一种逐次逼近法，其基本思想是将隐式方 程归结为一组显示的计算公式，就是说，迭代过程实质上是 一个逐步显示的过程。 k xk k xk k xk 0 1 2 1.5 1.35721 1.33086 3 4 5 1.32588 1.32494 1.32476 6 7 8 1.32473 1.32472 1.32472 继续迭代下去已经没有必要，因为结果显然会越来越大， 不可能趋于某个极限。这种不收敛的迭代过程称作是发散的。 一个发散的迭代过程，纵使进行了千百次迭代，其结果也毫 无价值。因此，迭代格式形式不同，有的收敛，有的发散，只 有收敛的迭代过程才有意义，为此要研究不动点的存在性及迭 代法的收敛性。 二、不动点的存在性与迭代法的收敛性 k xk k xk 1 2 3 1.484248034 1.472705730 1.468817314 4 5 6 1.467047973 1.466243010 1.465876820 k xk k xk 1 2 3 0.1 0.089482908 0.090639135 4 5 6 0.090512616 0.090526468 0.090524951 三、局部收敛性与收敛阶 k xk 迭代法(1) 迭代法(2) 迭代法(3) 迭代法(4) 0 1 2 3 ? x0 x1 x2 x3 ? 2 3 9 87 ? 2 1.5 2 1.5 ? 2 1.75 1.73475 1.732631 ? 2 1.75 1.732143 1.732051 ? k xk | xk- xk-1| 0 1 2 3 1.5 1.481248 1.482671 1.482563 0.018752 0.001423 0.000108 7.3 牛顿法 一、牛顿法及其收敛性 k xk k xk 0 1 0.5 0.57102 2 3 0.56716 0.56714 二、牛顿法应用举例 k xk 0 1 2 3 4 10 10.750000 10.723837 10.723805 10.723805 三、简化牛顿法与牛顿下山法 k xk xk xk f(xk) 0 1 2 3 4 1.5 1.34783 1.32520 1.32472 0.6 17.9 发散 0.6 -1.384 1.140625 -0.656643 1.36181 0.1866 1.32628 0.00667 1.32472 0.0000086 四、重根情形 k xk (1) (2) (3) 0 1 2 3 x0 x1 x2 x3 1.5 1.458333333 1.436607143 1.425497619 1.5 1.416666667 1.414215686 1.414213562 1.5 1.411764706 1.414211438 1.414213562 7.4 弦截法 k xk k xk 0 1 2 0.5 0.6 0.56532 3 4 0.56709 0.56714 知 识 结 构 图 七 方 程 近 似 求 根 基本概念（单根、重根、有根区间、不动点、收敛阶） 求根方法 二分法及其收敛性 不动点迭代法及其收敛性定理 （不动点迭代法的加速技巧） 牛顿迭代法及其收敛性 插值型迭代法（多点迭代） 弦截法 抛物线法 第八章 矩阵特征值问题计算 内容提要 8.1 引言 8.2 幂法及反幂法 6.2 基本迭代 一、雅可比迭代法 二、高斯—塞德尔迭代法 SOR迭代法的计算公式:对k=0,1,…, 三、逐次超松驰(SOR)迭代法 说明: 1)ω=1,即为GS（高斯-赛德尔迭代法）; 2)ω1，称为超松驰法； ω1，称为低松驰法； 3) SOR方法每迭代一次主要运算量是计算一次矩阵 与向量的乘法。 例6-3 用SOR迭代法解线性代