李建军-半导体器件模拟及数值分析第二章 基本数学方法.pptVIP

§2-3 线性代数方程组的求解 线性代数方程组的求解 线性代数方程组的求解 线性代数方程组 求解算法: 直接法 迭代法 高斯消去法 LU分解法 高斯消去法 线性代数方程组 A X = B 其中 2 1 2 1 ) , , ( ) , , ( T n T n b b b B x x x X L L = = 为了方便讨论，我们将其写成如下形式 2 22 21 1 12 11 n n a a a a a a a a a A L M M M M L L = n1 n2 nn 线性代数方程组的求解 高斯消去法 线性代数方程组 A X = B 设 第一次消元：由第 个方程减去第1个方程乘以 ,则将上面方程组中的第 其中记 线性代数方程组的求解 高斯消去法 个方程中的第一个未知数 消去，得到如下方程组 其中 设 由第 个方程减去上面方程组中的第2个方程乘以 ，则将其中的 个方程 线性代数方程组的求解 高斯消去法 中第2个未知数 消去，得到如下方程组 其中 继续这个过程，经过n-1次消元后，方程变为 线性代数方程组的求解 高斯消去法 其中 线性代数方程组的求解 高斯消去法 于是可以从第n个方程开始逐次向上解方程。 线性代数方程组的求解 高斯消去法 以上的计算是假设 ，事实上，即使 ， 但如果 的值很小，消去过程也无法进行，需要进行列选主元素，这里只解释为什么 的值很小时不可以。 经过第一次消元，第2个方程减去第1个方程乘以 得 假设计算机可以保留10位有效数字，用消元法解下面的方程组。 线性代数方程组的求解 高斯消去法 其中 解得 而真解为 线性代数方程组的求解 高斯消去法 交换方程组的第1行和第2行，得 再用消去法，可得真解。 线性代数方程组的求解 2.LU分解法 线性代数方程组A X = B A分解为上三角阵U和下三角阵L，即A=LU。 则AX＝LUX＝B，设UX＝Y，那么LY＝B 线性代数方程组的求解 2.LU分解法 下面介绍L，U的求法。 首先，用L阵的第1行分别乘U阵的各列，算出U阵的第1行各元素。 然后，用L阵的各行分别去乘U阵的第1列，算出L阵的第1列 线性代数方程组的求解 2.LU分解法 再用L阵的第2行分别乘U阵的各列，算出U阵的第2行个元素，用L阵的各行分别去乘U阵的第2列，算出L阵的第2列。 用L阵的第r行分别乘U阵的j列 ，得 线性代数方程组的求解 按照此方法一直算下去，现在假设已经算出U阵的前r-1行元素，L阵的r-1列元素，下面来算U阵的第r行元素，L阵的第r列元素。 2.LU分解法 所以，得U阵的第r行元素 再用L阵的第i行 分别乘U阵的第r列，得 ， 所以，得L阵的第r列元素 线性代数方程组的求解 2.LU分解法 由于AX＝LUX＝B，设UX＝Y，那么LY＝B Ax＝B，等价为 UX＝Y LY＝B 逐次用向后代入过程先解LY＝B得 线性代数方程组的求解 2.LU分解法 由于AX＝LUX＝B，设UX＝Y，那么LY＝B Ax＝B，等价为 UX＝Y LY＝B 然后再用逐次向前回带过程解UX=Y得 线性代数方程组的求解 3.迭代法 直接法一般工作量较大，本节介绍迭代法，将AX=B系数矩阵