Loading [MathJax]/jax/output/HTML-CSS/jax.js

数学物理学报, 2021, 41(2): 326-335 doi:

论文

常系数无穷维Hamilton系统的二阶循环算子结构

耿万鹏, 任文秀,, 程意苏,

Second Order Recursion Operator Structure of Constant Coefficient Infinite Dimension Hamiltonian System

Geng Wanpeng, Ren Wenxiu,, Cheng Yisu,

通讯作者: 任文秀, E-mail: renwenxiu2003@hotmail.com

收稿日期: 2020-01-9  

Received: 2020-01-9  

作者简介 About authors

程意苏,E-mail:1070448132@qq.com , E-mail:1070448132@qq.com

Abstract

By virtue of limited and formal constant coefficient Hamiltonian operator, it applies the method of general system recursion operator to Hamiltonian canonical system of infinite dimensional form. As to the result, the general structure of the recursion operator allowed by the next-order constant coefficient Hamiltonian operator under constraint condition and specific form of its coefficient are obtained. And then, it verifies the correctness and convenience of the conclusion by means of calculating example.

Keywords: Constant coefficient Hamiltonian operator ; Infinite dimensional linear Hamiltonian canonical system ; Recursion operator

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

本文引用格式

耿万鹏, 任文秀, 程意苏. 常系数无穷维Hamilton系统的二阶循环算子结构. 数学物理学报[J], 2021, 41(2): 326-335 doi:

Geng Wanpeng, Ren Wenxiu, Cheng Yisu. Second Order Recursion Operator Structure of Constant Coefficient Infinite Dimension Hamiltonian System. Acta Mathematica Scientia[J], 2021, 41(2): 326-335 doi:

1 引言

Olver[1]于1977年首次给出了循环算子的概念, 并指出循环算子就是把一个对称映射成另一个对称的线性算子. 之后的研究成果表明, 循环算子是可积系统的一个重要性质, 它与无穷多守恒律、Hamilton结构和方程的谱系等有着密切关系. 因此, 许多学者开始对循环算子的构造、性质及应用进行研究[2-8], 如Bluman, Gürses等. 通过文献查阅, 我们发现关于循环算子的研究成果多数是在一般体系下展开的工作, 而在Hamilton体系下的研究却很少, 且研究思路采用Hamilton算子为四个子块的形式, 因此我们的想法是对该体系下的循环算子及其形式在其它路径下重新进行剖析.

在微分方程中, 循环算子通常用有限和形式表示. 受其启发, 不同于其它文献的设定, 我们把Hamilton算子也写成类似的求和形式, 再结合文献[1]获得循环算子的方法, 将线性微分算子与循环算子之间的关系, 移植到无穷维线性Hamilton正则系统下, 使我们在处理Hamilton算子所允许循环算子的问题上可能实现统一化. 本文根据以上思想, 主要围绕以下两个方面进行展开: Hamilton算子与循环算子结构上的联系; 常系数Hamilton体系下所允许循环算子的统一形式. 具体工作是通过约化确定方程组, 得到了约束条件下一阶常系数Hamilton算子所允许循环算子的一般结构. 再利用Hamilton算子固有元素分情况进行讨论, 得到了一阶常系数Hamilton算子所允许循环算子的具体形式, 并对二阶常系数Hamilton算子所允许循环算子的结构进行探讨. 同时, 给出了几个算例.

2 研究思路

2.1 研究对象

考虑无穷维线性Hamilton正则系统具有以下形式

Δ=(IDxH)w=0,
(2.1)

其中Hamilton算子

H=(ABCA),
(2.2)

状态变量w=(w1,w2,,w2n)T, 且子块满足B=B,C=C. 利用文献[9] 中的引理2.1及注释2.3可知, (2.1)式中的H为常系数Hamilton算子时, 可得到(2.1)式允许的循环算子和线性微分算子的关系

[DΔ,R]=[IDxH,R]=0.
(2.3)

我们受循环算子有限和写法的启发, 将Hamilton算子(2.2)式重新约定为微分算子求和形式, 即

H=Pi(x,t,u)Dit=PiDit.

这不同于其它文献对于Hamilton算子(2.2)式的设定. 下面我们试图在这种形式下得到关于循环算子更多有效的信息.

由于Hamilton算子阶数的升高和变量的增多, 循环算子的计算过程越来越困难, 所以本文以一阶常系数Hamilton算子

H=P1Dt+P0,(P10)
(2.4)

为例, 讨论其所对应的方程(2.1)n=1时, 记w=(w1,w2)T, 所允许循环算子的结构及其子块系数. 在此约定下, 显然(2.3) 式中的DΔ可展开为

DΔ=IDxH=IDxP0P1Dt.

利用文献[10], 相应的系数Pi(i=0,1)的形式为

P0=(b(0)2b(0)1b(0)0b(0)2),P1=(b(1)200b(1)2),

其中Pi中各个元素都是常数.同时, 记β0=(b(0)2)2+b(0)1b(0)0,β1=(b(1)2)2.

2.2 解决思路

经查阅文献, 发现方程允许的循环算子结构有一些固定的特点, 如主对角相等、负对角相等以及呈现上(下) 三角矩阵等形式, 故我们约定二阶循环算子的一般形式如下

R=(R1R2R3R4),
(2.5)

其中Ri=2m=0a(i)mDmx (i=1,,4),且全部系数a(i)m是关于x,t的函数. 根据(2.3) 式, 可得循环算子满足的确定方程为

[DΔ,R]w=(IDxP0P1Dt)RwR(IDxP0P1Dt)w=((Dx00Dx)(b(0)2b(0)1b(0)0b(0)2)(b(1)200b(1)2)Dt)(a(1)2D2x+a(1)1Dx+a(1)0a(2)2D2x+a(2)1Dx+a(2)0a(3)2D2x+a(3)1Dx+a(3)0a(4)2D2x+a(4)1Dx+a(4)0)(w1w2)(a(1)2D2x+a(1)1Dx+a(1)0a(2)2D2x+a(2)1Dx+a(2)0a(3)2D2x+a(3)1Dx+a(3)0a(4)2D2x+a(4)1Dx+a(4)0)((Dx00Dx)(b(0)2b(0)1b(0)0b(0)2)(b(1)200b(1)2)Dt)(w1w2)=0.

对上述式子进行展开, 可得以下方程组

2m=0((Dxa(1)mb(0)1a(3)m+b(0)0a(2)mb(1)2Dta(1)m)w1mx+(Dxa(2)m2b(0)2a(2)m+b(0)1(a(1)ma(4)m)b(1)2Dta(2)m)w2mx)=0;2m=0((Dxa(3)m+2b(0)2a(3)m+b(0)0(a(4)ma(1)m)b(1)2Dta(3)m)w1mx+(Dxa(4)mb(0)0a(2)m+b(0)1a(3)mb(1)2Dta(4)m)w2mx)=0,
(2.6)

其中, 我们约定wi2x=wixx,wi1x=wix,wi0x=wi (i=1,2). 再用所考虑的正则方程及其求导方程, 并统一变量, 将以上方程组约化为关于w1,w2及其各阶导数的微分方程组. 经求解, 可得循环算子的结构以及全部子块系数a(i)m. 下面是根据一阶常系数Hamilton算子与循环算子之间的联系得到的一些结论.

3 主要工作

3.1 约束条件下一阶系统的循环算子的一般结构

本小节主要讨论无穷维线性Hamilton正则系统下, 对循环算子中元素的自变量做限制, 即子块Ri中的系数a(i)m仅关于自变量xt的函数(二者计算方法平行). 下面我们以其系数仅含自变量x的函数为例, 得到约束条件下一阶常系数Hamilton算子所允许循环算子的信息, 具体结论如下.

命题3.1   若H为形如(2.4) 式的一阶常系数Hamilton算子, 则方程(2.1)n=1时, 允许循环算子的结构为

R=(R1R2R3R1+2m=0cmDmx),
(3.1)

其中cm为任意常数.

   由H为形如(2.4) 式的一阶常系数Hamilton算子可知P10,b(1)20.用相应的正则方程及其求导方程进行约化, 并统一变量, 然后整理方程组(2.6), 得到关于wi2t,wit,wi(i=1,2)的多项式. 根据它们的线性无关性, 令其系数为零, 可得12个方程. 为了书写方便, 下式中的a(4)m用符号a(0)m(m=0,1,2)替代. 当i=1,2时, 方程如下

β1(Dxa(2i)2+(1)i(b(0)1a(3)2b(0)0a(2)2))=0;(1)i+12b(1)2b(0)2(Dxa(2i)2+(1)i(b(0)1a(3)2b(0)0a(2)2))+b(1)2(Dxa(2i)1+(1)i(b(0)1a(3)1b(0)0a(2)1))+2b(1)2b(0)i1(Dxa(i+1)2+(1)i(2b(0)2a(i+1)2b(0)2i(a(1)2a(0)2)))=0;β0(Dxa(2i)2+(1)i(b(0)1a(3)2b(0)0a(2)2))+(1)i+1b(0)2(Dxa(2i)1+(1)i(b(0)1a(3)1b(0)0a(2)1))+b(0)i1(Dxa(i+1)1+(1)i(2b(0)2a(i+1)1b(0)2i(a(1)1a(0)1)))+(Dxa(2i)0+(1)i(b(0)1a(3)0b(0)0a(2)0))=0;β1(Dxa(4i)2+(1)i(b(0)i1(a(1)2a(0)2)2b(0)2a(4i)2))=0;(1)i+12b(1)2b(0)2(Dxa(4i)2+(1)i(b(0)i1(a(1)2a(0)2)2b(0)2a(4i)2))+b(1)2(Dxa(4i)1+(1)i(b(0)i1(a(1)1a(0)1)2b(0)2a(4i)1))+2b(1)2b(0)i1(Dxa(i1)2+(1)i(b(0)0a(2)2b(0)1a(3)2))=0;β0(Dxa(4i)2+(1)i(b(0)i1(a(1)2a(0)2)2b(0)2a(4i)2))+(1)i+1b(0)2(Dxa(4i)1+(1)i(b(0)i1(a(1)1a(0)1)2b(0)2a(4i)1))+b(0)i1(Dxa(i1)1+(1)i(b(0)0a(2)1b(0)1a(3)1))+(Dxa(4i)0+(1)i(b(0)i1(a(1)0a(0)0)2b(0)2a(4i)0))=0.
(3.2)

上述方程组(3.2) 的后6个方程在计算R2R3系数之间明确的关系时存在困难, 所以只给出约化后的前6个方程如下

Dxa(2i)m+(1)ib(0)1a(3)m+(1)i+1b(0)0a(2)m=0(i=1,2;m=0,1,2).
(3.3)

为了得到循环算子子块系数之间的关系, 我们将这6个方程分为三组(m=0,1,2).i=1,2时, 固定m=2, 将对应(3.3) 式得到的两个方程相加, 得

Dx(a(1)2+a(4)2)=0,

a(1)2+a(4)2=c2.

同理, 我们可以分别固定m=0,1, 由相应的方程(3.3)式也可推出如下关系

a(1)m+a(4)m=cm,

R1+R4=2m=0cmDmx.

其中cm为任意常数, 即得命题3.1中循环算子的结构.

值得注意的是, 工作初期如果假定循环算子子块系数是关于x,t的函数, 那么通过上述计算过程, 可得

(Dxb(1)2Dt)(a(1)m+a(4)m)=0(m=0,1,2).

解之, 有

a(1)m+a(4)m=φm(b(1)2x+t).

R1+R4=2m=0φm(b(1)2x+t)Dmx,

其中φm是一个关于x,t的任意函数, 其物理意义值得我们进一步探讨.

3.2 约束条件下一阶系统的循环算子的具体形式

本小节根据命题3.1中循环算子的结构, 从而得到一阶常系数Hamilton系统下循环算子的统一形式.

命题3.2   若H为形如(2.4) 式的一阶常系数Hamilton算子, 对于方程的任意一组解w, 则方程(2.1)n=1时, 允许的形如(3.1) 式的循环算子子块系数(m=0,1,2)如下.

情形1  若b(0)10,

(ⅰ) 当β0>0时, 有

a(1)m=1b(0)1((j=0,3(b(0)2+(1)j+1β120)dm+j+1e(1)j2β120x)+b(0)2dm+7)+cm2;a(2)m=dm+1e2β120x+dm+4e2β120x+dm+7;a(3)m=1b(0)1((j=0,3((1)j2b(0)2β1202β0b(0)1+b(0)0)dm+j+1e(1)j2β120x)+b(0)0dm+7);

(ⅱ) 当β0<0时, 有

a(1)m=1b(0)1((j=0,3(b(0)2dm+j+1+(1)j+1(β0)12dmj+4)cos(jπ2+2(β0)12x))+b(0)2dm+7)+cm2;a(2)m=dm+1cos2(β0)12x+dm+4sin2(β0)12x+dm+7;a(3)m=1b(0)1((j=0,3((b(0)02β0b(0)1)dm+j+1+(1)j2b(0)2(β0)12dmj+4b(0)1)×cos(jπ2+2(β0)12x))+b(0)0dm+7);

(ⅲ) 当β0=0时, 有

a(1)m=12b(0)1(2b(0)2dm+1x2+2(b(0)2dm+4dm+1)xdm+4+2b(0)2dm+7)+cm2;a(2)m=dm+1x2+dm+4x+dm+7;a(3)m=1b(0)1(b(0)0dm+1x2+(b(0)0dm+4+2b(0)2dm+1b(0)1)x(dm+1b(0)2dm+4)b(0)1+b(0)0dm+7);

情形2  若b(0)1=0,

(ⅰ) 当b(0)20时, 有

a(1)m=b(0)02b(0)2dm+1e2b(0)2x+dm+4;a(2)m=dm+1e2b(0)2x;

a(3)m=(b(0)0)24(b(0)2)2dm+1e2b(0)2x+dm+7e2b(0)2x+b(0)0b(0)2dm+4b(0)0cm2b(0)2;

(ⅱ) 当b(0)2=0时,

a(1)m=b(0)0dm+1x+dm+4;a(2)m=dm+1;

a(3)m=(b(0)0)2dm+1x2+b(0)0(2dm+4cm)x+dm+7;

其中di (i=1,,9);cm为任意常数.

   由命题3.1可知循环算子的结构, 展开(3.2)式中的12个方程, 将其中的a(0)m (实际上为a(4)m)统一替换成cma(1)m(m=0,1,2), 再次整理约化得到关于wi2t,wit,wi(i=1,2)的多项式. 根据它们的线性无关性, 令其系数为零, 从而可得9个方程. 当m=0,1,2时, 方程组如下

Dxa(1)mb(0)1a(3)m+b(0)0a(2)m=0;Dxa(4i)m+(1)i(2b(0)2a(4i)m+2b(0)i1a(1)mb(0)i1cm)=0 (i=1,2).
(3.4)

为了得到循环算子子块系数, 我们将上述9个方程分为三组(m=0,1,2).i=1,2时, 固定m=2, 相应的(3.4) 式可得

Dxa(1)2b(0)1a(3)2+b(0)0a(2)2=0;
(3.5)

Dxa(3)2+2b(0)2a(3)22b(0)0a(1)2+b(0)0c2=0;
(3.6)

Dxa(2)22b(0)2a(2)2+2b(0)1a(1)2b(0)1c2=0.
(3.7)

下面根据(3.5)–(3.7) 式可以求出a(i)2 (i=1,2,3). 由于求解过程不唯一(可分b(0)1是否为零或b(0)0是否为零两种情况), 下面我们只给出b(0)10 (=0)的求解过程(后者类同它, 且所得到的循环算子子块系数与前者类似, 区别是某些项前系数及符号发生改变).

首先将(3.5)式代入到(3.7)式的求导方程中, 可得

D2xa(2)22b(0)2Dxa(2)2+2(b(0)1)2a(3)22b(0)1b(0)0a(2)2=0.
(3.8)

其次我们分两种情况进行讨论.

情形1  若b(0)10时, 由(3.7)和(3.8) 式反解, 有

a(1)2=12b(0)1(Dxa(2)2+2b(0)2a(2)2+b(0)1c2);
(3.9)

a(3)2=12(b(0)1)2(D2xa(2)2+2b(0)2Dxa(2)2+2b(0)1b(0)0a(2)2).
(3.10)

再将上述得到的a(1)2,a(3)2代入(3.6)式, 可得

D3xa(2)24β0Dxa(2)2=0,

其特征方程为

λ34β0λ=λ(λ24β0)=0.

这样, 分β0>0 (<0=0)三种情况求出对应a(2)2的表达式. 将其代入(3.9)和(3.10) 式可得a(1)2,a(3)2.

情形2   若b(0)1=0时, 由(3.5)和(3.7) 式可得

Dxa(1)2+b(0)0a(2)2=0;
(3.11)

Dxa(2)22b(0)2a(2)2=0.
(3.12)

利用(3.12)式求出a(2)2的通解, 再由(3.6)和(3.11) 式求出a(1)2,a(3)2的具体形式.

同理类似上面的计算过程, 可以分别固定m=0,1, 通过(3.4) 式得到的相应方程组, 就能求出其余系数a(i)m. 综上我们以b(0)10 (=0)两种情形为例, 得到了循环算子子块系数的一组解, 即命题3.2.

3.3 算例

实际上, 命题3.2揭示了Hamilton算子与循环算子间的关系, 即只要找到一个方程拥有的形如(2.4)式的Hamilton结构, 就可依据原理直接写出其允许的循环算子. 下面我们用两个例子说明它的正确性及便捷性.

例1   讨论一类方程: auxx+2buxt+cutteu=0 (a,b,c,e为常数) 的循环算子.

根据文献[10]可知, 其Hamilton正则形式为

x(w1w2)=((0 1ae 0)+(ba 00 ba)Dt+(0 0acb2a 0)D2t)(w1w2),

其中(w1,w2)T=(u,aux+but)T.acb2=0, 该类方程对应的Hamilton算子为一阶常系数算子. 下面我们以两种特殊形式为例.

特别地, 当方程中a=1,b=1,c=1,e=1时, 其Hamilton算子为

H=(0110)+(1001)Dt.

由上式可知: b(0)10,β0<0, 属于命题3.2中情形1(ⅱ), 从而得形如(3.1) 式的子块如下

R1=2m=0(dm+1sin2xdm+4cos2x+cm2)Dmx;

R2,3=2m=0(dm+1cos2x+dm+4sin2x±dm+7)Dmx.

特别地, 当方程中a=1,b=1,c=1,e=0时, 其Hamilton算子为

H=(0100)+(1001)Dt.

由上式可知: b(0)10,β0=0, 属于命题3.2中情形1(ⅲ), 即得子块为

R1=2m=0(dm+1xdm+42+cm2)Dmx;

R2=2m=0(dm+1x2+dm+4x+dm+7)Dmx;R3=2m=0(dm+1)Dmx.

例2   考虑方程: uxx+uxt+ut=0的循环算子.

根据文献[11] 可知, 其Hamilton正则形式中的Hamilton算子为

H=(1011)+(1001)Dt.

由上式可知: b(0)1=0,b(0)20,属于命题3.2中情形2(i), 故所允许循环算子的子块为

R1=2m=0(dm+12e2x+dm+4)Dmx;R2=2m=0(dm+1e2x)Dmx;

R3=2m=0(dm+14e2x+dm+7e2xdm+4+cm2)Dmx.

以上两个例子表明: 在这一途径下获得循环算子更为直观, 无需计算过程, 只需提供形如(2.4) 式的Hamilton结构, 即可写出.

4 结语

本文不同与已有文献中Hamilton算子形如(2.2) 式四个子块的设定, 将其写成有限和的形式, 使所含待定参数增多, 导致确定循环算子的方程组更为直观、具体. 但是, 在求解该方程组上存在着一定的困难. 为了克服困难, 我们对循环算子子块系数的自变量做限定, 从而得到以上结论.

实际上, 这一思路可以延续, 我们根据(2.3)和(2.5) 式还可以对更高阶Hamilton算子继续深入. 当二阶常系数Hamilton算子采用形式

H=P2D2t+P1Dt+P0,(P2=(b(2)2 b(2)1b(2)0 b(2)2)0)

时, 研究所允许循环算子的结构发现, 其所得到的28个方程繁琐且复杂, 且求解存在更多困难. 因此在约化确定方程组时, 我们将上述Hamilton算子中的个别待定参数做具体限制, 可以推测出所允许的候选循环算子类型如下.

型1   R=(R1 00 R1), 其中R1为常系数二阶算子.

型2   R=(R1 R2kR2 R1), 其中R1为常系数二阶算子, R2的系数为关于x的二阶算子且k为常数.

型3   (ⅰ) R=(R1 a(2)0R3 R4), (ⅱ) R=(R1 R2a(3)0 R4), 其中R1,R4的二阶项系数相等且为常数, 即a(1)2=a(4)2, 其余系数待定.

这里不给出证明过程, 只给出上述类型的参数限制条件, 如限定Hamilton算子中的元素b(2)i=0 (i=1,2)(b(1)2)2b(2)0b(0)10时, 允许循环算子的结构为型1 (a(3)2=0) 或型3(i) (a(3)20); 如限定b(2)i=0(i=0,2)(b(1)2)2b(2)1b(0)00时, 允许循环算子的结构为型1 (a(2)2=0) 或型3(ii) (a(2)20); 再如限定b(2)2=0,b(2)1b(2)00b(1)20,b(0)1b(0)00时, 允许循环算子的结构为型2. 对于其它类型的出现有待我们进一步研究. 下面给出算例加以说明.

例3   在例1中, 若acb20, 该类方程对应的Hamilton算子为二阶常系数算子. 特别地, 当a=2,b=12,c=2,e=0时, 其Hamilton算子为

H=(0 120 0)+(14 00 14)Dt+(0 0158 0)D2t.

由上述限定条件可知, 所允许循环算子的结构可能为型1或型3(ⅰ).

这种寻找候选循环算子类型的方法具有局限性, 且获得类型有限, 限制条件太强. 因此, 我们寄希望于寻找新的途径, 以得到二阶常系数Hamilton算子所允许循环算子的一般结构.

参考文献

Olver P J .

Evolution equation possessing infinitely many symmetries

Journal of Mathematical Physics, 1977, 18 (6): 1212- 1215

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

Bluman G W , Kumei S . Symmetries and Differential Equations. New York: Springer-Verlag, 1989

[本文引用: 1]

Gürses M .

On construction of recursion operators from lax representation

Department of Mathematies, 1999, 40 (12): 6473- 6490

URL    

Stampp M , Karasu A .

Recursion operators and Hamiltonian structures in sato's theory

Letters in Mathematical Physics, 1990, 20, 195- 210

DOI:10.1007/BF00398363     

Wang J P .

A list of 1+1 dimensional integrable equations and their properties

Journal of Nonlinear Mathematical Physics, 2002, 9 (1): 213- 233

URL    

Baldwin D E , Hereman W .

A symbolic algorithm for computing recursion operators of nonlinear PEDs

International Journal of Computer Mathematics, 2010, 87 (5): 1094- 1119

DOI:10.1080/00207160903111592     

Fokas A S , Fuchssteiner B .

On the structure of symplectic operators and hereditary ymmetries

Lettere al Nuovo Cimento, 1980, 28 (8): 299- 303

DOI:10.1007/BF02798794     

Wang Y , Li B , An H L .

Dark Sharma-Tasso-Olver equations and their recursion operators

Chinese Physics Letters, 2018, 35 (1): 010201

DOI:10.1088/0256-307X/35/1/010201      [本文引用: 1]

李静. 无穷维Hamilton系统下循环算子及反问题的研究[D]. 呼和浩特: 内蒙古工业大学, 2017

[本文引用: 1]

Li J. Study about recursion operators and inverse problem of infinite dimensional Hamiltonian system[D]. Huhehot: Inner Mongol University of Technology, 2017

[本文引用: 1]

霍晓霞, 任文秀.

二阶偏微分方程的Hamilton正则形式化的分类讨论

内蒙古工业大学学报, 2019, 38 (1): 1- 7

DOI:10.3969/j.issn.1001-5167.2019.01.001      [本文引用: 2]

Huo X X , Ren W X .

A discussion on the classification of Hamilton canonical formalization of second order partial differential equations

Journal of Inner Mongol University of Technology, 2019, 38 (1): 1- 7

DOI:10.3969/j.issn.1001-5167.2019.01.001      [本文引用: 2]

许晶. 获得Hamilton系统下循环算子、守恒律的方法及应用[D]. 呼和浩特: 内蒙古工业大学, 2013

[本文引用: 1]

Xu J. Obtaining methods and applications of recursion operators and conservation laws in the Hamiltonian system[D]. Huhehot: Inner Mongol University of Technology, 2013

[本文引用: 1]

/