[论文解读] An application of Bell polynomials in numerical solving of nonlinear differential equations

Josef Rebenda|arXiv (Cornell University)|Jan 27, 2019

Iterative Methods for Nonlinear Equations被引用 4

一句话总结

本文提出了一种偏普通贝尔多项式的新应用，通过微分变换法简化非线性微分方程的数值求解。通过使用贝尔多项式重构法卡·布鲁诺公式，该方法实现了无需符号导数即可高效计算非线性项，仅通过递归算术运算即可获得精确的级数解，如在两个初值问题中所示，通过模式识别确定了解析解。

ABSTRACT

Partial ordinary Bell polynomials are used to formulate and prove a version of the Fa\\`{a} di Bruno's formula which is convenient for handling nonlinear terms in the differential transformation. Applicability of the result is shown in two examples of solving the initial value problem for differential equations which are nonlinear with respect to the dependent variable.

研究动机与目标

• 为解决微分变换法中处理涉及未知解复合函数的非线性项时的计算复杂性问题。
• 开发一种系统且计算高效的非线性项变换方法，利用贝尔多项式的组合性质。
• 为现有半解析方法（如阿多米安分解法）提供一种实用替代方案，后者需要符号导数或复杂的积分计算。
• 证明所提出的方法可通过递归算术运算生成精确的级数解，并具备识别闭式解的潜力。

提出的方法

• 本文引入了基于偏普通贝尔多项式的法卡·布鲁诺公式改进版本，用于表达复合函数的高阶导数。
• 将非线性函数 $ f(u(t)) $ 的微分变换表述为涉及变换系数 $ U(k) $ 的贝尔多项式之和，从而避免对 $ f $ 进行符号导数运算。
• 该方法依赖于通过先前计算的 $ U(i) $（$ i < k $）递归计算 $ U(k) $，确保计算效率。
• 该算法使用逆微分变换将解重构为类似泰勒级数的形式，其系数由贝尔多项式展开导出。
• 该方法避免了数值积分、微分或迭代求解器，仅依赖基本算术运算。
• 该方法在两个初值问题上进行了验证，通过可识别的系数模式显示出收敛于精确解的特性。

实验结果

研究问题

• RQ1 能否利用偏贝尔多项式在无需符号微分的情况下，简化常微分方程中非线性项的微分变换？
• RQ2 如何利用贝尔多项式重构法卡·布鲁诺公式，以提高半解析方法中的计算效率？
• RQ3 当该方法应用于非线性初值问题时，其在多大程度上能产生闭式精确解？
• RQ4 与阿多米安分解法或同伦摄动法等现有方法相比，该方法在所需运算量方面具有何种计算优势？

主要发现

• 非线性项 $ f(u(t)) $ 的微分变换成功地通过偏普通贝尔多项式表达，使得可直接从变换系数 $ U(k) $ 计算，无需符号导数。
• 在第一个例子中，解的系数遵循模式 $ U(k)[0] = 0 $（当 $ k $ 为偶数时），以及 $ U(k)[0] = (-1)^{(k-1)/2} \frac{2^k}{k!} $（当 $ k $ 为奇数时），从而导出精确解 $ u(t) = \sin 2t $。
• 在第二个例子中，该方法生成的系数与 $ \sin 2t $ 的泰勒级数完全匹配，通过模式识别确认了解的精确性。
• 该算法计算高效，仅依赖算术运算和递归计算，避免了积分、微分或初始猜测的需求。
• 当系数模式与已知函数匹配时，该方法可识别出闭式精确解，其精度超越标准数值近似。
• 该方法具有通用性，可推广至其他问题（包括两点边值问题），并提供一种保持展开点附近解析性的函数解析近似。