Processing math: 25%

数学物理学报, 2024, 44(3): 670-686

广义b方程的孤立波解及周期波解

杨佼朋,1, 梁勇,2,*

1.广东外语外贸大学数学与统计学院 广州 510006

2.华南理工大学数学学院 广州 510640

Solitary and Periodic Solutions of the Generalized b-equation

Yang Jiaopeng,1, Liang Yong,2,*

1. School of Mathematics and Statistics, Guangdong University of Foreign Studies, Guangzhou 510006

2. School of Mathematics, South China University of Technology, Guangzhou 510640

通讯作者: *梁勇,E-mail:dyliang@scut.edu.cn

收稿日期: 2023-11-21   修回日期: 2024-01-4  

基金资助: 国家自然科学基金(12201136)
广东省青年创新人才项目(2021KQNCX019)
广州市基础与应用基础研究项目(202201010278)

Received: 2023-11-21   Revised: 2024-01-4  

Fund supported: NSFC(12201136)
Foundation for Innovative Young Talents of GuangDong(2021KQNCX019)
GuangZhou Basic and Applied Basic Research Foundation(202201010278)

作者简介 About authors

杨佼朋,E-mail:jpyang@oamail.gdufs.edu.cn

摘要

对于广义 b 方程的研究主要集中在 b0 的情况, 该文利用分支方法研究了 b=3 这类特殊广义 b 方程的分支、非线性波解及动力学特征. 在一定参数条件下, 得到了该方程的分支相图, 还发现了不同于 b>0 情况的新现象, 在行波系统中有无限多周期轨穿过奇异直线 φ=c. 同时, 给出了光滑孤立波解和光滑周期波解的存在性及其精确表达式, 共获得了 15 个非线性波解的显式表达式.

关键词: 广义 b 方程; 定性理论; 分支方法; 孤立波解; 周期波解

Abstract

The research of generalized b-equation mainly focuses on the case of b0. This paper uses the bifurcation method to investigate the bifurcation, nonlinear wave solutions and dynamical characteristics of generalized b-equation with b=3. Under certain parameter conditions, one obtains the bifurcation phase diagram of the equation. Meanwhile, a new phenomenon is found different from the case of b>0, in which an infinite number of periodic trajectories in the traveling wave system cross the singular line φ=c. The existence of solitary and periodic solutions is given, and the 15 exact expressions for nonlinear wave solutions are obtained.

Keywords: Generalized b-equations; Qualitative theory; Bifurcation method; Solitary solutions; Periodic solutions

PDF (1417KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

杨佼朋, 梁勇. 广义b方程的孤立波解及周期波解[J]. 数学物理学报, 2024, 44(3): 670-686

Yang Jiaopeng, Liang Yong. Solitary and Periodic Solutions of the Generalized b-equation[J]. Acta Mathematica Scientia, 2024, 44(3): 670-686

1 引言

考虑非线性波方程

utuxxt2u2ux+3uxuxxuuxxx=0,
(1.1)

方程(1.1)为 n 次广义 b 方程

ut+2kuxuxxt+(b+1)unuxbuxuxxuuxxx=0
(1.2)

的特殊情况, 即方程(1.2)中参数 b=3, k=0, n=2.

b=2n=1 时, 方程(1.2)转化为 Camassa 和 Holm[1] 在 1993 年提出的著名浅水波方程 CH 方程.当 k=0 时, Camassa 和 Holm[1]提出 CH 方程具有可积性并且有孤立尖波解.Cooper 和 Shepard[2]研究了 CH 方程的近似孤立波解.Constantin 和 Strauss[3]给出了 CH 方程孤立波的存在性及其数学解释, 并证明了 k=0 时 CH 方程孤立尖波解的稳定性.当 b=3n=1 时, 方程(1.2)可简化为 Degasperis 和 Procesi[4]提出的 DP 方程.Lundmark 和 Szmigielski[5]利用反散射法计算出 DP 方程 3 次孤立尖波解的精确表达式. Chen 和 Tang[6]得到了 DP 方程的扭波解.

CH 方程和 DP 方程引起了学者们广泛的兴趣, 在这两个方程的研究基础上,针对不同的参数 b, n 和波速 c 有更多深入的工作. 例如,当 b>1n=1 时, Guo 和 Liu[7]证明了方程(1.2)具有周期性尖波解,Li[8]讨论了方程(1.2)Cauchy 问题的非均匀依赖性,Freire[9]给出了方程(1.2)在 Sobolev 空间中强解的全局存在性.对于给定 b=2n=2,Tian 和 Song[10]给出了方程(1.2)的一些物理解释,Liu 和 Ouyang[11]得到了 c=2 情况下的孤立尖波解及其显式表达式,当 c=1 时, Khuri[12]给出了由三角函数组成的奇异波解,Daros 和 Arruda[13]研究了椭圆行波解的不稳定性.当 b=2, n=3 时, Zhou 等[14]得到了三次 CH 方程多种行波解的精确表达式.当 b=3n=2 时, 又有如下的研究结果:Liu 和 Ouyang[11]得到了 c=52 情况下孤立尖波解的存在性及其表达式.Liu 和 Tang[15]研究了周期波解及其分支.当 b2c=b+22, Wazwaz[16]得到了两种孤立波解及其表达式.当 b0,2,3 时, Chen 等[17]研究了显式周期波及其极限形式.当 b=0, n=2 时, Yang 等[18]研究了方程(1.2)的孤立波解及周期波解, 并给出了解之间的相关联系.当 k=0, b>1n 是任意正整数时, Yang 等[19]考虑了方程(1.2)孤立尖波解的存在性及其分支.当 k=0, b=0, n=2 时, Li 和 Liu[21]研究了方程(1.2)的精确解及其分支,Li 和 Liu[21]得到了方程五种类型的爆破波解并发现了它们之间的关系.

本文利用定性分析及分支方法研究了方程(1.1)孤立波和周期波的存在性及其显式表达式, 并发现了不同于 b>0 时的新现象, 具体安排如下. 在第 2 节中, 利用行波变换得到了方程(1.1)的行波系统及其动力学特征, 并给出了一定参数条件下行波系统的分支相图. 在第 3 节中, 基于分支相图, 研究了方程孤立波的存在性及其显式表达式, 并给出了对应孤立波解的数值模拟. 在第 4 节中, 研究了周期波的存在性及其显式表达式, 并给出了对应周期解的数值模拟. 第 5 节对论文做出简要总结.

2 行波系统和分支相图

在本节中, 先利用行波变换推导出方程(1.1)的行波系统, 再研究了该系统的奇点性质及分支曲线, 最后得到了系统的分支相图.

2.1 行波系统及其动力学分析

假设波速 c 为常数, 将行波变换 u=φ(ξ), ξ=xct 代入方程(1.1), 由此得出

cφ+cφ2φ2φ+3φφφφ=0.
(2.1)

对方程(2.1)关于 φ 积分一次, 可得

(φc)φ=gcφ23φ3+2(φ)2,

其中 g 是积分常数. 令 y=φ, 可以得到如下行波系统

{dφdξ=y,dydξ=gcφ23φ3+2y2φc.
(2.2)

行波系统(2.2)具有一条奇异直线 φ=c, 这使得系统产生了奇异性, 因此该系统也被称为方程(1.1)的奇异行波系统. 为了消除系统的奇异性, 进而对系统(2.2)做时间尺度变换 dτ=dξφc, 奇异行波系统则变换为以下正则行波系统

{dφdτ=(φc)y,dydτ=gcφ23φ3+2y2.
(2.3)

显而易见, 行波系统 (2.2) 和 (2.3) 具有相同的首次积分

H(φ,y)=h=(φc)4{y243[(φc)3+ˉa22(φc)2+ˉa13(φc)+ˉa04]},
(2.4)

其中

ˉa2=3c,  ˉa1=32(2c2+c),  ˉa0=12(2c3+3c23g),

并且 h 是另一个积分常数.

相同的首次积分 H(φ,y) 意味着系统 (2.2) 和 (2.3) 除了奇异直线 φ=c 外, 两个系统具有相同的拓扑相图. 因此, 我们可以通过正则行波系统 (2.3) 的动力学来理解奇异行波系统 (2.2) 的动力学, 进而分析方程(1.1)的非线性波解.

根据系统 (2.3) 的第一个方程, 易得该系统的奇点在 φ 轴或奇直线 φ=c 上. 假设

f(φ)=gcφ23φ3,  f0(φ)=cφ23φ3.

由方程

f0(φ)=c2φ2

易知, f0(φ) 有两个零点 (±φ0,0), 其中 φ0=c2, 并且 (±φ0,f(±φ0))f(φ) 的极值点.

由系统 (2.3) 的奇点个数及性质, 可得该系统有四条分支曲线 g1(c), g1(c), g2(c)g3(c), 具体表达式,

g1(c)=16(2c)32, g2(c)=23c3+c2,g3(c)=124(12c28c3+9c2c24c(c24c)32).

分支曲线 g2(c)g1(c) 有两个交点 (c1,g1(c1))(0,0),曲线 g2(c)g1(c) 有两个交点 (c2,g2(c2))(0,0),曲线 g3(c)g1(c) 有两个交点 (c1,g1(c1))(0,0),曲线 g3(c)g2(c) 有两个交点 (c1,g1(c1))(0,0),其中 c1=12,c2=2.

为了说明系统 (2.3) 的奇点情况, 首先计算在 φ 轴上的奇点, 即 f(φ) 的零点, 分别记为 φ1, φ2, φ3 (其中 φ3φ2φ1),

1) 当 g>g1(c)g<g1(c) 时, f(φ)=0 有一个实根和两个虚根 (虚根省略),

φ1=(6gΔ)13+(6g+Δ)132,

其中 Δ=36g2+8c3.

2) 当 g=g1(c) 时, f(φ)=0 有三个实根, 其中一个为二重根,

φ1=3gc, φ2=φ3=3g2c.

3) 当 g=g1(c) 时, f(φ)=0 有三个实根, 其中一个为二重根,

φ3=3gc, φ1=φ2=3g2c.

4) 当 g1(c)<g<g1(c) 时, f(φ)=0 有三个不等实根,

φ1=2ccosθ3,φ2=c2(cosθ33sinθ3),φ3=c2(cosθ3+3sinθ3),

其中 θ=arccos(3g2(c)32).

事实上, 由上述奇点分析可得 f(φ) 的函数图像 (如图1 所示).

图1

图1   f(φ) 的函数图像


此外, 当 gg2(c) 时, 在奇直线 φ=c 上系统 (2.3) 有奇点 (c,±y0), 其中

y0=13c3+12c212g.

因此, (φi,0)(i=1,2,3)(c,±y0) 都是系统 (2.3) 的奇点.

下面研究系统 (2.3) 的奇点性质. 系统 (2.3) 的线性化系统有如下 Jacobi 矩阵

J(φ,y)=(yφcf(φ)4y),

其对应的特征方程为

|JλI|=|yλφcf(φ)4yλ|=λ2y2(φc)f(φ)=0.

因此, 在 φ 轴上的奇点 (φi,0) 处, 系统 (2.3) 的线性化系统存在特征值

λ±(φi,0)=±(φic)f(φi)(i=1,2,3).

同时, 在奇直线 φ=c 上的奇点 (c,±y0) 处, 系统 (2.3) 的线性化系统分别有相对应特征值

λ1(c,y0)=y0, λ2(c,y0)=4y0,λ1(c,y0)=y0, λ2(c,y0)=4y0.

通过系统 (2.3) 奇点的特征值分析, 由平面定性理论可知, 当 (φic)f(φi)>0 时, 奇点 (φi,0) 为鞍点; 当 (φic)f(φi)=0 时, 奇点 (φi,0) 为退化奇点; 当 (φic)f(φi)<0 时, 奇点(φi,0) 为中心. 同时, 奇点 (c,y0) 为不稳定结点, 奇点 (c,y0) 为稳定结点.

2.2 系统 (2.3) 的分支相图

由 2.1 小节中系统 (2.3) 的动力学分析, 本小节给出了系统奇点的具体分类, 并进一步得到了在 cg 平面上的分支相图. 根据平面系统的定性分析, 系统 (2.3) 奇点有以下特性.

(1) 当 g>g1(c) 时, f(φ) 只有一个零点 φ1 且满足下面不等式 φ1>c, 这意味着 (φ1,0) 是系统 (2.3) 的一个中心.

(2) 当 g=g1(c) 时, f(φ) 有两个零点 φ1φ2=φ3, 并且满足下面关系

 0<c<c1,  φ3=φ2<c<φ1, c=c1,  φ3=φ2=c<φ1, c>c1,  c<φ3=φ2<φ1.

所以, 在系统 (2.3)中 (φ1,0) 是一个中心, (φ2,0) 是退化奇点.

(3) 当 g1(c)<g<g1(c)g>g2(c) 时, f(φ) 有三个零点 φ1, φ2φ3, 它们存在以下关系

 0>c>c1,  φ3<φ2<c<φ1, c<c1,  φ1<φ2<φ3<c.

因此, 若 0>c>c1, 在系统 (2.3) 中 (φ1,0)(φ2,0) 是两个中心, (φ3,0) 是一个鞍点. 若 c<c1, 则 (φ1,0)(φ3,0) 是两个中心, (φ2,0) 是一个鞍点.

(4) 当 g1(c)<g<g1(c)g=g2(c) 时, f(φ) 有三个零点 φ1, φ2φ3, 满足

 c1<c<0,  φ3<φ2=c<φ1, c2<c<c1,  φ3=c<φ2<φ1.

所以, 我们可以得到结论: 若 c1<c<0, (φ1,0) 是一个中心, (φ2,0) 是一个三重奇点, (φ3,0) 是一个鞍点. 若c2<c<c1, (φ1,0) 是一个中心, (φ2,0) 是一个鞍点, (φ3,0) 是一个三重奇点.

(5) 当 g3(c)<g<g2(c)c1<c<0 时, f(φ) 有三个零点 φ1, φ2φ3, 它们存在以下关系

φ3<c<φ2<φ1,

这意味着, 在系统 (2.3)中 (φ1,0) 是一个中心, (φ2,0)(φ3,0) 是两个鞍点.

(6) 当 g=g3(c)c1<c<0 时, f(φ) 有三个零点 φ1, φ2φ3, 它们存在关系φ3<c<φ2<φ1, 即系统 (2.3)中 (φ1,0) 是一个中心, (φ2,0)(φ3,0) 是两个鞍点.

g=g3(c)c<c1 时, f(φ) 有一个零点 φ1 且满足下面不等式 φ1>c, 这意味着 (φ1,0) 是系统 (2.3) 的一个中心.

(7) 当 g1(c)<g<g3(c)g<g2(c) 时, f(φ) 有三个零点 φ1, φ2φ3, 满足

 c2<c<0,  φ1<φ2<c<φ3.

这意味着, 在系统 (2.3)中 (φ1,0) 是一个中心, (φ2,0)(φ3,0) 是两个鞍点.

(8) 当 g=g1(c) 时, f(φ) 有两个零点 φ1=φ2φ3, 满足

 c2<c<0,  φ1<c<φ2=φ1, c=c2,  φ3=c<φ2=φ1, c<c2,  c<φ3<φ2=φ1.

通过不等式和对应特征值, 可知若 c2<c<0, (φ1,0) 是一个退化奇点, (φ3,0) 是一个鞍点.若 c=c2, (φ1,0) 是一个退化奇点, (φ3,0) 是一个三重奇点.若 c<c2, (φ1,0) 是一个退化奇点, (φ3,0) 是一个中心.

(9) 当 g<g1(c) 时, f(φ) 只有唯一零点 φ1, 它满足以下不等式

 g<g2(c),  φ1<c, g=g2(c),  φ1=c, g>g2(c),  φ1>c.

因此, 奇点 (φ1,0) 具有以下性质:若 g<g1(c)g<g2(c), (φ1,0) 是一个鞍点;若 g=g2(c)<g1(c), (φ1,0) 是一个三重奇点;若 g2(c)<g<g1(c), (\varphi_1,0) 是一个中心.

(10) 当 g\le g_2(c) 时, 系统 (2.3) 有奇点 (c, \pm y_0), 且满足 (c, y_0) 为不稳定结点, (c, -y_0) 为稳定结点.

根据系统 (2.3) 的奇点性质 (1)-(10) 以及分支曲线 g_1(c),-g_1(c)g_2(c), 可得系统在 c-g 平面的分支相图, 如图2 所示. 从奇点性质 (4)-(10) 和分支相图, 不难看出系统 (2.3) 有无限多闭轨穿过奇异直线 \varphi=c.

图 2

图 2   系统 (2.3) 的分支相图


3 孤立波的存在性与表达式

本节讨论系统 (2.3) 中孤立波的存在性, 并给出了存在的参数条件及表达式. 由分支相图以及动力系统分支理论中同宿轨对应方程(1.1)的孤立波, 可以得到以下命题.

命题 3.1 对于给定波速 c (c<0) 和积分常数 g, 方程 (1.1) 的峰型孤立波和谷型孤立波存在参数条件如下.

(i) 在以下四种情况中存在峰型孤立波

情形 A g_2(c)<g<g_1(c), c_1<c<0.

情形 B g=g_2(c), c_1<c<0.

情形 C g_3(c)<g<g_2(c), c_1<c<0.

情形 D -g_1(c)<g<g_3(c)g<g_2(c), c_2<c<0.

(ii) 在以下四种情况中存在谷型孤立波

情形 E g_2(c)<g<g_1(c)g>-g_1(c), c<c_1.

情形 F g=g_2(c), c_2<c<c_1.

情形 G g=g_2(c), c=c_2.

情形 H -g_1(c)<g<g_3(c)g<g_2(c), c_2<c<0.

不难看出, 当 -g_1(c)<g<g_3(c)g<g_2(c), c_2<c<0 时方程 (1.1) 同时存在一个峰型孤立波和一个谷型孤立波. 下面由命题3.2给出情形 A-H 中孤立波的精确表达式, 为了方便说明该命题, 引入以下记号

\begin{align*}& H_1=H(\varphi_1,0),\ H_2=H(\varphi_2,0),\ H_3=H(\varphi_3,0),\\& w_1=\sqrt{4-6\bar{a}_2H_3-18c^2H_3^2-12H_3\varphi_3-18H_3^2\varphi_3^2+12cH_3+36cH_3^2\varphi_3^2},\\& w_2=\sqrt{4-6\bar{a}_2H_2-18c^2H_2^2-12H_2\varphi_2-18H_2^2\varphi_2^2+12cH_2+36cH_2^2\varphi_2^2},\\&\alpha_1=\alpha_4=\frac{-2+6cH_3-3H_3\varphi_3+w_1}{3H_3},\ \alpha_2=\alpha_5=\frac{-2+6cH_3-3H_3\varphi_3-w_1}{3H_3},\\ & \alpha_3=\frac{-4+9cH_3-6H_3\varphi_3}{3H_3},\ \alpha_6=\alpha_8=\frac{-2+6cH_2-3H_2\varphi_2+w_2}{3H_2}, \\& \alpha_7=\alpha_9=\frac{-2+6cH_3-3H_2\varphi_2-w_2}{3H_2},\ \alpha_{10}=\frac{-4+9cH_2-6H_2\varphi_2}{3H_2}.\end{align*}

命题 3.2 对于不同的波速 c (c<0) 和积分常数 g, 方程 (1.1) 具有以下精确的孤立波解

(1) 当 g_2(c)<g<g_1(c), c_1<c<0 时, 有如下峰型孤立波解的表达式

\begin{equation}u_1(\xi)=\varphi_3+\frac{4\mu_1 a_1 {\rm e}^{-\eta_1 \xi}}{\mu_1^2 {\rm e}^{-2\eta_1 \xi}+b_1^2-2\mu_1 b_1 {\rm e}^{-\eta_1 \xi}-4a_1},\end{equation}
(3.1)

其中

\begin{align*}&a_1=(\alpha_1-\varphi_3)(\alpha_2-\varphi_3), \ b_1=2\varphi_3-\alpha_1-\alpha_2,\\&\eta_1=\sqrt{a_1H_3},\ \chi_1=(\alpha_2-\varphi_3)^2+b_1(\alpha_2-\varphi_3)+a_1, \\&\mu_1=\frac{2\sqrt{a_ 1\chi_1}+b_1(\alpha_2-\varphi_3)+2a_1}{\alpha_2-\varphi_3}.\end{align*}

(2) 当 g=g_2(c), c_1<c<0 时, 有一个峰型孤立波解, 其表达式如下

\begin{equation}u_2(\xi)=\varphi_3+\frac{4\mu_2 a_2 {\rm e}^{-\eta_2 \xi}}{\mu_2^2 {\rm e}^{-2\eta_2 \xi}+b_2^2-2\mu_2 b_2 {\rm e}^{-\eta_2 \xi}-4a_2},\end{equation}
(3.2)

其中

\begin{align*}&a_2=(c-\varphi_3)(\alpha_3-\varphi_3), \ b_2=2\varphi_3-c-\alpha_3,\\&\eta_2=\sqrt{a_2H_3},\ \chi_2=(c-\varphi_3)^2+b_2(c-\varphi_3)+a_2, \\&\mu_2=\frac{2\sqrt{a_ 2\chi_2}+b_2(c-\varphi_3)+2a_2}{c-\varphi_3}.\end{align*}

(3) 当 g_3(c)<g<g_2(c), c_1<c<0 时, 有一个峰型孤立波解, 其表达式如下

\begin{equation}u_3(\xi)=\varphi_3+\frac{4\mu_3 a_3 {\rm e}^{-\eta_3 \xi}}{\mu_3^2 {\rm e}^{-2\eta_3 \xi}+b_3^2-2\mu_3 b_3 {\rm e}^{-\eta_3 \xi}-4a_3},\end{equation}
(3.3)

其中

\begin{align*}&a_3=(\varphi_3-\alpha_4)(\varphi_3-\alpha_5), \ b_3=2\varphi_3-\alpha_4-\alpha_5,\\&\eta_3=\sqrt{a_3H_3},\ \chi_3=(\alpha_5-\varphi_3)^2+b_3(\alpha_5-\varphi_3)+a_3, \\&\mu_3=\frac{2\sqrt{a_ 3\chi_3}+b_3(\alpha_5-\varphi_3)+2a_3}{\alpha_5-\varphi_3}.\end{align*}

(4) 当 -g_1(c)<g<g_3(c)g<g_2(c), c_2<c<0 时, 有如下峰型孤立波解的表达式

\begin{equation}u_4(\xi)=\varphi_2+\frac{4\mu_4 a_4 {\rm e}^{\eta_4 \xi}}{\mu_4^2 {\rm e}^{2\eta_4 \xi}+b_4^2+4a_4-2\mu_4 b_4 {\rm e}^{\eta_4 \xi}},\end{equation}
(3.4)

其中

\begin{align*}&a_4=(\alpha_6-\varphi_2)(\varphi_2-\alpha_7), \ b_4=\alpha_6+\alpha_7-2\varphi_2,\\&\eta_4=\sqrt{-a_4H_2},\ \chi_4=-(\alpha_6-\varphi_2)^2+b_4(\alpha_6-\varphi_2)+a_4, \\&\mu_4=\frac{2\sqrt{a_4\chi_4}+b_4(\alpha_6-\varphi_2)+2a_4}{\alpha_6-\varphi_2}.\end{align*}

(5) 当 g_2(c)<g<g_1(c)g>-g_1(c), c<c_1 时, 有如下谷型孤立波解的表达式

\begin{equation}u_5(\xi)=\varphi_2-\frac{4\mu_5 a_5 {\rm e}^{\eta_5 \xi}}{\mu_5^2 {\rm e}^{2\eta_5 \xi}+b_5^2-2\mu_5 b_5 {\rm e}^{\eta_5 \xi}-4a_5},\end{equation}
(3.5)

其中

\begin{align*}&a_5=(\varphi_2-\alpha_8)(\varphi_2-\alpha_9), \ b_5=\alpha_8+\alpha_9-2\varphi_2,\\&\eta_5=\sqrt{a_5H_2},\ \chi_5=(\varphi_2-\alpha_9)^2+b_5(\varphi_2-\alpha_9)+a_5, \\&\mu_5=\frac{2\sqrt{a_5\chi_5}+b_5(\varphi_2-\alpha_9)+2a_5}{\varphi_2-\alpha_9}.\end{align*}

(6) 当 g=g_2(c), c_2<c<c_1 时, 有一个谷型孤立波解, 其表达式如下

\begin{equation}u_6(\xi)=\varphi_2-\frac{4\mu_6 a_6 {\rm e}^{\eta_6 \xi}}{\mu_6^2 {\rm e}^{2\eta_6 \xi}+b_6^2-2\mu_6 b_6 {\rm e}^{\eta_6 \xi}-4a_6},\end{equation}
(3.6)

其中

\begin{align*}&a_6=(\varphi_2-c)(\varphi_2-\alpha_{10}), \ b_6=c+\alpha_{10}-2\varphi_2,\\&\eta_6=\sqrt{a_6 H_2},\ \chi_6=(\varphi_2-c)^2+b_6(\varphi_2-c)+a_6, \\&\mu_6=\frac{2\sqrt{a_6\chi_6}+b_6(\varphi_2-c)+2a_6}{\varphi_2-c}.\end{align*}

(7) 当 g=g_2(c), c=c_2=-2 时, 有一个谷型孤立波解, 其表达式如下

\begin{equation}u_7(\xi)=1-\frac{12}{4-9H_1 \xi^2}.\end{equation}
(3.7)

(8) 当 -g_1(c)<g<g_2(c), c_2<c<0 时, 有如下谷型孤立波解的表达式

\begin{equation}u_8(\xi)=\varphi_2-\frac{4\mu_8 a_8 {\rm e}^{\eta_8 \xi}}{\mu_8^2 {\rm e}^{2\eta_8 \xi}+b_8^2+4a_8-2\mu_8b_8 {\rm e}^{\eta_8 \xi}},\end{equation}
(3.8)

其中

\begin{align*}&a_8=(\varphi_2-\alpha_6)(\alpha_7-\varphi_2), \ b_8=2\varphi_2-\alpha_6-\alpha_7,\\&\eta_8=\sqrt{-a_8H_2},\ \chi_8=(\varphi_2-\alpha_7)^2+b_8(\varphi_2-\alpha_7)+a_8, \\&\mu_8=\frac{2\sqrt{a_8\chi_8}+b_8(\varphi_2-\alpha_7)+2a_8}{\varphi_2-\alpha_7}.\end{align*}

为了方便给出命题3.2中孤立波解的推导过程, 我们做出情形 A-H 下行波系统 (2.3) 的同宿轨道相图, 如图3.

图3

图3   系统 (2.3) 的同宿轨道相图


(1) 当 g_2(c)<g<g_1(c), c_1<c<0 时, 系统 (2.3) 存在一个同宿轨道 \Gamma_1 (如图3A), \Gamma_1 具有表达方式

\begin{equation*}y^2=H(\varphi_3,0)(\alpha_1-\varphi)(\alpha_2-\varphi)(\varphi-\varphi_3)^2,\ \mbox{且}\ \varphi_3<\varphi<\alpha_2<\alpha_1.\end{equation*}

易得 y 的表达式并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿同宿轨道 \Gamma_1 积分

\begin{equation}\int_{\alpha_2}^{\varphi}\frac{\mathrm{d}s}{(s-\varphi_3)\sqrt{(\alpha_1-s)(\alpha_2-s)}}=\sqrt{H(\varphi_3,0)}|\xi|.\end{equation}
(3.9)

在(3.9)式中, 计算积分并关于 \varphi 求解方程, 可得

\begin{equation*}\varphi=\varphi_3+\frac{4\mu_1 a_1 {\rm e}^{-\eta_1 \xi}}{\mu_1^2 {\rm e}^{-2\eta_1 \xi}+b_1^2-2\mu_1 b_1 {\rm e}^{-\eta_1 \xi}-4a_1},\end{equation*}

u=\varphi(\xi), 得到峰型孤立波解 u_1(\xi) 如方程(3.1).

(2) 当 g=g_2(c), c_1<c<0 时, 系统 (2.3) 存在一个同宿轨道 \Gamma_2 (如图3B), \Gamma_2 具有表达式

\begin{equation*}y^2=H(\varphi_3,0)(\alpha_3-\varphi)(c-\varphi)(\varphi-\varphi_3)^2,\ \mbox{且}\ \varphi_3<\varphi<c<\alpha_3.\end{equation*}

计算可得 y 的表达式并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿同宿轨道 \Gamma_2 积分

\begin{equation}\int_{c}^{\varphi}\frac{\mathrm{d}s}{(s-\varphi_3)\sqrt{(\alpha_3-s)(c-s)}}=\sqrt{H(\varphi_3,0)}|\xi|.\end{equation}
(3.10)

在 (3.10) 式中, 计算积分并关于 \varphi 求解方程, 得峰型孤立波解 u_2(\xi) 如方程 (3.2).

(3) 当 g_3(c)<g<g_2(c), c_1<c<0 时, 系统 (2.3) 存在一个同宿轨道 \Gamma_3 (如图3C), \Gamma_3 具有表达式

\begin{equation*}y^2=H(\varphi_3,0)(\alpha_4-\varphi)(\alpha_5-\varphi)(\varphi-\varphi_3)^2,\ \mbox{且}\ \varphi_3<\varphi<\alpha_5<\alpha_4.\end{equation*}

由上式可知 y 的表达式并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿同宿轨道 \Gamma_3 积分

\begin{equation}\int_{\alpha_5}^{\varphi}\frac{\mathrm{d}s}{(s-\varphi_3)\sqrt{(\alpha_4-s)(\alpha_5-s)}}=\sqrt{H(\varphi_3,0)}|\xi|.\end{equation}
(3.11)

在(3.11)式中, 计算积分并关于 \varphi 求解方程, 得峰型孤立波解 u_3(\xi) 如方程(3.3).

(4) 当-g_1(c)<g<g_3(c)g<g_2(c), c_2<c<0 时, 系统 (2.3) 存在两个同宿轨道 \Gamma_4 (如图3D) 和\Gamma_8 (如图3H), 分别有如下表达式

\begin{equation*}\Gamma_4:\ y^2=-H(\varphi_2,0)(\alpha_6-\varphi)(\varphi-\alpha_7)(\varphi-\varphi_2)^2,\ \mbox{且}\ \alpha_7<\varphi_2<\varphi<\alpha_6,\end{equation*}
\begin{equation*}\Gamma_8:\ y^2=-H(\varphi_2,0)(\varphi_2-\varphi)^2(\alpha_6-\varphi)(\varphi-\alpha_7),\ \mbox{且}\ \alpha_7<\varphi<\varphi_2<\alpha_6.\end{equation*}

同理, 分别将 y 的表达式代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 再沿着同宿轨道 \Gamma_4\Gamma_8 积分

\begin{equation}\int_{\varphi}^{\alpha_6}\frac{\mathrm{d}s}{(s-\varphi_2)\sqrt{(\alpha_6-s)(s-\alpha_7)}}=\sqrt{-H(\varphi_2,0)}|\xi|,\end{equation}
(3.12)
\begin{equation}\int_{\alpha_7}^{\varphi}\frac{\mathrm{d}s}{(\varphi_2-s)\sqrt{(\alpha_6-s)(s-\alpha_7)}}=\sqrt{-H(\varphi_2,0)}|\xi|,\end{equation}
(3.13)

关于 \varphi 分别求解方程 (3.12) 和 (3.13), 得孤立波解 u_4(\xi)u_8(\xi) 的表达式, 如方程 (3.4) 和 (3.8).

(5) 当 g_2(c)<g<g_1(c)g>-g_1(c), c<c_1 时, 系统 (2.3) 存在一个同宿轨道 \Gamma_5 (如图3E), \Gamma_5 具有表达式

\begin{equation*}y^2=H(\varphi_2,0)(\varphi-\varphi_2)^2(\varphi-\alpha_8)(\varphi-\alpha_9),\ \mbox{且}\ \alpha_9<\alpha_8<\varphi<\varphi_3.\end{equation*}

易得 y 的表达式并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿同宿轨道 \Gamma_5 积分

\begin{equation}\int_{\alpha_8}^{\varphi}\frac{\mathrm{d}s}{(\varphi_2-s)\sqrt{(s-\alpha_8)(s-\alpha_9)}}=\sqrt{H(\varphi_2,0)}|\xi|.\end{equation}
(3.14)

在 (3.14) 式中, 计算积分并关于 \varphi 求解方程, 得到谷型孤立波解 u_5(\xi) 如方程 (3.5).

(6) 当 g=g_2(c), c_2<c<c_1 时, 系统 (2.3) 存在一个同宿轨道 \Gamma_6 (如图3F), \Gamma_6 具有表达式

\begin{equation*}y^2=H(\varphi_2,0)(\varphi_2-\varphi)^2(\varphi-c)(\varphi-\alpha_{10}),\ \mbox{且}\ \alpha_{10}<c<\varphi<\varphi_2.\end{equation*}

计算上式可得 y 的表达式并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿同宿轨道 \Gamma_6 积分

\begin{equation}\int_{c}^{\varphi}\frac{\mathrm{d}s}{(\varphi_2-s)\sqrt{(s-c)(s-\alpha_{10})}}=\sqrt{H(\varphi_2,0)}|\xi|.\end{equation}
(3.15)

在 (3.15) 式中, 计算积分并关于 \varphi 求解方程, 可得谷型孤立波解 u_6(\xi) 如方程 (3.6).

(7) 当 g=g_2(c), c=c_2=-2 时, 系统 (2.3) 存在一个同宿轨道 \Gamma_7 (如图3G), \Gamma_7 具有表达式

\begin{equation*}y^2=-H(\varphi_1,0)(\varphi_1-\varphi)^2(\varphi-c)(1-\varphi),\ \mbox{且}\ c<\varphi<\varphi_1=1.\end{equation*}

易得 y 的表达式并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿同宿轨道 \Gamma_7 积分

\begin{equation}\int_{-2}^{\varphi}\frac{\mathrm{d}s}{(1-s)\sqrt{-s^2-s+2}}=\sqrt{-H(\varphi_1,0)}|\xi|.\end{equation}
(3.16)

在 (3.16) 式中, 计算积分并关于 \varphi 求解方程, 有谷型孤立波解 u_7(\xi) 如方程 (3.7).

因此, 完成了命题 3.2 的推导.

对于给定的波速 c (c<0) 和积分常数 g, 可以利用软件数值模拟出孤立波解 u_1(\xi)-u_8(\xi) 的波形图. 例如, 情形 A 的u_1(\xi) 中取 c=-0.1, g=0.01;情形 B 的 u_2(\xi) 中取 c=-0.4, g=g_2(c);情形 C 的 u_3(\xi) 中取 c=-0.2, g=0.02;情形 D 的 u_4(\xi) 中取 c=-1.5, g=-0.5;情形 E 的 u_5(\xi) 中取 c=-1, g=0.4;情形 F 的 u_6(\xi) 中取 c=-1, g=g_2(c);情形 G 的 u_7(\xi) 中取 c=-2, g=g_2(c);情形 H 的 u_8(\xi) 中取 c=-1, g=0, 具体波形图见图4.

图4

图4   方程 (1.1) 孤立波解 u_1(\xi)-u_8(\xi) 的波形图


4 周期波的表达式

本节研究了系统 (2.3) 在情形 A-G 中周期波的存在性及其显式表达式. 根据动力系统分支理论中封闭轨道对应方程 (1.1) 的光滑周期波, 可以得到命题 4.1. 为了方便说明命题 4.1, 给出以下记号

\begin{align*}&w_3=\sqrt{4-6\bar{a}_2H_1-18c^2H_1^2-12H_1\varphi_1-18H_1^2\varphi_1^2+12cH_1+36cH_1^2\varphi_1^2},\\ &\alpha_{11}=\alpha_{14}=\alpha_{16}=\alpha_{18}=\frac{-2+6cH_1-3H_1\varphi_1-w_3}{3H_1},\\ &\alpha_{12}=\alpha_{15}=\alpha_{17}=\alpha_{19}=\frac{-2+6cH_1-3H_1\varphi_1+w_3}{3H_1},\\ &\alpha_{13}=\alpha_{20}=\frac{-4+9cH_1-6H_1\varphi_1}{3H_1},\ H=H(\alpha_{21},0),\\ \ & \rho=-216H-32,\ \alpha_{21}=\frac{1}{18H}\left((-8-36H)\frac{2^{\frac{8}{3}}(27H+4)}{\rho^{\frac{1}{3}}}+2^{\frac{4}{3}}\rho^{\frac{1}{3}} \right).\end{align*}

命题 4.1 对于不同的波速 c (c<0) 和积分常数 g, 方程 (1.1) 具有以下精确的光滑周期波解.

(1) 当 g_2(c)<g<g_1(c), c_1<c<0 时, 有如下周期波解的表达式

\begin{equation}u_9(\xi)=\varphi_1-\frac{2a_9}{\sqrt{q_1} \sin{(\mu_9-\eta_9 \xi)}-b_9},\end{equation}
(4.1)

其中

\begin{align*}&a_9=(\alpha_{11}-\varphi_1)(\varphi_1-\alpha_{12}), \ b_9=2\varphi_1-\alpha_{11}-\alpha_{12},\\&\eta_9=\sqrt{a_9H_1},\ q_1=b_9^2+4a_9, \\&\mu_9=\arcsin\left(\frac{b_9(\varphi_1-\alpha_{11})+2a_9}{\sqrt{q_1}(\varphi_1-\alpha_{11})}\right).\end{align*}

(2) 当 g=g_2(c), c_1<c<0 时, 有一个周期波解的表达式如下

\begin{equation}u_{10}(\xi)=\varphi_1-\frac{2a_{10}}{\sqrt{q_2} \sin{(\mu_{10}-\eta_{10} \xi)}-b_{10}},\end{equation}
(4.2)

其中

\begin{align*}&a_{10}=(c-\varphi_1)(\varphi_1-\alpha_{13}), \ b_{10}=2\varphi_1-c-\alpha_{13},\\&\eta_{10}=\sqrt{a_{10}H_1},\ q_2=b_{10}^2+4a_{10}, \\&\mu_{10}=\arcsin\left(\frac{b_{10}(\varphi_1-c)+2a_{10}}{\sqrt{q_2}(\varphi_1-c)}\right).\end{align*}

(3) 当 g_3(c)<g<g_2(c), c_1<c<0 时, 有一个周期波解的表达式如下

\begin{equation}u_{11}(\xi)=\varphi_1-\frac{2a_{11}}{\sqrt{q_3} \sin{(\mu_{11}-\eta_{11} \xi)}-b_{11}},\end{equation}
(4.3)

其中

\begin{align*}&a_{11}=(\alpha_{15}-\varphi_1)(\varphi_1-\alpha_{14}), \ b_{11}=2\varphi_1-\alpha_{14}-\alpha_{15},\\&\eta_{11}=\sqrt{a_{11} H_1},\ q_3=b_{11}^2+4a_{11}, \\&\mu_{11}=\arcsin\left(\frac{b_{11}(\varphi_1-\alpha_{14})+2a_{11}}{\sqrt{q_2}(\varphi_1-\alpha_{14})}\right).\end{align*}

(4) 当 -g_1(c)<g<g_3(c)g<g_2(c), c_2<c<0 时, 有如下周期波解的表达式

\begin{equation}u_{12}(\xi)=\varphi_1-\frac{2a_{12}}{\sqrt{q_4} \sin{(\mu_{12}-\eta_{12} \xi)}-b_{12}},\end{equation}
(4.4)

其中

\begin{align*}&a_{12}=(\alpha_{16}-\varphi_1)(\varphi_1-\alpha_{17}), \ b_{12}=2\varphi_1-\alpha_{16}-\alpha_{17},\\&\eta_{12}=\sqrt{a_{12}H_1},\ q_4=b_{12}^2+4a_{12}, \\&\mu_{12}=\arcsin\left(\frac{b_{12}(\varphi_1-\alpha_{16})+2a_{12}}{\sqrt{q_4}(\varphi_1-\alpha_{16})}\right).\end{align*}

(5) 当 g_2(c)<g<g_1(c)g>-g_1(c), c<c_1 时, 有如下周期波解的表达式

\begin{equation}u_{13}(\xi)=\varphi_1-\frac{2a_{13}}{\sqrt{q_5} \sin{(\mu_{13}-\eta_{13} \xi)}-b_{13}},\end{equation}
(4.5)

其中

\begin{align*}&a_{13}=(\alpha_{18}-\varphi_1)(\varphi_1-\alpha_{19}), \ b_{13}=2\varphi_1-\alpha_{18}-\alpha_{19},\\&\eta_{13}=\sqrt{a_{13}H_1},\ q_4=b_{13}^2+4a_{13}, \\&\mu_{13}=\arcsin\left(\frac{b_{13}(\varphi_1-\alpha_{19})+2a_{13}}{\sqrt{q_5}(\varphi_1-\alpha_{19})}\right).\end{align*}

(6) 当 g=g_2(c), c_2<c<c_1 时, 有一个周期波解, 其表达式如下

\begin{equation}u_{14}(\xi)=\varphi_1-\frac{2a_{14}}{\sqrt{q_5} \sin{(\mu_{14}-\eta_{14} \xi)}-b_{14}},\end{equation}
(4.6)

其中

\begin{align*}&a_{14}=(c-\varphi_1)(\varphi_1-\alpha_{20}), \ b_{14}=2\varphi_1-c-\alpha_{20},\\&\eta_{14}=\sqrt{a_{14}H_1},\ q_5=b_{14}^2+4a_{14}, \\&\mu_{14}=\arcsin\left(\frac{b_{14}(\varphi_1-c)+2a_{14}}{\sqrt{q_4}(\varphi_1-c)}\right).\end{align*}

(7) 当 g=g_2(c), c=c_2 时, 有一个周期波解, 其表达式如下

\begin{equation}u_{15}(\xi)=\frac{\alpha_{21}B-2A-(\alpha_{21}B+2A) \mathrm{cn}(\sqrt{-ABH},l)}{A+B+(A-B) \mathrm{cn}(\sqrt{-ABH},l)},\end{equation}
(4.7)

其中

\begin{align*}&H=H(\alpha_{21},0),\\&l^2=\frac{(\alpha_{21}-c)^2-(A-B)^2}{4 A B},\\&a_1^2=-\frac{(\beta_1-\beta_2)^2}{4},\ b_1=\frac{\beta_1+\beta_2}{2},\\&A^2=(\alpha_{21}-b_1)^2+a_1^2, \ B^2=(c-b_1)^2+a_1^2, \\\end{align*}

并且 \alpha_{21}, \beta_1\beta_2 由方程

\begin{equation*}(\varphi-c)^4+\frac{4(\varphi-c)^3}{3H}+\frac{4\bar{a}_2(\varphi-c)^2}{3H}+\frac{9\bar{a}_1(\varphi-c)}{9H}+\frac{\bar{a}_0}{3H}=(\varphi-\alpha_{21})(\varphi-c)(\varphi-\beta_1)(\varphi-\beta_2)\end{equation*}

所确定, 计算可得 \beta_1\beta_2 的表达式如下

\begin{align*}\beta_1=&\frac{729\cdot2^{\frac{1}{3}}H^3-18\left(3\cdot2^{\frac{2}{3}}\rho^{\frac{1}{3}}+2\rho^{\frac{2}{3}}-18\cdot2^{\frac{1}{3}}\right)H^2-\left(27\cdot2^{\frac{1}{3}}\rho_1+2^{\frac{11}{3}}\rho^{\frac{1}{3}}+8\rho^{\frac{2}{3}}-2^{\frac{16}{3}}\right)H}{18H^2\rho^{\frac{2}{3}}}\\&+\frac{\sqrt{-6\rho_1H\left(19683\cdot2^{\frac{2}{3}}H^3+2916\left(\rho^{\frac{1}{3}}+3\cdot2^{\frac{2}{3}}\right)H^2-27\rho_2H-4\rho_3\right)}}{18H^2\rho^{\frac{2}{3}}},\\\beta_2=&\frac{729\cdot2^{\frac{1}{3}}H^3-18\left(3\cdot2^{\frac{2}{3}}\rho^{\frac{1}{3}}+2\rho^{\frac{2}{3}}-18\cdot2^{\frac{1}{3}}\right)H^2-\left(27\cdot2^{\frac{1}{3}}\rho_1+2^{\frac{11}{3}}\rho^{\frac{1}{3}}+8\rho^{\frac{2}{3}}-2^{\frac{16}{3}}\right)H}{18H^2\rho^{\frac{2}{3}}}\\&-\frac{\sqrt{-6\rho_1H\left(19683\cdot2^{\frac{2}{3}}H^3+2916\left(\rho^{\frac{1}{3}}+3\cdot2^{\frac{2}{3}}\right)H^2-27\rho_2H-4\rho_3\right)}}{18H^2\rho^{\frac{2}{3}}},\\\end{align*}

其中

\begin{align*} & \rho=-216H-32, \ \rho_1=27H^2+4H, \\ & \rho_2=27\cdot2^{\frac{2}{3}}\rho_1-48\rho^{\frac{1}{3}}-2^{\frac{7}{3}}\rho^{\frac{2}{3}}-2^{\frac{20}{3}},\\ & \rho_3=27\rho^{\frac{1}{3}}\rho_1+54\cdot2^{\frac{2}{3}}\rho_1-32\rho^{\frac{1}{3}}-2^{\frac{7}{3}}\rho^{\frac{2}{3}}-2^{\frac{17}{3}}.\end{align*}

下面给出命题 4.1 的推导过程. 为了方便求出命题中周期波解, 给出七种情形下行波系统 (2.3) 的封闭轨道相图, 如图5.

图5

图5   系统 (2.3) 的封闭轨道相图


(1) 当 g_2(c)<g<g_1(c), c_1<c<0 时, 系统 (2.3) 存在一个封闭轨道 \Gamma_9 (如图5A), \Gamma_9 具有表达式

\begin{equation*}y^2=-H(\varphi_1,0)(\varphi_1-\varphi)^2(\alpha_{11}-\varphi)(\varphi-\alpha_{12}),\ \mbox{且}\ \alpha_{12}<\varphi<\alpha_{11}<\varphi_1.\end{equation*}

由上式可解得 y 的表达式, 并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿封闭轨道 \Gamma_9 积分

\begin{equation}\int_{\alpha_{11}}^{\varphi}\frac{\mathrm{d}s}{(\varphi_1-s)\sqrt{(\alpha_{11}-s)(s-\alpha_{12})}}=\sqrt{-H(\varphi_1,0)}|\xi|.\end{equation}
(4.8)

在 (4.8) 式中, 计算积分并关于 \varphi 求解方程, 可得周期波解 u_9(\xi)=\varphi, 如方程 (4.1).

(2) 当 g=g_2(c), c_1<c<0 时, 系统 (2.3) 存在一个封闭轨道 \Gamma_{10} (如图5B), \Gamma_{10} 具有表达式

\begin{equation*}y^2=-H(\varphi_1,0)(\varphi_1-\varphi)^2(c-\varphi)(\varphi-\alpha_{13}),\ \mbox{且}\ \alpha_{13}<\varphi<c<\varphi_1.\end{equation*}
(4.9)

计算可得 y 的表达式并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿封闭轨道 \Gamma_{10} 积分

\begin{equation}\int_{c}^{\varphi}\frac{\mathrm{d}s}{(\varphi_1-s)\sqrt{(c-s)(s-\alpha_{13})}}=\sqrt{-H(\varphi_1,0)}|\xi|.\end{equation}

在 (4.9) 式中, 计算积分并关于 \varphi 求解方程, 可得周期波解 u_{10}(\xi) 如方程(4.2).

(3) 当 g_3(c)<g<g_2(c), c_1<c<0 时, 系统 (2.3) 存在一个封闭轨道 \Gamma_{11} (如图5C), \Gamma_{11} 具有表达式

\begin{equation*}y^2=-H(\varphi_1,0)(\varphi_1-\varphi)^2(\alpha_{14}-\varphi)(\varphi-\alpha_{15}),\ \mbox{且}\ \alpha_{15}<\varphi<\alpha_{14}<\varphi_1.\end{equation*}

易得 y 的表达式并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿封闭轨道 \Gamma_{11} 积分

\begin{equation}\int_{\alpha_{14}}^{\varphi}\frac{\mathrm{d}s}{(\varphi_1-s)\sqrt{(\alpha_{14}-s)(s-\alpha_{15})}}=\sqrt{-H(\varphi_1,0)}|\xi|.\end{equation}
(4.10)

在 (4.10) 式中, 计算积分并关于 \varphi 求解方程, 可得周期波解 u_{11}(\xi) 如方程 (4.3).

(4) 当 -g_1(c)<g<g_2(c), c_2<c<0 时, 系统 (2.3) 存在一个封闭轨道 \Gamma_{12} (如图5D), \Gamma_{12} 具有表达式

\begin{equation*}y^2=-H(\varphi_1,0)(\varphi_1-\varphi)^2(\alpha_{16}-\varphi)(\varphi-\alpha_{17}),\ \mbox{且}\ \alpha_{17}<\varphi<\alpha_{16}<\varphi_1.\end{equation*}

计算可得 y 的表达式并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿封闭轨道 \Gamma_{12} 积分

\begin{equation}\int_{\varphi}^{\alpha_{16}}\frac{\mathrm{d}s}{(\varphi_1-s)\sqrt{(\alpha_{16}-s)(s-\alpha_{17})}}=\sqrt{-H(\varphi_1,0)}|\xi|.\end{equation}
(4.11)

在 (4.11) 式中, 计算积分并关于 \varphi 求解方程, 可得周期波解 u_{12}(\xi) 如方程 (4.4).

(5) 当 g_2(c)<g<g_1(c)g>-g_1(c), c<c_1 时, 系统 (2.3) 存在一个封闭轨道 \Gamma_{13} (如图5E), \Gamma_{13} 具有表达式

\begin{equation*}y^2=-H(\varphi_1,0)(\varphi_1-\varphi)^2(\alpha_{18}-\varphi)(\varphi-\alpha_{19}),\ \mbox{且}\ \alpha_{19}<\varphi<\alpha_{18}<\varphi_1.\end{equation*}

由上式可得 y 的表达式并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿封闭轨道 \Gamma_{13} 积分

\begin{equation}\int_{\alpha_{19}}^{\varphi}\frac{\mathrm{d}s}{(\varphi_1-s)\sqrt{(\alpha_{18}-s)(s-\alpha_{19})}}=\sqrt{-H(\varphi_1,0)}|\xi|.\end{equation}
(4.12)

在 (4.13) 式中, 计算积分并关于 \varphi 求解方程, 可得周期波解 u_{13}(\xi) 如方程 (4.5).

(6) 当 g=g_2(c), c_2<c<c_1 时, 系统 (2.3) 存在一个封闭轨道 \Gamma_{14} (如图5F), \Gamma_{14} 具有表达式

\begin{equation*}y^2=-H(\varphi_1,0)(\varphi_1-\varphi)^2(\alpha_{20}-\varphi)(\varphi-c),\ \mbox{且}\ c<\varphi<\alpha_{20}<\varphi_1.\end{equation*}

易得 y 的表达式并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿封闭轨道 \Gamma_{14} 积分

\begin{equation}\int_{c}^{\varphi}\frac{\mathrm{d}s}{(\varphi_1-s)\sqrt{(\alpha_{20}-s)(s-c)}}=\sqrt{-H(\varphi_1,0)}|\xi|.\end{equation}
(4.13)

在 (4.13) 式中, 计算积分并关于 \varphi 求解方程, 可得周期波解 u_{14}(\xi) 如方程 (4.6).

(7) 当 g=g_2(c), c=c_2 时, 系统 (2.3) 存在一个封闭轨道 \Gamma_{15} (如图5G), \Gamma_{15} 具有表达式

\begin{equation*}y^2=-H(\alpha_{21},0)(\alpha_{21}-\varphi)(\varphi-c)(\varphi-\beta_1)(\varphi-\beta_2),\ \mbox{且}\ c<\varphi<\alpha_{21},\ \beta_1 \mbox{和} \beta_2 \mbox{为虚根}.\end{equation*}

同理可得 y 的表达式并代入方程 \frac{\mathrm{d}\varphi}{y}=\mathrm{d}\xi, 然后将所得方程沿封闭轨道 \Gamma_{15} 积分

\begin{equation}\int_{c}^{\varphi}\frac{\mathrm{d}s}{(\alpha_{21}-s)(s-c)(s-\beta_1)(s-\beta_2))}=\sqrt{-H(\alpha_{21},0)}|\xi|.\end{equation}
(4.14)

在 (4.14) 式中, 计算积分并关于 \varphi 求解方程, 可得周期波解 u_{15}(\xi) 如方程 (4.7).

因此, 完成了命题 4.1 的推导.

对于给定的波速 c (c<0) 和积分常数 g, 可以数值模拟出周期波解 u_9(\xi)-u_{15}(\xi) 的波形图. 例如,情形 A 的 u_9(\xi) 中取 c=-0.2, g=0.04;情形 B 的 u_{10}(\xi) 中取 c=-0.1, g=g_2(c);情形 C 的 u_{11}(\xi) 中取 c=-0.2, g=0.02;情形 D 的 u_{12}(\xi) 中取 c=-1, g=0;情形 E 的 u_{13}(\xi) 中取 c=-1, g=0.4;情形 F 的 u_{14}(\xi) 中取 c=-1, g=g_2(c);情形 G 的 u_{15}(\xi) 中取 c=-2, g=g_2(c), H=-0.5,具体波形图见图6.

图6

图6   方程 (1.1) 周期波解 u_9(\xi)-u_{15}(\xi) 的波形图


5 总结

本文利用动力系统分支方法研究了方程 (1.1) 的分支、非线性波解及动力学特征, 并给出了孤立波解和周期波解的存在性及其精确表达式.由于方程 (1.1) 对应的奇异行波系统 (2.2) 具有一条奇异直线 \varphi=c, 并且该直线在正则行波系统 (2.3) 中有两个结点,因此方程 (1.1) 不存在孤立尖波解和紧孤立波解.此外, 还发现了在行波系统中有无限多周期轨穿过奇异直线 \varphi=c, 该现象在 b>0 时不存在.

参考文献

Camassa R, Holm D D.

An integrable shallow water equation with peaked solitons

Phys Rev Lett, 1993, 71: 1661-1664

DOI:10.1103/PhysRevLett.71.1661      PMID:10054466      [本文引用: 2]

Cooper F, Shepard H.

Solitons in the Camassa-Holm shallow water equation

Phys Lett A, 1994, 194: 246-250

[本文引用: 1]

Constantin A, Strauss W A.

Stability of peakons

Commun Pur Appl Math, 2000, 53: 603-610

[本文引用: 1]

Degasperis A, Procesi M.

Asymptotic integrability

Symmetry and Perturbation Theory, 1999, 1: 23-37

[本文引用: 1]

Lundmark H, Szmigielski J.

Multi-peakon solutions of the Degasperis-Procesi equation

Inverse Prob, 2003, 19: 1241-1245

[本文引用: 1]

Chen C, Tang M Y.

A new type of bounded waves for Degasperis-Procesi equations

Chaos Solit Fract, 2006, 27: 698-704

[本文引用: 1]

Guo B L, Liu Z R.

Periodic cusp wave solutions and single-solitons for the b-equation

Chaos Solit Fract, 2005, 23: 1451-1463

[本文引用: 1]

Li Y.

Nonunfiorm dependence for the Cauchy problem of the general b-equation

J Math Phys, 2001, 52: 033101

[本文引用: 1]

Freire I L.

Remarks on strong global solutions of the b-equation

Appl Math Lett, 2023, 146: 108820

[本文引用: 1]

Tian L X, Song X Y.

New peaked solitary wave solutions of the generalized Camassa-Holm equation

Chaos Solit Fract, 2004, 21: 621-637

[本文引用: 1]

Liu Z R, Ouyang Z Y.

A note on solitary waves for modified forms of Camassa-Holm and Degasperis-Procesi equations

Phys Lett A, 2007, 366: 377-381

[本文引用: 2]

Khuri S A.

New ansatz for obtaining wave solutions of the generalized Camassa-Holm equation

Chaos Solit Fract. 2005, 25: 705-710

[本文引用: 1]

Daros A, Arruda L K.

On the instability of elliptic traveling wave solutions of the modified Camassa-Holm equation

J Differ Equations, 2019, 266: 1946-1968

[本文引用: 1]

Zhou Y Q, Chen G R, Li J B.

Bifurcations, exact peakon, periodic peakons and solitary wave solutions of the cubic Camassa-Holm type equation

Int J Bifurcat Chaos, 2023, 33: 2350014

[本文引用: 1]

Liu Z R, Tang H.

Explicit periodic wave solutions and their bifurcations for generalized Camassa-Holm equation

Int J Bifurcat Chaos, 2010, 20: 2507-2519

[本文引用: 1]

Wazwaz A M.

New solitary wave solutions to the modified forms of Degasperis-Procesi and Camassa-Holm equations

Appl Math Comput, 2007, 186: 130-141

[本文引用: 1]

Chen Y R, Ye W B, Liu R.

The explicit periodic wave solutions and their limit forms for a generalized b-equation

Acta Math Appl Sin, 2016, 32: 513-528

[本文引用: 1]

Yang J P, Liu R, Chen Y R.

Bifurcations of solitary waves of a simple equation

Int J Bifurcat Chaos, 2020, 30: 2050138

[本文引用: 1]

Yang J P, Li Z G, Liu Z R.

The existence and bifurcation of peakon and anti-peakon to the n-degree b-equation

Int J Bifurcat Chaos, 2018, 28: 1850014

[本文引用: 1]

Li Z G, Liu R.

Bifurcations and exact solutions in a nonlinear wave equation

Int J Bifurcat Chaos, 2019, 29: 1950098

Li Z G, Liu R.

Blow-up solutions for a case of b-family equations

Acta Math Sci, 2020, 40B(4): 910-920

[本文引用: 2]

/