基于matlab非线性薛定谔方程数值算法的研究.pdfVIP

基于MATLAB的非线性薛定谔方程的数值算法研究 雷超 四川农业大学农业工程系，四川雅安 （625014） E­mail（leichaocool123@） 摘 要：本文简单介绍了一种基于MATLAB 来求解非线性薛定谔方程的分步傅立叶数值算 法。非线性薛定谔方程是由线性算符和非线性算符两部分组成，依据分步傅立叶数值算法的 思想，当光场每通过一段微小距离的时候，先发生线性算符代表的色散数学方程，再计算非 线性算符代表的非线性效应数学方程，从而得到近似的数值结果。在计算存在高阶偏微分项 时，在时域中计算非常不便，可利用傅立叶变换，把偏微分方程变换为代数方程来运算，进 而求出其相应的数值结果。从算法的编写角度出发，该算法在MATLAB中易于实现，而且 比较简单清晰，能够广泛的应用在非线性薛定谔方程的数值求解之中。 关键词：非线性薛定谔方程；傅立叶变换；分步傅立叶数值算法；MATLAB 中图分类号：O415 1. 引 言 非线性薛定谔方程是研究光脉冲在光纤中传输的基本方程，大多数文献都是直接引用非 线性薛定谔方程，采用分步傅里叶法数值求解非线性薛定谔方程。本文从通过MATLAB编 程来求解非线性薛定谔方程，并给出了分步傅里叶法数值求解算法程序。 MATLAB是美国MathWorks公司自20世纪80年代中期推出的数学软件，优秀的数值计算 能力和卓越的数据可视化能力使其很快在数学软件中脱颖而出。随着版本的不断升级，它在 数值计算及符号计算功能上得到了进一步完善。到目前为止，MATLAB已经发展成为多学 科、多种工作平台的功能强大的大型软件。 2. 非线性薛定谔方程 非线性薛定谔方程，简称NLS 方程，是一个非线性偏微分方程，通常情况下无法求出 解析解，只能求出它的数值解。在含有非线性色散介质地脉冲传输问题中应用非常广泛。 非线性薛定谔方程的基本形式为: 2 iu  u  2|u| u （2­1） t xx 其中u是未知的复值函数. 3. 分步傅立叶数值算法 目前，采用分步傅立叶算法(SplitstepFourierMethod)求解非线性薛定谔方程的数值解应 用比较多。分步傅立叶方法最早是在1937年开始应用的，这种方法己经被证明是相同精度下 数值求解非线性薛定愕方程最快的方法，部分原因是它采用了快速傅立叶变换算法(Fast FourierTransformAlgorithm)。基于MATLAB科学计算软件以及MATLAB强大的符号计算功 能，完全可以实现分步傅立叶数值算法来对脉冲形状和频谱进行仿真。 ­ 1 ­ 一般情况下，光脉冲沿光纤传播时受到色散和非线性效应的共同作用，假设当传输距离 很小的时候，两者相互独立作用，那么，根据这种思想可建立如下分步傅立叶数值算法的数 学模型： 把待求解的非线性薛定谔方程写成以下形式: U ˆ ˆ  (D N)U （3­1） z ˆ ˆ D 其中 是线性算符，代表介质的色散和损耗，N 是非线性算符，它决定了脉冲传输过程中 光纤的非线性效应。 一般来讲，沿光纤的长度方向，色散和非线性是同时作用的。分步傅立叶法假设在传输 过程中，光场每通过一小段距离h，色散和非线性效应可以分别作用，得到近似结果。也就 z 是说脉冲从 到z h的传输过程中分两步进行。第一步，只有非线性作用，方程(2)式中的 ˆ ˆ