2012数值分析 -lec3--高斯消去法x.pptVIP

数值分析 朱立永 北京航空航天大学 数学与系统科学学院 Email: numerical_analysis@ Password:beihang2012 答疑时间：星期四下午14:00－17:00 答疑地点：主216 第三讲线性方程组解法 问题: In Scientific Computing ↓ Large Linear Systems Ax=b as sub-problems/ as intermediate steps Gauss-Seidel method Jacobi method SOR method Conjugate Gradient method for symmetric systems Gaussian elimination LU factorization Cholesky factorization GMRES GCR Bi-CG CGS Bi-CGSTAB Bi-CGSTAB2 GPBi-CG Bi-CGSTAB(L) 例：用单精度计算 的根 精确解为 算法1：利用求根公式 在计算机内，109存为0.1?1010，1存为0.1?101。做加法时，两加数的指数先向大指数对齐，再将浮点部分相加。即1 的指数部分须变为1010，则：1 = 0.0000000001 ? 1010，取单精度时就成为： 109+1=01010+0?1010=0?1010 大数吃小数导致严重后果！ 算法2：先解出 再利用 注：求和时从小到大相加，可使和的误差减小。 例：按从小到大、以及从大到小的顺序分别计算 1 + 2 + 3 + … + 40 + 109 其它问题 奇异矩阵的求解：高斯消去法；迭代法 高斯消去法针对的是一般的方程组： 对称、非对称、正定、不定的都可以求解 作业 教材第45页，习题1. 掌握（或者自己实现）列主元高斯消去法的程序 课后阅读：《C数值算法》第二章 数值分析