15章 求矩阵特征值和特征向量.docx

1、15章 求矩阵特征值和特征向量第15章 求矩阵特征值和特征向量 幂 法 幂法规范化算法1. 输入矩阵A、初始向量u(0)，误差eps2. k13. 计算V(k) Au(k-1)4. mk max(V(k), mk-1 max(V(k-1)5. uk V(k)/mk 6. 如果| mk - mk-1|m,m2=xk;m=m1, k,1,Lengthx; m2v=a.u;m0=fmaxu;m1=fmaxv;t=Absm1-m0/N;k=0;Whileteps&k=nmax,Print迭代超限说明：本程序用于求矩阵A按模最大的特征值及其相应特征向量。程序执行后，先通过键盘输入矩阵A、迭代初值向量u(

2、0)、精度控制eps和迭代允许最大次数nmax，程序即可给出每次迭代的次数和对应的迭代特征值、特征向量及误差序列，它们都按10位有效数输出。其中最后输出的结果即为所求的特征值和特征向量序。如果迭代超出nmax次还没有求出满足精度的根则输出迭代超限提示，此时可以根据输出序列判别收敛情况。程序中变量说明a：存放矩阵Au：初始向量u(0)和迭代过程中的向量u（k）及所求特征向量v: 存放迭代过程中的向量V（k）m1: 存放所求特征值和迭代过程中的近似特征值nmax:存放迭代允许的最大次数eps:存放误差精度fmaxx: 给出向量x中绝对值最大的分量k:记录迭代次数t1:临时变量注：迭代最大次数可以修

3、改为其他数字。 例题与实验例1 用幂法求矩阵的按模最大的特征值及其相应特征向量，要求误差10-4。解：执行幂法程序后在输入的4个窗口中按提示分别输入133,6,135，44,5,46，-88,-6,-90，1,1,1，0.0001，20每次输入后用鼠标点击窗口的“OK”按扭，得如下输出结果。 k=1 特征值=44.42335766 误差=229.5766423 特征向量=1., 0.3467153285, -0.6715328467k=2 特征值=44.92343082 误差=0.5000731606 特征向量=1., 0.3341275058, -0.6672691423k=3 特征值=44

4、.99546459 误差=0.07203376236 特征向量=1., 0.3333729572, -0.6667020234k=4 特征值=44.99977337 误差=0.004308781874 特征向量=1., 0.3333351894, -0.6666684279k=5 特征值=44.99998937 误差=0.0002160020115 特征向量=1., 0.3333334179, -0.6666667492k=6 特征值=44.99999952 误差=0.0000101441501 特征向量=1., 0.3333333371, -0.6666666704此结果说明迭代6次，求得误

5、差为err=0.0000101441501的按模最大的特征值=44.99999952 及其对应的一个特征向量=1., 0.3333333371, -0.6666666704。本题矩阵A的3个特征值为45., 2., 1.，可见所求结果很好。但如果执行幂法程序后在输入的4个窗口中按提示分别输入133,6,135，44,5,46，-88,-6,-90，1,1,-1，0.0001，20每次输入后用鼠标点击窗口的“OK”按扭，得如下输出结果。k=1 特征值=2.5 误差=1.5 特征向量=1., 0.75, -1.k=2 特征值=2.2 误差=0.3 特征向量=1., 0.7, -1.k=3 特征值=

6、2.090909091 误差=0.1090909091 特征向量=1., 0.6818181818, -1.k=4 特征值=2.043478261 误差.0474* 特征向量=1., 0.6739130435, -1.k=5 特征值=2.021276596 误差.022* 特征向量=1., 0.670212766, -1.k=6 特征值=2.010526316 误差.010* 特征向量=1., 0.6684210526, -1.k=7 特征值=2.005235602 误差=0.005290713695 特征向量=1., 0.667539267, -1.k=8 特征值=2.002610966 误

7、差=0.002624636037 特征向量=1., 0.6671018277, -1.k=9 特征值=2.001303781 误差=0.001307185093 特征向量=1., 0.6668839635, -1.k=10 特征值=2.000651466 误差=0.0006523151668 特征向量=1., 0.6667752443, -1.k=11 特征值=2.000325627 误差=0.0003258389664 特征向量=1., 0.6667209378, -1.k=12 特征值=2.000162787 误差=0.0001628399197 特征向量=1., 0.6666937978

8、, -1.k=13 特征值=2.000081387 误差=0.00008140008032 特征向量=1., 0.6666802311, -1. 此结果说明迭代13次，求得误差为err=0.00008140008032的按模最大的特征值=2.000081387及其对应的一个特征向量=1., 0.6666802311, -1.。 选用不同的迭代初值获得两个不同结果，显然第二个特征值=2.000081387不是模最大的特征值。 上面实验说明使用幂法依赖于迭代初值的选取且有时得到的结果不是模最大的特征值（知道是什么原因吗？）。不过一般情况下，幂法是可以求出按模最大的特征值的。如果不放心，可以选用两个

9、不同的初值迭代计算，通过计算结果可以马上确定按模最大的特征值。例2 用幂法求矩阵的按模最大的特征值及其相应特征向量，要求误差 10-5。解：执行幂法程序后在输入的4个窗口中按提示分别输入2.,-2.,3., 1,1.,1, 1.,3,-1，1,0,1，0.00001，20每次输入后用鼠标点击窗口的“OK”按扭，得如下输出结果。 .k=9 特征值=2.990381958 误差.0454* 特征向量=1., 0.9810402963, 0.952738931k=10 特征值=2.974823934 误差=0.01555802398 特征向量=0.9684837058, 0.9810717388,

10、1.k=11 特征值=2.99573741 误差.020* 特征向量=1., 0.9915058875, 0.9787802529k=12 特征值=2.988679144 误差=0.007058265982 特征向量=0.985843744, 0.9915041722, 1.k=13 特征值=2.998102573 误差=0.009423429932 特征向量=1., 0.996208617, 0.9905232775k=14 特征值=2.994943938 误差=0.003158635286 特征向量=0.993679344, 0.9962073749, 1.k=15 特征值=2.99915

11、5514 误差=0.004211575789 特征向量=1., 0.9983114144, 0.9957787292k=16 特征值=2.997748176 误差=0.001407337687 特征向量=0.9971851559, 0.9983110678, 1.k=17 特征值=2.999624371 误差=0.001876194866 特征向量=1., 0.9992487853, 0.9981219847k=18 特征值=2.998998284 误差=0.0006260875486 特征向量=0.9987478473, 0.9992487055, 1.k=19 特征值=2.99983298

12、7 误差=0.0008347031267 特征向量=1., 0.9996659782, 0.9991649479k=20 特征值=2.999554616 误差=0.0002783707007 特征向量=0.9994432692, 0.9996659612, 1.迭代超限此结果说明迭代20次后还没有得到满足要求的解，但观察特征值序列发现其是收敛的，因此可以增大迭代次数以求得满足要求的解。本题将最大迭代次数设定为100后得出在迭代第30次时的满足要求的解为 k=30 特征值=2.999992274 误差=4.828768325 10-6 特征向量=0.9999903425, 0.999994205

13、5, 1.注意到本题按模最大的特征值为3，因此求解效果较满意。反 幂 法 反幂法规范化算法1. 输入矩阵A、初始向量u(0)，误差eps2. k13. 解方程AV(k) =u(k-1) 求出解V(k)4. mk max(V(k), mk-1 max(V(k-1)5. uk V(k)/mk 6. 如果| mk - mk-1|m,m2=xk;m=m1, k,1,Lengthx; m2;v=a.u;a1=Inversea;m0=fmaxu;m1=fmaxv;t=Absm1-m0/N;k=0;Whileteps&k=nmax,Print迭代超限说明：本程序用于求矩阵A按模最小的特征值及其相应特征向量。

14、程序执行后，先通过键盘输入矩阵A、迭代初值向量u(0)、精度控制eps和迭代允许最大次数nmax，程序即可给出每次迭代的次数和对应的迭代特征值、特征向量及误差序列，它们都按10位有效数输出。其中最后输出的结果即为所求的特征值和特征向量序。如果迭代超出nmax次还没有求出满足精度的根则输出迭代超限提示，此时可以根据输出序列判别收敛情况。程序中变量说明：a：存放矩阵Au初始向量u(0)和迭代过程中的向量u（k）及所求特征向量v: 存放迭代过程中的向量V（k）a1：存放逆矩阵A-1m1: 存放所求特征值和迭代过程中的近似特征值nmax:存放迭代允许的最大次数eps: 存放误差精度fmaxx: 给出向

15、量x中绝对值最大的分量k: 记录迭代次数t1 :临时变量注：迭代最大次数可以修改为其他数字。 例题与实验例3用反幂法求矩阵的按模最小的特征值及其相应特征向量，要求误差eps&km,m=Absai,j;p=i;q=j, i,1,n,j,i+1,n;mu=ap,p-aq,q;Ifmu=0,thi=Pi/4,thi=ArcTan2*ap,q/mu/2;s=Sinthi/N;c=Sqrt1-s2;a1=bbp;bbp=c*bbp+s*bbq;bbq=-s*a1+c*bbq;pp=ap,p*c*c+aq,q*s*s+2ap,q*s*c;qq=ap,p*s*s+aq,q*c*c-2ap,q*s*c;Doa

16、1=ap,j; ap,j=c*ap,j+s*aq,j; aj,p=ap,j; aq,j=c*aq,j-s*a1; aj,q=aq,j, j,1,n;ap,p=pp;aq,q=qq;ap,q=0;aq,p=0;ea=Sumai,j2,i,1,n,j,1,n-Sumai,i2,i,1,n/N;k=k+1;Print误差=,ea;Print相似矩阵A=;PrintMatrixForma;Print特征向量J;PrintMatrixFormTransposebb ;Ifk=nmax,Print“迭代超限” 说明：本程序用于求对称矩阵A的所有特征值及其相应特征向量。程序执行后，先通过键盘输入矩阵A、矩阵

17、阶数n、精度控制eps和迭代允许最大次数nmax，程序即可给出每次迭代的次数和对应的迭代特征值、特征向量及误差序列。其中最后输出的结果即为所求的特征值和特征向量。如果迭代超出nmax次还没有求出满足精度的根则输出迭代超限提示。此外，输出的特征值矩阵可以不是真正的对角矩阵，但它们的主对角元素就是满足要求的所有特征值。程序中变量说明a：存放矩阵A及其相似变换过程中的Akbb: 存放特征向量矩阵J的转置v: 存放迭代过程中的向量V（k）m1: 存放所求特征值和迭代过程中的近似特征值nmax:存放迭代允许的最大次数eps: 存放误差精度k: 记录迭代次数t1,mu,s,c,ea,p,q,m: 临时变量

18、a1,qq,pp：临时向量 例题与实验例4用Jacobi方法求矩阵的所有特征值及其相应特征向量，要求误差10-6。解：执行Jacobi方法程序后在输入的4个窗口中按提示分别输入2,-1,0，-1,2,-1，0,-1,2，3，0.000001，30每次输入后用鼠标点击窗口的“OK”按扭，得如下输出结。迭代次数k=0误差=2.相似矩阵A=1. 0 -0.7071070 3. -0.707107-0.707107 -0.707107 2特征向量J0.707107 -0.707107 00.707107 0.707107 00 0 1迭代次数k=1误差=1.相似矩阵A=0.633975 -0.3250

19、58 0-0.325058 3. -0.6279630 -0.627963 2.36603特征向量J0.627963 -0.707107 -0.3250580.627963 0.707107 -0.3250580.459701 0 0.888074迭代次数k=2误差=0.211325相似矩阵A=0.633975 -0.276837 -0.170364-0.276837 3.38645 0-0.170364 0 1.97958特征向量J0.627963 -0.431846 -0.6474340.627963 0.772574 0.09376140.459701 -0.465444 0.75633

20、2.迭代次数k=5误差=8.3089110-6相似矩阵A=0.585788 0.00203811 -0.00002415420.00203811 3.41421 0-0.0000241542 0 2.特征向量J0.499652 -0.50036 -0.7070980.707616 0.706597 0.00001206840.499628 -0.50036 0.707115迭代次数k=6误差=1.1668510-9相似矩阵A=0.585786 0 -0.00002415420 3.41421 -1.7405 10-8-0.0000241542 -1.7405 10-8 2.特征向量J0.500

21、012 -0.5 -0.7070980.707107 0.707107 0.00001206840.499988 -0.5 0.707115 此结果说明迭代6次，求得误差为1.1668510-9的矩阵A的所有特征值1=0.585786，2=3.41421，3=2 及对应的特征向量0.500012，0.707107，0.499988，-0.5，0.707107，-0.5，-0.707098，0.0000120684，0.707115本题矩阵A的3个特征值为2-sqrt2,2+sqrt2, 2。线性代数告诉我们，不同特征值对应的特征向量是正交的，用如上求出的特征向量验证，可以得到0.500012，0.707107，0.499988.-0.5，0.707107，-0.5= -5.5511210-17，0.500012，0.707107，0.499988.-0.707098，0.0000120684，0.707115=-5.5511210-17-0.5，0.707107，-0.5