计算方法列主元素消去法求线性方程组的解

最新推荐文章于 2024-04-20 11:02:39 发布

lf144090

最新推荐文章于 2024-04-20 11:02:39 发布

阅读量2.7k

点赞数

CC 4.0 BY-SA版权

分类专栏：计算方法文章标签： function 算法

本文链接：https://blog.csdn.net/lf144090/article/details/2261453

计算方法专栏收录该内容

14 篇文章

订阅专栏

本文介绍了一种使用列主元消去法解决线性方程组的方法，并提供了详细的MATLAB代码实现过程。该方法适用于求解唯一解的满秩方阵，通过列主元的选择提高数值稳定性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

function x=leizhuyuansu(A,b)
% 用列主元素消去法求线性方程组的解
% A为系数矩阵(要求解惟一，则必须为满秩方阵 n*n),b为右值 n*1
% x为所求的根 n*1

% 作者：野渡无人
% 最后修改日期：2008.3.27

% 如：
% A=[2 4 -2 -2;1 2 4 -3;-3 -3 8 -2;-1 1 6 -3],b=[-4 5 7 7]'
% 结果为：

% x=leizhuyuansu(A,b)
% AA =
%
%     2     4    -2    -2    -4
%     1     2     4    -3     5
%    -3    -3     8    -2     7
%    -1     1     6    -3     7
%
%
% AA =
%
%   -3.0000   -3.0000    8.0000   -2.0000    7.0000
%         0    2.0000    3.3333   -3.3333    0.6667
%         0         0    5.0000   -2.0000    7.0000
%         0         0         0    1.0000    4.0000
% x =
%
%    1.0000
%    2.0000
%    3.0000
%    4.0000

[m,n]=size(A);
if m~=n
error('请输入方阵！');
end
rankA=rank(A);
if rankA~=m
error('输入的矩阵不是满秩矩阵，请重新输入！');
end

AA=[A,b]
n=n+1;
for i=1:m-1
    % 找出列主元素值，以及相应的下标号flag
    max=abs(AA(i,i));
    flag=i;
    for j=i+1:m
        if max < abs(AA(j,i))
            flag=j;
            max=abs(AA(j,i));
        end
    end

    %如果列主元素为0，则报错
    if max==0
        error('本列元素已全0，不能继续！');
    end

     % 交换第i行与第flag行的元素
    if flag~=i
        temp=AA(i,:);
        AA(i,:)=AA(flag,:);
        AA(flag,:)=temp;
    end

    % 以下的算法类同Guass-Jordan 法，注释不再重述
    for k=i+1:m
        AA(k,i:n)=AA(k,i:n)-(AA(k,i)/AA(i,i))*AA(i,i:n);
    end
end
AA

% 求根
x=zeros(m,1);
x(m)=AA(m,n)/AA(m,m);
for i=m-1:-1:1
x(i)=(AA(i,n)-AA(i,i+1:m)*x(i+1:m,1))/AA(i,i);
end

% A=[11 -3 -2;-23 11 1;1 -2 2]
% b=[3 0 -1]'
% x=leizhuyuansu(A,b)
% 结果：
% AA =
%
%    11    -3    -2     3
%   -23    11     1     0
%     1    -2     2    -1
%
%
% AA =
%
% -23.0000   11.0000    1.0000         0
%         0    2.2609   -1.5217    3.0000
%         0         0    1.0192    1.0192
%
%
% x =
%
%    1.0000
%    2.0000
%    1.0000