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

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