波动问题的高精度重心有理插值配点法

提出一种数值求解波动问题的高精度重心有理插值配点法。对于给定的时间和空间上的计算节点,采用重心有理插值近似未知函数,建立未知函数关于时间和空间变量导数的微分矩阵。将未知函数的重心有理插值近似函数代入波动问题的控制方程,得到波动问题方程和定解条件的离散代数方程组。利用微分矩阵的记号,将离散后的代数方程组写成简洁的矩阵形式。通过置换法施加边界条件和初始条件,求解代数方程组,得到波动问题在计算节点处的位移值。数值算例表明,重心有理插值配点法具有计算公式简单、计算节点适应性好、程序实施方便和计算精度高的优点。     

极坐标下薄板弯曲问题的重心有理插值法

利用重心有理插值配点法(BRICM)研究了极坐标下薄板的弯曲问题,该方法是以重心有理插值近似未知函数强迫微分方程在离散节点处成立,得到微分方程的离散代数方程组,进而采用重心有理插值的微分矩阵将离散代数方程组表达为矩阵的形式。利用置换法施加边界条件,求解微分方程组。数值算例结果表明,该方法在解决极坐标下薄板弯曲问题上公式简单,程序实施方便且计算精度高。     

极坐标系下弹性问题的重心插值配点法

针对岩土工程中的孔洞及曲梁问题,提出一种在极坐标系下求解二维弹性问题的重心插值配点法.该方法分别在r和θ方向分别布置m和n个节点,生成求解区域上的节点.以一维重心Lagrange插值的张量积插值形式近似二维弹性问题的位移函数,代入位移表达的平衡方程和边界条件,平衡方程和边界条件分别在所有的计算节点和边界节点上精确成立,得到极坐标下弹性力学平衡方程和边界条件的离散代数表达式.利用一维重心Lagrange插值微分矩阵,将离散的平衡方程和边界条件表达为矩阵形式.利用置换法施加边界条件,求得在计算节点处的位移,进而通过微分矩阵直接求得计算节点处的应力.数值算例表明:极坐标下重心插值配点法具有计算格式简单、程序实施容易和计算精度高的特点.     

重心有理插值配点法分析压杆稳定问题

采用重心有理插值近似未知函数,得到未知函数的各阶微分矩阵．利用微分矩阵将压杆控制微分方程离散为代数方程组．将离散的边界条件采用置换法施加到代数方程组中,得到关于压杆屈曲载荷的特征方程．求解特征值问题得到压杆的屈曲载荷．算例表明,重心有理插值配点法具有数值稳定性好、计算精度极高,程序实施容易等优点．     

求解非线性伪抛物方程的重心Lagrange插值配点法

提出了重心Lagrange插值配点法求解一类非线性伪抛物方程。首先,介绍了重心Lagrange插值并给出了微分矩阵表达式。其次,构造了求解非线性伪抛物方程的直接线性化迭代格式、部分线性化迭代格式、Newton线性化迭代格式。再次,未知函数和初边值条件利用重心Lagrange插值函数来近似,利用配点法得到离散方程,获得了方程的矩阵表达式。最后,数值算例表明,重心Lagrange插值配点法具有高精度和高效率的优点。     

重心插值配点法求解Allen-Cahn方程

利用重心插值配点法(重心Lagrange插值配点法和重心有理插值配点法)构造包含时间、空间变量的近似函数,给定Chebyshev点族;将重心插值配点法代入Allen-Cahn方程及定解条件,得到离散方程组,并采用Newton迭代格式求解方程组.数值算例表明:文中的配点法具有较高精度;利用2种配点法的能量函数满足能量递减规律.     

重心插值配点法求解Volterra积分方程

文章提出了求解Volterra积分方程的一种高精度数值方法：重心插值配点法（包括重心Lagrange插值配点法和重心有理插值配点法）。该方法分为两步：首先对Volterra积分方程采用两种重心插值配点法进行离散，构造出Volterra积分方程的数值求解格式；然后，依次选取第二类Chebyshev节点和等距节点进行数值计算。文章主要研究积分项中含有未知函数的一阶导函数的Volterra积分方程的离散格式构造及数值实现。数值实验结果表明：在使用第二类Chebyshev节点时，用重心Lagrange插值配点法较好；在使用等距节点时，使用重心有理插值配点法较好。     

时间分数阶Black-Scholes方程的重心Lagrange插值配点法

针对欧式期权定价的时间分数阶Black-Scholes模型,设计一种重心Lagrange插值配点法格式.首先,采用Laplace变换近似Caputo型分数阶导数,将分数阶方程转化为整数阶方程;然后,在时-空方向上均采用重心Lagrange插值配点法进行离散,构造重心Lagrange插值配点法格式.结果表明：时间分数阶Black-Scholes方程的重心Lagrange插值配点法具有高精度和有效性.     

基于重心型插值的数值计算方法

对以重心型插值作为近似函数,数值求解微分方程初值问题和边值问题的数值计算方法作了介绍。给出了重心Lagrange插值和重心有理插值的计算公式和插值节点类型。归纳了求解微分方程初边值问题的重心插值配点法、重心插值Galerkin法和重心插值单元法的计算公式、边界条件/初始条件的离散和施加方法。     

重心插值配点法求解二维非线性椭圆型方程

首先利用重心插值配点法离散二维非线性椭圆型方程和边界条件,其次采用完全线性化迭代和Newton-Raphson迭代求出方程的近似解.实验结果表明:重心插值配点法理论简单,计算精度高; Newton-Raphson迭代法无论是在计算效率上,还是在计算精度上,都优于完全线性化迭代.     

Legendre-Gauss-Lobatto节点的一个注记

以Legendre-Gauss-Lobatto点为节点的Lagrange插值基函数,构造N阶插值多项式P_N(x)。对P_N(x)分别求一阶和二阶导数,得到一阶和二阶微分矩阵。利用Legendre-Gauss-Lobatto点的性质导出一阶和二阶微分矩阵的关系,由此可利用Lagrange插值多项式数值求解微分方程。     

重心插值配点法求解Cahn-Hilliard方程

对Cahn-Hilliard方程中的时、空方向均采用重心插值配点格式(重心Lagrange插值配点格式和重心有理插值配点格式)进行离散,非线性项采用一般迭代法,导出离散的线性代数方程组,并给出重心Lagrange插值的逼近误差估计.数值算例表明:两种重心插值配点格式均具有高精度,且满足能量递减规律.     

三元Lagrange插值正则性问题研究

针对当今许多科研领域中(如曲面拼接、散乱数据插值与拟合等)经常涉及到的三元Lagrange插值问题进行了研究。提出了沿空间代数曲线插值的基本概念,同时通过使用代数几何中的若干理论,得到了构造沿空间代数曲线及代数曲面插值正则结点组的迭加构造方法,该方法推广了文献[1-2]中的某些主要研究结果。     

球面上Lagrange插值问题研究

以代数几何中某些理论方法为工具,对球面上Lagrange插值问题进行了研究和探讨.将文献[1]中所给出构造关于球面插值适定结点组的添加平面法推广到了添加圆锥曲面的情形,该方法是以迭加过程来实现的,因此便于在计算机上实现其构造过程.     

SAV方法求解Allen-Cahn方程的数值比较

该文研究基于标量辅助变量(SAV)格式下求解Allen-Cahn方程的数值比较.首先给出1维Allen-Cahn方程的SAV格式; 然后,对方程的时间方向采用2阶向后差分(BDF2)格式和Crank-Nicolson(CN)格式离散,对方程的空间方向采用重心Lagrange插值配点法和2阶中心差分法离散,用离散正弦变换(DST)、快速傅里叶变换(FFT)求解差分导出的线性代数方程组; 最后,通过数值算例验证重心Lagrange插值配点法是指数收敛,与差分格式比较,配点格式用较少的点就能达到较高的精度且耗时少,并进一步验证几种SAV离散格式都满足能量递减规律.     

构造矩阵有理插值函数的方法

熟知的构造矩阵值有理插值函数的方法,是基于矩阵的古典逆或Samelson逆,利用连分式给出的,其算法可行性不易预知。借助构造向量值有理插值的方法,引入多个参数,定义一对多项式:代数多项式和矩阵值多项式,并利用两多项式相等的充分必要条件,通过求解方程组确定参数,并由此给出类似于多项式插值的矩阵值有理插值公式;该公式简单,便于实际应用。