幂法计算特征值--亮

clear all

clc

format long

a=[3 -4 3;-4 6 3;3 3 1]; v0=[1;1;1];

i=1;

while i<100

vi=a*v0;

t=max(abs(vi));

u=vi/t;

v0=u;

i=i+1;

end

disp(u)

disp(t)

计算结果：

U=

-0.604362874184780

1.000000000000000

0.150816554524373

t=

8.869901160312239