a=[3, -2 ,-4;-2 ,6 ,-2;-4, -2, 3];
v=[1,1,1];
for i=1:100
tma=max(abs(v));
v=v/tma;
v=v*a;
tem=max(abs(v));
if(abs(tma-tem)<=10^(-4))
break;
end
if(i==100)
disp('the i is too small');
end
end
tem=max(abs(v));
v=v/tem;
'迭代次数',disp(i);
'特征值',disp(tem);
'特征向量',disp(v);
% while (abs(max(abs(vv))-max(abs(v)))<10)
% vv=v*a;
% tem=max(abs(vv));
% vv=v/tem;
% end
评论0