数值分析.doc

《数值分析》 数 值 实 验 班级：101001 姓名：龙云祥 学号：101001113 指导老师：任林源 完成日期：2012.7.4 实验一：求解线性方程组 一．实验名称： 用高斯消去法和选列主元消去法求解线性方程组。 二．实验目的： 掌握高斯消去法和选列主元消去法的基本思想，能够熟练运用两种方法求解线性方程组，并进行比较两种方法的特点，能够根据不同的题目而选择最佳的求解方法。 运用计算机Matlab软件编程实现算法的程序设计，验证算法的准确性。 三．实验内容： Consider system of equations 0.4096x1+0.1234x2+0.3678x3+0.2943x4=0.4063 0.2246x1+0.3872x2+0.4015x3+0.1129x4=0.1550 0.3645x1+0.1926x2+0.3781x3+0.0643x4=0.4240 0.1784x1+0.4002x2+0.2786x3+0.39274x4=-0.2557 (a)Use the Gaussian elimination to solve the system of equations. (b)Use columnelimination with pivoting to resolve the system of equations and compare theresult with(a). 其中文翻译为： 考虑方程组 0.4096x1+0.1234x2+0.3678x3+0.2943x4=0.4063 0.2246x1+0.3872x2+0.4015x3+0.1129x4=0.1550 0.3645x1+0.1926x2+0.3781x3+0.0642小=0.4240 0.1784x1+0.4002x2+0.2786x3+0.3927x4=-0.2557 ( a ) 用高斯消元法解线性方程组 （b）用选列主元消去法解线性方程组 四．实验过程： （一）用高斯消元法求解线性方程组 1.算法思想（原理）： 首先逐步消去变元的系数，把原方程组化为系数矩阵为上三角的等价方程组，然后解系数矩阵为上三角行的方程组，即可得到原方程组的解，逐步消去变元的系数。将原方程组化为系数矩阵为上三角形等价方程组的过程为消元过程；把求系数矩阵为上三角方程组的解的过程为回代过程，求解线性方程组的高斯消元法，就是由消元过程和回代过程组成。 2.算法步骤： （1）消元过程 第一步，对第一列进行消元，将a21，a31，…an1消为零； 第二步，对第二列进行消元，将a32，a42，…an2消为零。 设akk!=0(k=1,2,…,n-1).做完第k-1步消元后，系数矩阵化为 A11 a21 … ak1 … an1 0 a22 … ak2 … an2 A= 0 0 … ak2… an2 一般对于第k次消元，采用如下公式： l= a/a(i=k+1,k+2…,n) a=a-la(i=k+1,k+2,…n;j=k.k+1,…,n) 直到做完n-1步，相应的方程组化为与原方程组等价的方程组为： Ax=y （2）回代过程.回代过程就是解方程组Ax=y，采用如下公式： x=y/a,x= y-ax/ a(k=n-1,n-2,…,2,1) （二）列主元消去法解线性方程组： ①算法思想： 列主元消去法有足够的精确度。只要detA≠0，从理论上讲，列主元消去法总可以进行下去。因为经选取主元后仍有a=0，则 必有矩阵第k列中对角线元素起以下均为0.故此列已消元完毕，可继续对第k+1列进行消元，最后得到的三角矩阵U，其对角线上有0元存在。 ②算法设计： 作第k步消元时，从A的子矩阵的第一列中选取|ak|=|a|,以a作为主元素，第i行称为主行。将主行与第k行交换，然后再作这一步的消元。 五．两种算法的分析比较： 从前述高斯消去法的消元过程有其不足之处，即在第k步消元是一定要假设a！=0，否则在消元过程中就会出现“主元素” a=0的情况，尽管det（A）！=0，此时消元过程也将无法进行下去。另外如果“主元素” a很小，由于计算机字长有限必然有舍入误差等因素的影响，其本身常常有较大的相对误差，用它做除数会导致其他因素量级严重增长和舍入误差的扩散，这样所得解极不准确，稳定，甚至可能产生溢出停机。而选列主元消去法解决了“主元素” a小或为零的情况，并且，列主元消去法有足够的精确度，而且，在选列主元消去法中，