偏微分方程数值解.ppt

②变分问题与前述二阶椭圆边值问题相应的变分问题是求，使其中③变分原理（变分问题与边值问题的等价性）对第一边值问题，无论齐次或非齐次边界条件，泛函是一样的，只是边界条件要作为强加边值条件加在所取的函数类上。设，是二阶椭圆边值问题的解，则使J(u)达到极小值；反之，若使J(u)达到极小值，则是二阶椭圆边值问题的解。其中对第二、三类边值问题，无论齐次或非齐次边界条件，二次泛函形式相对于第一边值问题有所改变，但函数类的选取与边界条件无关。（2）虚功原理问题其中设，以v乘方程两端后在G上积分，并利用Green公式，得变分方程虚功原理在力学里，表示虚功设是边值问题的解，则对任意，满足变分方程。反之，若，且对任意满足变分方程，则为边值问题的解。与极小位能原理类似，第一类边界条件为强加边界条件，第二、三类边界条件为自然边界条件。虚功原理比极小位能原理应用更广。目的：求解相应的变分问题或相应的变分方程。Ritz方法是近似求解变分问题(即二次泛函极小值)的算法。Galerkin方法是近似求解变分方程的算法，这两种算法统称为Ritz-Galerkin方法。Ritz-Galerkin方法的基本思想以下用V表示等Sobolev空间，L表示微分算子，?(u,v)为由L及边值条件决定的双线性泛函。4Ritz-Galerkin方法用有限维空间的函数代替变分问题(或变分方程)中无限维空间的函数，从而在有限维函数空间中求变分问题(或变分方程)的近似解，并要求当有限维空间的维数不断增加时，有限维近似解逼近原变分问题(或变分方程)的解。由极小位能原理得出的变分问题为：Ritz方法：求变分问题的近似解。（1）Ritz方法求，使其中，设是V的n维子空间，是的一组基底(称为基函数)。中任一元素可表示为即选择适当的，使取极小值。求，使Ritz方法：展开令则满足解出代入，则得Ritz方法步骤为：根据最小位能原理构造相应于微分方程或物理问题的变分问题；取作为的一组基底，即用近似代替无穷维空间V；根据二次函数取极值的必要条件，得到中所满足的方程组：求解关于的线性代数方程组。由虚功原理得出的变分方程为：Galerkin方法：求变分方程的近似解。（2）Galerkin方法设是V的n维子空间，是的一组基底(称为基函数)。中任一元素可表示为即选择适当的，使取极小值。Galerkin方法：求，使对，满足(2).跳点格式首先将网格点（j,n）按j+n等于偶数或奇数分成两组，分别称为偶数网点和奇数网点。从到的计算过程中，先在偶数网格点上用古典显式格式计算，再在奇数网点上用古典隐格式计算，即：三双曲型方程（a）一阶常系数线性双曲型方程（b）二阶常系数线性双曲型方程（波动方程）其中a为常数主要对象为这些方程的定解条件，可仅有初始条件，也可以有初始