重心插值配点法及其应用.docVIP

重心插值配点法及其应用摘要：重心Lagrange插值具有数值稳定性好、计算精度高的优点。采用重心Lagrange插值近似未知函数，建立未知函数各阶导数的微分矩阵。采用微分矩阵近似未知函数的导数，利用配点法将控制方程和边界条件离散为代数方程组，通过求解代数方程组，从而可求解偏微分方程。数值算例表明，重心插值配点法具有原理简单，易于程序实现和数值计算精度高的优点。 关键词：重心Lagrange插值；微分矩阵；配点法 中图分类号：0241文献标识码：A Barycentric Lagrange interpolation collocation Method and its Application Abstract：Barycentric Lagrange interpolation collocation method has excellent numerical stability and high accuracy. this paper presents the Barycentric Lagrange interpolation collocation method to get the differentiation Matrix of unknown function. So the control equation can be expressed as linear systems by the collocation method. According to those formulas, differential equations can be soluted. The principle of this method is simple and easy to programming. The accuracy and the numerical stability are very excellent. Key words: Barycentric Lagrange interpolation, differentiation Matrix, collocation method, 中图分类号：O241.3 文献标识码：A文章编号：2095-2104（2012）01-0020-02 0 引言 具有初、边值条件的常、偏微分方程的解析解常无法通过理论推导获得，一种有效的途径是采用数值求解方法获得具有一定数值精度的近似解，这些数值方法包括：有限差分法、有限单元法、边界单元法、无网格法及一些杂交使用的方法等数值求解方法。其中，有限差分法、有限单元法这两种方法要对求解区域划分单元，计算精度依赖于单元的大小。采用配点法求解边值问题不需要划分单元，公式简单，不需要积分，易于编程。目前用于求解工程中的常微分方程边值问题的配点法主要有拟谱法和微分求积法。拟谱法是根据谱方法发展出来的一种方法，虽然这种方法的理论研究已经有进一步的发展，但是工程技术人员对这种方法不是很了解。微分求积法的基本原理是将未知函数在区间上所有离散点的函数值的加权和来逼近该函数在某一离散点的偏导数或者积分，这种方 法中的权系数的确定通常是根据Lagrange多项式在网点处的导数值给出。但是这种方法的局限性是离散点不能取得太多，否则Lagrange多项式表示的曲线随多项式次数的升高而出现Runge现象，从而产生计算的不稳定性。重心Lagrange插值具有极好的数值稳定性和极高的近似精度，同时重心Lagrange插值公式具有紧凑的各阶导数的计算公式。 本文所采用的重心插值配点法就是用重心Lagrange插值多项式求出某一函数在各个离散点的微分矩阵，从而可以通过矩阵的运算来求出微分方程的解。而许多结构力学问题的求解最终也归结为在一定的边界条件和初始条件下的(偏)微分方程(组)的求解，所以，可以把重心插值配点法作为用于结构力学问题的求解的一种数值方法。本文把重心插值配点法用于矩形薄板的挠度问题的分析中，精度较高，得到了满意的效果。 1微分矩阵 设为一待求函数，它在区间[0，L]上连续可微，现沿轴设置个节点，并以节点函数值作为基本未知量，且令=。在全域内采用重心Lagrange高阶多项式插值逼近，则有： （1） 其中，为重心权， （2） 令：（3） 因为所求解的问题是线性的，因此可以定义一个的矩阵。则（3）式可以写成如下的形式： （4） 其中为未知函数一阶导数值构成的列向量，为节点函数值的列向量。 就是对求一次导数的微分矩阵。由（1）、（2）、（3）式可得： （5） （6） 同理，可以推得求次导数的微分矩阵。 （7）