一类广义浅水波KdV方程的可积性研究

一类广义浅水波KdV方程的可积性研究

郝晓红^,¹, 程智龙^,²

The Integrability of the KdV-Shallow Water Waves Equation

Hao Xiaohong^,¹, Cheng Zhilong^,²

通讯作者: 程智龙, E-mail: zhilong0793@sina.cn

收稿日期: 2017-04-10

基金资助:

安徽省自然科学研究项目. KJ2016A071

Received: 2017-04-10

Fund supported:

the Natural Science Foundation of Anhui Province. KJ2016A071

作者简介 About authors

郝晓红,haoxiaohong200866@163.com , E-mail：haoxiaohong200866@163.com

摘要

该文应用双Bell多项式，系统研究了一类广义浅水波KdV方程的可积性.先构造出双线性表达式、Bäklund变换，再通过Bäklund变换线性化得到孤子解与Lax对.最后通过级数展开式代入得到无穷守恒律，从而证明此方程具有可积性.

关键词： Bäklund变换 ; Lax对 ; 无穷守恒律

Abstract

In this paper, the binary Bell polynomials to construct bilinear forma, bilinear Bäcklund transformation, Lax pair of the KdV-shallow water waves equation. Through bilinear Bäcklund transformation, some soliton solutions are presented. Moreover, the infinite conservation laws are also derived by Bell polynomials, all conserved densities and fluxes are given with explicit recursion formulas.

Keywords： Bäcklund transformation ; Lax pair ; Infinite conservation laws

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

本文引用格式

郝晓红, 程智龙. 一类广义浅水波KdV方程的可积性研究. 数学物理学报[J], 2019, 39(3): 451-460 doi:

Hao Xiaohong, Cheng Zhilong. The Integrability of the KdV-Shallow Water Waves Equation. Acta Mathematica Scientia[J], 2019, 39(3): 451-460 doi:

1 引言

近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼 $\theta$ 函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注.

1996年, Lambert, Gilson和Nimmo建立了Bell多项式与Hirota双线性算子之间的关系^[14],通过转换关系得到双线性Bäcklund变换,这个方法很有效地避免了求Bäcklund变换过程中使用交换公式繁琐的计算,简洁实用.并且直接对其做线性化还可以得到方程的Lax对^{[9, 12]},最后将级数展开式代入计算求解得到方程的无穷守恒律.

本文重点应用Bell多项式研究一类广义浅水波KdV方程

$\begin{equation} r_{t}-\alpha r_{xxt}-4\alpha r r_{t}-2 \alpha r_{x}\int_x^{\infty}r_{t}{\rm d}x+ \alpha r_x+\beta r_{xxx}+6\beta rr_x=0, \end{equation}$

(1.1)

其中 $\alpha$ 与 $\beta$ 是任意的非零常数.方程(1.1)包涵一些特殊的物理方程:

(ⅰ)当 $\alpha$ =0, $\beta$ =1时

$\begin{equation} r_{t}+r_{xxx}+6rr_x=0, \end{equation}$

(1.2)

即为KdV方程^[1-3];

(ⅱ)当 $\alpha$ =1, $\beta$ =0时

$\begin{equation} r_{t}-r_{xxt}-4r r_{t}-2 r_{x}\int_x^{\infty}r_{t}{\rm d}x+r_x=0, \end{equation}$

(1.3)

即为浅水波方程^{[1, 4]}. Ablowitz等用反散射法研究了此类方程^[13]. Gilson, Nimmo和Willox考虑其Wronskian形式的孤子解^[18]. Shang和Hong通过其次平衡法与扩展双曲函数法得到大量的精确解与Bäcklund变换形式^[19], Wazwaz则应用简化形式的Hirota双线性法求解出 $N$ -孤子解^[4].

本文重点研究方程(1.1)的可积性,即为双线性表达式、Bäcklund变换、Lax对与无穷守恒律.本文结构如下:第二部分给出必要Bell多项式的定义与性质;第三部分应用Bell多项式法研究其可积性;最后部分给出结论与参考文献.

2 多维双Bell多项式

首先我们先给出Bell多项式的定义以及性质:

定义2.1 设 $f=f(x_{1}, x_{2}, \cdots, x_{n})$ 是具有 $n$ 个变量的 ${\cal C}^{\infty}$ 函数,则称

$\begin{equation}Y_{n_{1}x_{1}, \cdots, n_{l}x_{l}}(f)\equiv Y_{n_{1}, \cdots, n_{l}}(f_{r_{1}x_{1}, \cdots, r_{l}x_{l}})= e^{-f}\partial_{x_{1}}^{n_{1}}\cdots\partial_{x_{l}}^{n^{l}}e^{f} \end{equation}$

(2.1)

为Bell多项式^{[9, 14]}.其中 $f_{r_{1}x_{1}, \cdots, r_{l}x_{l}}=\partial_{x_{1}}^{n_{1}}\cdots\partial_{x_{l}}^{n^{l}}f (r_{1}=0, \cdots, n_{1};\cdots; r_{l}=0, \cdots, n_{l}.)$ .当 $f=f(x, t)$ 时,对应的2-维Bell多项式,举例为

$\begin{array}{l}&&Y_{x}{(f)}=f_{x}, \qquad Y_{2x}{(f)}=f_{2x}+f_{x}^2, \qquadY_{x, t}{(f)}=f_{x, t}+f_{x}f_{t}, \\&&Y_{2x, t}{(f)}=f_{2x, t}+f_{2x}f_{t}+2f_{x, t}f_{x}+f_{x}^{2}f_{t}, \\&&Y_{3x}{(f)}=f_{3x}+3f_{2x}f_{x}+f_{x}^{3}, \qquad\cdots .\end{array}$

(2.2)

定义2.2 基于上述Bell多项式定义,多维的双Bell多项式定义为

${\cal Y}_{n_{1}x_{1}, \cdots, n_{l}x_{l}}(v, w) = Y_{n_{1}, \cdots, n_{l}}(f)|_{f_{r_{1}x_{1}, \cdots, r_{l}x_{l}}}\\ = \left\{ \begin{array}{ll}v_{r_{1}x_{1}, \cdots, r_{l}x_{l}}, & r_{1}+r_{2}+ \cdots +r_{l}\ \mbox{为奇数, } \\w_{r_{1}x_{1}, \cdots, r_{l}x_{l}}, & n_{1}+n_{2}+\cdots+n_{l}\ \mbox{为偶数, }\end{array} \right.$

(2.3)

则双Bell多项式可以表示为函数具有 $v$ 和 $w$ 的形式,举例为

$\begin{array}{l}&&{\cal Y}_{x}(v)=v_{x}, {\cal Y}_{2x}(v, w)=w_{2x}+v_{x}^{2}, {\cal Y}_{x, t}(v, w)=w_{xt}+v_{x}v_{t}, \\&&{\cal Y}_{2x, t}(v, w)=v_{2x, t}+w_{2x}v_{t}+2w_{xt}v_x+v_x^2v_t, \\&&{\cal Y}_{3x}(v, w)=v_{3x}+3v_{x}w_{2x}+v_{x}^{3}, \cdots .\end{array}$

(2.4)

定义2.3 ${\cal Y}$ 多项式和Hirota双线性 $D$ 算子之间的关系可由如下等式得出

$\begin{equation}{\cal Y}_{n_{1}x_{1}, \cdots, n_{l}x_{l}}(v=\ln F/G, w=\lnFG)=(FG)^{-1}D_{x_{1}}^{n_{1}}\cdots D_{x_{l}}^{n_{l}}F\cdot G, \end{equation}$

(2.5)

其中 $n_{1}+n_{2}+\cdots+n_{l}\geq 1$ .算子 $D_{x_1}, \cdots, D_{x_l}$ 为Hirota双线性算子定义为

$D_{x_1}^{n_1}\cdots D_{x_l}^{n_l}F\cdot G=(\partial_{x_1}-\partial_{x_1' })^{n_1} \cdots(\partial_{x_l}-\partial_{x_l' })^{n_l}F(x_1, \cdots, x_l)\times G(x_1', \cdots, x_l' )\mid_{x_1' =x_1, \cdots, x_l' =x_l}.$

特别的,当 $F=G$ 时, $(2.5)$ 式则化为

$\begin{array}{ll}(F)^{-2}D_{x_{1}}^{n_{1}}\cdots D_{x_{l}}^{n_{l}}F\cdot F = {\cal Y}_{n_{1}x_{1}, \cdots, n_{l}x_{l}}(0, q=2\ln F)\\ = \left\{\begin{array}{ll} 0, &n_{1}+n_{2}+ \cdots+n_{l} \ \mbox{为奇数, } \\ P_{n_{1}x_{1}, \cdots, n_{l}x_{l}}(q), & n_{1}+n_{2}+ \cdots+n_{l}\ \mbox{为偶数.} \end{array}\right.\end{array}$

(2.6)

称之为 $P$ 多项式,此类多项式具有易于被识别的偶数部分划分结构的特点,如

$\begin{array}{l}&&P_{2x}(q)=q_{2x}, \qquad P_{x, t}(q)=q_{xt}, \\&&P_{4x}(q)=q_{4x}+3q_{2x}^{2}, \qquadP_{6x}(q)=q_{6x}+15q_{2x}q_{4x}+15q_{2x}^3 \cdots .\end{array}$

(2.7)

定义2.4 双Bell多项式 ${\cal Y}_{n_{1}x_{1}, \cdots, n_{l}x_{l}}(v, w)$ 可分离成 $P$ 多项式和 $Y$ 多项式的组合形式

$\begin{array}{l} (FG)^{-1}D_{x_{1}}^{n_{1}}\cdots D_{x_{l}}^{n_{l}}F\cdot G\\ = {\cal Y}_{n_{1}x_{1}, \cdots, n_{l}x_{l}}(v, w)|_{v=\ln F/G, w=\ln FG} ={\cal Y}_{n_{1}x_{1}, \cdots, n_{l}x_{l}}(v, v+q)|_{v=\ln F/G, q=2\ln G}\\ = \sum\limits_{n_{1}+\cdots+n_{l}=even}^{}\sum\limits_{r_{1} =0}^{n_{1}}\cdots\sum\limits_{r_{l}=0}^{n_{l}}\prod\limits_{i=1}^{l}(n_{i}, r_{i})' P_{r_{1}x_{1}, \cdots, r_{l}x_{l}}(q)Y_{(n_{1}-r_{1})x_{1}, \cdots, (n_{l}-r_{l})x_{l}}(v), \end{array}$

(2.8)

注意

$\begin{equation}Y_{n_{1}x_{1}, \cdots, n_{l}x_{l}}(v)|_{v=\ln\psi}=\frac{\psi_{n_{1}x_{1}, \cdots, n_{l}x_{l}}}{\psi}, \end{equation}$

(2.9)

这意味着一个Bell多项式 ${\cal Y}_{n_{1}x_{1}, \cdots, n_{l}x_{l}}(v, w)$ 可以通过Hopf-Cole变换 $v=\ln\psi$ (即 $\psi=F/G)$ 进行线性化.性质(2.8)、(2.9)在推导孤子方程的Lax对时起到重要作用.

接着应用Bell多项式求解无穷守恒律时,只需作一个新的变量变换

$\begin{equation}\eta=\frac{(q_{x_{k}}' -q_{x_{k}})}{2}, \end{equation}$

(2.10)

这里 $q$ 和 $q'$ 满足 $q=w-v$ , $q' =w+v$ , $x_{k}$ 是 $x_{1}, \cdots, x_{l}$ 中任意随机变量.则由条件

$\begin{equation}C(q', q)=E(q' )-E(q)=0\end{equation}$

(2.11)

出发可以得到一组约束条件,只要把这一组约束条件化为Bell多项式以及关于某个变量的微分的形式,即一个Riccati方程

$\begin{equation}\eta_{x_{k}}+f(\eta)=0\end{equation}$

(2.12)

和一个离散型的方程

$\begin{equation}\partial_{x_{1}}F_{1}(\eta)+\cdots+\partial_{x_{l}}F_{l}(\eta)=0.\end{equation}$

(2.13)

然后,通过将级数展开式代入计算就可以从中导出非线性方程的守恒律.

3 一类广义浅水波KdV方程

令 $r(x, t)=u_x(x, t)$ ,方程(1.1)化成

$\begin{equation}u_{xt}-\alpha u_{xxxt}-4\alpha u_xu_{xt}-2\alpha u_{xx}u_{t}+ \alpha u_{xx}+\beta u_{xxxx}+6 \beta u_xu_{xx}=0, \end{equation}$

(3.1)

在本部分,通过Bell多项式的应用,给出方程(3.1)双线性表达式, Bäklund变换, Lax对和无穷守恒律.

3.1 双线性表达式

定理3.1 作 $u=q_x,$ 变换,方程(3.1)线性化为

$\begin{array}{l}&&(D_{x}^4+D_{x}D_{z})G\cdot G=0, \\&&[D_{x}D_{t}+\alpha(-\frac{2}{3}D_{t}D_{x}^{3}+\frac{1}{3}D_{t}D_{z}+D_{x}^{2})+\betaD_{x}^{4}]G\cdot G=0, \end{array}$

(3.2)

其中 $z$ 为辅助变量.

证为了得到方程(3.1)的双线性化表达式,首先引入一个辅助变量 $q$ ,且令

$\begin{equation}u=cq_{x}, \end{equation}$

(3.3)

其中 $c$ 是一个待选择的可以与方程(3.1) $P$ -多项式相关联的常数.将方程(3.3)代入方程(3.1),得到

$\begin{equation}q_{2x, t}+\alpha(-\frac{2}{3}q_{4x, t}-2cq_{3x}q_{xt}-4cq_{2x}q_{2x, t}-\frac{1}{3}q_{4x, t}+q_{3x})+\beta(q_{5x}+6cq_{2x}q_{3x})=0, \end{equation}$

(3.4)

并将其对 $x$ 积分一次可得

$\begin{equation}E(q)\equivq_{xt}+\alpha\big(-\frac{2}{3}(q_{3x, t}+3cq_{2x}q_{xt})-\frac{1}{3}\partial_x^{-1}\partial_t(q_{4x}+3cq_{2x}^2)+q_{2x}\big)+\beta(q_{4x}+3cq_{2x}^2)=0, \end{equation}$

(3.5)

应用性质(2.7),得 $c=1$ ,则方程(3.5)可化为

$\begin{equation}E(q)\equivq_{xt}+\alpha\big(-\frac{2}{3}(q_{3x, t}+3q_{2x}q_{xt})-\frac{1}{3}\partial_x^{-1}\partial_t(q_{4x}+3q_{2x}^2)+q_{2x}\big)+\beta(q_{4x}+3q_{2x}^2)=0.\end{equation}$

(3.6)

为了将方程(3.6)写成双线性形式,需要消除 $\partial_x^{-1}$ 积分项,为此引入一个辅助变量 $z$ ,并附加一个约束条件

$\begin{equation}q_{4x}+3q_{2x}^2=-q_{xz}, \end{equation}$

(3.7)

那么,方程(3.6)可化为

$\begin{equation}E(q)\equivq_{xt}+\alpha\big(-\frac{2}{3}(q_{3x, t}+3q_{2x}q_{xt})+\frac{1}{3}q_{tz}+q_{2x}\big)+\beta(q_{4x}+3q_{2x}^2)=0.\end{equation}$

(3.8)

依据性质(2.5),方程(3.7)与方程(3.8)可被转化成如下 $P$ -多项式组合形式

$\begin{array}{l}&&P_{4x}(q)+P_{xz}(q)=0, \\&&P_{xt}(q)+\alpha(-\frac{2}{3}P_{3x, t}(q)+\frac{1}{3}P_{tz}(q)+P_{2x}(q))+\beta P_{4x}(q)=0.\end{array}$

(3.9)

最后,由性质(2.6),在因变量变换

$\begin{equation}q=2\ln G\Longleftrightarrow u=cq_x=2(\ln G)_x\end{equation}$

(3.10)

作用下, (3.9)式可得方程(3.1)的双线性形式:

$\begin{array}{l}&&(D_{x}^4+D_{x}D_{z})G\cdot G=0, \\&&[D_{x}D_{t}+\alpha(-\frac{2}{3}D_{t}D_{x}^{3}+\frac{1}{3}D_{t}D_{z}+D_{x}^{2})+\betaD_{x}^{4}]G\cdot G=0. \end{array}$

(3.11)

证毕.

3.2 Bäcklund变换,孤子解与Lax对

定理3.2 假设 $F$ 为双线性方程(3.11)的一个解,如果 $G$ 满足

$\begin{array}{l}&&(D_{x}^2-\lambda)F\cdot G=0, \\&&[(1-3\alpha\lambda)D_{t}+(\alpha+3\beta\lambda)D_{x}-\alphaD_{x}^{2}D_{t}+\beta D_{x}^{3}]F\cdotG=0, \end{array}$

(3.12)

则 $G$ 即为方程(3.11)的另解,其中 $\lambda$ 为任意参数.

证令 $q$ 和 $q'$ 为方程(3.6)的两个不同的解

$\begin{equation} q=2\ln F, \qquad q' =2\ln G, \end{equation}$

(3.13)

相应地,引入两个新的变量

$\begin{equation}w=\frac{q' +q}{2}=\ln (FG), \qquadv=\frac{q' -q}{2}=\ln(\frac{F}{G}), \end{equation}$

(3.14)

则二场条件为

$\begin{array}{l}E(q' )-E(q) = E(w+v)-E(w-v)\\ = 2v_{xt}+\alpha[-2v_{3x, t}-4w_{2x}v_{x, t}-4w_{x, t}v_{2x} \\ -4\partial_x^{-1}(w_{2x}v_{2x, t}+w_{2x, t}v_{2x})+2v_{xx}]+\beta(2v_{4x}+12w_{xx}v_{xx})\\ = 2\partial_x[{\cal Y}_t(v)-\alpha{\cal Y}_{2x, t}{(v, w)}+\beta {\cal Y}_{3x}{(v, w)}]+R{(v, w)}=0, \end{array}$

(3.15)

其中

$\begin{eqnarray*}R{(v, w)} = \alpha\big(2\partial_x[(w_{2x}+v_x^2)v_t]-4w_{2x}v_{xt}+4w_{2x, t}v_x-4\partial_x^{-1}(w_{2x}v_{2x, t}+w_{2x, t}v_{2x})+2v_{xx}\big)\\ +6\beta(w_{2x}v_{xx}-v_xw_{3x}-v_x^2v_{2x}).\end{eqnarray*}$

这个二场条件可以认为是在适当限制条件下便于求得Bäcklund变换.

为了将二场条件(3.15)写成一对限制条件,可加入一个限制条件,则 $R(v, w)$ 可以表示成 ${\cal Y}$ -多项式的 $x$ -导数结合形式.可选

$\begin{equation}{\cal Y}_{2x}{(v, w)}=w_{2x}+v_x^2=\lambda, \end{equation}$

(3.16)

其中 $\lambda$ 为任意参数.由(3.16)式, $R(v, w)$ 可化为

$\begin{array}{l}R{(v, w)} = \alpha[2\lambdav_{xt}-4w_{2x}v_{xt}+4w_{2x, t}v_x-4\partial_x^{-1}(w_{2x}v_{2x, t}+w_{2x, t}v_{2x})+2v_{xx}]\\+6\beta(w_{2x}v_{xx}-v_xw_{3x}-v_x^2v_{2x})\\ = -6\alpha\lambda v_{xt}+2\alpha v_{2x}+6\beta\lambda v_{2x}, \end{array}$

(3.17)

此处应用 $w_{2x, t}=-2v_xv_{xt}$ 与 $w_{2x}=\lambda-v_x^2.$

结合(3.15)-(3.17)式,可得到 ${\cal Y}$ -多项式系统

$\begin{array}{l}&&{\cal Y}_{2x}(v, w)-\lambda=0, \\&&\partial_x[{(1-3\alpha\lambda)}{\cal Y}_{t}{(v)}+(\alpha+3\beta\lambda){\cal Y}_{x}{(v)}-\alpha{\cal Y}_{2x, t}{(v, w)}+\beta{\cal Y}_{3x}{(v, w)}]=0, \end{array}$

(3.18)

其中第二个方程保留导数形式,这在之后的求解方程的守恒律时是非常重要.应用性质(2.5),从方程(3.15)可立即得到双线性Bäcklund变换

$\begin{array}{l}&&(D_{x}^2-\lambda)F\cdot G=0, \\&&[(1-3\alpha\lambda)D_{t}+(\alpha+3\beta\lambda)D_{x}-\alphaD_{x}^{2}D_{t}+\beta D_{x}^{3}]F\cdotG=0, \end{array}$

(3.19)

其中我们对系统(3.18)中第二个方程 $x$ 进行积分.

通过此Bäcklund变换,我们可以很容易的求出其孤子解.接下来我们以一孤子解与二孤子解为例.

从平凡解 $u(x, t)=0$ 出发,选择 $F=1$ ,代入方程(3.19),可得方程(3.19)的另解,令 $\lambda=\frac{k_1^2}{4}$ ,解得

$\begin{equation}G_1=e^{\frac{\xi_1}{2}}+e^{-\frac{\xi_1}{2}}, \quad\xi_1=k_1x-\frac{\alpha k_1+\beta k_1^3}{1-\alphak_1^2}t+\xi_1^{(0)}, \end{equation}$

(3.20)

其中 $k_1, \xi_1^{(0)}$ 为任意实数.方程(3.20)化为

$\begin{equation}G_1=e^{-\frac{\xi_1}{2}}(1+e^{\xi_1}), \end{equation}$

(3.21)

依据(3.13)式有 $q'=2\ln G$ , (3.3)式有 $u=q_x$ , (3.1)式有 $r=u_x$ 得一孤子解为

$\begin{equation}r=2\ln[1+e^{k_1x-\frac{\alpha k_1+\beta k_1^3}{1-\alphak_1^2}t+\xi_1^{(0)}}]_{xx}, \end{equation}$

(3.22)

再令 $F=e^{\frac{\xi_1}{2}}+e^{-\frac{\xi_1}{2}}$ ,可得解

$\begin{equation}G_2=(k_1-k_2)(e^{\frac{\xi_1+\xi_2}{2}}+e^{-\frac{\xi_1+\xi_2}{2}})-(k_1+k_2)(e^{\frac{\xi_1-\xi_2}{2}}+e^{-\frac{\xi_1-\xi_2}{2}}), \quad\lambda=\frac{k_2^2}{4}, \end{equation}$

(3.23)

其中 $\xi_j=k_j x-\frac{\alpha k_j+\beta k_j^3}{1-\alpha k_j^2}t+\xi_j^{(0)}, j=1, 2.$

作变换 $\xi_j=\eta_j+\frac{1}{2}A_{12}$ 其中 $e^{A_{12}}=\frac{(k_1-k_2)^2}{(k_1+k_2)^2}, j=1, 2$ ,则(3.23)式可化为

$\begin{equation}G_2=1+e^{\eta_1}+e^{\eta_2}+e^{\eta_1+\eta_2+A_{12}}, \end{equation}$

(3.24)

其对应方程的二孤子解

$\begin{equation}r=2\ln [1+e^{\eta_1}+e^{\eta_2}+e^{\eta_1+\eta_2+A_{12}}]_{xx}.\end{equation}$

(3.25)

接下来求解三孤子解,将 $F=(k_1-k_2)(e^{\frac{\xi_1+\xi_2}{2}}+e^{-\frac{\xi_1+\xi_2}{2}})-(k_1+k_2)(e^{\frac{\xi_1-\xi_2}{2}}+e^{-\frac{\xi_1-\xi_2}{2}})$ 代入方程组(3.19),可得

$\begin{array}{l}G_3 = a(e^{\frac{\xi_1+\xi_2+\xi_3}{2}}+e^{-(\frac{\xi_1+\xi_2+\xi_3}{2}})+b(e^{\frac{-\xi_1+\xi_2+\xi_3}{2}}+e^{-\frac{-\xi_1+\xi_2+\xi_3}{2}})\\+c(e^{\frac{\xi_1-\xi_2+\xi_3}{2}}+e^{-\frac{\xi_1-\xi_2+\xi_3}{2}})+d(e^{\frac{\xi_1+\xi_2-\xi_3}{2}}+e^{-\frac{\xi_1+\xi_2-\xi_3}{2}}), \end{array}$

(3.26)

其中

$\begin{eqnarray*}\label{a1}&&a=(k_1-k_2)(k_1-k_3)(k_2-k_3), \quad b =(k_1+k_2)(k_1+k_3)(k_2-k_3), \\&&c=(k_1+k_2)(k_1-k_3)(k_2+k_3), \quad d=(k_1-k_2)(k_1+k_3)(k_2+k_3)\end{eqnarray*}$

和 $\xi_j=k_j x-\frac{\alpha k_j+\beta k_j^3}{1-\alpha k_j^2}t+\xi_j^{(0)}\ (j=1, 2, 3)$ .

以此类推,得出 $N$ -孤子解为

$\begin{array}{l}&&G_n=\sum\limits_{\epsilon={\pm1}}\prod \limits_{1\leqj <l}^{n}\epsilon_l(\epsilon_jk_j-\epsilon_lk_l)e^{\frac{1}{2}\sum \limits _{j=1}^{n}\epsilon_j\xi_j}, \\&&\xi_j=k_j x-\frac{\alpha k_j+\beta k_j^3}{1-\alphak_j^2}t+\xi_j^{(0)}\ (j=1, 2, \cdots, n), \quad\lambda=\frac{k_n^2}{4}, \end{array}$

(3.27)

其中 $\sum \limits_{\epsilon={\pm1}}$ 表示的是 $\epsilon _j=1, -1\ (j=1, 2, \cdots, n)$ 的所有可能形式.

接下来,我们求解方程(3.1)的Lax对.

定理3.3 方程(3.1) Lax对为

$\begin{array}{l}&&L\psi=\psi_{xx}+(u_{x}-\lambda)\psi=0, \\&&M\psi=(1-3\alpha\lambda-u_{x})\psi_{t}+(\alpha+3\beta\lambda+3\betau_{x}-2\alpha u_{t})\psi_{x}-\alpha\psi_{xxt}+\beta\psi_{xxx}=0, \end{array}$

(3.28)

其中 $\lambda$ 任意参数.

证由Hopf-Cole变换 $v=\ln \psi$ ,由公式(2.8)和(2.9)可得

$\begin{array}{l}&&{\cal Y}_{x}(v)=\frac{\psi_{x}}{\psi}, \qquad{\cal Y}_{t}=\frac{\psi_{t}}{\psi}, \qquad{\cal Y}_{2x}(v, w)=q_{2x}+\frac{\psi_{2x}}{\psi}, \\&&{\cal Y}_{3x}(v, w)=3q_{2x}\frac{\psi_{x}}{\psi}+\frac{\psi_{3x}}{\psi}, \qquad{\cal Y}_{2x, t}(v, w)=2q_{xt}\frac{\psi_{x}}{\psi}+q_{2x}\frac{\psi_{t}}{\psi}+\frac{\psi_{2x, t}}{\psi}.\end{array}$

(3.29)

结合(3.29)式, (3.19)式可被线性化为带有参数 $\lambda$ 的系统

$\begin{array}{l}L\psi=0, M\psi=0, \end{array}$

(3.30)

其中

$\begin{array}{l}&&L=\partial_{x}^2+q_{2x}-\lambda, \\&&M=(1-3\alpha\lambda-q_{2x})\partial_{t}+(\alpha+3\beta\lambda+3\betaq_{2x}-2\alpha q_{xt})\partial_{x}-\alpha\partial_t\partial_{x}^2+\beta \partial_{x}^3, \end{array}$

(3.31)

用 $q_{x}$ 代替 $u$

(3.32)

得出Lax对,且容易验证相容性条件

$\begin{equation} [L, M]=q_{2x, t}+\alpha(-q_{4x, t}-4 q_{2x}q_{2x, t}-2 q_{3x}q_{xt}+q_{3x})+\beta(q_{5x}+6q_{2x}q_{3x})=0. \end{equation}$

(3.33)

证毕.

3.3 无穷守恒律

定理3.4 方程(3.1)具有如下无穷守恒律

$\begin{equation}I_{n, t}+F_{n, x}=0, n=1, 2, 3, \cdots , \end{equation}$

(3.34)

其中守恒密度 $I_n$ 递推公式为

$$\begin{array}{l}&&I_{1}=-\frac{1}{2}u_{x}, \qquadI_{2}=-\frac{1}{2}I_{1, x}=\frac{1}{4}u_{2x}, \\& &\cdots\cdots\\&&I_{n}=-\frac{1}{2}\bigg[I_{n-1, x}+\sum\limits_{k=1}^nI_{k}I_{n-1-k}\bigg], \ n=2, 3, 4, \cdots. \end{array}

(3.35)

连带流 $F_n$ 递推关系为

$\begin{array}{l}&&F_{1}=-\frac{1}{2}[\alpha(u_{x}-u_{xxt}-2\partial_x^{-1}(u_xu_{xt})-2u_tu_x)+\beta(u_{3x}+3u_x^2)], \\&&\cdots\cdots\\&&F_{n}=\alpha\bigg[I_n- I_{n, xt}+4\partial_x^{-1}\bigg(\sum\limits_{k=1}^{n}I_{k}I_{n+1-k, t}\bigg)+4\partial_x^{-1}\bigg(\sum\limits_{k=1}^nI_{k, t}\bigg)I_{n+1-k}\\&& +4\partial_x^{-1}\bigg(\sum\limits_{i+j+k=n}I_{i}I_{j, t}I_{k}\bigg)\bigg]\\&& +\beta\bigg(I_{n, 2x}-2\sum\limits_{i+j+k=n}I_{i}I_{j}I_{k}-6\sum\limits_{k=1}^nI_kI_{n+1-k}\bigg), \ n=2, 3, 4, \cdots .\end{array}$

(3.36)

证首先从公式(3.19)出发,由关系 $\partial_x{\cal Y}_t(v)=\partial_t{\cal Y}_x(v)=v_{xt}$ ,则(3.19)式可化为

$\begin{array}{l}&&{\cal Y}_{2x}(v, w)-\lambda=0, \\&&\partial_t(1-3\alpha\lambda){\cal Y}_{x}+\partial_x[{(\alpha+3\beta\lambda){\cal Y}_{x}{(v)}-\alpha{\cal Y}_{2x, t}{(v, w)}+\beta{\cal Y}_{3x}{(v, w)}}]=0.\end{array}$

(3.37)

引入一个势函数

$\begin{equation}\eta=\frac{q_{x}' -q_{x}}{2}, \end{equation}$

(3.38)

那么从关系式(3.15)可得到

$\begin{equation}v_{x}=\eta, w_{x}=q_{x}+\eta, \end{equation}$

(3.39)

将(3.39)式代入(3.37)式,可得到一个Riccati -类型的方程

$\begin{equation}\eta_{x}+\eta^{2}+q_{2x}=\lambda\end{equation}$

(3.40)

和一个离散型方程

$\begin{equation}(1-4\alpha\lambda)\eta_{t}+\partial_{x}\{(\alpha+6\beta\lambda)\eta-\alpha\eta_{xt}+4\alpha[\partial_x^{-1}(\eta\eta_t)]\eta+\beta\eta_{2x}-2\beta\eta^3\}=0, \end{equation}$

(3.41)

其中我们取 $\lambda=\varepsilon^2$ .

进一步,在 $\eta=\varepsilon+\tilde{\eta}$ 变换下, (3.40)式和(3.41)式可化为

$\begin{equation}\tilde{\eta}_x +\tilde{\eta}^2+2\varepsilon\tilde{\eta}+q_{xx}=0\end{equation}$

(3.42)

和

$\begin{equation}\tilde{\eta}_t+\partial_x\{\alpha\tilde{\eta}-\alpha\tilde{\eta}_{xt}+4\alpha\partial_x^{-1}(\tilde{\eta}\tilde{\eta}_t)\varepsilon+4\alpha\partial_x^{-1}(\varepsilon\tilde{\eta}_t+\tilde{\eta}\tilde{\eta}_t)\tilde{\eta}+\beta\tilde{\eta}_{2x}-2\beta\tilde{\eta}^3-6\beta\varepsilon\tilde{\eta}^2\}=0.\end{equation}$

(3.43)

将级数展开式

$\begin{equation}\tilde{\eta}=\sum\limits_{n=1}^{\infty}I_{n}(q, q_{x}, \cdots)\varepsilon^{-n}\end{equation}$

(3.44)

代入(3.42)式且令上式中 $\varepsilon$ 的幂次系数相等,就可以得到守恒密度 $I_{n}$ 的递推关系式

$\begin{array}{l}&&I_{1}=-\frac{1}{2}u_{x}, \qquadI_{2}=-\frac{1}{2}I_{1, x}=\frac{1}{4}u_{2x}, \\&&I_{3}=-\frac{1}{2}(I_{2, x}+I_{1}^{2})=-\frac{1}{8}(u_{xxx}+u_x^2), \\&&I_{4}=-\frac{1}{2}(I_{3, x}+2I_{1}I_{2})=\frac{1}{16}(u_{4x}+4u_xu_{2x}), \\&&\cdots\cdots\\&&I_{n}=-\frac{1}{2}[I_{n-1, x}+\sum\limits_{k=1}^nI_{k}I_{n-1-k}], n=2, 3, 4, \cdots .\end{array}$

(3.45)

再将展开式(3.44)代入(3.43)式,则有

$\begin{array}{l}&&\sum\limits_{n=1}^{\infty}I_{n, t}\varepsilon^{-n}+\partial_{x}\bigg[\alpha\sum\limits_{n=1}^{\infty}I_{n}\varepsilon^{-n}-\alpha\sum\limits_{n=1}^{\infty}I_{n, xt}\varepsilon^{-n}+4\alpha\partial_x^{-1}\bigg(\sum\limits_{n=1}^{\infty}I_{n}\varepsilon^{-n}\sum\limits_{n=1}^{\infty}I_{n, t}\varepsilon^{-n}\bigg)\cdot\varepsilon\\&&+4\alpha\partial_x^{-1}\bigg(\varepsilon\cdot\sum\limits_{n=1}^{\infty}I_{n, t}\varepsilon^{-n}+\sum\limits_{n=1}^{\infty}I_{n}\varepsilon^{-n}\sum\limits_{n=1}^{\infty}I_{n, t}\varepsilon^{-n}\bigg)\sum\limits_{n=1}^{\infty}I_{n}\varepsilon^{-n}+\beta\sum\limits_{n=1}^{\infty}I_{n, 2x}\varepsilon^{-n}\\&&-2\beta\bigg(\sum\limits_{n=1}^{\infty}I_{n}\varepsilon^{-n}\bigg)^{3}-6\beta\varepsilon\bigg(\sum\limits_{n=1}^{\infty}I_{n}\varepsilon^{-n}\bigg)^2\bigg]=0, \end{array}$

(3.46)

即为无穷守恒律

$\begin{equation}I_{n, t}+F_{n, x}=0, n=1, 2, 3, \cdots, \end{equation}$

(3.47)

在无穷守恒律(3.47)中, (3.45)式给出守恒密度 $I_{n}$ ,连带流 $F_n$ 即为

$\begin{array}{l}&&F_{1}=-\frac{1}{2}[\alpha(u_{x}-u_{xxt}-2\partial_x^{-1}(u_xu_{xt})-2u_tu_x)+\beta(u_{3x}+3u_x^2)], \\&&F_{2}=\frac{1}{4}[\alpha(u_{2x}-u_{3x, t}-4(u_xu_{xt})-2u_tu_{2x})+\beta(u_{4x}+6u_xu_{2x})], \\&&\cdots\cdots\\&&F_{n}=\alpha\bigg[I_n- I_{n, xt}+4\partial_x^{-1}\bigg(\sum\limits_{k=1}^{n}I_{k}I_{n+1-k, t}\bigg)+4\partial_x^{-1}\bigg(\sum\limits_{k=1}^nI_{k, t}\bigg)I_{n+1-k}\\&& +4\partial_x^{-1}\bigg(\sum\limits_{i+j+k=n}I_{i}I_{j, t}I_{k}\bigg)\bigg]\\&& +\beta\bigg(I_{n, 2x}-2\sum\limits_{i+j+k=n}I_{i}I_{j}I_{k}-6\sum\limits_{k=1}^nI_kI_{n+1-k}\bigg), \ n=2, 3, 4, \cdots .\end{array}$

(3.48)

据此可以验证(3.47)式的第一个方程即为方程(3.1).

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

Hirota

The Direct Method in Soliton Theory. Cambridge: Cambridge University Press, 2004

[本文引用: 3]

[2]

Guo

Y C

Nonlinear Partial Differential Equations. Beijing: Tsinghua University Press, 2008

[3]

W X

, Wu

H Y

Solving the Korteweg-De Vries equation by its bilinear form: Wronskian solutions

Transactions of the American Mathematical Society, 2004, 375 (5): 1753- 1778

[本文引用: 2]

[4]

Wazwaz

A M

$N$ -soliton solutions for shallow water waves equations in $(1+1$ ) and $(2+1)$ -dimensions

Appiedl Mathematics Computation, 2011, 217: 8840- 8845

DOI:10.1016/j.amc.2011.03.048 [本文引用: 3]

[5]

Zhang

, Cheng

Z L

Exact soliton solutions and quasi-periodic wave solutions to the forced vc-KdV equation

Intern J Mod Phys B, 2012, 26: 1250072

DOI:10.1142/S0217979212500725 [本文引用: 2]

[6]

Zhang

, Cheng

D Y

Bäcklund transformation and soliton solutions for the shallow water waves equation

Chaos, Solitons and Fractals, 2004, 20: 343- 351

DOI:10.1016/S0960-0779(03)00394-1 [本文引用: 1]

[7]

Cheng

Z L

, Hao

X H

The periodic wave solutions for a $(2 + 1)$ -dimensional AKNS equation

Appl Math Comput, 2014, 234: 118- 126

[8]

Zhang

, Cheng

Z L

, Hao

X H

Riemann theta functions periodic wave solutions for the variable-coefficient mKdV equation

Chin Phys B, 2012, 21: 120203

DOI:10.1088/1674-1056/21/12/120203 [本文引用: 1]

[9]

郝晓红, 程智龙.

$(2+1)$ 维AKNS方程的可积性研究

动力学与控制学报, 2018, 16 (3): 201- 205

[本文引用: 4]

Hao

X H

, Cheng

Z L

Research on integrability of a $(2 + 1)$ -dimensional AKNS equation

Journal of Dynamics and Control, 2018, 16A (3): 201- 205

[本文引用: 4]

[10]

Hirota

A new form of Bäcklund transformations and its relation to the inverse scattering problem

Prog Theor Phys, 1974, 52: 1498- 1512

DOI:10.1143/PTP.52.1498 [本文引用: 1]

[11]

W X

Bilinear equations and resonant solutions characterized by Bell polynomials

Reports on Mathematical Physics, 2013, 72: 41- 56

DOI:10.1016/S0034-4877(14)60003-3

[12]

Fan

E G

, Chow

K W

Darboux covariant Lax pairs and infinite conservation laws of the $(2+1)$ -dimensional breaking soliton equation

J Math Phys, 2011, 52: 023504

DOI:10.1063/1.3545804 [本文引用: 2]

[13]

Ablowitz

M J

, Kaup

D J

, Newell

A C

, Segur

The inverse scattering transform: fourier analysis nonlinear problems

Stud Appl Math, 1974, 53: 249- 315

DOI:10.1002/sapm1974534249 [本文引用: 1]

[14]

Gilson

C R

, Lambert

, Nimmo

J J C

, Willox

On the combinatorics of the Hirota $D$ -operators

Proc R Soc Lond A, 1996, 452: 223- 234

DOI:10.1098/rspa.1996.0013 [本文引用: 2]

[15]

Cheng

Z L

, Hao

X H

The travelling wave solutions for a $(2 + 1)$ -dimensional AKNS equation

Chinese Quarterly Journal of Mathematics, 2015, 3 (30): 323- 329

[本文引用: 1]

[16]

W X

, Lee

J H

A transformed rational function method and exact solutions to the 3+1 dimensional Jimbo-Miwa equation

Chaos, Solitons and Fractals, 2009, 42: 1356- 1363

DOI:10.1016/j.chaos.2009.03.043

[17]

W X

Lump solutions to the Kadomtsev-Petviashvili equation

Phys Lett A, 2015, 379 (36): 1975- 1978

DOI:10.1016/j.physleta.2015.06.061 [本文引用: 1]

[18]

Gilson

C R

, Nimmo

J J C

, Willox

R A

$(2+1)$ -dimensional generalization of the AKNS shallow water wave equation

Phys Lett A, 1993, 180: 337- 345

DOI:10.1016/0375-9601(93)91187-A [本文引用: 1]

[19]

Shang

Y D

, Huang

Bäcklund transformations and abundant explicit exact solutions for the AKNS-SWW equation

Commun Nonlinear Sci Numer Simulat, 2011, 16: 2445- 2455

DOI:10.1016/j.cnsns.2010.09.024 [本文引用: 3]

2004

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

... 即为KdV方程^[1-3]; ...

... 即为浅水波方程^{[1, 4]}. Ablowitz等用反散射法研究了此类方程^[13]. Gilson, Nimmo和Willox考虑其Wronskian形式的孤子解^[18]. Shang和Hong通过其次平衡法与扩展双曲函数法得到大量的精确解与Bäcklund变换形式^[19], Wazwaz则应用简化形式的Hirota双线性法求解出

$N$

-孤子解^[4]. ...

2008

Solving the Korteweg-De Vries equation by its bilinear form: Wronskian solutions

2004

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

... 即为KdV方程^[1-3]; ...

$N$

-soliton solutions for shallow water waves equations in

$(1+1$

) and

$(2+1)$

-dimensions

2011

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

$N$

-孤子解^[4]. ...

... [4]. ...

Exact soliton solutions and quasi-periodic wave solutions to the forced vc-KdV equation

2012

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

... [5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

B?cklund transformation and soliton solutions for the shallow water waves equation

2004

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

The periodic wave solutions for a

$(2 + 1)$

-dimensional AKNS equation

2014

Riemann theta functions periodic wave solutions for the variable-coefficient mKdV equation

2012

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

$(2+1)$

维AKNS方程的可积性研究

2018

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

... -9, 15-17, 19]越来越受到一些专家学者的关注. ...

... 1996年, Lambert, Gilson和Nimmo建立了Bell多项式与Hirota双线性算子之间的关系^[14],通过转换关系得到双线性Bäcklund变换,这个方法很有效地避免了求Bäcklund变换过程中使用交换公式繁琐的计算,简洁实用.并且直接对其做线性化还可以得到方程的Lax对^{[9, 12]},最后将级数展开式代入计算求解得到方程的无穷守恒律. ...

... 为Bell多项式^{[9, 14]}.其中

$f_{r_{1}x_{1}, \cdots, r_{l}x_{l}}=\partial_{x_{1}}^{n_{1}}\cdots\partial_{x_{l}}^{n^{l}}f (r_{1}=0, \cdots, n_{1};\cdots; r_{l}=0, \cdots, n_{l}.)$

.当

$f=f(x, t)$

时,对应的2-维Bell多项式,举例为 ...

$(2+1)$

维AKNS方程的可积性研究

2018

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

... -9, 15-17, 19]越来越受到一些专家学者的关注. ...

... 为Bell多项式^{[9, 14]}.其中

$f_{r_{1}x_{1}, \cdots, r_{l}x_{l}}=\partial_{x_{1}}^{n_{1}}\cdots\partial_{x_{l}}^{n^{l}}f (r_{1}=0, \cdots, n_{1};\cdots; r_{l}=0, \cdots, n_{l}.)$

.当

$f=f(x, t)$

时,对应的2-维Bell多项式,举例为 ...

A new form of B?cklund transformations and its relation to the inverse scattering problem

1974

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

Bilinear equations and resonant solutions characterized by Bell polynomials

2013

Darboux covariant Lax pairs and infinite conservation laws of the

$(2+1)$

-dimensional breaking soliton equation

2011

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

The inverse scattering transform: fourier analysis nonlinear problems

1974

$N$

-孤子解^[4]. ...

On the combinatorics of the Hirota

$D$

-operators

1996

... 为Bell多项式^{[9, 14]}.其中

$f_{r_{1}x_{1}, \cdots, r_{l}x_{l}}=\partial_{x_{1}}^{n_{1}}\cdots\partial_{x_{l}}^{n^{l}}f (r_{1}=0, \cdots, n_{1};\cdots; r_{l}=0, \cdots, n_{l}.)$

.当

$f=f(x, t)$

时,对应的2-维Bell多项式,举例为 ...

The travelling wave solutions for a

$(2 + 1)$

-dimensional AKNS equation

2015

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

A transformed rational function method and exact solutions to the 3+1 dimensional Jimbo-Miwa equation

2009

Lump solutions to the Kadomtsev-Petviashvili equation

2015

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

$(2+1)$

-dimensional generalization of the AKNS shallow water wave equation

1993

$N$

-孤子解^[4]. ...

B?cklund transformations and abundant explicit exact solutions for the AKNS-SWW equation

2011

... 近几十年以来,通过双线性方法^{[1-3, 5]}, Bäcklund变换^{[6-10, 19]},黎曼

$\theta$

函数法^[5-8]等研究非线性偏微分方程中可积性^{[9, 12]}与精确解^{[4-9, 15-17, 19]}越来越受到一些专家学者的关注. ...

... , 19]越来越受到一些专家学者的关注. ...

$N$

-孤子解^[4]. ...

〈

〉