求解任意方程组的代码

function x=equ_group(A,b)

%该函数用于求解线性方程组

%输入参数,A:方程组的系数矩阵，b:方程组的常数向量（列向量）

%输出参数,x:方程组的解

%%使用实例

% A=[1,1,-3,-1;

% 3,-1,-3,4;

% 1,5,-9,-8];

% b=[1;4;0];

B=[A,b];%增广矩阵

[m,n]=size(B);

num1=0;

for i=1:m

num1=num1+i;

end

C=zeros(1,n);

i=1;

j=1;

n1=n-1;%系数矩阵或增广矩阵的列数

flag1=0;

while j<=n1

if B(i,j)~=0

B(i,:)=B(i,:)/B(i,j);

for k=i+1:m

B(k,:)=B(k,:)-B(k,j)*B(i,:);

end

C(1,j)=i;

if sum(C)~=num1&&j==n1&&flag1==0%系数矩阵在消元过程中，若出现对角线及其一下元素均为0时，将n1变为增广矩阵的列数。

n1=n1+1;%在j等于系数矩阵的列时，n1增加1,变为增广矩阵的列。

flag1=1;%flag1保证if内的语句，只执行1次。

end

i=i+1;

j=j+1;

else

flag=0;

k=i+1;

while k<=m

if B(k,j)~=0

tt=B(i,:);

B(i,:)=B(k,:);

B(k,:)=tt;

flag=1;

if flag==1

break;

end

end

k=k+1;

end

if flag==0

j=j+1;

end

end

end

j=n-1;

while j>=1

i=C(1,j);

if i~=0

k=i-1;

while k>=1

B(k,:)=B(k,:)-B(k,j)*B(i,:);

k=k-1;

end

end

j=j-1;

end

for i1=m:-1:1%i1:增广矩阵最后1列，非0行的行数

if B(i1,n)~=0

break;

end

end

for i2=m:-1:1%i1:系数矩阵最后1列，非0行的行数

if B(i2,n-1)~=0

break;

end

end

if i1>i2

disp('方程无解!');

x=nan;

elseif i2==m

disp('方程有唯一解!');

x=B(:,n);

else

disp('方程有无限多解!');

disp('方程增广矩阵的行最简形为:');

x=B;

end

end