数值分析幂法和反幂法.docxVIP

数值分析幂法和反幂法 数值分析幂法和反幂法 数值分析幂法和反幂法 幂法和反幂法求矩阵特点值 目 具 随机 生一 称矩 ， 不一样的原点位移和初 (起码取 3 个)分 使用 体 法求 算矩 的主特点 及主特点向量，用反 法求 算矩 的按模最小特点内 及特点向量，并比 不一样的原点位移和初 明收 。 容 1. 真 ，认识 的数学原形； 要2. 适合 求解的数 算方法； 求3. 程序并 行 算； 4. 果 行解 明； 于 法和反 法求解矩 特点 和特点向量的 将从 剖析，算法 和流程 ，理 依照，程序及 果 行 述 。 一．问题的剖析： 求 n 方 A 的特点 和特点向量，是 算中经常遇到的 ，如： 机械、 构或 磁振 中的固有 等。 于 n 矩 A，若存在数 和 采 n 向量 x 足 Ax= x （1） 用 称 矩 A 的特点 ， x 相 的特点向量。 方 由高等代数知 可知，特点 是代数方程 法 | I-A|= n +a1 n 1 +?+a n 1 +a n =0 （2） 及 的根。从表面上看，矩 特点 与特点向量的求解 仿佛很 ，只要 求解方程（ 2）的根，就能获得特点 ，再解 次方程 果 （ I-A ） x=0 （3） 的解，便可获得相 的特点向量。 明 上述方法 于 n 很小 是能够的。但当 n 稍大 ， 算工作量将以惊 人的速度增大， 并且因为 算 有 差， 方程（2）未必是精准的特点方程， 自然就不用 求解方程（ 2）与（ 3）的困 了。 法是一种 算矩 主特 征 （矩 按模最大的特点 ）及 特点向量的迭代方法，特 是用于 大型稀少矩 。反 法是 算海森伯格 或三角 的 一个 定近似特 征 的特点向量的有效方法之一。 二．算法设计及流程图 1、 法算法 （1）取初始向量 u ( 0 ) （比如取 u( 0) =(1,1, ?1) T ）, 置精度要求 ，置 k=1. （2） 算 v ( k ) =Au( k 1) ,m k =max(v( k ) ), u (k ) = v ( k) / m k （3）若 | m k = mk 1 | ， 停止 算（ mk 作 最大特点 1 ，u ( k) 作 相 的特点向量）否 置 k=k+1， （ 2） 2、反 法算法 （1）取初始向量 u ( 0 ) （比如取 u( 0) =(1,1, ?1) T ）, 置精度要求 ，置 k=1. （2） A 作 LU分解，即 A=LU （3）解 性方程 Ly (k ) =u( k 1) ,Uv ( k ) =y (k ) （4） 算 m k =max(v(k ) ), u ( k ) = v ( k ) / m k （5）若 |m k =mk 1 | ， 停止 算（ 1/m k 作 最小特点 n ，u ( k) 作 相 的特点向量） ; 否 置 k=k+1， （ 3）. 法流程 ： 开始 输入 A；[m,u,index] =pow(A,1e-6) k=0;m1= v=A*u [vmax,i]=max(abs(v)) m=v(i);u=v/m abs(m-m1) 1e-6 index=1;break; 输出： m，u，index 结束  m1=m;k=k+1 反幂法流程图 开始 输入 A；[m ,u,index] =pow_inv(A,1e-6) k=0;m1=0 v=invA*u m1=m;k=k+1 [vmax,i]=max(abs(v)) m=v(i);u=v/m abs(m-m1) 1e-6 index=1;break; 输出： m，u，index 输入 A； 结束 三、算法的理论依照及其推导 （一）幂法算法的理论依照及推导 幂法是用来确立矩阵的主特点值的一种迭代方法，也即，绝对值最大的特 征 。略微改正 方法，也能够用来确立其余特点 。 法的一个很实用的特 性是它不 能够生成特点 ，并且能够生成相 的特点向量。 上， 法 常用来求通 其余方法确立的特点 的特点向量。 1、 法的迭代格式与收 性 n 矩 A 的特点 1， 2,?, n 是按 大小 号的， x i (i=1,2, ?,n)i 的特点向量，且 1 根，即 | 1|| 2 | ≥?≥ | n | 算最大特点 与特点向量的迭代格式 v (k ) =Au( k 1) ,m k =max(v( k) ), u (k ) = v ( k ) / m k （ 1） 此中 max(v( k) ) 表示向量 v ( k) 的最大重量。 2、 于 法的定理 按式（ 1） 算出 mk 和 u (k ) 足 lim