一阶偏微分方程求解方法.ppt

1.偏微分方程求解－－有限元法的原理（加权余量法和变分法） 2.数值求解方法 2.数值求解方法 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 3.电磁场位函数偏微分方程的数值求解方法－加权余量法 3. 加权余量法－－例 3. 加权余量法－－例 3. 加权余量法－－例 3. 加权余量法－－例 3. 加权余量法－－例 3. 加权余量法－－例 3. 加权余量法－－例 4. 加权余量法求解一般化偏微分方程的归纳 4. 加权余量法求解一般化偏微分方程的归纳 4. 加权余量法求解一般化偏微分方程的归纳 5. 加权余量法的进一步优化（边界条件的处理） 5. 加权余量法求解一般化方法的进一步优化 5. 加权余量法求解一般化方法的进一步优化 5. 加权余量法求解一般化方法的进一步优化 5. 加权余量法求解一般化方法的进一步优化 5. 加权余量法求解一般化方法的进一步优化 5. 加权余量法求解一般化方法的进一步优化 6. 简化后加权余量法 例2 6. 简化后加权余量法 例2 6. 简化后加权余量法 例2 6. 简化后加权余量法 例2 6. 简化后加权余量法 例2 6. 简化后加权余量法 例2 3.代公式计算矩阵元素 （边界矩阵b为0） 4. 封装矩阵： 5. 求解矩阵，得近似解： * * 解析法 应用范围有限，适用于理论求解，但有强烈的物理含义（常系数微分方程） 某些复杂问题，很考虑根本找不到解析解 2. 数值法 工程实际中应用广泛，复杂场域问题，但物理含义不很清楚。任何问题总可以找到数值解（数学方法） 2/4 1. 基本思想： 以偏微分方程的近似解来代替其真解，只要近似解与真解足够接近，就可以近似解作为问题的解，并满足足够的精度。 2. 基本方法： 假设一个近似解，该解一组（形式上）简单函数 的线性组合来表示，线性组合的系数就是一组待定系数 然后建立一种考虑了微分方程和边界条件的关于真解 和近似解间误差的目标函数 F 用适当的算法使得该目标函数最小化――最小化的过程就确定了待定系数，从而也就得到了问题的近似解。 尝试函数，基函数，形函数 2/4 目标函数最小化的目的：一方面，使得近似解最大程度接近真解；另一方面，求得构成近似解的待定系数。 数学上，构成目标函数的方法很多，不同的构成方法就形成了不同的数值解法，电磁场中就常见的是：加权余量法和变分法。 电磁场问题总可以用位函数的偏微分方程和相应的边界条件表述 两个偏微分方程形式相同，故以电位方程的求解过程为例。磁位矢量的方程可以分解到个分量上变为标量方程。 在求解场域内，偏微分方程的真解为 ，近似解为 它由一组简单函数 的线性组合表达，表达中有待定系数 即： 加权余量法 简单函数，一般选用简单形式的函数，一旦选定就是已知的了 待定系数是真正的求解目标 问题的自由度 近似解 加权余量法就是一种定义近似解与真解之间误差（即余数），并设法使其最小的方法。 加权余量法误差（即余数）的定义： 注意：一般余数并不表示近似解与真解间的差（场域内），加权余量法的采用拉普拉斯算子作用后的差别（即余数），来代表近似解接近偏微分方程真解的程度。 问题的自由度 当余数小于要求的精度时，就可以认为近似解就是偏微分方程的解。 要减少余数，我们可以通过寻求适当的待定系数来实现。 为有效表达减小余数的效果，还选取适当的加权函数，以使余数和该加权函数的积分为0。－－“加权余量法”的来由。 加权余数的定义： 加权函数的选取方法很多：如点重合、子域重合、最小二乘法、迦辽金法。 效果较好的、运用较多的是迦辽金法： 即：迦辽金法选取尝试函数本身为加权函数 由此构建加权量法的目标函数： 上述过程中，已经将偏微分方程转化为j个代数方程组，便于计算机求解。 关于函数是函数，称为：泛函数，或泛函 例1.两极电容板内部电场分布问题： 根据问题特点将3维问题简化为2维， 进一步简化为1维。 该问题是静态电场问题， 偏微分方程和边界条件： 加权余量法求解： 1.选取尝试函数、构造近似解： 2.结合问题，写出余数表达式： 理论上任意选取，操作中越简单越好 2.结合问题，写出余数表达式： 3. 加权余数表达式： 3. 加权余数表达式： 4. 求解上述两个代数方程组