偏微分方程和边界条件： 加权余量法求解.PPT

第三讲 1.偏微分方程求解－－有限元法的原理（加权余量法和变分法） 2.数值求解方法 2.数值求解方法 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 3. 加权余量法－－例1 3. 加权余量法－－例1 3. 加权余量法－－例1 3. 加权余量法－－例1 3. 加权余量法－－例1 3. 加权余量法－－例1 3. 加权余量法－－例1 4. 加权余量法求解一般化偏微分方程的归纳 4. 加权余量法求解一般化偏微分方程的归纳 4. 加权余量法求解一般化偏微分方程的归纳 5. 加权余量法的进一步优化（边界条件的处理） 5. 加权余量法求解一般化方法的进一步优化 5. 加权余量法求解一般化方法的进一步优化 5. 加权余量法求解一般化方法的进一步优化 5. 加权余量法求解一般化方法的进一步优化 5. 加权余量法求解一般化方法的进一步优化 5. 加权余量法求解一般化方法的进一步优化 6. 简化后加权余量法 例2 6. 简化后加权余量法 例2 6. 简化后加权余量法 例2 6. 简化后加权余量法 例2 6. 简化后加权余量法 例2 6. 简化后加权余量法 例2 7. 简化后加权余量法 求解一般化的微分方程 例3 7. 简化后加权余量法 求解一般化的微分方程 例3 7. 简化后加权余量法 求解一般化的微分方程 例3 7. 简化后加权余量法 求解一般化的微分方程 例3 7. 简化后加权余量法 求解一般化的微分方程 例3 7. 简化后加权余量法 求解一般化的微分方程 例3 7. 简化后加权余量法 求解一般化的微分方程 例3 8. 归纳加权余量求解偏微分方程步骤 8. 归纳加权余量求解偏微分方程步骤 9. 变分法简介 9. 变分法简介－－拉普拉斯方程 9. 变分法简介－－拉普拉斯方程 9. 变分法简介—帕松方程 9. 变分法简介—帕松方程 9. 变分法简介—赫姆霍兹和一般化偏微分方程（省略） 9. 变分法简介—赫姆霍兹和一般化偏微分方程（省略） 3.代公式计算矩阵元素 （边界矩阵b为0） 4. 封装矩阵： 5. 求解矩阵，得近似解： 该有源静态电场问题的真解（解析解：） 真解与近似解相同是由于尝试函数选择的刚好，通常有差别。如例3 偏微分方程描述的问题如下： 加权余量法求解： 1.初选尝试函数、构造近似解： 利用问题及其边界条件，对尝试函数进行优化（使近似解满足边界条件） 通过尝试函数的选取，近似解满足1类边界条件，使得1类边界条件在方程中消失 两个方程，两个独立未知数，消a1、a2，重定尝试函数，边界条件自动满足，简化求解过程 2. 修正尝试函数，以满足1、2类边界条件： 余数为： 结合问题，余数的具体表达式为： 问题的加权余数（目标泛函）为： 4. j=2,3时得代数方程： 5. 求解矩阵，得近似解： 5. 求解矩阵，得待定系数和近似解： 真解（解析解：） 加权余量法求解流程： 1.初步选取尝试函数、构造近似解 2.结合问题的边界条件对尝试函数进行修正，以简化求解 3.写出余数表达式 3. 写出加权余数表达式（迦辽金方法选取加权函数） 4. 令权余数表达式在各尝试函数下为0，得到代数方程组，解之得到待定系数，从而确定近似解 加权余数法求解一般性偏微分方程的方法： 方程的近似解被表示为一系列独立的尝试函数的线性组合，其中包括未知的待定系数。 通常用迦辽金原理选取加权函数，（即令加权函数等于尝试函数本身），从而完成对加权余数的定义，(尝试函数的选取满足边界条件) 通过对加权函数在区域内和在边界上的积分使其平均值为零，也就是说，使近似解与精确解之间的差别在某种指标下达到最小化。 如此可以形成一个矩阵形式的代数方程组，求解该矩阵方程可以确定待定系数．从而得到偏微分方程的唯一近似解。 另外一种求解偏微分方程的一般方法，即变分法。 变分法与加权余数法类似，近似解也用一系列线性独立的尝试函数表示．包括未知的待定系数。 与加权余数法不同的是，变分法用另外的方法来形成求解待定系数的矩阵方程。在变分法中，首先要构成一个近似解的函数，称为泛函。从广义来说，加权余数积分(即平均值)也是一种泛函。 然后使该泛函最小化，从而减小近似解的误差。一般说来，要找到一个适合于偏微分方程及边界条件的泛函是一项难度很大的工作。由于前人已做了许多研究工作，已找到了适合于许多常