数值方法线性方程组的高斯消元法

最新推荐文章于 2021-04-09 01:44:11 发布

perfyJoe

最新推荐文章于 2021-04-09 01:44:11 发布

阅读量1k

点赞数

文章标签： matlab

本文链接：https://blog.csdn.net/perfyJoe/article/details/11678525

版权

该博客介绍了如何使用高斯消元法解决线性方程组的问题，提供了一个MATLAB函数`uptrbk`，用于处理非奇异矩阵A和矩阵B，输出解矩阵X。在函数中，进行了部分列交换以优化过程，并通过消除过程逐步求解。如果矩阵A为奇异矩阵，则提示无唯一解。

摘要由CSDN通过智能技术生成

function X=uptrbk(A,B)
%Input     -A is an N x N nonsingular matrix
%          -B is an N x 1 matrix
%Output    -X is an N x 1 matrix containing the solution to AX=B
%Initialize X and the temporary storage matrix C
[N N]=size(A);
X = zeros(N,1);
C = zeros(1,N+1);
%Form the augmented matrix : Aug=[A|B]
Aug=[A B];
for p=1:N-1
    % Partial povoting for column p
    [Y,j]=max(abs(Aug(p:N,p)));
    %Interchange row p and j
    C=Aug(p,:);
    Aug(p,:)=Aug(j+p-1,:);
    Aug(j+p-1,:)=Aug(j+p-1,:);
    Aug(j+p-1,:)=C;
if Aug(p,p)==0
     'A was singular.   No unique solution';
end
%Elimination process for column p
for k=p+1:N
     m=Aug(k,p)/Aug(p,p);
     Aug(k,p:N+1)=Aug(k,p:N+1)-m*Aug(p,p:N+1);
end
end