一类广义不稳定时空分数阶薛定谔方程的近似解

时间:2023-08-18 11:10:02 来源:网友投稿

洪宝剑

(南京工程学院 数理学院,江苏 南京 211167)

众所周知,许多自然科学和工程问题都可以归结为非线性方程的求解,尤其是分数阶方程几乎应用于社会生活的每一个领域.因此,研究非线性方程的解就成了国内外学者需要关注的首要问题.迄今为止,人们提出了寻求非线性方程精确解及近似解的各种方法,如贝克隆变换法、约束方程法、同伦分析法、Adomian展开法、变分迭代法等[1-7].近年来,Laplace变换法被广泛应用于各类非线性系统的求解[8-10].论文将借助Adomian多项式和同伦摄动的理论,将Laplace变换法应用于一类广义分数阶薛定谔方程,得到较好的结果.

考虑下列广义不稳定时空分数阶薛定谔方程

(1)

其中:u=u(x,t),a,γ为非零常数,v(x)为实函数.

为了方便后面计算,关于分数阶导数和Laplace变换,给出几个定义[7,23].

定义1一个实函数f(t),当μ∈,t>0时,若存在实数p>μ,使得f(t)=tpf1(t),其中:f1(t)∈C[0,∞),则称f(t)∈Cμ,若f(n)(t)∈Cμ,n∈,则称

定义2函数u(x,t)关于t的Caputo分数阶导数定义为[13,18,23]

(2)

Caputo分数阶导数具有下列性质

DαC=0,C为常数,

Dαf[u(x)]=fuu(α)(x).

(3)

关于Laplace变换迭代法的具体步骤可参阅文献[8-10,24-27],对(1)进行简单变换得

(4)

(5)

(6)

对(6)两边作Laplace逆变换

(7)

构造同伦方程

(8)

显然当p:0→1时,u:u0→u,若令

(9)

则(1)的精确解为

(10)

比较p的同次幂,构造迭代式

(11)

其中:An是(8)中非线性项N(u)=|u2r|u通过Adomian多项式线性化后的第n项[4,13,21].

(12)

(13)

当r≠1时,若u为实函数,则按牛顿二项公式展开,若为复函数,则按三角形式展开,均可得到An.再由(10),(11)式就可以得到方程(1)的各级近似解或精确解,下面讨论(1).

情形1当r=1,v(x)=-cos2x,u0=sinx时,有

(14)

由(11),(13)式,有

u0=sinx.

(15)

(16)

[a2sin(x+2πβ)+4a(2γ+1)sin2xsin(x+πβ)-2asin(x+πβ)+

(17)

(18)

方程(14)的一个n级近似解为

(19)

(20)

情形2当v(x)=2λ,u0=eix时,有

(21)

构造迭代式

(22)

u0=eix.

(23)

(24)

(25)

(26)

(27)

方程(21)的一个n级近似解为

(28)

所以方程(21)的一个精确解为

(29)

注2当γ=a,r=1,λ=0时,u1,u2,u3与文献[23]中的结果一致.当γ=1,λ=0时,Un表达式比文献[13]中的结果(3.18)标准.当γ=a=0.5,r=1,λ=0时,(29)式转化为文献[28]中的(49)式.

情形3当v(x)=k(1-x2),u0=x2时,有

(30)

接下来尝试对迭代式的构成方式进行修正.考虑到u(x,t)=φ(x,t)eiψ(x,t),|u(x,t)|=|φ(x,t)|,利用文献[23]中的修正思想,构造迭代式

(31)

u0=x2.

(32)

(33)

k2x2-2k2x4+k2x6],

(34)

4γk2x4r+2+8γk(γ-k)x4r+4-4γk(2γ-k)x4r+6+4γ2kx8r+2+

4γ2(2γ-k)x8r+4+k2(2γ-k)x8-k2(4γ-3k)x6+k2(2γ-3k)x4+k3x2].

(35)

方程(30)的近似解为

Un=u0+u1+u2+u3+….

借助Mathematica软件,表1,2和图1,2分别对方程(1)的n级近似解(28)式和精确解(29)式进行对比和数值模拟,分别讨论当α=1,β=1,γ=1,a=1,r=1,v(x)=0时的整数阶方程情形和当α=0.5,β=0.8,γ=1,a=1,r=1,v(x)=0时的分数阶方程情形.图3,4对方程(30)当a=k=r=1,γ=0.5时,在α,β不同参数取值的近似解的实部和虚部进行了模拟.

表1 当α=1,β=1,γ=1,a=1,r=1,v(x)=0时,(28)和(29)虚部的数值结果

注3Absolute Error精确位数取为小数点后20位,迭代到50次后误差几乎消失了.

表2 当α=0.5,β=0.8,γ=1,a=1,r=1,v(x)=0时,(28)和(29)实部的数值结果

部分数值模拟结果如图1~4所示.

图1 当α=1,β=1,γ=1,a=1,r=1,v(x)=0时,精确解与近似解迭代5,25次虚部图

图2 当α=0.5,β=0.8,γ=1,a=1,r=1,v(x)=0时,近似解迭代5,25,50次实部图

图3 不同参数下各级近似解实部变化 图4 不同参数下各级近似解虚部变化

图1,2的数值结果表明,U25,U50已经和精确解靠得很近.在图3,4中随着α,β的增大,实部和虚部函数稳定地向α=β=0.05及α=β=1靠拢,因而该迭代算法有效.

论文通过将Laplace变换和同伦摄动法相结合,求解一种不稳定时空分数阶薛定谔方程,得到了方程的各级近似解,并进行了误差分析和数值模拟.研究结果表明,这一方法对于薛定谔类复偏微分方程近似解的处理规范有效.

猜你喜欢虚部薛定谔实部复数知识核心考点综合演练中学生数理化·高一版(2022年3期)2022-04-15拟相对论薛定谔方程基态解的存在性与爆破行为数学物理学报(2022年1期)2022-03-16Chern-Simons-Higgs薛定谔方程组解的存在性数学物理学报(2021年6期)2021-12-21两类特殊多项式的复根虚部估计中等数学(2021年6期)2021-08-14例谈复数应用中的计算两次方法数学学习与研究(2020年23期)2020-01-11一类相对非线性薛定谔方程解的存在性数学物理学报(2019年1期)2019-03-21薛定谔的馅幽默大师(2019年6期)2019-01-14浅谈正Γ型匹配网络的设计卷宗(2016年8期)2016-11-15一种基于电涡流和实部互阻抗检测的金属温度监测方法电测与仪表(2016年2期)2016-04-12

推荐访问:广义 不稳定 近似