三阶两点边值问题的数值解与重心有理插值的搭配方法

文摘

数值解的一种三阶边值问题进行了探讨。重心有理插值的搭配方法,矩阵形式的三阶两点边值问题,得到了收敛性和误差分析。此外,据报道,一些数值例子证实了理论分析。

1。介绍

微分方程可以充分发挥数学不同学科的优势。微分方程的理论与实际问题相结合可以建立实际问题的模型。许多工程和物理问题可以转化为微分方程的初边值问题。在这些问题中,只有几个简单的案例分析可以解决,最需要解决的工程问题数值方法。与多项式插值相比,有理函数插值具有较高的插值精度,可以有效地克服插值的不稳定(1- - - - - -4]。重心有理插值不仅具有特殊的分布式节点上插值精度高但也有等距节点的插值精度高(5- - - - - -7]。这种方法已被用于解决某些问题,如沃尔泰拉积分方程(2,8,9),延迟沃尔泰拉积分微分的方程(10,11[],平面弹性问题12[],非线性问题13],热传导方程[14],等等15- - - - - -17]。

三阶微分方程已经广泛应用在很多科学领域和重要的理论价值,应用数学和物理学等。因此,三阶边值问题已经受到了许多学者的广泛关注(18- - - - - -20.]。在本文中,我们考虑三阶两点边值问题的数值解, 重心有理插值的搭配方法。

重心有理插值搭配方法意味着使用重心插值多项式找到函数的微分矩阵每个离散点;因此,微分方程的解可以通过矩阵运算。重心有理插值具有良好的数值稳定性和近似精度高,和重心有理插值公式有一个紧凑的计算公式的衍生品。因此,搭配重心有理插值方法是一种有效的方法求解微分方程的边值问题。

2。公式的重心插值搭配方法

离散化区间 成统一的部分 ,假设 是一个未知函数的函数值在离散的节点 。

对于任何 , ,插值函数在点吗 ;然后,我们有 ,和 在哪里

通过改变多项式拉格朗日插值形式 并结合(4)和(5)在一起,我们得到的 在哪里 。

然后,我们得到 它的基函数在哪里吗

等距点的权函数

第二类切比雪夫点, 权函数是

由公式(8),米th阶的导数在节点 可以表示为 然后(12)可以写成矩阵形式 在哪里 和 。

通过使用重心插值函数 方程(1)可以用数值形式

通过使用的符号微分矩阵,(15)也可以表示为 或者简单的矩阵形式

边界条件(2)可分为 在哪里 ,

3所示。收敛性和误差分析

在本节中,我们将考虑误差的等距问题插值节点:

让的解决方案(1);对于任何 ,假设 ,重心插值函数在点 ;然后,我们有 ,和 在哪里

然后,误差函数被定义为 和 在哪里 。

以数值的形式, 并结合(24)和(1),我们得到 在哪里 。

引理1。为中定义的(23),我们有

让的解决方案(1),是数值解;然后,我们有 和

结果可以在[1]。

根据引理1,我们可以得到下面的定理。

定理1。让 , ,和 ;然后,我们有

证明。让 添加第二列,第三列,…,列n列1, 然后,我们有 与 , ,在哪里
通过 这意味着 在哪里矩阵的元素是什么 ,我们有 完成证明。

4所示。数值例子

作为一个例子,我们考虑的两点边值问题:

对于这个问题,我们可以找到一个函数这样的分析解决方案 在哪里是一个自由选择的参数。

用(38)(36),我们得到

为不同的值和不同数量的节点,我们可以计算相应的相对误差和收敛速度;一些数据如表所示1和2。

在表1,用不同的等距节点的收敛速度是 ;在表2的收敛速度与不同的第二类切比雪夫点是 。

为不同的值和不同数量的节点,我们可以计算相应的相对误差;一些数据如表所示3和4。

从表3和4,我们可以发现,不同的价值观,收敛速度可以达到 对等距和nonequidistant节点。

5。结论

摘要重心合理搭配方法求解三阶两点边值方程,和误差函数的收敛速度也获得了。常系数和变系数的两点边值方程,数值结果表明,收敛速度可以达到对等距节点和第二类切比雪夫点 所以重心合理搭配方法是一种有效的方法。与其他方法相比,这种方法的优点是,矩阵方程可以很容易地获得,程序简单,计算精度高,可以通过使用一些点。

数据可用性

的数据支持本研究的发现可以从相应的作者在合理的请求。

的利益冲突

作者宣称没有利益冲突有关的出版。

确认

工业大学合作的支持合作教育项目(201801123024)。