《Symmetry》：Investigating Bifurcation, Chaos, Multistability, and Localized Interaction Structures in the Nurshuak–Tolknay–Myrzakulov (NTM-III) Equation
编辑推荐：
本文采用动力系统理论与精确解技术对非线性Nurshuak–Tolkunay–Myrkakulov（NTM-III）方程进行了分析。通过适当的约化，将NTM-III模型化为二阶非线性常微分方程，进而写成平面哈密顿动力系统。研究人员利用雅可比矩阵（Jacob
本文采用动力系统理论与精确解技术对非线性Nurshuak–Tolkunay–Myrkakulov（NTM-III）方程进行了分析。通过适当的约化，将NTM-III模型化为二阶非线性常微分方程，进而写成平面哈密顿动力系统。研究人员利用雅可比矩阵（Jacobian matrices）和本征值理论，对平衡点、分岔结构和稳定性性质进行了解析研究。不同的参数组合以各种相图呈现，其中存在鞍点和中心平衡态并标明了其稳定性。为了研究复杂的非线性动力学，研究人员向系统引入外部周期微扰。借助相图、时间序列分析、回归映射（return maps）、李雅普诺夫指数（Lyapunov exponents）、灵敏度分析和多稳定性诊断，对所得强迫动力学模型进行了研究。数值模拟揭示了混沌行为、对初值的强依赖以及同一参数组下多个共存吸引子的现象。此外，研究人员采用多元广义指数有理积分函数（Multivariate Generalized Exponential Rational Integral Function，MGERIF）方法，获得了一组精确解析解，包括三角函数、双曲函数和指数波结构。所得解展现了有趣的非线性波相互作用、多峰子（multi-peakon）形成和局域化传播模式。研究结果表明NTM-III方程中分岔、混沌、多稳定性与非线性波传播之间存在复杂关系，并为该方程的数学与物理性质提供了新的认识。
Nurshuak–Tolknay–Myrzakulov（NTM-III）方程是NTM可积系统家族中的第三类约化形式。该方程描述了两个主要复场和两个辅助场之间的时空耦合动力学，包含了振幅相关势项、对流通量导数、线性时间演化项以及高阶混合色散项，能够刻画多组分非线性波现象，如多峰子形成、弹性孤子碰撞和局域lump型相互作用，在非线性光学、等离子体物理和流体动力学中具有重要的数学框架意义。
非线性微分方程（NDEs）在数学、物理学和应用科学中有着广泛的建模作用。自1967年逆散射变换（IST）的提出为现代孤子理论奠定基础以来，可积非线性系统的研究持续受到关注。经典的孤子模型如Korteweg–de Vries方程、非线性薛定谔方程、sine-Gordon方程和Heisenberg自旋模型已被广泛用于描述局域稳定波结构。近年来，可积系统的多维推广成为研究热点。NTM方程族作为重要的可积框架，其完整可积性和几何解释已通过Lax对形式得以建立。然而，对于NTM-III方程而言，系统性的分岔结构分析、混沌动力学诊断、多稳定性表征以及精确解析解的构造仍亟待深入。为此，研究人员采用非线性动力系统理论和解析求解方法，对NTM-III方程开展了系统性的研究，旨在揭示其分岔、混沌、多稳定性与局域波结构之间的内在关联。
研究人员首先将NTM-III方程通过行波变换约化为非线性常微分方程，再进一步写成平面哈密顿动力系统。通过雅可比矩阵和本征值分析，对平衡点的类型（鞍点、中心或尖点）和稳定性进行了判定。在约化过程中采用了适当的可积约化条件，将NTM系统分别简化为NTM-I、NTM-II和NTM-III三类形式，其中NTM-III在实函数约化下获得。为研究外部扰动对系统的影响，研究人员将周期微扰项引入哈密顿系统，通过相图、时间序列、回归映射、李雅普诺夫指数、多稳定性诊断和灵敏度分析等多种工具系统地表征了系统的复杂动力学行为。该研究主要依赖以下关键技术方法：
（1）行波约化方法：将偏微分方程化为常微分方程；（2）平面动力系统与哈密顿系统理论：构造约化模型的Hamilton函数并分析相图结构；（3）Jacobian矩阵与特征值分析法：用于平衡点的稳定性判定；（4）周期摄动方法：将系统扩展为外力驱动系统以激发混沌；（5）定性混沌诊断工具：包括相图、时间序列分析、回归映射、Lyapunov指数计算和灵敏度分析；（6）MGERIF方法：用于构造精确解析解，涵盖三角函数、双曲函数和指数函数形式。数值模拟均采用四阶Runge–Kutta（RK4）方法以保证结果的准确性和稳定性。研究基于理论建模和数值仿真，没有临床样本队列来源。
研究结果分为以下多个方面：
**分岔结构与稳定性交换。** 研究人员将约化后的NTM-III模型转化为平面哈密顿系统，得到两个平衡点，参数M和N分别控制线性恢复力和非线性二次相互作用。当M>0且N>0时，点表现为鞍点，点表现为中心；当M<0时，二者角色互换。分岔边界为M=0，此时两个平衡点合并且特征值消失，构成两种相平面状态的退化极限情形。相图（图1）展示了鞍-中心构型随参数演化产生的稳定性交换。
**混沌动力学检测。** 在加入外部周期摄动项后，系统（29）呈现丰富的混沌行为。当摄动振幅较大（F=5.2）而频率极低（ω=0.001）时，系统表现出强混沌特征，相轨迹高度离散（图2）；将振幅降至F=0.05且频率升至ω=2.5时，混沌强度减弱，运动相对有组织（图3）；当振幅增至F=0.15时，非线性相互作用增强，混沌重新加剧（图4）；当F=4.5、ω=3.5时，系统进入更深的混沌区域，相图完全无序（图5）。外部摄动的振幅是控制混沌强度的关键参数，增大振幅使不规则振荡增强，而频率变化主要调节振荡结构和转迁模式。
**多稳定性分析。** 在参数完全相同的条件下，选择不同初始条件（如图6中的(0.3,0.2)和(0.1,0.02)）可使系统的长期状态收敛到不同的吸引子。相图和时间序列呈现截然不同的几何结构和振荡特性。图8进一步表明，即使初始条件非常接近（如(0.3,0.1)与(0.2,0.1)），轨迹仍会发散并收敛至不同吸引子。多稳定性的产生与二次恢复项和外部强迫项之间的非线性竞争密切相关，这导致相空间被划分为多个吸引域。
**回归映射分析。** 在四组不同的参数组合下，将连续采样值对绘制成回归映射图形，得到的结果均呈现不规则且非周期的散点分布，而不是收敛到固定点或简单光滑闭合曲线，定性证明了系统的混沌特征。
**Lyapunov指数定量验证。** 对于参数(M=2.2,N=?2,F=3.5,ω=1.8)的系统，最大Lyapunov指数接近正值1.42；对于(M=3.3,N=?2,F=4.8,ω=2.5)的参数组合，最大Lyapunov指数同样接近正值1.08。两个正值都表明邻近轨迹以指数速率分离，定量确认了混沌动力学的存在。
**灵敏度分析。** 在无扰动系统（31）中，研究人员保持所有系统参数不变，仅略微改变初始条件。例如，将初始条件从(0.4,0)变为(0.7,0)（初始差值为0.3），轨迹出现明显分化；将初始条件(0.04,0)与(0.9,0)比较，波型差异更加显著。即使初始差小至0.2，轨迹仍然呈现显著不同的演化特征。三组初始条件(0.4,0)、(0.7,0)和(0.9,0)的比较进一步确认了系统对初值条件的强烈依赖。
**MGERIF方法精确解。** 利用多元广义指数有理积分函数方法，通过齐次平衡原理确定N=2，研究人员构造了包含任意参数的精确解析解。其中余弦函数表达式的参数设置产生了三角函数形式的精确解，如解；正弦函数形式产生了涉及正弦函数的波结构；指数函数形式给出指数型孤波解；双曲余弦形式给出了双曲函数解。这些解经过与行波变换结合，得到了NTM-III方程的多个精确解析波解，分别表示为三角函数、双曲函数和指数函数等不同形式。数值可视化显示这些解可以表现为具有lump-soliton相互作用和多峰子结构的局域波模式。
**讨论和结论。** 在讨论部分，研究人员分析了参数对系统动力行为的影响。外部摄动产生的混沌行为源于共振波相互作用机制，这与光纤、等离子体约束装置、流体输运系统、大气波和机械振子等实际系统中的不规则波传播密切相关。多稳定性的存在表明同一运行条件下系统可以支持多个允许的波传播状态，微小的环境扰动可能引起动力状态之间的阶跃转变，这对光开关器件、神经网络、激光系统、气候模型和非线性通信系统具有重要意义。灵敏度分析表明该系统具有混沌系统典型的初值敏感性，长期预测能力受到本质限制。MGERIF方法导出的精确解为局域能量输运提供了数学模型。
研究结论部分指出，通过合理的约化和行波变换，将NTM-III方程转换为非线性哈密顿系统；利用Jacobian矩阵和本征值分析确定了平衡点并研究了其稳定性，分岔相图证明了鞍-中心构型的存在并揭示了参数变化引起的稳定性交换；外部周期强迫效应的研究表明扰动系统出现混沌振荡、初值敏感性和多稳定性等现象；这些结果通过相图、时间响应、Lyapunov指数等多种手段得到验证，多吸引子在相同参数下被计算出来以证明系统的多稳定特性；MGERIF方法还获得了多个精确解析解，涵盖三角函数、双曲函数和指数解，产生了周期波、局域波、多峰子和lump型孤子等多种非线性波结构。所有这些解析和数值结果都证实了所采用的方法在研究非线性波动方程方面的有效性，并加深了对NTM-III模型动力学的理解。该研究可为可积系统和非可积系统中非线性波传播、分岔和混沌行为的后续研究提供有价值的参考。该研究成果发表在《Symmetry》期刊上。
