计算方法大作业.docxVIP

B1334091 1133409004 孙容瑾 XXXX XXXX XXXX PAGE 16 计算方法大作业 班级 XXXXXX 学号 XXXXXXX 任课老师 贺力平 姓名 XXX 2013年 12月 PAGE 1 实验一 幂法与矩阵特征值 幂法求主特征值思路 幂法的主要思想就是对假设的任意初始列向量作用n次A矩阵（左乘A矩阵）后，初始向量就接近A矩阵的主特征值对应的特征向量。由于左乘n次A矩阵有可能会造成计算量溢出，所以每次都对列向量作归一化处理。 开始 开始 输入矩阵 输入矩阵A,v(0) v(1)=Av(0)，i=0 v(1)=Av(0)，i=0 u(i)=v(i)/max(v(i)) u(i)=v(i)/max(v(i)) i=i+1 i=i+1 v(i+1)=Au(i) v(i+1)=Au(i) max(v(i+1))-max(v(i)) max(v(i+1))-max(v(i))ε? no no yes yes u=v(i)/max(v(i)) u=v(i)/max(v(i)) 结束输出矩阵u,λλ=u 结束 输出矩阵u,λ λ=uTAu/uTu 程序代码 function [ ] = mifa( A,v ) %UNTITLED Summary of this function goes here % Detailed explanation goes here A =[ 1 21 2 3 34 5 2 1 54 2 6 4 2 2 59 0]; v=[1 1 1 1]; u(:,1)=v(:,1); fori=1:100 v(:,i+1)=A*u(:,i); if abs(max(v(:,i+1))-max(v(:,i)))10^-4 break end u(:,i+1)=v(:,i+1)/max(v(:,i+1)); end disp(u(:,i)); disp(max(v(:,i))); k=u(:,i)*A*u(:,i)/(u(:,i)*u(:,i)); disp(k); [x,c]=eig(A); disp(c); disp(x); disp(i); end 结果比较和结论 初始矩阵 A =[ 1 21 2 3 34 5 2 1 54 2 6 4 2 2 59 0]; 初始向量 v=[1 1 1 1]; 幂法求得特征向量 12.5147 14.9140 25.6935 39.6330 归一化后特征向量 0.3158 0.3763 0.6483 1.0000 列向量最大值近似主特征值 39.6330 Rayleigh商求出主特征值 39.6330 用eig()函数算出的特征值和特征向量 39.6331 0 0 0 0 -18.2401 + 7.4985i 0 0 0 0 -18.2401 - 7.4985i 0 0 0 0 8.8471 -0.2450 -0.0471 + 0.0965i -0.0471 - 0.0965i -0.0574 -0.2919 0.1147 - 0.0939i 0.1147 + 0.0939i -0.1747 -0.5029 0.2862 - 0.1187i 0.2862 + 0.1187i 0.1535 -0.7758 -0.9330 -0.9330 0.9709 达到精度要求所需次数：17 结论： 可以看出初始列向量经过多次迭代后，用幂法求出的特征值和用eig()函数求出的A的特征值，满足计算精度在，并且特征向量也具有数乘关系。另外发现取