高斯消元法求解线性方程组 - 详细代码解析 - 常规

高斯消元法求解线性方程组 - 详细代码解析

本文将详细解析使用高斯消元法求解线性方程组的 MATLAB 代码，并解释其原理和步骤。

function x = gauss512(A, b)
    A = [10 -7 0 1;
         -3 2.099999 6 2;
         5 -1 5 -1;
         2 1 0 2];
    b = [8;
         5.900001;
         5;
         1];
    [n,n] = size(A);
    x = zeros(n, 1);
    P = [1:n];
    Aug = [A,b]; % 增广矩阵
    for k = 1:n-1
        [piv,r] = max (abs (Aug(k:n, k))); % 找列主元所在子矩阵的行
        rr = r + k - 1; % 列主元所在大矩阵的行
        if r > k
            temp = Aug(k, :);
            Aug(k, :) = Aug(r, :);
            Aug(r, :) = temp;
        end
        if Aug(k, k) == 0
            error('对角元出现 0')
        end

        % 把增广矩阵消元成为上三角
        for p = k+1:n
            Aug(p, :) = Aug(p, :) - Aug(k, :) * Aug(p, k) / Aug(k, k);
        end
    end

    % 解上三角方程组
    A = Aug(:,1:n); b = Aug(:,n+1);
    x(n) = b(n) / A(n, n);
    for k = n-1:-1:1
        x(k) = b(k);
        for p = n:-1:k+1
            x(k) = x(k) - A(k, p) * x(p);
        end
        x(k) = x(k) / A(k, k);
    end
    P
x
detA = det(A)
end

代码逐句解析

定义函数 gauss512，输入参数为矩阵 A 和列向量 b。

该函数用于求解线性方程组 Ax = b，其中 A 为系数矩阵，b 为常数向量，x 为未知向量。
初始化矩阵 A 和列向量 b。

这里给出了一个具体的线性方程组的系数矩阵 A 和常数向量 b，用户可以根据实际情况修改。
获取矩阵 A 的行数和列数，同时初始化解向量 x 为全零向量。

size(A) 用于获取矩阵 A 的行数和列数，并将其赋值给变量 n。解向量 x 被初始化为全零向量，其维度与 A 的列数相同。
初始化置换向量 P，将其初始化为 $1$ 到 $n`。

置换向量 P 用于记录高斯消元过程中行交换的操作，初始情况下 P 的元素为 $1$ 到 $n$，表示没有进行行交换。
将矩阵 A 和列向量 b 拼接成增广矩阵 Aug。

增广矩阵 Aug 由系数矩阵 A 和常数向量 b 拼接而成，它包含了线性方程组的所有信息。
进行 $n-1$ 次高斯消元操作。

高斯消元法的核心思想是将增广矩阵通过一系列行变换，将系数矩阵化为上三角矩阵，然后通过回代法求解方程组。
- 找到列主元所在子矩阵的行。
  
  [piv,r] = max (abs (Aug(k:n, k))) 用于找到从第 k 行到第 n 行的第 k 列元素中绝对值最大的元素，并将其索引 r 记录下来。
- 交换列主元所在行与第 k 行。
  
  如果 r 不等于 k，则需要将第 r 行与第 k 行交换位置，以保证第 k 行的第 k 列元素不为零。
- 判断对角元是否为零。
  
  如果 Aug(k, k) 为零，则说明该方程组不存在唯一解，程序将抛出异常。
- 利用第 k 行将第 k+1 到第 n 行消元。
  
  对于第 k+1 到第 n 行，使用第 k 行将它们消元，使其第 k 列元素变为零。
解上三角方程组。

当增广矩阵的系数矩阵被化为上三角矩阵后，可以使用回代法求解方程组。
- 将矩阵 A 和列向量 b 从增广矩阵中分离出来。
- 从最后一行开始，依次求解每个未知量 x_i 的值，利用已知的 x_{i+1} 至 x_n 的值进行回代计算。
- 输出置换向量 P、解向量 x 和矩阵 A 的行列式。

高斯消元法的原理

高斯消元法是一种求解线性方程组的经典方法，它利用一系列行变换将系数矩阵化为上三角矩阵，然后通过回代法求解方程组。

高斯消元法的步骤

构建增广矩阵。
将系数矩阵化为上三角矩阵。
- 找到列主元。
- 交换列主元所在行与当前行。
- 利用当前行将其他行消元。
解上三角方程组。
- 从最后一行开始，依次求解每个未知量。
- 利用已知的未知量进行回代计算。

总结

本文详细解析了使用高斯消元法求解线性方程组的 MATLAB 代码，包含代码逐句解释，并解释了高斯消元法的原理和步骤。希望本文能帮助读者更好地理解高斯消元法及其应用。