数值分析报告.docxVIP

数值分析报告 名称 学号姓名 任课老师数值分析航空航天学院1124139027孙诚骁贺 力平 2012年12月院 部分一 幕法求主特征值思路 幕法的主要思想就是对假设的任意初始列向量作用n 次A矩阵后，初始向量就接近A矩阵的主特征值对应的特 征向量。由于左乘n次A矩阵有可能会造成计算量溢出， 所以每次都对列向量作归一化处理 1 程序代码 function[]=mifa(A, v) %UNTITLEDSummaryoft hisfunctio ngoeshere %Detailed explanatio ngoeshere A=[12123 345 21 54264 225 90]; v=[lll 1]; u(:, l)=v(:, 1); fo ri=l:100 v(:, i+1)二A*u (:, i); if abs(max(v(:, i+1))-ma x(v(:, i)))break end u(:, i + l)=v (:, i + l)/ma x(v(:, i + 1)); end disp(u(:, i)); di sp(max (v (:, i))); k=u (:, i) *A*u (:, i) / (u (:, i) *u (:, i)); disp(k); [x, c]=ei g(A)； disp(c); disp(x); disp(i); e nd 结果比较和结