求矩阵特征值和特征向量课件1

《求矩阵特征值和特征向量课件1》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、求矩阵特征值和特征向量的方法本章探讨求矩阵特征值及特征向量的常用数值方法的构造和原理，主要介绍在计算机上常用的求矩阵特征值和特征向量的的常用方法和有关知识。重点论述幂法的构造内容。6.1实际案例旅游地选择问题通过层次分析法可以转化为求成对比较矩阵的绝对值最大的特征值及其对应的特征向量的问题。求矩阵的特征值及特征向量的问题在实际的科研和工程问题中经常遇到，在这些问题中解出矩阵（特别是高阶矩阵）特征值或特征向量成为解决问题的关键。求矩阵的特征值及特征向量的计算机解法也称为代数特征问题的计算方法。6.2问题的描述与基本概念定义6

2、.1设矩阵，称关于变量的行列式函数为矩阵A的特征多项式，称方程为特征方程。定义6.2若存在某个实数或复数及非零向量满足，则称是矩阵A的一个特征值，而称为对应的一个特征向量。是关于的n次多项式，矩阵A的特征值就是的零点。在线性代数中，有求解矩阵A的特征值和特征向量的解法，该解法理论很严密，但由于将特征多项式化为一个n次多项式很复杂且特征方程对舍入误差很敏感，特别当n较大时，这些问题更突出。由于这些原因，实用中在求解代数特征值问题时一般不用如上的线性代数的方法，而采用本章介绍的迭代加变换的计算机求解方法，这些方法具有编程简单，

3、对舍入误差不敏感等优点。6.3幂法幂法---把最大特征值直接从矩阵乘出来！幂法是求矩阵按模最大的特征值及其相应特征向量的方法。基本思想利用矩阵的特征值与特征向量的关系构造迭代向量序列来求矩阵按模最大的特征值及其相应特征向量。1、构造原理设方阵,是A的n个线性无关的特征向量，其对应的特征值为，任取一个非零向量,则有用A左乘，并利用有记，可得假设，，因为,有令V(k)的第i个分量为，考虑分量比，有当k充分大时，有，是对应的一个近似特征向量。用如上求矩阵按模最大的近似特征值及其相应特征向量的方法称为幂法。2.分析当时，导致的计算

4、出现上益错误。定理设方阵，是A的n个线性无关的特征向量，是对应的特征值，任取一个非零向量,按构造规范化向量序列，其中表示的绝对值最大的分量，则有证明由式（6.3）有一般的有记，由，有，再由有利用定理可以写出规范化幂法算法1.输入矩阵A、初始向量，误差eps，实用中一般取；2.kÜ13.计算V(k)ÜAu(k-1)4.mkÜmax(V(k)),mk-1Ümax(V(k-1))5.u(k)ÜV(k)/mk6.如果

5、mk-mk-1

6、

7、满足怎样用幂法求按模最小的特征值及相应特征向量？设特征值为对应的特征向量为，有因为A的n个特征值都不为零，故A可逆，有这说明是A-1特征值，x(k)是对应的特征向量。由，有于是，求A按模最小的特征值相当于求A-1按模最大的特征值，此时，只要将幂法中的A换为A-1即可。用幂法求出A-1按模最大的特征值后取其倒数就得到A按模最小的特征值，相应特征向量不变。用如上方法求矩阵按模最小的特征值及其相应特征向量的称为反幂法。由于求逆是很费时的，在反幂法迭代公式V(k)=A-1u(k-1)常用解线性方程组AV(k)=u(k-1)的方法求

8、得V(k)。数值实验案例编写幂法的通用程序，并用该程序求矩阵按模最大的特征值及其特征向量，要求误差<10-4。观察选择不同初值计算的结果。幂法规范化算法1.输入矩阵A、初始向量u(0)，误差eps2.kÜ13.计算V(k)ÜAu(k-1)4.mkÜmax(V(k)),mk-1Ümax(V(k-1))5.ukÜV(k)/mk6.如果

9、mk-mk-1

10、

11、迭代向量u(0)="];n=Length[u];eps=Input["误差精度eps="];nmax=Input[“迭代允许最大次数nmax=”];fmax[x_]:=Module[{m=0,m1,m2},Do[m1=Abs[x[[k]]];If[m1>m,m2=x[[k]];m=m1],{k,1,Length[x]}];m2]v=a.u;m0=fmax[u];m1=fmax[v];t=Abs[m1-m0]//N;k=0;While[t>eps&&k

12、x[v];t=Abs[m1-m0]//N；Print["k=",k,"特征值=",N[m1,10],"误差=",N[t,10]];Print["特征向量=",N[u,10]]];If[k>=nmax,Print["迭代超限"]]说明：本程序用于求矩阵A按模最大的特征值及其相应特征向量。程序执行后，先通过