ranchos.m

来自「实现兰乔斯迭代」· M 代码 · 共 17 行

%200612009佟跃
%兰乔斯迭代
function Q=ranchos(A)
m=length(A);b=rand(m,1);t=zeros(m-1,1);d=zeros(m-1,1);
Q=zeros(m);Q(:,1)=b/norm(b);
Q(:,2)=A*Q(:,1)-Q(:,1)'*A*Q(:,1)*Q(:,1)/norm(A*Q(:,1)-Q(:,1)'*A*Q(:,1)*Q(:,1));
for n=2:m-1
    v=zeros(m,1);
    v=A*Q(:,n);
    d(n)=Q(:,n)'*v;
    v=v-t(n-1)*Q(:,n-1)-d(n)*Q(:,n);
    t(n)=norm(v);
    Q(:,n+1)=v/t(n);
end