摘要

我们构造了(+1)-用Exp-function方法对KdV-Zakharov-Kuznetsev方程进行了量纲修正。该方法有效地获得了大量任意参数的精确行波解。结果表明,对于求解具有任意参数的高维非线性偏微分方程的解析解,Exp-function方法是一种有效而直观的数学工具。

1.简介

非线性偏微分方程(NLPDEs)在应用科学的不同分支中发挥着重要作用。近年来,许多研究人员研究了非线性偏微分方程的精确行波解,这些行波解对揭示复杂物理现象至关重要。在过去的几十年中,各种有效且强大的方法,如变分迭代法[1——],tanh-coth方法[4],同伦摄动法[5——7],Fan子方程法[8],射影Riccati方程法[9],微分变换法[10],直接代数法[11],第一积分法[12],Hirota双线性方法[13],改进的扩展直接代数方法[14],扩展tanh方法[15],Backlund变换[16],分叉法[17],Cole-Hopf变换方法[18],sech-tanh方法[19],(𝐺/𝐺)-膨胀法[20——22],已修改(𝐺/𝐺)-膨胀法[23],多波法[24],已扩展(𝐺/𝐺)-膨胀法[25,26]、和其他[27——33]用于寻求非线性发展方程(NLEE)的精确行波解。

最近,他和吴[34]提出了一种新的方法,称为Exp-function方法,用于搜索数学物理中出现的非线性发展方程的行波解。Exp-function方法被广泛应用于多种非线性偏微分方程,例如好的Boussinesq方程[35],非线性微分方程[36],高阶边值问题[37],非线性问题[38],Calogero-Degasperis-Fokas方程[39],非线性反应扩散方程[40],2D Bratu型方程[41],非线性晶格微分方程[42],广义Zakharov方程[43],(3+1)维Jimbo-Miwa方程[44],修正的Zakharov-Kuznetsov方程[45],布鲁塞尔反应扩散模型[46],非线性热方程[47]以及其他重要的NLPDE[48——51].

在本文中,我们应用Exp-function方法[34]获得非线性偏微分方程的解析解,即(3+1)维修正KdV-Zakharov-Kuznetsev方程。

2.Exp函数方法的描述

考虑一般非线性偏微分方程𝑃𝑢,𝑢𝑡,𝑢𝑥,𝑢𝑦,𝑢𝑧,𝑢𝑡𝑡,𝑢𝑥𝑡,𝑢𝑥𝑥,𝑢𝑥𝑦,𝑢𝑦𝑦,𝑢𝑦𝑡,𝑢𝑧𝑧,𝑢𝑧𝑡,𝑢𝑧𝑥,𝑢𝑧𝑦,=0.(2.1)Exp-function方法的主要步骤[34]如下所示。

步骤1。将复杂变量视为𝑢(𝑥,𝑦,𝑧,𝑡)=𝑢(美国),美国=𝑥+𝑦+𝑧𝑉𝑡.(2.2)现在使用(2.2), (2.1)转换为非线性常微分方程𝑢(美国)𝑄𝑢,𝑢,𝑢,𝑢,=0,(2.)其中素数表示关于美国.

第2步。我们假设(2.3)可以表示为[34]𝑢(美国)=𝑑𝑛=𝑐𝑎𝑛e(电子)x个第页(𝑛美国)𝑞𝑚=𝑝𝑏𝑚=𝑎e(电子)x个第页(𝑚美国)𝑐e(电子)x个第页(𝑐美国)++𝑎𝑑e(电子)x个第页(𝑑美国)𝑏𝑝e(电子)x个第页(𝑝美国)++𝑏𝑞e(电子)x个第页(𝑞美国),(2.4)哪里𝑐,𝑑,𝑝、和𝑞是稍后确定的正整数,以及𝑎𝑛𝑏𝑚是未知常数。方程式(2.4)可以用以下等效形式重写:𝑎𝑢(美国)=𝑐e(电子)x个第页(𝑐美国)++𝑎𝑑e(电子)x个第页(𝑑美国)𝑏𝑝e(电子)x个第页(𝑝美国)++𝑏𝑞e(电子)x个第页(𝑞美国).(2.5)

步骤3。为了确定𝑐𝑝,我们将最高阶线性项与最高阶非线性项进行平衡,并确定𝑑𝑞,我们在中平衡最低阶线性项和最低阶非线性项(2.3). 因此,我们得到了𝑐,𝑑,𝑝、和𝑞.

步骤4。替换的值𝑐,𝑑,𝑝、和𝑞到(2.5),然后替换(2.5)进入(2.3)简化后,我们得到𝑖美国航空航天局𝑖e(电子)x个第页(±𝑖美国)=0,𝑖=0,1,2,,.(2.6)然后每个系数美国航空航天局𝑖=0是设置,生成一个代数方程组𝑎𝑐𝑠𝑏𝑝𝑠.

步骤5。我们假设未知𝑎𝑐𝑠𝑏𝑝𝑠可以通过求解步骤中获得的代数方程组来确定4.将这些值放入(2.5),我们得到了精确的行波解(2.1).

3.方法的应用

在本节中,我们应用该方法构造了(3+1)维修正KdV-Zakharov-Kuznetsev方程的行波解。获得的解决方案将显示在图中1,2,,4,5、和6使用Maple 13软件。

我们考虑(3+1)维修正KdV-Zakharov-Kuznetsev方程𝑢𝑡+𝛼𝑢2𝑢𝑥+𝑢𝑥𝑥𝑥+𝑢𝑥𝑦𝑦+𝑢𝑥𝑧𝑧=0,(.1)哪里𝛼是一个非零常数。

扎耶德[52]已解决(3.1)使用(𝐺/𝐺)-扩展方法。稍后,在文章中[53],他用广义的(𝐺/𝐺)-膨胀法。

在这里,我们将用Exp-function方法求解这个方程。

现在,我们使用转换(2.2)进入(3.1),它产生𝑉𝑢+𝛼𝑢2𝑢+𝑢=0,(.2)其中素数表示关于美国.

根据步骤2,解决方案(3.2)可以写为(2.5). 要确定𝑐𝑝,根据步骤,我们平衡𝑢具有最高阶非线性项𝑢2𝑢英寸(3.2)也就是说,𝑢𝑢2𝑢因此,我们有以下几点:𝑐𝑢=1[]e(电子)x个第页(𝑝+𝑐)美国+𝑐2[],𝑢e(电子)x个第页4𝑝美国+2𝑐𝑢=[(]e(电子)x个第页𝑝+𝑐)美国+𝑐4[],e(电子)x个第页4𝑝美国+(.)哪里𝑐𝑗仅为简单起见的系数;来自(3.3),我们得到𝑝+𝑐=𝑝+𝑐,w个小时c小时e(电子)d日𝑝=𝑐.(.4)要确定𝑑𝑞,我们平衡的最低阶线性项𝑢具有最低阶非线性项𝑢2𝑢英寸(3.2). 我们有𝑢=+𝑑1[]e(电子)x个第页(𝑑𝑞)美国+𝑑2[],𝑢e(电子)x个第页4𝑞美国2𝑢=+𝑑[]e(电子)x个第页(𝑑𝑞)美国+𝑑4[],e(电子)x个第页4𝑞美国(.5)哪里𝑑𝑗仅为了简单而确定系数;来自(3.5),我们获得(𝑑𝑞)=(𝑑𝑞),w个小时c小时e(电子)d日𝑞=𝑑.(.6)可以考虑任何实际值𝑐𝑑,因为它们是自由参数。但最终的解决方案(3.1)不依赖于选择𝑐𝑑.

案例1。我们设置了𝑝=𝑐=1𝑞=𝑑=1.
对于这种情况,试验解决方案(2.5)减少到𝑎𝑢(美国)=1𝑒美国+𝑎0+𝑎1𝑒美国𝑏1𝑒美国+𝑏0+𝑏1𝑒美国.(.7)自,𝑏10, (3.7)可以简化𝑎𝑢(美国)=1𝑒美国+𝑎0+𝑎1𝑒美国𝑒美国+𝑏0+𝑏1𝑒美国.(.8)通过替换(3.8)进入(3.2)并将e(电子)x个第页(±𝑛美国),𝑛=0,1,2,,借助Maple 13,我们得到了一组代数方程𝑎1,𝑎0,𝑎1,𝑏1,𝑏0、和𝑉1𝐴美国航空航天局𝑒美国+美国航空航天局2𝑒2美国+美国航空航天局1𝑒美国+美国航空航天局0+美国航空航天局1𝑒美国+美国航空航天局2𝑒2美国+美国航空航天局𝑒美国=0.(.9)并且,设置e(电子)x个第页(±𝑛美国),𝑛=0,1,2,,,到零,我们得到美国航空航天局=0,美国航空航天局2=0,美国航空航天局1=0,美国航空航天局0=0,美国航空航天局1=0,美国航空航天局2=0,美国航空航天局=0.(.10)为了确定未知数,我们求解得到的代数系统(3.10)借助Maple 13,我们得到了四组不同的解。设置1。我们得到了𝑏1=𝑏1,𝑎1=6𝑏12𝛼,𝑎0=0,𝑎16=±2𝛼,𝑏0=0,𝑉=6,(.11)哪里𝑏1是自由参数。设置2。我们得到了𝑎0=𝑎0,𝑏0=𝑏0,𝑎11=122𝛼2𝛼𝑎20+9𝑏20,𝑎1=±2𝛼,𝑏1=118𝛼𝑎20+14𝑏20,𝑉=2,(.12)哪里𝑎0𝑏0是自由参数。设置3。我们得到了𝑎1=𝑎1,𝑏0=𝑏0,𝑎1=𝑏202𝛼𝑎21+98𝛼𝑎1,𝑎0=𝑏0𝛼𝑎21+9𝛼𝑎1,𝑏1=𝑏202𝛼𝑎21+98𝛼𝑎21,𝑉=+𝛼𝑎21,(.1)哪里𝑎1𝑏0是自由参数。设置4。我们得到了𝑎0=𝑎0,𝑎1=0,𝑎1=0,𝑏1=172𝛼𝑎20,𝑏0=0,𝑉=,(.14)哪里𝑎0是自由参数。
现在,替换(3.11)进入(3.8),我们得到了行波解𝑢(美国)=±6𝑒美国6𝑏1𝑒美国𝑒2𝛼美国+𝑏1𝑒美国.(.15)方程式(3.15)可以简化为𝑢(美国)=±62𝛼12𝑏1(co(o)小时美国n个小时美国)1+𝑏1co(o)小时美国+1𝑏1n个小时美国,(.16)哪里美国=𝑥+𝑦+𝑧+6𝑡.
如果𝑏1=1来自(3.16),我们获得𝑢(美国)=±6𝑖2𝛼t吨n个小时美国.(.17)替换(3.12)进入(3.8)简化后,我们得到了行波解𝑢(美国)=±2𝛼1+12±𝑎02𝛼+𝑏022𝛼𝑎20+9𝑏20(co(o)小时美国n个小时美国)6+2𝛼𝑎20+9𝑏20co(o)小时美国+62𝛼𝑎209𝑏20n个小时美国+6𝑏0,(.18)哪里美国=𝑥+𝑦+𝑧+(/2)𝑡.
如果𝛼是负的,也就是说,𝛼=𝛽,𝛽>0,𝑏0=2𝑎0=0,然后从(3.18),我们获得𝑢(美国)=±美国2𝛽t吨n个小时2.(.19)替换(3.13)进入(3.8)简化后,我们得到𝑢(美国)=𝑎11+72𝑏08𝛼𝑎21+𝑏202𝛼𝑎21+9co(o)小时美国+8𝛼𝑎21𝑏202𝛼𝑎21+9n个小时美国+8𝛼𝑎21𝑏0,(.20)哪里美国=𝑥+𝑦+𝑧(+𝛼𝑎21)𝑡.
如果𝑏0=1,𝛼=6、和𝑎1=1/2, (3.20)成为1𝑢(美国)=2+1+2co(o)小时美国.(.21)替换(3.14)进入(3.8)简化后,我们得到𝑢(美国)=72𝑎072+𝛼𝑎20co(o)小时美国+72𝛼𝑎20n个小时美国,(.22)哪里美国=𝑥+𝑦+𝑧𝑡.
如果𝑎0=𝛼=8, (3.22)成为𝑢(美国)=2e(电子)c美国.(.2)

案例2。我们设置了𝑝=𝑐=2𝑞=𝑑=1.
对于这种情况,试验解决方案(2.5)减少到𝑎𝑢(美国)=2𝑒2美国+𝑎1𝑒美国+𝑎0+𝑎1𝑒美国𝑏2𝑒2美国+𝑏1𝑒美国+𝑏0+𝑏1𝑒美国.(.24)因为,在(3.24)为了简单起见,我们可以考虑𝑏2=1𝑏1=0。然后是解决方案(3.24)简化为𝑎𝑢(美国)=2𝑒2美国+𝑎1𝑒美国+𝑎0+𝑎1𝑒美国𝑒2美国+𝑏1𝑒美国+𝑏0.(.25)执行案例中描述的相同程序1,我们得到了四组解。

设置1。我们得到了𝑏0=𝑏0,𝑎1=0,𝑎0=6𝑏02𝛼,𝑎1=0,𝑎26=±2𝛼,𝑏1=0,𝑉=6,(.26)哪里𝑏0是自由参数。

设置2。我们得到了𝑎1=𝑎1,𝑏1=𝑏1,𝑎1=0,𝑎01=122𝛼2𝛼𝑎21+9𝑏21,𝑎2=±2𝛼,𝑏0=118𝛼𝑎21+14𝑏21,𝑉=2,(.27)哪里𝑎1𝑏1是自由参数。

设置3。我们得到了𝑎2=𝑎2,𝑏1=𝑏1,𝑎1=0,𝑎0=𝑏212𝛼𝑎22+98𝛼𝑎2,𝑎1=𝑏1𝛼𝑎22+9𝛼𝑎2,𝑏0=𝑏212𝛼𝑎22+98𝛼𝑎22,𝑉=𝛼𝑎22+,(.28)哪里𝑎2𝑏1是自由参数。

设置4。我们得到了𝑎1=𝑎1,𝑎1=0,𝑎0=0,𝑎2=0,𝑏0=𝛼𝑎2172,𝑏1=0,𝑉=,(.29)哪里𝑎1是一个自由参数。
使用(3.26)进入(3.25)简化后,我们得到𝑢(美国)=±62𝛼12𝑏0(co(o)小时美国n个小时美国)1+𝑏0co(o)小时美国+1𝑏0n个小时美国.(.0)如果𝑏0=1,来自(3.30),我们得到𝑢(美国)=±6𝑖2𝛼t吨n个小时美国,(.1)哪里美国=𝑥+𝑦+𝑧+6𝑡.
替换(3.27)进入(3.25)简化后,我们得到𝑢(美国)=±2𝛼1+12±𝑎12𝛼+𝑏122𝛼𝑎21+9𝑏21(co(o)小时美国n个小时美国)6+2𝛼𝑎21+9𝑏21co(o)小时美国+62𝛼𝑎219𝑏21n个小时美国+6𝑏1.(.2)如果𝛼是负的,也就是说,𝛼=𝛽,𝛽>0,𝑏1=2𝑎1=0, (3.32)可以简化为𝑢(美国)=±美国2𝛽t吨n个小时2,(.)哪里美国=𝑥+𝑦+𝑧+(/2)𝑡.
替换(3.28)进入(3.25)简化后,我们得到𝑢(美国)=𝑎21+72𝑏18𝛼𝑎22+𝑏212𝛼𝑎22+9co(o)小时美国+8𝛼𝑎22𝑏212𝛼𝑎22+9n个小时美国+8𝛼𝑎22𝑏1.(.4)如果𝑏1=1,𝛼=6、和𝑎2=1/2, (3.34)成为1𝑢(美国)=2+1+2co(o)小时美国,(.5)哪里美国=𝑥+𝑦+𝑧(+𝛼𝑎22)𝑡.
使用(3.29)进入(3.25)简化后,我们得到𝑢(美国)=72𝑎172+𝛼𝑎21co(o)小时美国+72𝛼𝑎21n个小时美国.(.6)如果𝑎1=、和𝛼=8, (3.36)成为𝑢(美国)=2e(电子)c美国,(.7)哪里美国=𝑥+𝑦+𝑧𝑡.

案例3。我们设置了𝑝=𝑐=2𝑞=𝑑=2.
对于这种情况,试验解决方案(2.5)减少到𝑎𝑢(美国)=2𝑒2美国+𝑎1𝑒美国+𝑎0+𝑎1𝑒美国+𝑎2𝑒2美国𝑏2𝑒2美国+𝑏1𝑒美国+𝑏0+𝑏1𝑒美国+𝑏2𝑒2美国.(.8)因为,在(3.38),我们可以考虑𝑏2=1,𝑎2=0,𝑏2=0、和𝑏1=0因此(3.38)减少到(3.25). 这表明案件等同于案例2.方程式(3.38)可以重写为𝑎𝑢(美国)=2𝑒美国+𝑎1+𝑎0𝑒美国+𝑎1𝑒2美国+𝑎2𝑒美国𝑏2𝑒美国+𝑏1+𝑏0𝑒美国+𝑏1𝑒2美国+𝑏2𝑒美国.(.9)
如果我们把𝑎2=0,𝑎1=0,𝑏2=1,𝑏2=0,𝑏1=0到(3.39),我们得到的解形式为(3.8). 这意味着本案等同于案例1.
此外,如果我们考虑𝑝=𝑐=𝑞=𝑑=可以看出,本案例也与案例等效12.
因此,我们认为没有必要再次寻找解决方案。

值得注意的是,解决方案(3.17)和(3.31)是相同的,解决方案(3.19)和(3.33)是相同的,解决方案(3.21)和(3.35)是相同的,并且解决方案(3.23)和(3.37)都是相同的。

超出表格范围1,扎耶德[52]得到另一个解(3.39)。但是,我们得到了另外两个新的解决方案(3.21)和(3.23).

解决方案的图形表示
图中借助Maple 13显示了上述解决方案。

4.结论

利用Exp-function方法,借助符号计算软件Maple 13,构造了(3+1)维修正KdV-Zakharov-Kuznetsev方程的新的精确行波解。重要的是,获得的一些解决方案与文献中可用的解决方案相同,还有一些是新的。这些解可以用来描述对复杂物理现象的洞察力。

确认

本文得到了USM短期拨款304/PMATHS/6310072的支持,作者对USM数学科学学院提供的相关研究设施表示感谢。作者还感谢裁判的宝贵意见和建议。