线性方程组的数值解法
——高斯列主元消去法解线性方程组的MATLAB实现
班级:MATH 20XX 学号:xxxxxxxxxx 姓名:LI Wen
高斯列主元消去法解线性方程组的MATLAB实现
摘 要
在自然科学和工程技术中许多问题的解决常常归结为解线性代数方程组,并且在实际计算中很多线性方程组是无法求其精确解的,所以只能求他的数值解。而解线性方程组的数值解法一般有两类:直接法和迭代法。其中直接法就是经过有限的算术运算,求得线性方程解的方法。但在实际计算中由于舍入误差的存在和影响,这种方法也只能求的线性方程组的近似解。本论文将阐述这类方法中的最基本的高斯消去法及其变形而来的列主元消去法,并运用MATLAB程序来实现它。由于高斯消去法在计算的过程中,会出现某个主元等于零的情况,从而无法计算;或主元绝对值很小的情况,也使得结果与真实结果相差甚远。为了避免高斯消去法消元过程中出现的上述两个问题,一般在实际计算中采用所谓的选择主元法。目前计算机上常用的是按列选主元的方法。因此本文中做的也是列选主元高斯消去法。下面本文将详细介绍高斯消去法和列主元消去法的计算过程,并给出列主元消去法的程序算法。
关键词:线性方程组,消去法,列主元,MATLAB,数值计算
引言
在自然科学和工程技术中许多问题的解决常常归结为解线性代数方程组,例如电学中的网络问题,船体数学放样中建立三次样条函数问题,用最小二乘法求实验数据的曲线拟合问题,解非线性方程组问题,用差分法或者有限元法解常微分方程,偏微分方程边值问题等都导致求解线性方程组,而且后面几种情况常常归结为求解大型线性方程组。
关于解线性方程组的数值解法一般有两类:直接法和迭代法。直接法就是经过有限的算术运算,可求的线性方程组精确解的方法(若计算过程中没有舍入误差)。但在实际计算中由于舍入误差的存在和影响,这种方法也只能求的线性方程组的近似解。本论文将阐述这类方法中的最基本的高斯消去法及其变形而来的列主元消去法。
一、问题的提出
我们都知道,高斯列主元素消去法是计算机上常用来求解线性方程组的一种直接的方法。就是在不考虑舍入误差的情况下,经过有限步的四则运算可以得到线性方程组的准确解的一类方法。实际运算的时候因为只能有限小数去计算,因此只能得到近似值。在实际运算的时候,我们很多时候也常用高斯消去法。但是高斯消去法在计算机中运算的时候常会碰到两个问题。
1)一旦遇到某个主元等于0,消元过程便无法进行下去。
2)在长期使用中还发现,即使消元过程能进行下去,但是当某个主元的绝对值很小时,求解出的结果与真实结果相差甚远。
为了避免高斯消去法消元过程中出现的上述两个问题,一般采用所谓的选择主元法。其中又可以分为列选主元和全面选主元两种方法。目前计算机上常用的按列选主元的方法。因此我在这里做的也是按列选主元高斯消去法。
二、高斯消去法的基本思想
大家知道,如果一个线性方程组的系数矩阵是上三角矩阵时,即这种方程组我们称
之为上三角方程组,它是很容易求解的。我们只要把方程组的最下面的一个方程求解出来,在把求得的解带入倒数第二个方程,求出第二个解,依次往上回代求解。然而,现实中大多数线性方程组都不是上面所说的上三角方程组,所以我们有可以把不是上三角的方程通过一定的算法化成上三角方程组,由此我们可以很方便地求出方程组的解。高斯消元法的目的就是把一般线性方程组简化成上三角方程组。于是高斯消元法的基本思想是:通过逐次消元将所给的线性方程组化为上三角形方程组,继而通过回代过程求解线性方程组。
三、高斯消去法的计算过程
设有线性方程组
线性方程组的数值解法
——高斯列主元消去法解线性方程组的MATLAB实现
班级:MATH 20XX 学号:xxxxxxxxxx 姓名:LI Wen
高斯列主元消去法解线性方程组的MATLAB实现
摘 要
在自然科学和工程技术中许多问题的解决常常归结为解线性代数方程组,并且在实际计算中很多线性方程组是无法求其精确解的,所以只能求他的数值解。而解线性方程组的数值解法一般有两类:直接法和迭代法。其中直接法就是经过有限的算术运算,求得线性方程解的方法。但在实际计算中由于舍入误差的存在和影响,这种方法也只能求的线性方程组的近似解。本论文将阐述这类方法中的最基本的高斯消去法及其变形而来的列主元消去法,并运用MATLAB程序来实现它。由于高斯消去法在计算的过程中,会出现某个主元等于零的情况,从而无法计算;或主元绝对值很小的情况,也使得结果与真实结果相差甚远。为了避免高斯消去法消元过程中出现的上述两个问题,一般在实际计算中采用所谓的选择主元法。目前计算机上常用的是按列选主元的方法。因此本文中做的也是列选主元高斯消去法。下面本文将详细介绍高斯消去法和列主元消去法的计算过程,并给出列主元消去法的程序算法。
关键词:线性方程组,消去法,列主元,MATLAB,数值计算
引言
在自然科学和工程技术中许多问题的解决常常归结为解线性代数方程组,例如电学中的网络问题,船体数学放样中建立三次样条函数问题,用最小二乘法求实验数据的曲线拟合问题,解非线性方程组问题,用差分法或者有限元法解常微分方程,偏微分方程边值问题等都导致求解线性方程组,而且后面几种情况常常归结为求解大型线性方程组。
关于解线性方程组的数值解法一般有两类:直接法和迭代法。直接法就是经过有限的算术运算,可求的线性方程组精确解的方法(若计算过程中没有舍入误差)。但在实际计算中由于舍入误差的存在和影响,这种方法也只能求的线性方程组的近似解。本论文将阐述这类方法中的最基本的高斯消去法及其变形而来的列主元消去法。
一、问题的提出
我们都知道,高斯列主元素消去法是计算机上常用来求解线性方程组的一种直接的方法。就是在不考虑舍入误差的情况下,经过有限步的四则运算可以得到线性方程组的准确解的一类方法。实际运算的时候因为只能有限小数去计算,因此只能得到近似值。在实际运算的时候,我们很多时候也常用高斯消去法。但是高斯消去法在计算机中运算的时候常会碰到两个问题。
1)一旦遇到某个主元等于0,消元过程便无法进行下去。
2)在长期使用中还发现,即使消元过程能进行下去,但是当某个主元的绝对值很小时,求解出的结果与真实结果相差甚远。
为了避免高斯消去法消元过程中出现的上述两个问题,一般采用所谓的选择主元法。其中又可以分为列选主元和全面选主元两种方法。目前计算机上常用的按列选主元的方法。因此我在这里做的也是按列选主元高斯消去法。
二、高斯消去法的基本思想
大家知道,如果一个线性方程组的系数矩阵是上三角矩阵时,即这种方程组我们称
之为上三角方程组,它是很容易求解的。我们只要把方程组的最下面的一个方程求解出来,在把求得的解带入倒数第二个方程,求出第二个解,依次往上回代求解。然而,现实中大多数线性方程组都不是上面所说的上三角方程组,所以我们有可以把不是上三角的方程通过一定的算法化成上三角方程组,由此我们可以很方便地求出方程组的解。高斯消元法的目的就是把一般线性方程组简化成上三角方程组。于是高斯消元法的基本思想是:通过逐次消元将所给的线性方程组化为上三角形方程组,继而通过回代过程求解线性方程组。
三、高斯消去法的计算过程
设有线性方程组