Processing math: 100%

数学物理学报, 2021, 41(2): 289-295 doi:

论文

环图上Sturm-Liouville算子的部分反谱问题

官声玉,1, 吴东杰1, MuratSat2, 杨传富1

A Partial Inverse Problem for the Sturm-Liouville Operator on Quantum Graphs with a Loop

Guan Shengyu,1, Wu Dongjie1, Murat Sat2, Yang Chuanfu1

通讯作者: Murat Sat杨传富

收稿日期: 2020-01-14  

基金资助: 国家自然科学基金.  11871031
江苏省自然科学基金.  BK20201303

Received: 2020-01-14  

Fund supported: the NSFC.  11871031
the Jiangsu NCF.  BK20201303

作者简介 About authors

官声玉,E-mail:guanshengyu@njust.edu.cn , E-mail:guanshengyu@njust.edu.cn

Abstract

This deals with the Sturm-Liouville operator on the quantum graphs with a loop. Given the potential on a part of edges, we try to recover the remaining potential from the subspectrum. The uniqueness theorem and a constructive algorithm for the solution of this partial inverse problem are provided.

Keywords: Partial inverse spectral problem ; Sturm-Liouville operator ; Quantum graph ; Riesz-basis

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

本文引用格式

官声玉, 吴东杰, MuratSat, 杨传富. 环图上Sturm-Liouville算子的部分反谱问题. 数学物理学报[J], 2021, 41(2): 289-295 doi:

Guan Shengyu, Wu Dongjie, Murat Sat, Yang Chuanfu. A Partial Inverse Problem for the Sturm-Liouville Operator on Quantum Graphs with a Loop. Acta Mathematica Scientia[J], 2021, 41(2): 289-295 doi:

1 引言

量子图上微分算子是现代数学物理研究的一个重要方向. 该模型的重要性出现在碳纳米结构、声学和光子晶体理论研究中[1, 6-7, 11-12]. 本文主要是构造量子图上Sturm-Liouville算子中的未知势函数. 这个问题在星型图上已经解决[2, 8-9]. 本文考虑带环量子图:在假定两条边上的势函数是已知的情况下去重构剩余未知势函数. 杨和Bondarenko[10]考虑了带有一个环的套索图. 在此基础上, Bodarenko和Shieh又研究了具有一个环和n1条边的量子图, 利用Riesz基和Weyl函数得到了部分反问题的唯一性定理和构造解[3].

本文在量子环图上考虑除环外所有边的势函数都是已知的Sturm-Liouville算子. 讨论其特征值的渐近式, 唯一性定理及重构算法。

2 环图上微分算子

考虑如下由一环和两条边构成的量子图G. 各条边的长度分别为l1=1,l2=2,l3=3.引入参数xj[0,lj],j=¯1,3, 且在边界顶点处满足x1=x2=0, 在顶点处满足x1=1,x2=2,x3=0,x3=3.

图 1

图 1   量子环图


在图G上研究Sturm-Liouville边值问题L

yj(xj)+qj(xj)yj(xj)=λyj(xj),xj(0,lj),j=1,2,3,
(2.1)

y1(0)=y2(0)=0,
(2.2)

y1(1)=y2(2)=y3(3)=y3(0),
(2.3)

y1(1)+y2(2)+y3(3)=y3(0).
(2.4)

这里λ为谱参数, qjL2[0,lj]且是实的.

定义2.1  令B2,a为指数不大于a的Paley-Wienner函数. 符号κk,odd(ρ),κk,even(ρ)分别表示属于B2,a的奇、偶函数

κk,odd(ρ)=k0K(t)sinρtdt,    κk,even(ρ)=k0N(t)cosρtdt.

Cj(xj,λ),S(xj,λ)为等式(2.1) 满足初值条件: Cj(0,λ)=Sj(0,λ)=1,Cj(0,λ)=Sj(0,λ)=0下的解, 则由文献[4]得到如下关系式

{Cj(lj,ρ)=cosρlj+κlj,even(ρ),Sj(lj,ρ)=sinρljρ+κlj,odd(ρ)ρ,Sj(lj,ρ)=cosρlj+κlj,even(ρ).
(2.5)

3 特征值分布

研究边值问题L的特征值, 需研究示性函数的性质. 示性函数

Δ(λ)=S1(1,λ)S2(2,λ)[2C3(3,λ)S3(3,λ)][S1(1,λ)S2(2,λ)+S1(1,λ)S2(2,λ)]S3(3,λ).
(3.1)

下面我们考虑它的零点的分布情况, 基于儒歇定理([4], 定理1.1.3)及渐近分析有如下引理.

引理3.1  问题L特征值为{λnk}nZ,k=1,2{λ(i)n3}nZ,i=1,2{λ(i)n4}nN,i=1,2, 且满足

{ρnk=λnk=|2πn+αk|+κn,   nZ,k=1,2,ρ(i)n3=λ(i)n3=|2πn+α3|+κn,   nZ,i=1,2,ρ(i)n4=λ(i)n4=nπ+κn,   nN,i=1,2,
(3.2)

其中α1=arccos2+24,α2=arccos224,α3=2π3.

4 部分反问题唯一性定理

选择子谱Λ:={λnk}nZ,k=1,2, 记符号I={(n,j):ρnjΛ}. 本文假设

(1) Λ中所有特征值是不同的;

(2) Λ中所有特征值是正的;

(3) 对于λnjΛ, S1(1,λnj)S2(2,λnj)S3(3,λnj)0.

在假设(1)–(3)下, 研究如下反问题.

反问题: 给定势函数q2,q3, 和子谱Λ, 重构函数q1.

由关系式(3.2), 得到

{S1(1,ρ)=sinρρ+1ρ10K(t)sinρtdt,S1(1,ρ)=cosρ+10N(t)cosρtdt,
(4.1)

其中函数K,N, 是属于L2(0,1)的实函数. 将此式代入(3.1)式, 得到

10K(t)anjsinρnjtdt+10N(t)bnjcosρnjtdt=fnj,    (n,j)I,
(4.2)

anj=[2S2(2,λnj)S2(2,λnj)C3(3,λnj)S2(2,λnj)S3(3,λnj)S2(2,λnj)S3(3,λnj)]ρnj,bnj=S2(2,λnj)S3(3,λnj)ρ2nj,fnj=anjsinρnjbnjcosρnj.
(4.3)

f(t)=[K(t)N(t)],vnj(t)=[anjsinρnjtbnjcosρnjt],(n,j)I,
(4.4)

设向量函数g, h是属于Hilbert空间H, 则在空间H上内积为

(g,h)=10(g1(t)h1(t)+g2(t)h2(t))dt,g=[g1g2],h=[h1h2],g,hH.

向量函数vnk, f属于实Hilblet空间H:=L2(0,1)L2(0,1), (4.2)式可写为

(f,vnj)H=fnj,    (n,j)I.
(4.5)

引理4.1   向量函数系V={vnj}(n,j)IH中完备.

(反证法)   设存在w1,w2L2(0,1), 使得

10w1(t)anjsinρnjtdt+10w2(t)bnjcosρnjtdt=0,    (n,j)I,
(4.6)

由假设(3) 和(3.1)式, 可知

anj=S1(1,λnj)bnjρnjS1(1,λnj).

再将此式代入(4.6)式, 有

10(w1(t)S1(1,λnj)sinρnjtρnjw2(t)S1(1,λnj)cosρnjt)dt=0,
(4.7)

W(λ):=10(w1(t)S1(1,λ)sinρtρw2(t)S1(1,λ)cosρt)dt
(4.8)

的零点为Λ. 结合(4.1)式, 有估计式

W(λ):=O(e2|Imρ||ρ|),    |ρ|.
(4.9)

在假设(1)–(2)下构造无穷乘积

D(λ)=λnjΛ(1λλnj).
(4.10)

易知W(λ)D(λ)是整函数, 再由渐近式(3.2), 函数D(λ)可表示为

D(λ)=C2k=1(cosρcosαk)+κ2,even(ρ),
(4.11)

这里C是非零常数, 而且对于正数ε, ρ, 我们有估计式

|D(ρ2)|Ce2|Imρ|,    ε<argρ<πε,    |ρ|ρ,
(4.12)

因此结合(4.9)式得到

W(λ)D(λ)=o(1),    λ=ρ2,    ε<argρ<πε,    |ρ|ρ.
(4.13)

这里我们注意到若{ξn}nN, {ηn}nN{0}分别是S1(1,λ)S1(1,λ)的零点. 则将λ=ξnλ=ηn代入(4.8) 式容易得到

10w1(t)sinξntξndt=0 (nN),10w2(t)cosηndt=0 (nN{0}).

而又因为{sinξntξn}nN{cosηnt}nN{0}L2(0,1)里完备[5], 所以w1=0, w2=0, 因此VH中完备.

下面基于引理4.1证明反问题解的唯一性. 首先定义在边值问题L的基础上, 考虑同样形式的问题˜L, 但有不同的势函数q, 如果有某个符号γ表示与L相关的对象, 则对应的符号˜γ表示与˜L相关的类似对象.

定理4.2   设边值问题L˜L及其子谱Λ˜Λ定义如上, 且满足假设(1)–(3), q2(x)=˜q2(x),x(0,2);q3(x)=˜q3(x),x(0,3), 以及Λ=˜Λ, 则在(0,1)上有q1(x)=˜q1(x).

  因为q2(x)=˜q2(x), q3(x)=˜q3(x), 得到S2(2,λ)˜S2(2,λ),S2(2,λ)˜S2(2,λ),S3(3,λ)˜S3(3,λ),S3(3,λ)˜S3(3,λ). 再由(3.11)式和Λ=˜Λ, 观察到vnj=˜vnj; fnj=˜fnj, (n,j)I. 而由引理3.2又知系VH中完备, 因此由(3.13)式可知S1(1,λ)˜S1(1,λ), S1(1,λ)˜S1(1,λ). 所以由文献[6]得到在(0,1)上, 有q1(x)=˜q1(x).

5 部分反问题重构算法

引理5.1   向量函数系V={vnj}(n,j)IH中是Riese基.

  由(2.5)和(4.3)式, 有

anj=2sin2ρnj(1cos3ρnj)cos2ρnjsin3ρnj+κn,

bnj=sin2ρnjsin3ρnj+κn,    (n,j)I.

因此我们有{||vnjv0nj||H}(n,j)Il2, 其中

v0nk(t)=[[2sin2αk(1cos3αk)cos2αksin3αk]sin|2πn+αk|tsin2αksin3αkcos|2πn+αk|t].

注意到由(3.2)式可以看出2sin2αk(1cos3αk)cos2αksin3αk0,sin2αksin3αk0. 下证函数系V0:={v0nj}(n,j)IH中是Riesz基. 易知{sin(2πn+α1)t}nZ, {cos(2πn+α2)t}nZL2(0,1)中的Riese基. 引入线性算子A, B: HHvH, 满足

Av=A[v1v2]=[v1v2+sin2α1sin3α12sin2α1(1cos3α1)cos2α1sin3α1f(v1)],

Bv=B[v1v2]=[v1+2sin2α2(1cos3α2)cos2α2sin3α2sin2α2sin3α2g(v2)v2],

其中

f(u1)(t)=nZcn(u1)cos|2πn+α1|t,

u1(t)=nZcn(u1)sin|2πn+α1|t,

g(u2)(t)=nZcn(u2)sin|2πn+α2|t,

u2(t)=nZcn(u2)cos|2πn+α2|t,

这里cn(u1)cn(u2)分别是函数u1,u2在Riesz基{sin|2πn+α1|t}nZ, {cos(2πn+α2)t}nZ下的坐标. 从Riesz基性质可知存在正数C1, C2使得

C1||u1||L2||f(u1)||L2C2||u1||L2,

C1||u1||L2||g(u2)||L2C2||u1||L2.

由此可知算子A, B及其逆

A1v=A1[v1v2]=[v1v2sin2α1sin3α12sin2α1(1cos3α1)cos2α1sin3α1f(v1)],

B1v=B1[v1v2]=[v12sin2α2(1cos3α2)cos2α2sin3α2sin2α2sin3α2g(v2)v2]

H中有界. 且算子A, B将序列V0映射成Riezs基

Av0n1(t)=(2sin2α1(1cos3α1)cos2α1sin3α1)[sin|2πn+α1|t0],

Bv0n2(t)=sin2α2sin3α2[0cos|2πn+α2|t].

即函数系{Av0n1}{Bv0n2}H中是Riesz基, 因此V0H中也是Riesz基. 又由引理3.2可知系VH中完备且和Riese基V0是按l2逼近, 因此VH中是Riezs基.

重构算法: 给定函数q2,q3, 和特征值Λ.

步骤1   解初值问题

ljSj(xj,λnj)=λnjSj(xj,λnj),Sj(0,λnj)=0,Sj(0,λnj)=1,

ljCj(xj,λnj)=λnjCj(xj,λnj),Cj(0,λnj)=1,Cj(0,λnj)=0,j=2,3.

构造函数S1(1,λnj),S1(1,λnj), S2(1,λnj),S2(1,λnj), 其中(n,j)I.

步骤2   通过(4.2)和(4.5)式得到向量函数vnj, fnj.

步骤3  构造向量函数f, 通过它在Riesz基上的坐标, 即构造K(t), N(t).

步骤4   由(4.1)式计算得到S1(1,ρ), S1(1,ρ).

步骤5   参考文献[4], 由Weyl函数重构函数q1.

参考文献

Belishev M I .

Boundary spectral inverse problem on a class of graphs (trees) by the BC-method

Inverse Problems, 2004, 20 (3): 647- 672

DOI:10.1088/0266-5611/20/3/002      [本文引用: 1]

Bondarenko N P .

A partial inverse problem for the Sturm-Liouville operator on a star-shaped graph

Analysis and Mathematical Physics, 2018, 8 (1): 155- 168

DOI:10.1007/s13324-017-0172-x      [本文引用: 1]

Bondarenko N P , Shieh C T .

Partial inverse problems for quadratic differential pencils on a graph with a loop

Journal of Inverse and Ill-Posed Problems, 2020, 28 (3): 449- 463

DOI:10.1515/jiip-2018-0104      [本文引用: 1]

Freiling G , Yurko V A . Inverse Sturm-Liouville Problems and Their Applications. New York: Nova Science Publishers, 2001

[本文引用: 3]

He X , Volkmer H .

Riesz bases of solutions of Sturm-Liouville equations

Journal of Fourier Analysis and Applications, 2001, 7 (3): 297- 307

DOI:10.1007/BF02511815      [本文引用: 1]

Kuchment P .

Quantum graphs: I. Some basic structures

Waves in Random Media, 2004, 14 (1): S107- S128

DOI:10.1088/0959-7174/14/1/014      [本文引用: 2]

Kurasov P , Nowaczyk M .

Inverse spectral problem for quantum graphs

Journal of Physics A: Mathematical and General, 2005, 38 (22): 4901- 4915

DOI:10.1088/0305-4470/38/22/014      [本文引用: 1]

Wang Y P , Shieh C T .

Inverse problems for Sturm-Liouville operators on a star-shaped graph with mixed spectral data

Applicable Analysis, 2019, 99 (14): 2371- 2380

URL     [本文引用: 1]

Yang C F .

Inverse spectral problems for the Sturm-Liouville operator on a d-star graph

Journal of Mathematical Analysis and Applications, 2010, 365 (2): 742- 749

DOI:10.1016/j.jmaa.2009.12.016      [本文引用: 1]

Yang C F , Bondarenko N P .

A partial inverse problem for the Sturm-Liouville operator on the lasso-graph

Inverse Problems and Imaging, 2018, 13 (1): 69- 79

URL     [本文引用: 1]

Yurko V A .

Inverse spectral problems for Sturm-Liouville operators on graphs

Inverse Problems, 2005, 21 (3): 1075- 1086

DOI:10.1088/0266-5611/21/3/017      [本文引用: 1]

Yurko V A .

Inverse spectral problems for differential operators on arbitrary compact graphs

Journal of Inverse and Ill-Posed Problems, 2010, 18 (3): 245- 261

URL     [本文引用: 1]

/