矩阵征值计算ppt课件.ppt

《矩阵征值计算ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第八章矩阵特征值计算计算方法——幂法与反幂法1本章内容特征值基本性质幂法与反幂法正交变换与矩阵分解QR方法2本讲内容特征值基本性质幂法幂法的加速反幂法3特征值性质Ax=x(C,x0)性质(1)特征值与特征向量(2)(3)(4)若A对称，则存在正交矩阵Q，使得4圆盘定理定理：(Gerschgorin圆盘定理)设是A的特征值，则i=1,2,...,n设A=(aij)Rnn，记Gerschgorin圆盘若有m的圆盘互相连通，且与其它圆盘都不相连，则这m个圆盘内恰好包含m个特征值。5Ray

2、leigh商定理：设A是n阶实对称矩阵，其特征值为则对任意非零向量x，有且称为矩阵A关于x的Rayleigh商。6(1)任取一个非零向量v0，要求满足(x1,v0)0(2)对k=1,2,...，直到收敛，计算幂法计算矩阵的主特征值（按模最大）及其特征向量假设：(1)

3、1

4、>

5、2

6、…

7、n

8、0(2)对应的n个线性无关特征向量为：x1,x2,...,xn计算过程：幂法（乘幂法，幂迭代）7幂法的收敛性收敛性分析设越小，收敛越快8幂法的收敛性当k充分大时，有又(j=1,2,...,n)vk为

9、1的近似特征向量9幂法的收敛性定理：设A有n个线性无关的特征向量，其特征值满足则由幂法生成的向量满足注：幂法的收敛速度取决于的大小10幂法改进方法：规范化幂法中存在的问题11幂法1的计算12改进的幂法定理：设A有n个线性无关的特征向量，其特征值满足则由改进的幂法生成的向量满足(1)任取一个非零向量v0，要求满足(x1,v0)0(2)对k=1,2,...，直到收敛，计算改进的幂法13举例例：用改进的幂法计算下面矩阵的主特征值和对应的特征向量ex81.m14幂法的加速幂法的收敛速度取决于的大小

10、当r接近于1时，乘幂法收敛会很慢！幂法的加速：原点平移法令B=A–pI，则B的特征值为：i-p选择适当的p满足：(1)(j=2,...,n)(2)用幂法计算矩阵B的主特征值：1-p保持主特征值加快收敛速度带位移的幂法15举例例：用带位移的幂法计算下面矩阵的主特征值和对应的特征向量，取p=0.75ex82.m16反幂法计算矩阵的按模最小的特征值及其特征向量假设：(1)

11、1

12、

13、2

14、…

15、n-1

16、>

17、n

18、>0反幂法(2)对应的n个线性无关特征向量为：x1,x2,...,xnA-1的特征

19、值为：对应的特征向量仍然为x1,x2,...,xn反幂法：对矩阵A-1使用幂法17反幂法定理：设A有n个线性无关的特征向量，其特征值满足则由反幂法生成的向量满足(1)任取一个非零向量v0，要求满足(x1,v0)0(2)对k=1,2,...，直到收敛，计算反幂法18反幂法的加速反幂法的收敛速度取决于的大小当r接近于1时，反乘幂法收敛会很慢！可以使用原点平移法对反幂法进行加速问题：如何选择参数p？离n越近越好（但不能相等）19Rayleigh商加速Rayleigh商加速(1)任取一个非零向量v0

20、，要求满足(x1,v0)0(2)对k=1,2,...，直到收敛，计算20几点注记带位移的反幂法中需要计算带位移的反幂法可以用于计算任何一个特征值k将参数p取为k附近若已知特征值，计算特征向量时，可使用带位移的反幂法令p足够靠近k21作业教材276页，习题3(1)22