基于maple数学软件求schrodinger方程的精确解

《基于maple数学软件求schrodinger方程的精确解》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、基于Maple数学软件求Schrodinger方程的精确解摘要：指数函数法是求非线性偏微分方程精确解的简便、有效地方法之一，利用这种方法可以求出一大类非线性偏微分方程的精确解，是进一步获得具有物理背景的孤子解、周期解等的有效途径。本文用Maple数学软件实现了指数函数法的过程，并求得了非线性Schrodinger方程的一类精确解。关键词：Maple数学软件；指数函数法；非线性Schrodinger方程；精确解中图分类号：0241文献标识码：A文章编号:1007-9599(2011)23-0000-01RequirementsoftheExactSolutionofSchrodinge

2、rEquationBasedonMapleMathema廿calSoftwareXiaHongming(SchoolofMathematicsandStatisticsTianshuiNormalUniversity,Tianshui741001,China)Abstract:Theexponentialfunctionmethodistofindexactsolutionsofnonlinearpartialdifferentialequationissimple,effectivemethods,theuseofthismethodcanfindalargeclassofexac

3、tsolutionsofnonlinearpartialdifferentialequationsisobtainedwithaphysicalbackgroundforfurthersolitonsolutions,periodicsolutionsandothereffectiveway.Inthispaper,mathematicalsoftwareMapleexponentialfunctionmethodintheprocess,andobtainthenonlinearSchrodingerequationofaclassofexactsolutions.Keywords

4、:Maplemathematicalsoftware;Exponentialfunctionmethod;NonlinearSchrodingerequation;Exactsolution一、引言非线性偏微分方程的求解是各门学科中都会遇到的问题，一般来说，只有在非常特殊的情况下才能求得有显式表达式的精确解，如孤子解通常是偏微分方程的行波解，目前已有多种求显式行波解的方法，如反散射方法、Bäcklund变换、Hirota双线性方法、指数函数法等，其中指数函数法是近年来求非线性偏微分方程精确解的简便而有效地方法之一，用这种方法可以求出一大类非线性偏微分方程的精确解。二、指数函

5、数法简述对于一个给定的非线性偏微分方程(1)为求方程的行波解，令其中为待定常数，将(2)带入(1)约化成的常微分方程.(3)令,(4)其中和和是待定常数.和为待定常数.利用齐次平衡法确定与与的关系，并且通过解代数方程组能够确，就求出了偏微分方程的精确行波解.三、求非线性Srodinger方程精确解的Maple实现非线性薛定铐方程是数学物理中一类重要的非线性演化方程，在量子力学、非线性光学、电磁学等众多领域有着广泛的应用，故对该方程的求解进行研究有着重要的物理意义.本文使用指数函数法、借助于Maple求非线性薛定铐方程的精确解.通过变换(5)可将非线性薛定谴方程(6)化为•(7)下面是

6、求解方程(7)的Maple源代码：>restart；>eql：=-r(x)*omega+diff(r(x),x,x)+2*(l+B)*(r(x))A3：>r(x)：=(a*exp(alpha*x)A2+b*exp(alpha*x)+c)/(f*exp(alpha*x)A2+g*exp(alpha*x)+h)：>eql：二numer(eql)：>eql：=collect(eql,exp(alpha*x))：>sysl：=coeffs(eql,exp(alpha*x))：>sysl：={-omega*c*hA2+2*cA3+2*B*cA3,-omega*a*hA2+6*a*cA2+6*B

7、*a*cA2-omega*c*gA2-2*omega*c*f*h+alphaA2*c*gA2-alphaA2*b*g*h+6*B*bA2*c-4*alphaA2*c*f*h+4*alphaA2*a*hA2-2*omega*b*g*h+6*bA2*c,-omega*b*hA2-2*omega*c*g*h-alphaA2*c*g*h+6*b*cA2+6*B*b*cA2+alphaA2*b*hA2,-omega*b*gA2+3*alphaA2*c*f*g+2*B