jacobi.m

来自「数值线性代数的Matlab应用程序包 共13个程序函数」· M 代码 · 共 14 行

function A=Jacobi(A)
%标准行方式次序的雅可比算法
[m,n]=size(A);
for i=1:m
    for j=i+1:n
        p=atan(2*A(i,j)/(A(i+1,j)-A(i,j-1)))/2;
        temp=cos(p)*A(i,:)-sin(p)*A(j,:);
        A(j,:)=sin(p)*A(i,:)+cos(p)*A(j,:);
        A(i,:)=temp;
        temp=cos(p)*A(:,i)-sin(p)*A(:,j);
        A(:,j)=sin(p)*A(:,i)+cos(p)*A(:,j);
        A(:,i)=temp;
    end
end