matlab使用高斯消去法呵列主元高斯消去法解方程

MATLAB可以使用高斯消去法和列主元高斯消去法来解决线性方程组Ax=b的问题。其中，高斯消去法是一种基本的线性代数算法，它通过消元的方式将系数矩阵A转化为一个上三角矩阵，然后通过回代的方式求解出未知数向量x。而列主元高斯消去法则是在高斯消去法的基础上，每次选取主元时选择当前列中绝对值最大的元素作为主元，以提高数值稳定性。

下面是使用MATLAB进行高斯消去法和列主元高斯消去法求解线性方程组的示例代码：

% 高斯消去法求解线性方程组Ax=b
function x = gauss_elimination(A, b)
    n = size(A, 1);
    for k = 1:n-1
        for i = k+1:n
            factor = A(i,k) / A(k,k);
            A(i,k:n) = A(i,k:n) - factor * A(k,k:n);
            b(i) = b(i) - factor * b(k);
        end
    end
    x = zeros(n, 1);
    x(n) = b(n) / A(n,n);
    for i = n-1:-1:1
        x(i) = (b(i) - A(i,i+1:n)*x(i+1:n)) / A(i,i);
    end
end

% 列主元高斯消去法求解线性方程组Ax=b
function x = gauss_elimination_partial_pivot(A, b)
    n = size(A, 1);
    for k = 1:n-1
        [~, pivot] = max(abs(A(k:n,k)));
        pivot = pivot + k - 1;
        if pivot ~= k
            A([k,pivot],k:n) = A([pivot,k],k:n);
            b([k,pivot]) = b([pivot,k]);
        end
        for i = k+1:n
            factor = A(i,k) / A(k,k);
            A(i,k:n) = A(i,k:n) - factor * A(k,k:n);
            b(i) = b(i) - factor * b(k);
        end
    end
    x = zeros(n, 1);
    x(n) = b(n) / A(n,n);
    for i = n-1:-1:1
        x(i) = (b(i) - A(i,i+1:n)*x(i+1:n)) / A(i,i);
    end
end

相关问题：