高斯约当法矩阵求逆（Matlab实现）

数学模型到具体算法仍有一段距离，需要考虑计算机执行速度和“四舍五入”等问题，这里主要利用矩阵的初等变换进行矩阵求逆，也就是高斯约当法，代码如下：

主函数：main.m

%测试
A=[1 -2 3; 2 3 4; 3 5 7];
b=[-10;5;7];
%内置函数
x1=A\b
%高斯约当法
x=inverse(A,b)

函数：inverse.m

%求可逆矩阵（AX=b）
%高斯约当法
function x=inverse(A,b)
row=size(A,1);
%C用于存储A的逆矩阵
C=eye(row);
A=[A C];
row_add=size(A,2);
for i=1:row
    %求单列最大元素的行号
    [col_max,nmax]=max(abs(A(i:row,i)));%注意这里i:n
    A(nmax+i-1,:)=A(nmax+i-1,:)/A(nmax+i-1,i);%理解这里nmax+i-1，nmax为A(i:n,i)所在行号，i-1为第i行上面的行数
    %交换最大元素那一行和第i行
    temp=A(i,:);
    A(i,:)=A(nmax+i-1,:);
    A(nmax+i-1,:)=temp;
    for j=1:row
        if(j~=i)
           A(j,:)= A(j,:)-A(j,i)*A(i,:);
        end        
    end  
end
x=A(:,row+1:row_add)*b;%AX=b方程的解
C=A(:,row+1:row_add);%A的逆矩阵
end

结果输出：

 参考视频：陈博士编程序matlab第27期高斯约当法解线性方程组_哔哩哔哩_bilibilihttps://www.bilibili.com/video/BV1nq4y1J7eB?spm_id_from=333.337.search-card.all.click