利用Matlab求解线性方程组

2021-10-03 14:32:41 阅读：228 来源： 互联网

利用高斯消元法编写了一个能够计算线性方程组，无解，有唯一解，无穷多解情况的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

标签：方程,系数,end,求解,增广,线性方程组,矩阵,Matlab,n1
来源： https://www.cnblogs.com/siying-zhang/p/15364109.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

利用Matlab求解线性方程组