使用增强直接代数方法的具有适形分数阶导数的近轴波动方程的精确行波解

《Symmetry》：Exact Traveling Wave Solutions of the Paraxial Wave Equation with Conformable Fractional Derivative Using an Enhanced Direct Algebraic Method

　　

本文使用增强直接代数方法（Enhanced Direct Algebraic Method, EDAM）研究具有适形分数阶导数（Conformable Fractional Derivative）的近轴波动方程（Paraxial Wave Equation）的精确行波解。该方法因其算法简单、计算效率高以及在统一代数框架内生成多种解类型的独特能力而被选中。首先，通过行波变换将具有适形分数阶导数的近轴波动方程转化为常微分方程。然后，系统应用EDAM构建精确的行波解，包括扭曲孤子、钟形孤子、奇异孤子、魏尔斯特拉斯椭圆函数解（Weierstrass elliptic function solutions）和雅可比椭圆函数解（Jacobi elliptic function solutions），展示了该方法在处理复杂非线性结构方面相比标准展开技术的优越性。最后，基于Matlab软件绘制部分解的三维和二维图形，这些图形直观地展示了分数阶参数α对波形局域化和振幅的调节作用。具体而言，随着α的增加，孤子局域化减弱，波包展宽，为光纤中的色散管理提供了洞见。本文所用方法操作简单，解类型丰富，为研究分数阶非线性偏微分方程提供了一种有效途径。

论文《Exact Traveling Wave Solutions of the Paraxial Wave Equation with Conformable Fractional Derivative Using an Enhanced Direct Algebraic Method》发表在《Symmetry》期刊上。该研究聚焦于具有适形分数阶导数的近轴波动方程（Paraxial Wave Equation with Conformable Fractional Derivative）的精确行波解求解问题。分数阶导数（Fractional Derivative）概念起源于1695年，历经三个多世纪发展，已成为数学物理的重要分支。与整数阶模型相比，分数阶模型能更精确地描述复杂系统中的记忆效应、遗传特性和非局部响应。在非线性光学中，光脉冲在光纤、液晶等介质中的传播受衍射、色散和Kerr非线性效应影响，通常由分数阶非线性薛定谔方程描述。然而，分数阶导数存在多种定义（如Riemann–Liouville、Caputo、适形分数阶导数及M-截断导数），各具不同数学性质和物理解释。尽管对适形分数阶导数是否属于“真正”分数阶算子存在争议（因其局部性和满足标准链式法则），但近年来该导数被广泛用于构建各种分数阶非线性偏微分方程的行波解，且计算优势显著：链式法则、乘积法则和商法则均成立，使其特别适合代数求解方法。目前已有多种解析方法，如平面动力系统分析、-展开法、神经网络法、辅助方程法等，但多数研究仅获得一两种特定解（如亮/暗孤子、周期波解），缺乏能同时生成局域孤波、周期椭圆函数和高度局域化魏尔斯特拉斯椭圆结构等综合解谱的统一分析框架。此外，分数阶参数α对不同波形结构的系统调节机制仍有待阐明。基于此，研究人员采用增强直接代数方法（Enhanced Direct Algebraic Method, EDAM）对上述方程展开系统研究。

研究人员开展的研究内容如下：首先，通过行波变换将具有适形分数阶导数的近轴波动方程转化为常微分方程；然后，应用EDAM构建精确行波解，包括钟形孤子、扭结孤子、奇异孤子、雅可比椭圆函数解（Jacobi elliptic function solutions）和魏尔斯特拉斯椭圆函数解（Weierstrass elliptic function solutions）；最后，利用Matlab软件绘制部分解的三维和二维图形，直观展示分数阶参数α对波形局域化和振幅的调节作用。研究得出以下结论：随着α增大，孤子局域化减弱，波包展宽，对应光纤中增强的色散效应；对于周期解，α有效调节振荡频率和空间周期性。该研究的重要意义在于，EDAM首次在统一框架内为适形分数阶近轴波动方程生成多种解类型，弥补了先前研究仅获得孤立解类型的不足，并为其他分数阶非线性系统提供了高效通用方法。

为开展研究，作者使用了以下关键技术方法：增强直接代数方法（EDAM），其核心步骤包括：①通过行波变换将分数阶非线性偏微分方程（FNL PDE）简化为常微分方程（ODE）；②通过平衡最高阶导数项与非线性项确定平衡数N；③假设解结构为有理展开形式，其中辅助函数φ(ξ)满足一阶可积常微分方程（ODE）；④将假设解与辅助方程代入简化后的ODE，收集相同幂次项并令系数为零，得到超定代数方程组；⑤求解该方程组得到未知参数，再根据辅助方程已知解的分类构造各种精确行波解。该方法无需求解复杂辅助方程，仅需简单代数运算，且能自然适应适形分数阶导数的链式法则。

研究结果按原文小标题总结如下：

### 3.1 Case 1:
通过求解超定方程组得到参数关系，导出钟形孤子解φ_1,1和奇异孤子解φ_1,2，分别由sech函数和csch函数表示。钟形孤子描述Kerr介质中色散与非线性平衡下的局域波包，奇异孤子表征奇点附近的极端局域化现象。分数阶参数α控制脉冲包络陡峭度：α越小局域化越尖锐，α越大结构越平滑且非局域化。

### 3.2 Case 2:
Result 1：得到扭结波解φ_2,1和奇异孤子解φ_2,2，由tanh函数和coth函数表示。扭结波描述两个不同渐近态之间的平滑过渡，常见于等离子体物理和非线性光学中的畴壁或激波。α调节扭结过渡层陡峭度：α越大过渡层越宽越平滑，α越小过渡越尖锐。Result 2：得到组合解φ_2,3，由kink和奇异孤子混合构成。

### 3.3 Case 3:
Result 1：得到雅可比椭圆函数解φ_3,1（cn型）、φ_3,2（sn型）和φ_3,3（dn型）。当模数n→1时，这些解分别退化为Case 1和Case 2中的孤子解。Result 2：得到复合解φ_3,4（cn型与倒数组合）、φ_3,5（sn型与倒数组合）和φ_3,6（dn型与倒数组合）。雅可比椭圆函数解描述周期介质中的非线性波模式，如光子晶体或光晶格中的玻色-爱因斯坦凝聚。模数n控制非线性程度（n→0恢复三角函数线性波，n→1恢复局域孤子结构），α调制空间周期和振荡频率。

### 3.4 Case 4:
Result 1：得到魏尔斯特拉斯椭圆函数解φ_4,1和φ_4,2，分别对应不同参数条件。Result 2：得到另一组魏尔斯特拉斯解φ_4,3和φ_4,4。这些解描述高度局域化的周期结构，具有尖锐脉冲序列，适用于深水怪波或非线性光纤中的极端光脉冲等极端局域化现象。其双周期特性可产生复杂时空模式。

### 3.5 Case 5:
得到多种复合解φ_5,1至φ_5,4，由双曲函数与三角函数的组合构成，涵盖不同参数条件。

数值模拟部分（Section 4）通过Matlab绘制三维和二维图形，验证了上述解的存在性，并定量分析了α的影响。钟形孤子（图1）随α增大，峰值降低、波包展宽；扭结孤子（图2）过渡区间变宽、振幅减小；雅可比cn型解（图3）的周期性受α调节；sn型解（图4）在整数阶时展示稳定周期波；复合dn型解（图5）中α增大使奇异峰振幅降低、背景振荡更规则；魏尔斯特拉斯解（图6）呈现高度局域化尖锐脉冲结构；复合背景孤子（图7）在α变化下振幅缓慢演变。残差验证（表1）表明三种代表性解的最大残差均处于10^-2量级，确认解满足原方程，微小残差源于适形分数阶导数数值计算中的有限差分离散误差。EDAM的计算效率高，整个求解过程仅需求解7个多项式方程，无需数值迭代，CPU时间约3-5秒。

总结讨论部分：研究结论指出，本文成功将EDAM应用于具有适形分数阶导数的分数阶近轴波动方程，系统构建了多种精确行波解。与已有方法（如-模型展开法、-展开法、统一技术）相比，EDAM能在统一框架内同时生成所有五种解类型，其中雅可比椭圆、魏尔斯特拉斯椭圆及组合孤子解为首次获得，体现了方法的优越性和通用性。数值模拟系统分析了α对波形演化的影响：α增大削弱波包局域化、提高传播稳定性；对于周期解，α有效调制振荡频率和空间周期性。该研究不仅展示了EDAM求解复杂分数阶非线性偏微分方程（FNL PDE）的强大能力，还为更一般的分数阶非线性偏微分方程的行波解研究提供了有价值的理论参考。未来工作将致力于将该方法推广至高维分数阶系统，并探索符号计算与机器学习相结合以发现更复杂FNL PDE的行波解。