追赶法，高斯消元法，逆矩阵法.docVIP

追赶法，高斯消元法，逆矩阵法，迭代法 —— 解线性方程组 精仪学院 马金玉 1012202030 本文主要详细介绍了追赶法，高斯法，逆矩阵法的方法原理，运用这三种方法分别进行线性方程的求解举例，给出MATLAB相应程序，最后做结果分析，比较说明追赶法和高斯法的特点。最后对三种典型迭代方法Jacobi迭代,Gauss-Seidel迭代，SOR迭代进行简单的分析比较。 追赶法 1.1）.追赶法方法介绍 追赶法用于求解以下形式的方程组(三对角方程组) 其中 ，系数矩阵(三对角矩阵) 系数矩阵A的元素满足 第一步：实现A=LU的分解，按照递推公式 计算 ： 第二步：求解方程组LY=f,相应的递推公式 第三部：求解方程组UX=Y，相应的递推公式 求得x 因为计算 及 的过程是追赶的过程，结出结果X。 1.2）.追赶法解线性方程组的matlab实例 解线性方程组 第一步：编写M文件如下： function [x,y,beta]=zhuiganfa(a,b,c,f) %a,b,c是三对角阵的对角线上的元素,f是自由项. n=length(b); beta(1)=c(1)/b(1); for i=2:n beta(i)=c(i)/(b(i)-a(i)*beta(i-1)); end y(1)=f(1)/b(1); for i=2:n y(i)=(f(i)-a(i)*y(i-1))/(b(i)-a(i)*beta(i-1)); end x(n)=y(n); for i=n-1:-1:1 x(i)=y(i)-beta(i)*x(i+1); end disp(sprintf(k x(k) y(k) beta(k))); for i=0:n-1 disp(sprintf(%d %15.4f %15.4f %15.4f,i,x(i+1),y(i+1),beta(i+1))); end 追赶法M文件程序截图如图1所示 图1 追赶法M文件程序截图 第二步：根据所求方程，在命令窗口中输入如下命令，并按ENTER键确认。 a=[0,1,1,2];b=[2,3,1,1];c=[1,1,1,0];f=[1,2,2,0];[x,y,beta]=zhuiganfa(a,b,c,f) 追赶法命令窗口的运行结果如图2所示， 图2 追赶法运行结果截图 高斯法 2.1）高斯消元法方法介绍 高斯顺序消去法:设线性方程组AX＝b, 对增广矩阵[A|b ]顺序作初等行变换，使矩阵A化为上三角形矩阵，再回代，从而求得线性方程组的解。主要有消元和回代两个过程。 消元过程。设 ,对 计算 回代过程。 2.2）高斯消元法解线性方程组的matlab实例 解线性方程组 第一步：编写M文件如下： function [RA,RB,n,X]=gaus(A,b) B=[A b]; n=length(b);RA=rank(A); RB=rank(B);zhica=RB-RA; if zhica0, 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=RBn，所以此方程组有无穷多解.) end end 高斯消元法M文件程序截图如图3所示 图3 高斯消元法M文件程序截图 第二步：根据所求方程，在命令窗口中输入如下命令，并按ENTER键确认。 A=[2,1,0,0;1,3,1,0;0,1,1,1;0,0,2,1];b=[1,2,2,0]; [RA,RB,n,X] =gaus (A,b) 高斯法命令窗口的运行结果如图4所示， 图4 高斯法命令窗口的运行结果 3)逆矩阵法 本文用逆矩阵法进行检验。方法介绍：因为AX=b,得到 ，从而的 ，即可算出结果。