线性方程组AX=B的数值计算方法实验课案.docxVIP

《数值计算方法》实验报告 线性方程组AX=B的数值计算方法实验 学号：2013040207012 姓名：梁哲豪 实验描述 在自然科学和工程技术中很多问题的解决常常归结为解线性代数方程组。例如电学中的网络问题，船体数学放样中建立三次样条函数问题，用最小二乘法求实验数据的曲线拟合问题，解非线性方程组问题，用差分法或者有限元法解常微分方程，偏微分方程边值问题等都导致求解线性方程组，而且后面几种情况常常归结为求解大型线性方程组。线性代数方面的计算方法就是研究求解线性方程组的一些数值解法与研究计算矩阵的特征值及特征向量的数值方法。 关于线性方程组的数值解法一般有两类： 直接法：若在计算过程中没有舍入误差，经过有限步算术运算，可求得方程组的精确解的方法。 迭代法：用某种极限过程去逐步逼近线性方程组精确解的方法。迭代法具有占存储单元少，程序设计简单，原始系数矩阵在迭代过程中不变等优点，但存在收敛性及收敛速度等问题。 上三角线性方程组的求解： 基本算法： 高斯消元法：将原方程组化为三角形方阵的方程组： lik=aik/akk aij=aij-lik*akj （k=1,2,…,n-1; i=k+1,k+2, …,n ;j=k+1,k+2, …,n+1） 由回代过程求得原方程组的解： xn=ann+1/ann xk=(akn+1-akjxj)/akk LU分解法： 将系数矩阵A转化为A=L*U