[论文解读] Numerical solution of one-dimensional Sine--Gordon equation using Reproducing Kernel Hilbert Space Method

Ali Akgül, Inc, Mustafa|arXiv (Cornell University)|Apr 2, 2013

Numerical methods in engineering参考文献 28被引用 3

一句话总结

本文提出了一种再生核希尔伯特空间方法（RKHSM）来求解具初值和边界条件的一维 sine-Gordon 方程。该方法在再生核希尔伯特空间中将精确解表示为一个收敛级数，避免了空间和时间的离散化，从而实现了高精度的数值解，且计算量极小。数值例子验证了方法的均匀收敛性，误差低于 10^{-6}。

ABSTRACT

In this paper, we propose a reproducing kernel Hilbert space method (RKHSM) for solving the sine--Gordon (SG) equation with initial and boundary conditions based on the reproducing kernel theory. Its exact solution is represented in the form of series in the reproducing kernel Hilbert space. Some numerical examples have been studied to demonstrate the accuracy of the present method. The results obtained from the method are compared with the exact solutions and the earlier works. Results of numerical examples show that the presented method is simple and effective.

研究动机与目标

• 开发一种无网格的数值方法，用于求解具初值和边界条件的一维 sine-Gordon 方程。
• 利用再生核希尔伯特空间理论，提供精确解的收敛级数表示。
• 通过数值例子展示该方法的高精度和计算效率。
• 将 RKHSM 的结果与精确解及其他现有数值方法进行比较，表明其性能更优或具有竞争力。

提出的方法

• 该方法利用再生核希尔伯特空间理论构建解空间，使精确解以无穷级数形式表示。
• 在再生核空间中定义一个有界线性算子，将原偏微分方程转化为可解的算子方程。
• 通过截断核函数的级数迭代逼近解，确保近似解在再生核空间中一致收敛于精确解。
• 该方法避免了空间和时间的离散化，无需求解大型线性系统或更新雅可比矩阵。
• 算法直接从核函数计算近似解，显著降低计算成本并避免数值不稳定性。
• 收敛性分析证明，近似解在再生核空间中一致收敛于精确解。

实验结果

研究问题

• RQ1 再生核希尔伯特空间方法（RKHSM）是否能在不进行离散化的情况下，有效求解具初值和边界条件的一维 sine-Gordon 方程？
• RQ2 RKHSM 的精度与精确解及其他数值方法（如径向基函数配点法或有限差分格式）相比如何？
• RQ3 RKHSM 近似解的收敛行为如何？其是否一致收敛于精确解？
• RQ4 与需要求解大型线性系统的牛顿迭代法相比，RKHSM 的计算效率如何？
• RQ5 RKHSM 是否能够以高精度应用于 sine-Gordon 方程的线性和非线性形式？

主要发现

• RKHSM 实现了高精度，测试实例中绝对误差始终低于 5 × 10^{-7}，最佳情况下相对误差低于 3.6 × 10^{-7}。
• 在示例 5.1 中，最大绝对误差为 1.756 × 10^{-7}，相对误差为 2.987 × 10^{-7}，显示出极高的精度。
• 在示例 5.2 中，由于对称性精确且数值稳定性高，对称点（x = 0, ±0.4, ±0.8）处的绝对误差为零。
• 每个示例的 CPU 时间在 0.546 至 1.719 秒之间，表明尽管精度极高，计算成本仍很低。
• 该方法表现出近似解对精确解的一致收敛性，验证了理论收敛结果。
• 与文献 [37] 的结果比较表明，RKHSM 在所有测试点均产生了更小的绝对误差和相对误差，证实了其更高的精度。