利用Matlab求解线性方程组

 

利用高斯消元法编写了一个能够计算线性方程组，无解，有唯一解，无穷多解情况的matlab代码。

程序说明：变量n1表示系数矩阵或者增广矩阵的列数。当增广矩阵的秩与系数矩阵的秩相等时（方程有唯一解时），n1表示系数矩阵的列数。当方程组无解或者有无数多解时，n1表示增广矩阵的列数。

处理办法为：

if sum(C)~=num1&&j==n1&&flag1==0%系数矩阵在消元过程中，若出现对角线及其一下元素均为0时，将n1变为增广矩阵的列数。
            n1=n1+1;%在j等于系数矩阵的列时，n1增加1,变为增广矩阵的列。
            flag1=1;%flag1保证if内的语句，只执行1次。
        end
当j执行到系数矩阵的列n1，且sum(C)~=num1（即系数消元过程中，出现了对角线及其一下元素均为0，如图1所示）时，将n1+1.

 

图1

 

function x=liner_equ_v2(A,b)
%该函数用于求解线性方程组
%输入参数,A:方程组的系数矩阵，b:方程组的常数向量（列向量）
%输出参数,x:方程组的解
%时间，2021.10.3
%版权所有人，zsy

%%使用实例
% 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

　　

计算实例：

1、

A=[1,-2,2,-1;
2,-4,8,0;
-2,4,-2,3;
3,-6,0,-6];
b=[1;2;3;4];

x=liner_equ_v2(A,b)

方程无解!

x =

NaN

 

2、

A=[2,1,-3;
1,2,-2;
-1,3,2];
b=[1;2;-2];

x=liner_equ_v2(A,b)

方程有唯一解!

x =

-4.0000
0
-3.0000

 

3、

 

A=[1,1,-3,-1;
3,-1,-3,4;
1,5,-9,-8];
b=[1;4;0];

x=liner_equ_v2(A,b)

 

方程有无限多解!
方程增广矩阵的行最简形为:

x =

1.0000 0 -1.5000 0.7500 1.2500
0 1.0000 -1.5000 -1.7500 -0.2500
0 0 0 0 0