增广矩阵求逆矩阵MATLAB实现

月的高

已于 2023-10-20 10:46:05 修改

阅读量996

点赞数 1

文章标签：矩阵 matlab 线性代数

于 2022-10-22 13:17:05 首次发布

本文链接：https://blog.csdn.net/qq_43699298/article/details/127460728

版权

矩阵逆矩阵高斯消元法行变换主子式

关键词由CSDN通过智能技术生成

%输入矩阵A，矩阵阶数n，输出矩阵W为A的逆矩阵，注意，该代码需要保证矩阵A的各阶主子式都可逆

function W = matrixinv(A,n)
    
    WW=zeros(n,2*n);
    W=zeros(n,n);
    WW(:,1:n)=A;
    WW(:,n+1:2*n)=eye(n)
    for i=1:n-1
        for j=1:n-i
           WW(i+j,:)=WW(i,:)*(-WW(i+j,i)/WW(i,i))+WW(i+j,:); 
        end
    end
    for i=1:n-1
        for j=1:n-i
           WW(n+1-i-j,:)=WW(n+1-i,:)*(-WW(n+1-i-j,n+1-i)/WW(n+1-i,n+1-i))+WW(n+1-i-j,:); 
        end
    end
    for i=1:n
        WW(i,:)=WW(i,:)/WW(i,i)
    end
    W=WW(:,n+1:2*n);
 end