资源描述：
《线性代数方程组数值解法及matlab实现综述》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、.线性代数方程组数值解法及MATLAB实现综述廖淑芳20122090数计学院12计算机科学与技术1班（职教本科）一、分析课题随着科学技术的发展，提出了大量复杂的数值计算问题，在建立电子计算机成为数值计算的主要工具以后，它以数字计算机求解数学问题的理论和方法为研究对象。其数值计算中线性代数方程的求解问题就广泛应用于各种工程技术方面。因此在各种数据处理中，线性代数方程组的求解是最常见的问题之一。关于线性代数方程组的数值解法一般分为两大类：直接法和迭代法。直接法就是经过有限步算术运算，可求的线性方程组精确解的方法（若计算过程没有舍入误差），但实际
2、犹如舍入误差的存在和影响，这种方法也只能求得近似解，这类方法是解低阶稠密矩阵方程组级某些大型稀疏矩阵方程组的有效方法。直接法包括高斯消元法，矩阵三角分解法、追赶法、平方根法。迭代法就是利用某种极限过程去逐步逼近线性方程组精确解的方法。迭代法具有需要计算机的存储单元少，程序设计简单，原始系数矩阵在计算过程始终不变等优点，但存在收敛性级收敛速度问题。迭代法是解大型稀疏矩阵方程组（尤其是微分方程离散后得到的大型方程组）的重要方法。迭代法包括Jacobi法SOR法、SSOR法等多种方法。二、研究课题-线性代数方程组数值解法一、直接法1、Gauss消
3、元法通过一系列的加减消元运算，也就是代数中的加减消去法，以使A对角线以下的元素化为零，将方程组化为上三角矩阵；然后，再逐一回代求解出x向量。1.1消元过程1.高斯消元法(加减消元)：首先将A化为上三角阵，再回代求解。步骤如下：...第一步：第二步：类似的做下去，我们有：第k步：。n－1步以后，我们可以得到变换后的矩阵为：注意到，计算过程中处在被除的位置，因此整个计算过程要保证它不为0。所以，Gauss消元法的可行条件为：。就是要求A的所有顺序主子式均不为0，即因此，有些有解的问题，不能用Gauss消元求解。另外，如果某个很小的话，会引入大的
4、误差。...例用Gauss消去法解方程组：（1）（1）对增广矩阵进行初等变换得等价方程组回代得，，，。第一步：将(1)/3使x1的系数化为1,再将(2)、(3)式中x1的系数都化为零，即由(2)-2×(1)(1)得...由(3)-4×(1)(1)得第二步：将(2)(1)除以2/3，使x2系数化为1，得再将(3)(1)式中x2系数化为零，由(3)(1)-(-14/3)*(2)(2)，得第三步：将(3)(2)除以18/3，使x3系数化为1，得经消元后，得到如下三角代数方程组：1.2回代过程由(3)(3)得x3=1,将x3代入(2)(2)得x2=
5、-2,将x2、x3代入(1)(1)得x2=1,所以，本题解为[x]=[1,2,-1]T1.3高斯消元的公式综合以上讨论，不难看出，高斯消元法解方程组的公式为第一步，消元（1）令aij(1)=aij,（i,j=1,2,3,…,n）bi(1)=bi,（i=1,2,3,…,n）（2）对k=1到n-1，若akk(k)≠0，进行lik=aik(k)/akk(k),（i=k+1,k+2,…,n）aij(k+1)=aij(k)-lik*akj(k),（i,j=k+1,k+2,…,n）bi(k+1)=bi(k)-lik*bk(k),（i=k+1,k+2,…
6、,n）第二步，回代若ann(n)≠0xn=bn(n)/ann(n)xi=(bi(i)–sgm(aij(i)*xj)/-aii(i),（i=...n-1,n-2,…,1）,(j=i+1,i+2,…,n)2、LU分解法求解线性代数方程组除了高斯消元法外，还常用LU分解法（三角形分解法）。LU分解法的优点是当方程组左端系数矩阵不变，仅仅是方程组右端列向量改变，即外加激励信号变化时，能够方便地求解方程组。矩阵的三角分解法可分为直接三角分解法,列主元三角分解法,平方根法,三对角方程组的追赶法。下面讨论直接三角分解法。设n阶线性方程组Ax=b。假设能将
7、方程组左端系数矩阵A，分解成两个三角阵的乘积，即A=LU，式中，L为主对角线以上的元素均为零的下三角矩阵,且主对角线元素均为1的上三角矩阵；U为主对角线以下的元素均为零所以有，LUx=b令Ux=y，则Ly=b由A=LU,由矩阵的乘法公式：a1j=u1j,j=1,2,…,nai1=li1u11,i=1,2,…,n推出u1j=a1j,j=1,2,…,nli1=ai1/u11,i=1,2,…,n这样就定出了U的第一行元素和L的第一列元素。设已定出了U的前k-1行和L的前k-1列，现在确定U的第k行和L的第k列。由矩阵乘法：当r>k时，lkr=0,
8、且lkk=1，因为...所以，同理可推出计算L的第k列的公式：因此得到如下算法——杜利特（Doolittle）算法:(1)将矩阵分解为A=LU,对k=1,2,…,n；j=k,k+
