大连海事大学《数值分析》课件-第九讲（高老师）.pptVIP

方程形态的演变（k）方程形态的演变（n）高斯消去法的回代过程高斯消去法的计算量分析消元第k步的计算量为(n-k+1)(n-k+1)回代第k步的计算量为(k-1)高斯消去法总的乘除运算量为：高斯主元素消去法由高斯消去法知道，在消元的过程中可能出现的情况，这是消去法将无法进行即使主元素但很小时，用其作除数，会导致其他元素数量级的严重增长和舍入误差的扩散，最后也使得计算解不可靠。例2错误算法例2改进算法全主元高斯消去法第k步消元时选A(k)中绝对值最大的元素为主元，即①先选取全主元：ifik?kthen交换第k行和第ik行；ifjk?kthen交换第k列和第jk列；③消元列交换改变了xi的顺序，须记录交换次序，解完后再换回来。全主元高斯消去法具有很好的稳定性，但选全主元比较费时，故在实际计算中很少使用。定理定理1定理2假设方程组（4）是对角占优的，则（k=1，2，…，n）全不为0.假设方程组（4）对称并且是对角占优的，则（k=1，2，…，n）全是主元素.定理1假设方程组（4）是对角占优的，则（k=1，2，…，n）全不为0.证：定理1定理1定理1变换以后仍然是对角占优的，以此类推，则命题得证。例题2用列主元消去法求解下列方程组举例（续）SOR迭代格式取w=1.1，得x(1)=(0.5500,3.1350,-1.0257)Tx(7)=(2.0000,3.0000,-1.0000)T如何确定SOR迭代中的最优松弛因子是一件很困难的事。矩阵分裂法Jacobi迭代GS迭代SOR迭代A=M-NM=D,N=M–A=L+UM=D–L,N=U第五章线性方程组的迭代解法方程组系数矩阵的分类低阶稠密矩阵（例如，阶数不超过150）（一般用直接法来求解）大型稀疏矩阵（即矩阵阶数高且零元素较多）（一般用迭代法来求解）线性方程组的数值解法分类直接法经过有限步算术运算，可求得方程组精确解的方法。迭代法用某种极限过程去逐步逼近线性方程组精确解的方法。消去法的基本思想通过将一个方程乘或除以某个常数，以及将两个方程相加减这两种手续，逐步减少方程中变元的数目，最终使每个方程仅含一个变元，从而得出所求的解.约当（Jordan)消去法高斯（Gauss)消去法约当消去法特点它的每一步仅在一个方程中保留某个变元，而从其余的各个方程中消去该变元，这样经过反复消元后，所给方程组最终被加工成一个方程仅含一个变元的形式，从而得出所求的解.AX=b方程组的向量表示形式方程形态的演变（0）方程形态的演变（1）方程形态的演变（1）方程形态的演变（2）方程形态的演变（2）方程形态的演变（k-1）方程形态的演变（k）方程形态的演变（k）方程形态的演变（n）约当消去法第（1）步初始状态第（1）步约当消去法第（1）步运算过程约当消去法第（k）步第（k-1）步第（k）步约当消去法第（k）步运算过程约当消去法的计算量与存储空间计算量存储空间高斯消去法虽然高斯消去法是一个古老的求解线性方程组的方法（早在公元前250年我国就掌握了解方程组的消去法），但是由它改进、变形得到的选主元素消去法、三角分解法仍然是目前计算机上常用的有效方法。举例（一）解：例：直接法解线性方程组高斯消去法的基本思想考虑n阶线性方程组：矩阵形式=方程形态的演变（0）方程形态的演变（1）方程形态的演变（1）方程形态的演变（2）方程形态的演变（2）方程形态的演变（k-1）方程形态的演变（k）*关键的部分在于如何在计算机上实现这些算法，并使得计算复杂度和空间复杂度最好*通过分析每一步消元以后的状态，学习编写求解方程组的程序编写方法，思路清晰。*第一列的值不必计算*第一列的值不必计算*第一列的值不必计算*第一列的值不必计算*第一列的值不必计算*第一列的值不必计算*第一列的值不必计算*第一列的值不必计算*第一列的值不必计算*第一列的值不必计算第五章线性方程组的迭代法高云大连海事大学数值分析迭代法迭代法：从一个初始向量出发，按照一定的迭代格式，构造出一个趋向于真解的无穷序列迭代解法是目前求解大规模线性方程组的主要方法。只需存储系数矩阵中的非零元素运算量不超过O