计算方法第二单元报告.pptVIP

目录 一、 引言 二、 算法描述 三、 试验及分析 四、 自己的特色 五、 总结 一、引言 在我们的实际生活和生成当中会出现大量的方程组，在这些方程组需要求解，如果是要用人工的去求解的时候无疑是一个庞大的工作，而且很能避免错误的出现，这样会浪费大量的人力物力和时间。于是我们今天在前辈们关于方程组的求解总结出来的高斯消去法和直接三角分解法（LU分解）的基础之上借助编程过程中的循环来让计算机对方程组经行求解。这样可以节省大量的时间同时避免了错误的出现。下面将是我们本着这个目的去实现高斯消去法和直接三角分解法（LU分解）的过程。虽然有很多的难题还没有解决，但是这是我们小组共同努力的成果。 二、算法描述 1、高斯列主元消去法的原理 基本思想：一次用前面的方程消去后面的未知数，从而将方程 组化为等价形式，即消元 然后再回代，对于 ，若 依次计算 然后将其回代得到： 以上过程就是高斯消去。 说明： 高斯消去法在消元的过程中有可能会出现 的情况，这时消元就无法进行了，即使主元数 但是很小时，其做除数，也会导致其他元素数量级的严重增长和舍入误差的扩散。 因此，为了减少误差，每次消元选取系数矩阵的某列中绝对值最大的元素作为主元素。然后换行使之变到主元位置上，再进行消元计算。即高斯列主元消去法。 二、算法描述 2、直接三角分解法（LU分解）的原理 基本思想：先将矩阵A直接分解为则求解方程组的问题就等价于求解两个三角形方程组。直接利用矩阵乘法，得到矩阵的三角分解计算公式为： 由上面的式子得到矩阵A的LU分解后，求解Ux=y的计算公式为 以上就是LU分解法。 三、实验及分析 高斯消去 测试数据一：行数3 列数4 1 2 3 6 2 2 3 7 1 4 4 9 答案： 1 1 1 测试数据二：3 4 0 0 1 3 3 0 0 2 0 2 0 1 答案：0.667 0.5 3 LU分解算法： 实例1: 实例2: 实例3: 四、自己的特色(1) 1.用高斯算法解一般的方程组 测试数据一：3 5 10 200 300 4 233 22 22 44 1 22 1 20 30 0.4 233 测试数据二：4 6 100 34 345 334 11 552 12 0.2 0.12 12 23 344 122 53 12 43 23 45 123 234 0.3 34 345 34 测试数据三：行4 列5 21 345 67 15 456 0.12 35 02.4 234 3455 122 34 13