《北航数值分析第一次大作业(幂法反幂法)》由会员分享，可在线阅读，更多相关《北航数值分析第一次大作业(幂法反幂法)（18页珍藏版）》请在金锄头文库上搜索。
1、一、问题分析及算法描述1.问题的提出：（1）用幂法、反幂法求矩阵A=aij2020的按摸最大和最小特征值，并求出相应的特征向量。其中aij=sin0.5i+0.2j ij1.5cosi+1.2j i=j要求：迭代精度达到10-12。（2）用带双步位移的QR法求上述的全部特征值，并求出每一个实特征值相应的特征向量。2.算法的描述：(1)幂法幂法主要用于计算矩阵的按摸为最大的特征值和相应的特征向量。其迭代格式为：任取非零向量u0=h1(0)，hn(0)Thr(k-1)=max1jnhrk-1 yk-1=uk-1hrk-1 uk=Ayk-1=h1k，hnkT k=sgnhrk-1hrk (k=1，2
2、，)终止迭代的控制选用。幂法的使用条件为nn实矩阵A具有n个线性无关的特征向量x1，x2，xn，其相应的特征值1，2，n满足不等式123n或1=2=m1m+1m+2n幂法收敛速度与比值21或m+11有关，比值越小，收敛速度越快。(2)反幂法反幂法用于计算nn实矩阵A按摸最小的特征值，其迭代格式为：任取非零向量u0Rn k-1=uk-1Tuk-1 yk-1=uk-1k-1 Auk=yk-1 k=yk-1uk k=1，2， 每迭代一次都要求解一次线性方程组Auk=yk-1。当k足够大时，n1k，yk-1可近似的作为矩阵A的属于n的特征向量。比值nn-1越小，收敛的越快。反幂法要求矩阵A非奇异。(3
3、)带双步位移的QR分解法QR方法适用于计算一般实矩阵的全部特征值，尤其适用于计算中小型实矩阵的全部特征值。本算例中采用带双步位移的QR方法，可加速收敛，其迭代格式为：A1=A(n-1) Ak-s1(k)I=QkRk 对Ak-s1kI作QR分解 Ak+1=RkQk+s1kI Ak+1-s2kI=Qk+1Rk+1 (对Ak+1-s2kI作QR分解)Ak+2=Rk+1Qk+1+s2kI (k=1，3，5，) 二、计算结果及分析1.计算结果：(1)幂法：初始条件：最大迭代次数L=1000；向量u0=1，0，0计算结果：第1次迭代结果：最大特征值：0.00000e+000第2次迭代结果：最大特征值：2.
4、48910e+000 相对误差：1.00000e+000第3次迭代结果：最大特征值：1.67719e+000 相对误差：4.84085e-001第4次迭代结果：最大特征值：-2.10960e+000 相对误差：1.79503e+000第5次迭代结果：最大特征值：-6.13203e-001 相对误差：2.44030e+000第794次迭代结果：最大特征值：-1.97638e+000 相对误差：4.30074e-011第795次迭代结果：最大特征值：-1.97638e+000 相对误差：3.04354e-013*最终迭代结果*特征值：-1.97638e+000 相对误差：3.04354e-013迭
5、代次数：795(2)反幂法：初始条件：最大迭代次数L=1000；向量u0=1，0，0运行结果：第1次迭代结果：最大特征值：1.07542e+000第2次迭代结果：最大特征值：-3.66550e+000相对误差：1.29339e+000第3次迭代结果：最大特征值：1.22709e+001相对误差：1.29871e+000第4次迭代结果：最大特征值：-1.03421e+000相对误差：1.28650e+001第5次迭代结果：最大特征值：-5.46339e-001相对误差：8.92983e-001第995次迭代结果：最大特征值：-3.24922e-001相对误差：6.61387e-003第996次迭
6、代结果：最大特征值：-3.27176e-001相对误差：6.88964e-003第997次迭代结果：最大特征值：-3.29570e-001相对误差：7.26527e-003第998次迭代结果：最大特征值：-3.32147e-001相对误差：7.75893e-003第999次迭代结果：最大特征值：-3.34960e-001相对误差：8.39573e-003第1000次迭代结果：最大特征值：-3.38073e-001相对误差：9.20985e-003*超过最大设定迭代次数，迭代失败！(3)带双步位移的QR法：初始条件：最大迭代次数L=1000；向量u0=1，0，0运行结果：全部特征值：0.7477
7、38+i*0.000000-0.316674+i*0.025741-0.316674+i*-0.0257410.488580+i*0.1393760.488580+i*-0.1393761.045898+i*0.000000-0.630981+i*0.0000000.071427+i*0.5398270.071427+i*-0.5398271.265389+i*0.000000-1.459878+i*0.000000-1.521321+i*0.000000-1.412499+i*0.148349-1.412499+i*-0.148349-1.976376+i*0.0000001.810854+
8、i*0.0000001.362562+i*1.0661171.362562+i*-1.0661170.169262+i*1.9094240.169262+i*-1.909424特征向量（经谱范数归一化）：实特征值0.747738对应特征向量：-0.062705 -0.022368 0.304372 0.064466 0.521833 -0.157024 0.136942 -0.218108 0.250264 -0.043064 -0.228688 -0.184632 -0.072871 0.124721 0.029070 0.102566 -0.136358 0.167727 0.085747
9、 0.546165 实特征值1.045898对应特征向量：-0.018001 0.019652 0.273447 0.070528 0.274896 -0.144015 0.048385 0.376439 -0.583051 -0.054008 -0.168682 -0.113430 -0.034709 0.009204 0.472291 0.125664 -0.190617 0.113145 0.046278 0.059871 实特征值-0.630981对应特征向量：0.106861 0.087709 -0.024967 -0.020897 0.064302 0.034047 0.53514
10、3 0.046383 0.028832 0.003479 -0.097276 -0.383801 0.089445 -0.039560 -0.036928 -0.021330 0.014811 0.705836 -0.108904 0.082022 实特征值1.265389对应特征向量：-0.055201 0.003399 0.242191 0.102847 0.372470 -0.372826 0.113953 0.240659 -0.310401 -0.076590 -0.244632 -0.192549 -0.077259 0.263328 0.201662 0.154166 -0.40
11、7814 0.186782 0.094649 0.173302 实特征值-1.459878对应特征向量：0.427828 -0.546801 0.007822 -0.382580 0.025199 0.012788 0.033241 0.005389 -0.004065 0.043524 -0.032112 -0.044233 0.135395 -0.006564 0.001214 0.020165 0.011678 0.050001 -0.585765 0.013115 实特征值-1.521321对应特征向量：0.236032 -0.139250 -0.008143 0.638527 -0.
12、009049 -0.002911 -0.001307 0.003054 0.006515 -0.030134 0.012712 0.011368 -0.018792 -0.001753 -0.005749 -0.014290 -0.005292 -0.014591 0.717590 0.001369 实特征值-1.976376对应特征向量：-0.227404 -0.048154 0.022615 0.297305 0.070372 0.039927 0.078503 0.015822 -0.012182 0.605334 -0.083616 -0.106270 -0.573963 -0.019
13、907 0.003839 0.051362 0.036567 0.115613 0.332707 0.036954 实特征值1.810854对应特征向量：-0.027768 -0.051081 -0.159642 -0.054573 -0.084441 0.118378 0.029553 0.211088 0.203867 0.048627 0.075470 0.026824 -0.011103 -0.584846 -0.196009 -0.097668 0.673159 -0.036833 0.003394 0.0812442.结果分析以上三种方法中，幂法计算共进行了795次迭代才达到收敛，计算量较大，收敛性不好；反幂法计算结果未能收敛，通过进一步分析发现，这是因为反幂法迭代程序未考虑按模最小特征值为复数的情况，造成迭代失败。按双步位移的
