第十章 偏微分方程数值解.pptVIP

目的：求解相应的变分问题或相应的变分方程。 Ritz方法是近似求解变分问题(即二次泛函极小值)的算法。Galerkin方法是近似求解变分方程的算法，这两种算法统称为Ritz-Galerkin方法。 Ritz-Galerkin方法的基本思想 以下用V表示 等Sobolev空间，L表示微分算子，?(u,v)为由L及边值条件决定的双线性泛函。 4 Ritz-Galerkin方法 用有限维空间的函数代替变分问题(或变分方程)中无限维空间的函数，从而在有限维函数空间中求变分问题(或变分方程)的近似解，并要求当有限维空间的维数不断增加时，有限维近似解逼近原变分问题(或变分方程)的解。 由极小位能原理得出的变分问题为： Ritz方法：求变分问题的近似解。 （1）Ritz方法 求 ，使 其中 ， 设 是V 的n维子空间， 是 的一组基底(称为基函数) 。 中任一元素 可表示为 即选择适当的 ，使 取极小值。 求 ，使 Ritz方法： 展开 令 则 满足 解出 代入 ，则得 Ritz方法步骤为： 根据最小位能原理构造相应于微分方程或物 理问题的变分问题； 取 作为 的一组基底，即用 近似代替无穷维空间V； 根据二次函数取极值的必要条件，得到 中 所满足的方程组： 求解关于 的线性代数方程组。 由虚功原理得出的变分方程为： Galerkin方法：求变分方程的近似解。 （2） Galerkin方法 设 是V 的n维子空间， 是 的一组基底(称为基函数) 。 中任一元素 可表示为 即选择适当的 ，使 取极小值。 Galerkin方法： 求 ，使对 ， 满足 由 的任意性，取 作为v ，则得 将 代入变分方程，则 解出 代入 ，则得 Galerkin步骤为： 根据虚功原理构造相应于微分方程或物理问 题的变分方程； 取 作为 的一组基底，即用 近似代替无穷维空间V； 求解关于 的线性代数方程组。 取 作为v ，将 代 入变分方程，得到 满足的方程组： 有限元法广泛应用的原因 Ritz-Galerkin方法应用的困难 ① 基函数选取必须满足强加边界条件，因此 选取困难； ② 计算量、存储量巨大； ③ 方程组求解病态严重。 充分发挥了变分形式和Ritz-Galerkin方法的 优点； ② 摆脱了传统的基函数取法； ③ 各种问题的结构程序格式统一。 有限元方法基于变分原理， 又具有差分方法的一些特点，并且适于较复杂的区域和不同粗细的网格。 二 椭圆型方程的有限元方法 差分法解偏微分方程，解得的结果就是准确解u在节点上的近似值； Ritz-Galerkin方法得到近似的解析解，但对一般区域，却往往难以实现。 有限元方法与传统Ritz-Galerkin方法的差别在于有限维函数空间的构造方法。Ritz-Galerkin方法选用的基函数在整个定解区域上整体光滑，有限元则取分段或分片连续且局部非零的基函数。 考虑两点边值问题： 1 一维问题的线性元 将区间[a,b]分割为n个子区间 。 第i个单元记为 ，其长度