径向基函数的无网格Galerkin方法.docVIP

径向基函数的无网格Galerkin方法 摘要：首先我们把径向基函数的理论应用到了Galerkin方法解偏微分方程的领域中。在给了一个总的描述之后，我们展示了光滑问题在任意维当中的收敛性并作出了误差估计。 1引言 径向基函数插值在多元近似理论中已经成为了一个强有力的工具，特别对于紧支撑径向基函数出现后。这篇文章我们描述了径向基函数怎样被用来求解椭圆型偏微分方程的数值解。这里我们选择了同古典的有限元方法之中相同的Galerkin方法，得到的结果是可以同古典有限元方法比较的。 与有限元方法相比，使用径向基函数建立有限维子空间的结果与当前子空间的维数没有关系，那么原则上它能解决量子力学中的高维问题。其次，古典有限元方法关于网格的技术细节要花费很多时间，尤其是对于运动边界随时间变化的问题。网格的形成不仅要适应解的奇异性，同时要适应域的改变。无网格方法不需要处理像这样的问题，因为它们仅仅使用了无关的离散中心。最后，光滑解能同非光滑解一样简单的被建立。 第二部分描述了更多偏微分方程的细节及Galerkin方法。第三部分简单的概括了径向基函数插值理论。在第四部分我们展示了这个理论怎样被应用到Galerkin法和Rayleigh-Hitz法近似中来，并且在Sobolev空间得到了一种特殊的基函数。最后一部分，我们把这些结果推广到更一般的基函数中，即使我们在不知道精确解光滑性的情况下，也给出它的逼近阶。 2 PDE和Galerkin法 在有界域及其上考虑如下问题： （2.1） （2.2） 其中，，， 为边界上的单元标准向量。假设矩阵在上满足椭圆条件也就是说存在常数对所有都有 求（2.1）式的变分方程 任取，，乘以（2.1）两端并在上积分： （2.1.1） 对上式第一项由格林公式可得： （2.1.2） 由（2.2）式得： 则 我们可以得到连续的双线性形式： （2.3） 和连续的线性形式： 那么（2.1）~（2.2）的变分问题就是，求函数对于任意都有 （2.4） 成立，由Lax-Milgram定理知函数是唯一的。 这种方法可以在整个Sobolev空间上进行，边界条件本身融入了双线性形式和线性形式F。 为了得到（2.4）的数值解，构造解空间的有限维子空间，并离散化 则有，求函数对于任意都有 （2.5） (2.4)式得真解与数值解之间的误差是收敛的，可由Cea引理得 （2.6） C为一常数。 我们可以要求并且，为当前空间维数，事实上，根据Sobolev空间的嵌入定理，这就要求被逼近的函数至少是一个连续函数。 3 径向基函数 本文用有限维空间，，有如下形式： 函数至少是一阶连续的可导函数，为次数小于m的多项式空间，是由两两不同的点集组成。当为紧支撑函数且m=0时，这种情况下的刚度矩阵 是稀疏的，而且，且有一个单变量函数，刚度矩阵的很多项可以很容易的计算出来。 下面我们通过（2.6）式考虑逼近误差，因此我们需要调用径向基函数的一些理论。 定义3.1 函数是m阶条件正定的当且仅当对于任意互异的点集以及所满足的条件 不等式 是成立的，当m=0，称为正定函数。 如果知道u在一个两两互异的离散点上的值，其径向基函数插值就是构造如下的插值函数 P为次数小于m的多项式，通过插值，满足。因为条件，额外的自由度是有界的，p是多项式空间的基。这个表明始终存在一个满足所需的条件。 径向基函数 名称 m 薄板样条 Sobolev样条 MacDonald函数 0 紧支撑函数， p为多项式 0