matlab第10章 特征值与特征向量的计算.pptVIP

**第10章矩阵特征值与特征向量的计算 10.1幂法及反幂法 10.2Jacobi方法 10.3QR方法 10.4特征值与特征向量的MATLAB函数求解 10.5实例解析本章目标：计算矩阵的特征值及对应的特征向量一、幂法条件：A有特征根|?1||?2|?…?|?n|?0，对应n个线性无关的特征向量………|?i/?1|1当k充分大时，有这是A关于?1的近似特征向量思路：从任意出发，要求10.1幂法及反幂法?规范化为避免大数出现，需将迭代向量规范化，即每一步先保证，再代入下一步迭代。一般用。记：则有：一般地,不妨设：A为实方阵,有特征值|?1|?|?2|?…?|?n|,对应n个线性无关的特征向量当|?1|=|?2|时,需分情况加以讨论:(1)?1=?2;(2)?1=-?2;(3)?1,?2为共轭复数;?算法的一般化——实际计算中的幂法1°?1=-?2;|?1||?3|?…?|?n|从任意出发，不妨假定当k充分大时,有:同号同号所以可以证明，对应于?1的A的特征向量为：事实上，类似地，对应于?2的A的特征向量为：2°|?1|=|?2||?3|?…?|?n|此时,?1和?2有可能是共轭复数(也可能?1=?2,也可能是情况1°?1=-?2);|?1||?3|.不妨假设当k充分大时,有:Q:如何找到表示?1(?2)的较好的关系呢?消元?降次?不难验证:间近似地成立下述线性关系为求得?1和?2,可任取两组分量,并解下列方程组得p,q:其余分量是否也满足关系式?若满足即,?1和?2是方程?2+p?+q=0的两个根:显然:p24q,?1和?2是共轭复根;若p2=4q,则?1=?2;若p=0,则?1=-?2.不难验证:A的对应于?1的特征向量为：A的对应于?2的特征向量为：小结:根据以上讨论,用乘幂法进行计算:当k充分大时,检查是否出现下列三种情况之一:(1)趋近于某一常数;(2)趋近于某一常数;(3)的波动不规律,但相继三向量间满足关系若是,就可求出A的模最大的特征值和相应的特征向量.注:实际计算时,由于其他情况的存在,上述三种情况均不出现时,需考虑其他算法.因此乘幂法在某种意义上来说只能用来试算.?原点平移法决定收敛的速度，特别是|?2/?1|希望|?2/?1|越小越好。不妨设?1?2?…??n，且|?2||?n|。?1?2?nOp=(?2+?n)/2思路令B=A?pI，则有|?I?A|=|?I?(B+pI)|=|(??p)I?B|??A?p=?B。而，所以求B的特征根收敛快。p是假定的,p究竟是多少?p的选择或凭借于经验,或通过多次试算而得.二、反幂法若A有|?1|?|?2|?…|?n|,则A?1有对应同样一组特征向量。11111lll?…?-nnA?1的主特征根A的绝对值最小的特征根Q:Howmustwecomputeineverystep?A:SolvealinearsystemwithAfactorized.若知道某一特征根?i的大致位置p,即对任意j?i有|?i?p||?j?p|,并且如果(A?pI)?1存在，则可以用反幂法求(A?pI)?1的主特征根1/(?i?p)，收敛将非常快。思路(1)用反幂法求A的按模最小的特征根(2)利用反幂法将特征根的近似值精确化若?j’是A的特征值?j的近似值,且设?j是A的