数值分析41(免费阅读).ppt

反幂法 反幂法计算的主要步骤： 1．对A进行三角分解A=LU； 2．求整数r，使得 ，计算 3．解方程组 反幂法 例 用反幂法求矩阵 的按模最小的特征值和相应的特征向量。 取 ，要求误差不超过10-3. 解 矩阵A作LU分解，得 作归一化计算 反幂法 k y(k)（归一化向量） max(x(k)) 0 1.0000 1.0000 1.0000 1.0000 1 0.4348 1.0000 -0.4783 0.5652 2 0.1902 1.0000 -0.8834 0.9877 3 0.1843 1.0000 -0.9124 0.8245 4 0.1831 1.0000 -0.9129 0.8134 5 0.1832 1.0000 -0.9130 0.8134 按模最小特征值 小结 计算矩阵主特征值和主特征向量的幂法、反幂法 习题： 讲义第四章 1（选作一个小题） 课程设计三 构建具有10~20个节点的微型网络，编程，用幂法计算其PageRank. 1）至少设置两种不同的概率p，（如p=0.85,p=0.75），比较不同的概率p对页面的PageRank的影响； 2）分析影响某个网页的PageRank的因素，并在指向该网页的链接总数不变的前提下，找到提高其PageRank的方法。 3）还可进一步自行增加其它研究内容… 参考《电子科技大学学报》投稿格式 1. 题目和摘要（中英文） 2. 正文： a 引言; b 原理； c 数值算例； d 结论。 3. 附件：程序 第九周周三上课时交打印版和电子版（发邮件：xsyang@126.com。 每章约一个设计题目，每人总共做2个设计题目。 格式要求 幂法 A=[2,-1,0;0,2,-1;0,-1,2]; x=[0;0;1]; KM=15; eps=1e-3; k=0; r0=0; r1=norm(x,inf) y=x; while kKM abs(r1-r0)eps k=k+1; r0=r1; x=A*y; %计算下一个向量 r1=norm(x,inf); %向量无穷范数 y=x/r1; end 幂法 k y(k)（归一化向量） max(x(k)) 0 0.000000 0.000000 1.000000 1.000000 1 0.000000 -0.500000 1.000000 2.000000 2 0.200000 -0.800000 1.000000 2.500000 3 0.428571 -0.928571 1.000000 2.800000 4 0.609756 -0.975610 1.000000 2.928571 5 0.737705 -0.991803 1.000000 2.975610 6 0.824658 -0.997260 1.000000 2.991803 7 0.882998 -0.999086 1.000000 2.997260 8 0.921975 -0.999695 1.000000 2.999086 9 0.947978 -0.999898 1.000000 2.999695 幂法 2. 按模最大特征值λ1是r重根，即 矩阵A仍有n个线性无关的特征向量。 只要α1, …, αr不全为零，k充分大时，有 最大特征值 特征向量 注意：一般来说，用r个不同的初始向量，可以得到r个不同的特征向量。 幂法 3. ，矩阵A有n个线性无关的特征向量。 {x(k+1)}是个摆动序列，当k充分大时 幂法 又由 两个特征向量 和 {x(k+1)/x(k) }交错地收敛于两极限，利用这一特性，可以判断矩阵A是否有两个实的反号的主特征值。 幂法 例 求矩阵A的主特征值及主特征向量 解 x(k)的部分结果列于下表 k 0 0.2 0.4 0.6 23 24 25 幂法 x(24)和x(23)各分量之比0.8683, 1.006, 1.179，无规律，而x(25)和x(23)各分量之比是0.2509，0.2509，0.2514。非常接近。 取 小结:（1）对于给定的矩阵，无法事先判断其特征值属于何种情况。通常，采用先算下去，再观察特征值的分布规律，随时判断，分别对待；