微电子器件及工艺CAD第二章 基本数学方法.pptVIP

-1- 微分方程的有限差分解法 2.1.1 微分方程数值解的引入 微分方程的有限差分解法 微分方程的数值解法有许多种，有限差分法是微分方程数值解法的一种，是半导体器件数值分析中最常用的方法。后面的章节中，我们还将介绍另一种应用更普遍的离散数值分析法-----有限元法。 微分方程的有限差分解法 2.1.2.2 微分方程的有限差分法求解 微分方程的有限差分解法 常微分方程的边值问题又可以分为下列三种类型。设有二阶常微分方程y(x)，定义区间为[a,b]。 微分方程的有限差分解法 非线性微分方程组的求解 本节介绍一种非线性方程组的求解方法---- Newton迭代法，下一节将介绍线性方程组的求解方法。 非线性方程组的求解 线性代数方程组的求解 线性代数方程组的求解 线性代数方程组的求解 4. 追赶法 如果线性方程组的系数是三对角线的，则高斯消去法可以简化。 线性代数方程组的求解 作业 1. 说明Jacobi迭代和Gauss-Seidel迭代的区别。 2. 设矩阵A=LU，L为下三角矩阵，U为上三角矩阵，求L，U。 §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列 ，得