数学模型到具体算法仍有一段距离,需要考虑计算机执行速度和“四舍五入”等问题,这里主要利用矩阵的初等变换进行矩阵求逆,也就是高斯约当法,代码如下:
主函数: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
结果输出: