matlab中孤子解,Matlab在数学物理方法中应用

论 著 16 Matlab 在数学物理方法中应用 杨 红 黄勇刚 邓 科 王小云 (吉首大学物理与机电工程学院，湖南 吉首 416000) 摘 要：Matlab 在数学物理方法教学中起到了很重要的作用。本文利用 Matlab 中的快速傅里叶变换，研究 了与数学物理方法中的非线性薛定谔方程解的回归问题，并探讨如何培养本科生的科研创新能力。 关键词：Matlab；数学物理方法；非线性 DOI: 10.3969/j.issn.1671-6396.2015.07.005 1 引言 Matlab (Marrix laboratory) 初是由新墨西哥大学 Cleve Moler 博士教授开发的，用于辅助线性代数课程教学。后来其发展成一门高性能的计算语言，并有着强大的数值计算和工程运算能力，强大的科学数据可视化能力。作者在从事本科生的数学物理方法教学中发现，若能在学习的过程中引入 Matlab，将抽象的，难理解的函数用数值仿真，这将使得问题更加直观明了，更容易接受。并且，使用 Matlab 可以比较容易的解决一些较难的问题。本文利用Matlab 中的快速傅里叶变换，组织学生一起研究了与数学物理方法中的非线性薛定谔方程解的回归问题。希望这种方式，为学生营造一种科研氛围，较早接触科研，并能一起讨论，加深对问题的理解。 2 非线性薛定谔方程及解的回归 三维立方非线性薛定谔方程为： 2 0t xx yy zzi ψ α ψ β ψ γ ψ λ ψ ψ + + + + = 其中 (),,,x yzt ψ 为复函数， ,,, α βγλ 为实常 量。当 11 , 1, , 1 22 αβγλ = = − = =时，我们考虑关于以上方程的周期性边界条件解，并假设初始条件如下[1-3]： ( ) ( ) ( ) ( ){ } ( ) ( ) ( ){} 1/2 2 2 2 1 1 1 1/2 2 2 2 1 1 1 , , ,0 sec 5 sec 5 x y z h x x y y z z h x x y y z z ψ ⎡⎤ = + + + + + ⎣⎦ ⎡⎤ + − + − + − ⎣⎦ 其中 x1=1.1325，y1=1.1325，z1=1.1325。为了模拟这个方程解的运动情况，我们使用了分裂步长傅立叶法，其中的快速傅里叶变换直接可以用 Matlab 中的相关子程序。 首先将非线性薛定谔方程分解为线性和非线性两部分： 考虑时间间隔[tl，tl+1]，时间的步长为Δt，非线性部分可以精确求解： )),(exp(),(),(~ 2 1ttxitxtx lll Δ= + ψλ ψψ 得到的结果作为解线性部分的初始条件，线性部分用离散快速傅立叶变换(FFT) 来解。 概括的说为： 分裂步长傅立叶法之所以效率高是因为在线性部分的计算中采用了快速傅立叶变换(FFT)。傅立叶变换使得在实空间的求导等价于在傅立叶空间的乘以一个变量。快速傅立叶变换优点是将傅立叶变换中须要进行的大部分乘法计算转换为加法计算，而在计算机中进行一次乘法运算的时间是进行一次加法计算的许多倍，从而提高了运算速度[4-5]。从图1 中可以发现，两个孤子解在 t 为 T/4 时(T 为一个完整的演化周期)，分解为 4 个孤子解；当 t 为 T/2 时，解又变回了两个，但是它们的位置发生了改变，处在不同的对角线上；当 t 为 T 时，解又回归到原来的位置。 为了进一步研究此回归现象，我们定义三维振幅偏移量函数为： () () 111 1 ,,, ** NNN i j k i j k t x y z t NNN φ