平移反幂法

平移反幂法 求解距离p最近的特征值λ与单位特征向量x

clear
%% 程序说明
% 平移反幂法求解距p最近的特征值及其对应的特征向量

%% 输入端
A=input('请输入目标矩阵A: \n');
p=input('请输入平移单位p: \n');
n=length(A);
u=ones(n,1);
pI=p*eye(n);
k=0;
% 精度选择
disp('选择默认精度，请输入0；自定义请输入1');
chr=input("");
if (chr ~= 0)
err=input('请输入指定精度:');
else
    err=1e-2;
end
%迭代次数选择
disp('选择默认最大迭代次数，请输入0;自定义请输入1');
chr=input("");
if (chr ~= 0)
N=input('请指定最大迭代次数:');
else
    N=1e4;
end


%% 循坏
while (1)
    v=(A-pI)\u;
    w=v/max(v);
    if norm((u-v),2)<err
        break;
    end
    u=w;
    k=k+1;
    if k>N
        break;
    end
end

%% 输出端
fprintf('\n');
fprintf('\n 迭代步数k: %d ',k);
fprintf('\n 与P距离最近的特征值为：');
disp(p+1/max(v));
fprintf('\n 对应的特征向量为：\n');
disp(u);
fprintf('\n 转化为单位向量为：')
disp(u/norm((u),2));
     

迭代次数显示上还有点问题，正在改进…