[论文解读] Quintic Spline Solutions of Fourth Order Boundary-Value Problems

Shahid S. Siddiqi, Ghazala Akram|ArXiv.org|Jun 25, 2003

Advanced Numerical Analysis Techniques参考文献 11被引用 30

一句话总结

本论文提出了一种五次样条方法，用于高精度求解四阶线性两点边值问题。通过推导一致的端点条件，并基于五次样条插值求解线性方程组，该方法在逼近解及其导数时实现了六阶收敛（O(h⁶)），并通过逐步减小步长 h 的数值例子得到验证。

ABSTRACT

In this paper Quintic Spline is defined for the numerical solutions of the fourth order linear special case Boundary Value Problems. End conditions are also derived to complete the definition of spline.The algorithm developed approximates the solutions, and their higher order derivatives of differential equations. Numerical illustrations are tabulated to demonstrate the practical usefulness of method.

研究动机与目标

• 开发一种高阶数值方法，用于求解梁弯曲和弹性地基模型中出现的四阶线性两点边值问题。
• 推导在等距节点上使用五次样条插值的一致端点条件，以确保全局收敛性。
• 在逼近解及其高阶导数时实现六阶精度（O(h⁶)）。
• 提供一种稳定且精确的算法，可同时计算 y(x)、y′(x)、y′′(x)、y′′′(x) 和 y⁗(x) 的近似值。
• 通过逐步减小 h 的数值实验验证该方法，展示收敛速率和误差降低情况。

提出的方法

• 在区间 [a,b] 上定义五次样条 Q(x)，节点为 xi = a + ih，使用值 yi = Q(xi)，mi = Q′(xi)，Mi = Q′′(xi)，ni = Q′′′(xi)，Ni = Q⁗(xi)。
• 通过四次积分四阶导数构造样条，在每个区间 [xi−1, xi] 上得到一个五次多项式，其常数由连续性条件和边界条件确定。
• 利用恒等式 (2.9)–(2.19) 推导线性方程组，将每个节点处 yi、mi、Mi、ni、Ni 的值关联起来，以确保光滑性和一致性。
• 利用关系式 (2.14)–(2.19) 推导端点条件，以闭合系统并确保六阶收敛。
• 求解大型稀疏线性系统，以计算每个节点处的未知量 mi、Mi、ni、Ni、yi，从而实现对解及其导数的同步逼近。
• 该算法被数值实现，并在三个基准问题上进行测试，逐步减小 h 以验证收敛速率。

实验结果

研究问题

• RQ1 五次样条方法能否在具有通用边界条件的四阶边值问题中实现六阶收敛（O(h⁶)）？
• RQ2 在等距节点上使用五次样条插值时，应采用何种端点条件以确保一致收敛？
• RQ3 该五次样条方法在逼近解 y(x) 及其一阶、二阶、三阶和四阶导数时的精度如何？
• RQ4 当步长 h 减小时，y、y′、y′′、y′′′ 和 y⁗ 的最大绝对误差如何变化？
• RQ5 该方法能否在不同类型四阶边值问题（包括变系数和非线性项）中保持高精度？

主要发现

• 五次样条方法在逼近解及其导数时实现了六阶收敛（O(h⁶)），这通过 h 逐步减小带来的误差降低现象得到证实。
• 在例1中，y(x) 的最大误差从 h=1/8 时的 2.1×10⁻³ 降低到 h=1/1024 时的 1.07×10⁻⁷，表明接近 O(h⁶) 的收敛行为。
• 在例2中，y(x) 的误差从 h=1/8 时的 5.38×10⁻⁴ 降低到 h=1/1024 时的 6.37×10⁻⁸，支持 O(h⁶) 收敛。
• 在例3中，y(x) 的误差从 h=1/8 时的 1.4×10⁻³ 降低到 h=1/1024 时的 3.19×10⁻⁸，再次确认了 O(h⁶) 行为。
• 该方法对四阶以内的所有导数均提供了高度精确的逼近，y⁗(x) 的误差也随 h⁶ 比例减小。
• 推导出的端点条件 (2.14)–(2.19) 对实现完整收敛阶次和确保所有测试问题中的稳定性至关重要。