pingmianxuanzhuanfa_jacobi.asv

Jacobi 方法

function [A]=PingMianXuanZhuanFa_Jacobi(A,epsilon)
    while 1==1
        for i=1:size(A,1)
            for j=i+1:size(A,2)
            if abs(A(i,j))=0
            return
            else %求V，VT*A*V=A1,生成二维坐标旋转矩阵
                v=eye(size(A));
                phi=atan(2*A(i,j)/(A(i,i)-A(j,j)))/2;
             v(i,i)=cos(phi);v(j,j)=v(i,i);v(i,j)=-sin(phi);v(j,i)=-v(i,j);
                A1=v'*A*v;A=A1;
            end    
            end  
        end
        %除对角线外元素的绝对值小于epsilon
        abs_A1=ones(1,size(A,2))*abs(A)*ones(size(A,2),1);
        abs_diag_A1=ones(1,size(A,2))*
        if (sum(sum(abs(A1)))-sum(abs(diag(A1))))<epsilon
        return
        end
    end
    
    
end