计算物理学（数值分析）上机实验答案7、解线性方程组的直接法

实验七、解线性方程组的直接法

​ 解线性方程组的直接法是指经过有限步运算后能求得方程组精确解的方法。但由于实际计算中舍入误差是客观存在的，因而使用这类方法也只能得到近似解。目前较实用的直接法是古老的高斯消去法的变形，即主元素消去法及矩阵的三角 分解法。引进选主元的技巧是为了控制计算过程中舍入误差的增长，减少舍入误 差的影响。一般说来，列主元消去法及列主元三角分解法是数值稳定的算法，它 具有精确度较高、计算量不大和算法组织容易等优点，是目前计算机上解中、小 型稠密矩阵方程组可靠而有效的常用方法。

一、实验目的

1、了解求线性方程组的直接法的有关理论和方法；

2、会编制列主元消去法、LU 分解法的程序；

3、通过实际计算，进一步了解各种方法的优缺点，选择合适的数值方法。

二、算法实例

1.用高斯消元法解线性方程组 AX = b 的 MATLAB 程序

function [RA,RB,n,X]=gaus(A,b)
B=[A b]; n=length(b); RA=rank(A);
RB=rank(B);zhica=RB-RA;
if zhica>0,
 disp('请注意：因为RA~=RB，所以此方程组无解.')
    return
   end
if RA==RB
 if RA==n
        disp('请注意：因为RA=RB=n，所以此方程组有唯一解.')
        X=zeros(n,1); C=zeros(1,n+1);
        for p= 1:n-1
            for k=p+1:n
                m= B(k,p)/ B(p,p);
                B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1);
            end
        end
        b=B(1:n,n+1);A=B(1:n,1:n);
        X(n)=b(n)/A(n,n);
        for q=n-1:-1:1
            X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)))/A(q,q);
        end
    else
        disp('请注意：因为RA=RB<n，所以此方程组有无穷多解.')
    end
   end

例 7.1

​ 用高斯消元法和 MATLAB 程序求解下面的非齐次线性方程组，并且用逆矩阵解方程组的方法验证。
$$\{\begin{array}{l}{x}_1-x_2+x_3-3x_4=1\\ -x_2-x_3+x_4=0\\ 2x_1-\end{array}$$