谱方法和边界值法求解二维薛定谔方程

摘  要 薛定谔方程是物理系统中量子力学的基础方程，它可以清楚地说明量子在系统中随时间变化的规律。通过求解微观系统所对应的薛定谔方程，我们能够得到其波函数以及对应的能量，从而计算粒子的分布概率，进一步来了解其性质。在化学和物理等诸多科学研究领域当中，薛定谔方程求解的结果都与实际很相符。近年来，很多学者通过各种方法研究具有复杂势函数的薛定谔方程，解释了很多重要的物理现象，因此对薛定谔方程的求解具有相当重要的意义。 本文主要是用Galerkin-Chebyshev谱方法和边界值法求解二维薛定谔方程。首先运用Galerkin-Chebyshev谱方法来对空间导数进行近似，离散二维薛定谔方程，从而将原问题转化为复数域上的线性常微分方程组。然后用边界值法求解该方程组，所求得的数值解即为原问题的解，之后进行误差分析。最后利用Matlab进行数值模拟，给出数值解的图像以及误差曲面图像，结果显示此方法精度高且具有很好的稳定性。 关键词：薛定谔方程；Galerkin-Chebyshev谱方法；边界值法；数值解； 精度高；稳定 Abstract The Schrödinger equation is the basic equations of quantum mechanics in the physical system. It can clearly describe the regular of the quantum evolves over time. By solving the Schrödinger equation which the micro system correspond, we can get the wave function and energy, and thus calculate the probability distribution of the particles, further understand the nature of it. In chemistry, physics and other fields of scientific research, the  results of solving the Schrodinger equation are basically consistent with the actual. In recent years,  many researchers used a variety of methods to investigate the Schrödinger equation with complex potential function, and explained a lot of important phenomena. Thus solving the Schrödinger equation has very important significance. The main purpose of this paper is to solve the two dimensional Schrödinger equation through the Galerkin-Chebyshev spectral method and the boundary value method. First we use the spectral method to approximate the spatial derivation, discretize the two dimensional Schrödinger equation，and transform the original problem into a set of linear ordinary differential equations in  the complex number field. Then by using the boundary value method to solve the equations, that the numerical solutions is the solutions of the original problem, and then analyze the error. Finally we use Matlab to conduct the numerical simulation, and give the images of the numerical solutions and errors, which show that the methods have high precision and good stability. Keywords: Schrödinger equation, Galerkin-Chebyshev spectral method, boundary value method, numerical solutions, high precision, stability 目  录 摘  要    I Abstract    II 第1章 绪  论    1 1.1 课题研究的背景和意义    1 1.2 国内外研究现状    2 1.3 本文的主要研究内容    2 第2章 预备知识    4 2.1 克罗内克积的简介    4 2.2 Chebyshev多项式介绍及其性质    5 2.3 Chebyshev正交逼近的性质    6 2.4 投影算子的性质    7 2.5 本章小结    8 第3章 Galerkin-Chebyshev谱方法和边界值法    9 3.1 用Galerkin-Chebyshev谱方法求解椭圆型方程    9 3.2 用边界值法求解常微分方程    10 3.3本章小结    14 第4章 求解二维薛定谔方程    15 4.1 区域和边界条件的处理    15 4.1.1 区域的处理    15 4.1.2 边界条件的处理    17 4.2 二维薛定谔方程的求解    20 4.3 误差分析    21 4.4 本章小结    26 第5章 数值模拟    27 结  论    32 参考文献    33 哈尔滨工业大学学位 论文 政研论文下载论文大学下载论文大学下载关于长拳的论文浙大论文封面下载 原创性声明及使用授权说明    37 致  谢    38 第1章 绪  论 1.1 课题研究的背景和意义 薛定谔方程是一个偏微分方程，它可以清楚地说明量子在物理系统中随时间如何在变化，它是量子力学的一个基本的假设，也是量子力学的基础方程，由物理学家薛定谔提出而得名[1]。在经典力学和量子力学当中，人们分别是用牛顿第二定律和薛定谔方程来描述物体的运动的，这两者在物理系统当中具有相同的地位。薛定谔方程式可以描述任何的微观系统，通过求解该微观系统所对应的薛定谔方程，我们能够得到其波函数以及对应的能量，从而进一步来了解该微观系统的性质。 薛定谔方程可以分为与时间有关和与时间无关两种类型，其中量子系统的波函数随着时间的演化过程是通过与时间有关的薛定谔方程来描述的，而与时间无关的薛定谔方程则描述的是固定状态的量子系统的物理性质，方程的解即是该量子系统固定状态的波函数。 本文考虑的是二维与时间有关的薛定谔方程：       (1-1) 初始条件为： 边界值条件为：                           (1-2) 其中是任意的势函数，，是波函数，且在定义域内连续。 薛定谔方程是反应微观粒子随着时间变化的非相对论波动函数，它仅适用于速度比较缓慢的非相对论粒子。其中波函数可以很好地描述微观粒子的状态，在势函数中微观粒子运动的薛定谔方程即为方程(1-1)。我们可以通过给定的初始条件和边界值条件以及波函数所满足的条件，来求解出波函数，进而计算粒子的分布概率。 薛定谔方程被广泛地应用于化学和物理等领域中，如量子器件的建模[2]，光纤传播模型[3]，光电子器件的 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 [4]，电磁波的传播[5]，天体系统的量子化[6]，轴近似波动方程的水下声学[7]，量子动力学计算的应用[8,9]，化学核外电子的运动状态描述[10]等。它被应用到原子、核等诸多方面问题中，所得到的结果都与实际很相符。近年来，很多学者通过各种方法研究具有复杂势函数的薛定谔方程所描述的问题[11-14]，解释了很多重要的物理现象，因此对薛定谔方程的求解具有相当重要的意义。 1.2 国内外研究现状 到目前为止，对薛定谔方程(1-1)的求解已经有了很多种数值方法， 大多都是采用的有限差分法[15-17]，或者是用三角正交函数系或幂级数函数展开的谱方法[18,19]。Subasi给出了具有二阶精度的有限差分方法[20]，Kalita等人给出了一个隐式的半离散高阶紧凑方法[21]，Antonie等人给出了一个Crank-Nicolson隐格式方法[22]，Dehghan给出了不同的有限差分方法包括三个全隐式和两个全显示差分方法以及交替方向隐式法和Barakat和Clark类型的显示方法[23]，Dehghan和Shokri还给出了使用配置和薄板样条径向基函数的数值方法[24]，此外Dehghan和Mohebbi还给出了求解方程(1-1)的紧凑有限差分法[25]，Gao和xie还给出了紧凑的交替方向隐式有限差分法[26]，该方法在空间上具有四阶精度，在时间上具有二阶精度。Li等人还给出了多元二次（MQ）和薄板样条（TPS）径向基函数的MPS方法求解薛定谔方程，该方法类似于有限差分法[27]。Dehghan和Taleei还提出了一种紧凑的分布有限差分方法来求解薛定谔方程[28]，该方法通过使用四阶精度紧致差分格式，来提高分布有限差分方法的准确性，而且还具有无条件稳定的性质。 谱方法的思想起源于傅立叶分析，它是一种既古老又新兴的求解偏微分方程的方法。求解偏微分方程的三种最基本的方法分别是谱方法，有限差分方法和有限元方法。谱方法和另两种方法相比，具有“无穷阶”收敛的特点，即它的收敛速度会随着真解的光滑程度变高而变快，从而谱方法就能用限制自由度的方式来得到较高的精度[29]，另两种方法在这一点上是无法比拟的。 1.3 本文的主要研究内容 本文主要是用Galerkin-Chebyshev谱方法和边界值法求解二维Schrödinger方程，运用Galerkin-Chebyshev谱方法对空间导数进行近似,离散薛定谔方程(1-1)，从而将原问题转化为复数域上的线性常微分方程组，然后再用边界值法求解该方程组，所求得的数值解即为原问题的解，之后再进行误差分析，得到误差分析结果，最后再通过Matlab进行数值模拟，给出数值解的图像以及误差曲面图像。谱方法求解偏微分方程具有高精度的性质，边界值法求解常微分方程同样具有高精度和稳定的特点，这样问题即得到解决。 在第一章中我们阐述了薛定谔方程在当前科学研究中的应用， 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 明求解薛定谔方程具有很深远的意义，还介绍了现阶段求解该方程的主要方法，以及本文即将采用的方法。紧接着在第二章中，我们介绍了本论文所需要的一些预备的基础知识，为后面论文的顺利进行做好准备工作。 在第三章当中，我们采用Galerkin-Chebyshev谱方法求解椭圆型方程，以及用边界值法求解常微分方程，并给出求解特殊常微分方程组的求解格式，这两个方法求解微分方程都具有很高的精度和很好的稳定性。 第四章中，先对原问题进行区域映射处理，以及对边界条件进行齐次化处理以后，然后运用Galerkin-Chebyshev谱方法对二维薛定谔方程进行离散，将其转化成常微分方程组，然后对该微分方程组进行求解，得到数值解，接着对该方法进行误差分析，得到误差估计结果。 第五章进行数值模拟，根据前面的内容，编程得到问题的数值解，并和相应的精确解进行比较，分析其误差，画出误差曲面图像。     最后是本文的一个总结，以及研究此问题的意义和前景展望。 第2章 预备知识 2.1 克罗内克积的简介 定义2.1[30]：设是一个行列的矩阵，，是一个行列的矩阵，，克罗内克积可以表示成： 它是一个的分块矩阵。 克罗内克积具有如下的一些性质： 性质1：满足结合律与双线性的性质： 如果矩阵存在，则 ； 如果矩阵存在，则；     ，其中是常数； . 性质2：，，和是四个矩阵，如果矩阵乘积和存在，那么就有 性质3：是可逆的当且仅当和是可逆的，其逆矩阵是： 性质4：. 定义2.2：设是一个行列的矩阵，，那么把矩阵按列将后一列堆在前一列后面，形成的一个新的列的向量记为，即： 定理2.1：设是一个行列的矩阵，是一个行列的矩阵，是一个行列的矩阵，也是一个行列的矩阵，那么有： 证明：先将矩阵,,写成如下的形式： 其中,,分别是矩阵,,第列的列向量，， 则有： , 从而原题得证。 2.2 Chebyshev多项式介绍及其性质 定义2.3：在区间上的权函数以递归的形式定义的正交多项式称为Chebyshev多项式，它可写成：。 Chebyshev多项式具有如下的性质： 性质(1) [31]：正交性 性质(2) [31]：递推关系 性质(3) [31]：是阶多项式，是阶多项式，是阶多项式，满足：，，其中 。 定理2.2[32]：设 ，则: 2.3 Chebyshev正交逼近的性质 我们讨论Chebyshev逼近问题，需要借助带权的Sobolev空间，下面记以为权的阶空间为，它的内积和范数定义分别为 记。 设区间是一个非空集，且是Lebesgue可测的，记的范数为： 当时，。     接下来定义空间，设空间是有界的，且，有 在空间和上的全体次连续可微的函数所构成的集合分别记为和。记，，，其中是广义导算子，接下来定义弱导数。 定义2.4[33]：设，满足上面的式子，称是的阶弱导数，记为，如果 有时又称在弱的意义下。 下面定义空间，设区域是有界的，是非负整数，，，有： 其中为空间上面的范数。     在空间上的闭包记为，当时，。 定理2.3（Gronwall不等式）[34]：设和是上的非负的连续函数，并且在是可微的，如果存在常数满足，使得对任意的，都有： 或者等价的还有： 那么就有： 2.4 投影算子的性质 记是一个多项式空间，其最大自由度是，是到的正交投影算子，，是到的椭圆投影算子，则有如下的定义和性质定理： 定义2.5[35]：空间 中从到的正交投影算子为： 定义2.6[35]: 空间 中从到的椭圆投影算子为： 定理2.4[36]：对任意的非负整数，，都有下面的不等式： 定理2.5[37]：对任意的非负整数，，都有下面的不等式： 2.5 本章小结 本章给出了完成这篇论文所需要的一些必备的基础知识，首先介绍了克罗内克积的定义以及性质，然后介绍了切比雪夫多项式性质以及一些重要的关系定理，之后介绍了Chebyshev正交逼近的性质，其中包括内积和范数的定义，空间的定义和性质，最后还介绍了投影算子的定义和不等式性质，为论文的进行做好准备工作。 第3章 Galerkin-Chebyshev谱方法和边界值法 3.1 用Galerkin-Chebyshev谱方法求解椭圆型方程 考虑用Galerkin-Chebyshev谱方法来求解如下的椭圆型方程             (3-1) 边界条件是： 由Chebyshev多项式的定义和性质，设，，，则方程(3-1)的Galerkin-Chebyshev谱方法是求使得对任给的都满足 其中，，。 令,取 则有 将上式用矩阵表示即可写成 其中：，，和满足定理2.2中的条件关系， ，，且。   由定理2.1有         (3-2) 对方程(3-2)进行求解，就可以求出其数值解，从而得到方程(3-1)的数值解。 3.2 用边界值法求解常微分方程 边界值法是最近求解常微分方程数值解的常用方法，简称为BVMs，它是线性多步法的一个推广，和其他常微分初值问题的数值解法相比较，BVMs具有高精度和无条件稳定的特点[38-41]，是一个很好的方法。 考虑下面的初值问题               (3-3) 用步线性多步法离散上面的方程即可得到             (3-4) 其中，为系数。 由泰勒展开有： 从而令： 则有：         (3-5) 如果有次的连续微商，那么就可以选取和使得，即选取使其满足           (3-6) 此时就有 其中为截断误差，略去，就得到了线性多步法(3-4)，该方法的精度是阶的。 求解方程(3-4)需要个初始边界条件和个结尾边界条件，即我们需要和，初始边界条件可以由方程(3-3)得到。个初始边界条件和个结尾边界条件则来自于以下等式             (3-7) 和       (3-8) 其中系数和的选择，要满足使基于最初与最后的边界条件的方法的截断误差与基于 公式 小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载 (3-6)的方法具有相同的阶。 方程(3-4)—(3-8)用矩阵形式表示可以写为 其中. 用代替矩阵中的，即为矩阵，并且. 对进行划分 将第一列分离出来，可以得到的等价式                         (3-9) 其中是一个未知量，且有 在这里我们用四阶BVMs近似方程(3-3)，并取，由(3-6)得到 求解该方程组，得到其基础解系，取其中的三组解，一组代入到方程(3-4)，另两组分别代入到(3-7)和(3-8)，即可以得到下面的关系式           (3-10) 其中 。 其对应的初始边界条件和结尾边界条件分别为               (3-11) 其中 。           (3-12) 其中 。 把上面的三个式子化为等式(3-9)的形式，则其中的分别为 ,  如果我们考虑的是特殊的线性常微分方程组               (3-13) 其中是的矩阵,且 那么我们可以将(3-13)化为矩阵形式如下     (3-14) 其中是的单位矩阵，，且 如果线性常微分方程组可以化为下面的形式                 (3-15) 其中是的非奇异矩阵, 那么用四阶BVM法可以将方程(3-15)化为   (3-16) 3.3本章小结 本章给出了利用Galerkin-Chebyshev谱方法求解椭圆型方程的数值方法格式，利用该方法将偏微分方程进行离散以后得到常微分方程组，再利用常微分方程的解法对其进行求解即可达到将偏微分方程进行求解的目的，谱方法精度很高，稳定性也很好，对于求解偏微分方程是一个非常好的方法。还给出了利用边界值法求解微分方程组的过程，并给出了几种特殊形式的微分方程组的边界值法求解的数值格式，边界值法求解微分方程具有很高的精度，对于求解常微分方程组也是一个很好的方法。 第4章 求解二维薛定谔方程 4.1 区域和边界条件的处理 由于Galerkin-Chebyshev谱方法只能解决齐次边值条件的问题，故针对本文的二维Schrödinger方程问题，需要先进行区域的映射处理，对非齐次的边值问题进行齐次化处理，将其转化成方程的一般形式，再进行求解。 4.1.1 区域的处理 原问题中，我们在这里对其进行一定的变换处理，使区域变成。 令 ， 于是方程(1-1)就转化为               (4-1) 初始条件为： 边界值条件为：     对方程(4-1)进行简化，可以将其表示为     (4-2) 初始条件为： 边界值条件为： 其中 4.1.2 边界条件的处理 在本文中对边界条件的处理过程，就是一个对边界条件进行齐次化的过程。 由方程(4-2)，我们很容易得到 ，，， ，，， 于是令 ，则有 令 ，即可得到 同理可得 从而有 再令，则有 同理有 ，， 即 由可以得到 于是方程(4-2)就转化为下面的方程 即可写成               (4-3) 初始条件为： 边界值条件为： 其中 于是原问题就转化成了 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 问题。 接下来我们需要去复数域，令： , 带入方程(4-3)得到 根据复数的性质，即可以得到             (4-4)                     初始条件为： 边界值条件为： 4.2 二维薛定谔方程的求解 对上面的方程组(4-4)进行移项变换，得到         (4-5) 令 ， ， 令，由第三章的方法即可以得到解的弱形式   (4-6)         将(4-6)用矩阵表示即可写成                   (4-7) ， 其中 ，， 利用定理2.1，方程组(4-7)等价于           (4-8) 然后利用第三章的边界值法对上面的常微分方程组进行求解，所得的数值解即为原问题的数值解。 4.3 误差分析 由定理2.4可以得到 ， 那么问题(4-4)的弱形式是求使得           (4-9) Galerkin-Chebyshev谱方法是求,使得 (4-10) 令 ， ,，, 则有 ， 定理4.1：假设，，，满足上面的关系，那么对任意的都有 其中 证明：由定理2.5可以得到 由定理2.4可以得到 同理可得 接下来估计和，令 由方程组(4-9)，(4-10)可以得到     从而有 同理得到 由投影算子的性质，我们可以得到 即 同理有 分别取，分别带入到上面两个式子中，就有 由于，且 故得到 从而由Cauchy-Schwarz不等式有 得到 对上面的式子两边同时积分即得到 其中 从而得到 由于 所以有 即有 从而得到最后的误差分析的结果为 即有 其中 从而得证，得到最优误差估计。 4.4 本章小结 在本章中，首先是对二维薛定谔方程进行区域映射以及对边界条件进行齐次化处理，使其换化成便于Galerkin-Chebyshev谱方法进行求解的标准格式，接着利用Galerkin-Chebyshev谱方法离散空间变量得到常微分方程组，再利用边界值法对该方程组进行求解，即可以得到问题的数值解。最后再对本文所给出的方法进行误差分析，得到误差分析结果。 第5章 数值模拟 在这一章中，我们列举了两个实例进行数值模拟，通过Matlab求出其数值解，然后和精确解进行比较，来验证该方法的精度和稳定性。 对于微分方程(1-1)，数值模拟实例1当中，我们考虑的是，， ，相应的初始条件是： 此问题的精确解是：                 (5-1) 通过精确解(5-1)式即可很快得到此问题的边界值条件，即为 ， ，   针对上述问题，采用Galerkin-Chebyshev谱方法和边界值法对该问题进行求解，取，，，利用matlab编程在实部处求得该问题的数值解见下图5-1，以及得到绝对误差分析曲面图像见下图5-2，再取，，，利用matlab编程在虚部处求得该问题的数值解见下图5-3，以及得到绝对误差分析曲面图像见下图5-4。 图5-1：时刻实部的数值解图像 图5-2：时刻实部的绝对误差分析曲面图像 图5-3：时刻虚部的数值解图像 图5-4：时刻虚部的绝对误差分析曲面图像 对于微分方程(1-1)，数值模拟实例2当中，我们考虑的是，， ，相应的初始条件是： 此问题的精确解是：                           (5-2) 通过精确解(5-2)式即可很快得到此问题的边界值条件，即为 ， ，   针对上述问题，采用Galerkin-Chebyshev谱方法和边界值法对该问题进行求解，取，，，利用matlab编程在实部处求得该问题的数值解见下图5-5，以及得到绝对误差分析曲面图像见下图5-6，再取，，，利用matlab编程在虚部处求得该问题的数值解见下图5-7，以及得到绝对误差分析曲面图像见下图5-8。 图5-5：时刻实部的数值解图像 图5-6：时刻实部的绝对误差分析曲面图像 图5-7：时刻虚部的数值解图像 图5-8：时刻虚部的绝对误差分析曲面图像 上面两个实例的数值模拟，很好的表明了Galerkin-Chebyshev谱方法和边界值法求解问题(1-1)具有很高的精度和很好的稳定性。 结  论 本文主要研究的是用Galerkin-Chebyshev谱方法和边界值法求解二维薛定谔方程，先用Galerkin-Chebyshev谱方法求解椭圆型方程，以及用边界值法求解常微分方程，并给出求解格式，之后结合这两个方法对二维薛定谔方程进行求解，先对它进行区域映射处理，再对边界条件进行齐次化的处理来实现对二维薛定谔方程进行求解，得到其数值解，然后进行误差分析，在离散的加权的范数下我们得到半离散格式的谱精度，得到最优误差估计，最后再进行数值模拟，在数值模拟中所得到的结果显示我们的方法在离散空间具有高阶收敛性以及很好的稳定性。 薛定谔方程化学和物理等领域中都有非常广泛的应用，求解薛定谔方程具有相当重要的意义，本文所采用的方法所求得的数值解的误差很小，为二维薛定谔方程的提供了一个很好的解决方法，为物理和化学等领域的研究提供了一个很好的解决办法。 本文所给出的方法求得的数值解已经具备了较高的精度，文中所采用的边界值法给出的是定步长的格式，要是能实现变步长的边界值法来进行求解，效果应该能达到更好。 参考文献 [1] E.Schrödinger. An Undulatory Theory of the mechanics of atoms and molecules[J]. Physical Review, 1926, 28(6): 1049-1070. [2] Anton Arnold. Numerically absorbing boundary conditions for quantum evolution equations[J]. VLSI Design, 1998, 6(1): 313-319. [3] 李均，黄德修，张新亮. 光纤模型的数值计算研究[J]. 光电子技术与信息, 2003, 16(2):9-12. [4] Weiping Huang, Chenglin Xu, Sai-Tak Chu, Sujeet K.Chaudhuri. The finite-difference vector beam propagation method: analysis and assessment[J].  Journal of Ligntwave Technology, 1992, 10(3): 295-305. [5] Mireille Levy. Parabolic equation methods for electromagnetic wave propagation[M]. London: The Institution of Electrical Engineers, 2000: 57-61. [6] Floretin Smarandache, Vic Christianto. Schrödinger equation and the quantization of celestial systems[J]. Progress in Physics, 2006(2), 63-67. [7] Fred D.Tappert. The parabolic approximation method[J]. Wave Propagation and Underwater Acoustics Lecture Notes in Physics, 1977(70): 224-287. [8] Quan Minh Nguyen. Effects of randomness, dissipation and interaction on solitons of the cubic nonlinear Schrödinger equation and related nonlinear wave models[D]. Buffalo: State University of New York at Buffalo, 2011, 73(1): 351-468. [9] F Y Hajj. Solution of the Schrödinger equation in two and three dimensions[J]. Journal of Physics B:Atomic and Molecular Physics, 1985(18): 1-11. [10] 苑壮东, 考秀娟. 薛定谔方程在化学中的应用[J]. 济宁学院学报, 2008, 29(16): 42-43. [11] Simos T E. New high order multiderivative explicit four-step methods with vanished phase-lag and its derivatives for the approximate solution of the Schrödinger equation.Part I: construction and theoretical analysis[J]. Journal of Mathematical Chemistry, 2013, 51(1): 194-226. [12] 丁凌, 肖氏武, 姜海波. 非齐次边界条件下的具有复合级数非线性项的薛定谔方程[J]. 西南师范大学学报, 2011, 36(5): 30-34． [13] Hillel Tal-Ezer, Ronnie Kosloff, Ido Schaefer. New，highly accurate propagator for the linear and nonlinear Schrödinger equation[J]. Journal of Scientific Computing, 2012, 53(1): 211-221. [14] Othmar Koch, Christof Neuhauser, Mechthild Thalhammer. Error analysis of high-order splitting methods for nonlinear evolutionary Schrödinger equations and application to the MCTDHF equations in electron dynamics[J]. Mathematical Modeling and Numerrical Analysis, 2013, 47(5): 1265-1286. [15] Mehdi Dehghan, Ameneh Taleei. A compact split-step finite difference method for solving the nonlinear Schrödinger equations with constant and variable coefficients[J]. Computer Physics Communications, 2010, 181(1): 43-51. [16] Linghua Kong, Yali Auan, Lan Wang, Xiuling Yin, Yuanping Ma. Spectral-like resolution compact ADI finite difference method for the multi-dimensional Schrödinger equations[J]. Mathematical and Computer Modelling, 2012, 6(55): 1798-1812. [17] J.D.Cooper, A.Valavanis, Z. lkoni, P.Harrison, J.E.Cunningham. Finite difference method for solving the Schrödinger equation with band nonparabolicity in mid-infrared quantum cascade lasers[J]. Jonurnal of Applied Physics, 2010(108): 113109. [18] M.D Feit, J.A Fleck Jr, A Steiger. Solution of the Schrödinger equation by a spectral method[J]. Journal of Computational Physics, 1982, 47(3): 412-433. [19] Kosloff. A fourier method solution for the time dependent Schrödinger equation as a tool in molecular dynamics[J]. Journal of Computational Physics, 1983, 52(1): 35-53. [20] Murat Subasi. On the finite differences schemes for the numerical solution of two dimensional Schrödinger equation[J]. Numerical Methods for Partial Differential Equations, 2002, 18(6): 752-758. [21] Jiten C.Kalita, Puneet Chhabra, Sudhanshu Kumar. A semi-discrete higher order compact scheme for the unsteady two-dimensional Schrödinger equation[J]. Journal of Computational and Applied Mathematics, 2006, 197(1): 141-149. [22] Xavier Antoine, Christophe Besse, Vincent Mouysset. Numerical schemes for the simulation of the two-dimensional Schrödinger equation using non-reflecting boundary conditions[J]. Mathematics of Computation, 2004(73): 1779-1799. [23] Mehdi Dehghan. Finite difference procedures for solving a problem arising in modeling and design of certain optoelectronic devices[J]. Mathematics and Computers in Simulation, 2006, 71(1): 16-30. [24] Mehdi Dehghan, Ali Shoki. A numerical method for two-dimensional Schrödinger equation using collocation and radial basis functions[J]. Computers & Mathematics with Applications, 2007, 54(1): 136-146. [25] Akbar Mohebbi, Mehdi Dehnhan. The use of compact boundary value method for the solution of two-dimensional Schrödinger equation[J]. Journal of Computational and Applied Mathematics, 2009, 225(1): 124-134. [26] Zhen Gao, Shusen Xie. Fourth-order alternating direction implicit compact finite difference schemes for two-dimensional Schrödinger equations[J]. Applied Numerical Mathematics, 2011, 61(4): 593-614. [27] Tao Li, Guodong Wang, Ziw Jiang. A numerical method for two dimensional Schrödinger equation using MPS[J]. Information Computing and Applications, 2011(243): 44-51. [28] Mehdi Dehghan, Ameneh Taleei. A compact split-step finite difference method for solving the nonlinear Schrödinger equations with constant and variable coefficients[J]. Computer Physics Communications, 2010, 181(1): 43-51. [29] Mehdi Dehghan, Farhad Fakhar-lzadi. The spectral collocation method with three different bases for solving a nonlinear partial differential equation arising in modeling of nonlinear waves[J]. Mathematical and Computer Modelling. 2011, 53(9): 1865–1877. [30] Roger A. Horn, Charles R. Johnson. Topics in Matrix Analysis[M]. Cambridge : Cambridge University Press, 1991: 242-260. [31] Jie Shen, Tao Tang. Spectral and high-order methods with applications[M]. Beijing: Science Press, 2006: 15-18. [32] J.Shen. Efficient spectral-Galerkin method II. Direct solvers of second- and fourth-order equations using Chebyshev polynomials[J]. SIAM Journal on Scientific Computing, 1995, 16(1): 74-87. [33] 王明新. 索伯列夫空间[M]. 北京：高等教育出版社, 2013: 40-41. [34] 向新民. 谱方法的数值分析[M]. 北京：科学出版社, 2000: 44-45. [35] 黄明游. 发展方程数值计算方法[M]. 北京：科学出版社, 2004: 92-106. [36] Benyu Guo. Spectral methods and their applications[M]. Hong Kong: World Scientifici Singapore, 1998. [37] Jie Shen, Tao Tang, Lilian Wang. Spectral methods[M]. Berlin: Springer-Verlag Berlin and Heidelberg GmbH & Co. K. 2011. [38] Haiwei Sun, Jun Zhang. A high-order compact boundary value method for solving one-dimensional heat equations[J]. Numerical methods for partial differential equations, 2003, 19(6): 846-857. [39] L.Brugnano, D.Trigiante. Stability properties of some boundary value methods[J]. Applied Numerical Mathematics, 1993, 13(4): 291-304. [40] L.Brubnano, D.Trigiante. Solving differential problems by multistep initial and boundary value methods[M]. Amsterdam: Gordon and Breach Science Publishers, 1998. [41] Donald G Truhlar. Finite Finite difference boundary value method for solving one-dimensional eigenvalue equations[J]. Journal of Computational Physics, 1972, 10(1): 123-132. 文档已经阅读完毕，请返回上一页！