谱方法和边界值法求解二维薛定谔方程硕士学位论文1

1、硕士学位论文硕士学位论文 谱方法和边界值法求解二维薛定谔方程谱方法和边界值法求解二维薛定谔方程 摘摘 要要 薛定谔方程是物理系统中量子力学的基础方程，它可以清楚地说明量子在系 统中随时间变化的规律。通过求解微观系统所对应的薛定谔方程，我们能够得到 其波函数以及对应的能量，从而计算粒子的分布概率，进一步来了解其性质。在 化学和物理等诸多科学研究领域当中，薛定谔方程求解的结果都与实际很相符。 近年来，很多学者通过各种方法研究具有复杂势函数的薛定谔方程，解释了很多 重要的物理现象，因此对薛定谔方程的求解具有相当重要的意义。 本文主要是用 galerkin-chebyshev 谱方法和边界值法求解二维

2、薛定谔方程。 首先运用 galerkin-chebyshev 谱方法来对空间导数进行近似，离散二维薛定谔方 程，从而将原问题转化为复数域上的线性常微分方程组。然后用边界值法求解该 方程组，所求得的数值解即为原问题的解，之后进行误差分析。最后利用 matlab 进行数值模拟，给出数值解的图像以及误差曲面图像，结果显示此方法精度高且 具有很好的稳定性。 关键词：薛定谔方程；galerkin-chebyshev 谱方法；边界值法；数值解； 精度高；稳定 abstract the schrdinger equation is the basic equations of quantum mechani

3、cs in the physical system. it can clearly describe the regular of the quantum evolves over time. by solving the schrdinger 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 nat

4、ure 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 schrdinger equation with complex potential function, and

5、explained a lot of important phenomena. thus solving the schrdinger equation has very important significance. the main purpose of this paper is to solve the two dimensional schrdinger equation through the galerkin-chebyshev spectral method and the boundary value method. first we use the spectral met

6、hod to approximate the spatial derivation, discretize the two dimensional schrdinger 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 soluti

7、ons 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: schrdinger equation, galerkin-

8、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 章 galer

9、kin-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 课题研究的背景

10、和意义 薛定谔方程是一个偏微分方程，它可以清楚地说明量子在物理系统中随时间 如何在变化，它是量子力学的一个基本的假设，也是量子力学的基础方程，由物 理学家薛定谔提出而得名1。在经典力学和量子力学当中，人们分别是用牛顿第 二定律和薛定谔方程来描述物体的运动的，这两者在物理系统当中具有相同的地 位。薛定谔方程式可以描述任何的微观系统，通过求解该微观系统所对应的薛定 谔方程，我们能够得到其波函数以及对应的能量，从而进一步来了解该微观系统 的性质。 薛定谔方程可以分为与时间有关和与时间无关两种类型，其中量子系统的波 函数随着时间的演化过程是通过与时间有关的薛定谔方程来描述的，而与时间无 关的薛定谔方程

11、则描述的是固定状态的量子系统的物理性质，方程的解即是该量 子系统固定状态的波函数。 本文考虑的是二维与时间有关的薛定谔方程： (1-1) 22 22 ( , , )( , , )( , , )( , ) ( , , ) uuu ix y tx y tx y tw x y u x y t txy ( , , ),0,x y ta ba bt 初始条件为： ( , ,0)( , )u x yx y 边界值条件为： (1-2) 12 34 ( , , )( , ),( , , )( , ),0 ( , , )( , ),( , , )( , ),0 u x a tx tu x b tx tt u a

12、 y ty tu b y ty tt 其中是任意的势函数，是波函数，且在定义域内连续。( , )w x y1i ( , , )u x y t 薛定谔方程是反应微观粒子随着时间变化的非相对论波动函数，它仅适用于 速度比较缓慢的非相对论粒子。其中波函数可以很好地描述微观粒子的( , , )u x y t 状态，在势函数中微观粒子运动的薛定谔方程即为方程(1-1)。我们可以通( , )w x y 过给定的初始条件和边界值条件以及波函数所满足的条件，来求解出波函数 ，进而计算粒子的分布概率。( , , )u x y t 薛定谔方程被广泛地应用于化学和物理等领域中，如量子器件的建模2，光 纤传播模型3，

13、光电子器件的设计4，电磁波的传播5，天体系统的量子化6，轴 近似波动方程的水下声学7，量子动力学计算的应用8,9，化学核外电子的运动状 态描述10等。它被应用到原子、核等诸多方面问题中，所得到的结果都与实际很 相符。近年来，很多学者通过各种方法研究具有复杂势函数的薛定谔方程所描述 的问题11-14，解释了很多重要的物理现象，因此对薛定谔方程的求解具有相当重 要的意义。 1.2 国内外研究现状 到目前为止，对薛定谔方程(1-1)的求解已经有了很多种数值方法， 大多都 是采用的有限差分法15-17，或者是用三角正交函数系或幂级数函数展开的谱方法 18,19。subasi 给出了具有二阶精度的有限差

14、分方法20，kalita 等人给出了一个隐 式的半离散高阶紧凑方法21，antonie 等人给出了一个 crank-nicolson 隐格式方 法22，dehghan 给出了不同的有限差分方法包括三个全隐式和两个全显示差分方 法以及交替方向隐式法和 barakat 和 clark 类型的显示方法23，dehghan 和 shokri 还给出了使用配置和薄板样条径向基函数的数值方法24，此外 dehghan 和 mohebbi 还给出了求解方程(1-1)的紧凑有限差分法25，gao 和 xie 还给出了紧凑 的交替方向隐式有限差分法26，该方法在空间上具有四阶精度，在时间上具有二 阶精度。li

15、等人还给出了多元二次（mq）和薄板样条（tps）径向基函(3)m 数的 mps 方法求解薛定谔方程，该方法类似于有限差分法27。dehghan 和 taleei 还提出了一种紧凑的分布有限差分方法来求解薛定谔方程28，该方法通过使用四 阶精度紧致差分格式，来提高分布有限差分方法的准确性，而且还具有无条件稳 定的性质。 谱方法的思想起源于傅立叶分析，它是一种既古老又新兴的求解偏微分方程 的方法。求解偏微分方程的三种最基本的方法分别是谱方法，有限差分方法和有 限元方法。谱方法和另两种方法相比，具有“无穷阶”收敛的特点，即它的收敛 速度会随着真解的光滑程度变高而变快，从而谱方法就能用限制自由度的方式

16、来 得到较高的精度29，另两种方法在这一点上是无法比拟的。 1.3 本文的主要研究内容 本文主要是用 galerkin-chebyshev 谱方法和边界值法求解二维 schrdinger 方 程，运用 galerkin-chebyshev 谱方法对空间导数进行近似,离散薛定谔方程(1-1)， 从而将原问题转化为复数域上的线性常微分方程组，然后再用边界值法求解该方 程组，所求得的数值解即为原问题的解，之后再进行误差分析，得到误差分析结 果，最后再通过 matlab 进行数值模拟，给出数值解的图像以及误差曲面图像。谱 方法求解偏微分方程具有高精度的性质，边界值法求解常微分方程同样具有高精 度和稳定

17、的特点，这样问题即得到解决。 在第一章中我们阐述了薛定谔方程在当前科学研究中的应用，表明求解薛定 谔方程具有很深远的意义，还介绍了现阶段求解该方程的主要方法，以及本文即 将采用的方法。紧接着在第二章中，我们介绍了本论文所需要的一些预备的基础 知识，为后面论文的顺利进行做好准备工作。 在第三章当中，我们采用 galerkin-chebyshev 谱方法求解椭圆型方程，以及 用边界值法求解常微分方程，并给出求解特殊常微分方程组的求解格式，这两个 方法求解微分方程都具有很高的精度和很好的稳定性。 第四章中，先对原问题进行区域映射处理，以及对边界条件进行齐次化处理 以后，然后运用 galerkin-c

18、hebyshev 谱方法对二维薛定谔方程进行离散，将其转 化成常微分方程组，然后对该微分方程组进行求解，得到数值解，接着对该方法 进行误差分析，得到误差估计结果。 第五章进行数值模拟，根据前面的内容，编程得到问题的数值解，并和相应 的精确解进行比较，分析其误差，画出误差曲面图像。 最后是本文的一个总结，以及研究此问题的意义和前景展望。 第 2 章 预备知识 2.1 克罗内克积的简介 定义定义 2.130：设是一个行列的矩阵，是一个行列的矩阵，amn m n (a ) ij a bpq ，克罗内克积可以表示成：() ijp q bb ab 11121 21222 12 n n mmmn a b

19、a ba b a b a ba b ab a b aba b 它是一个的分块矩阵。mpnq 克罗内克积具有如下的一些性质： 性质 1：满足结合律与双线性的性质： 如果矩阵存在，则 ；bc()abcabac 如果矩阵存在，则；ab()abcacbc ，其中是常数；()()()kabakbk abk .()()abcabc 性质 2：，和是四个矩阵，如果矩阵乘积和存在，那么就有abcdabcd ()()ab cdacbd 性质 3：是可逆的当且仅当和是可逆的，其逆矩阵是：abab 111 ()abab 性质 4：.()t tt abab 定义定义 2.2：设是一个行列的矩阵，那么把矩阵按列将后一列

20、amn m n (a ) ij a a 堆在前一列后面，形成的一个新的列的向量记为，即：mn( )vec a 1112121231 ( )(a ,a ,a ,a ,a,a ,a)t nnmn vec a 定理定理 2.1：设是一个行列的矩阵，是一个行列的矩阵，是一个ammbnnu 行列的矩阵，也是一个行列的矩阵，那么有：mnfmn ()( )( ) t aubfba vec uvec f 证明：证明：先将矩阵,写成如下的形式：buf 121212 , , t nnn bb bbuu uuffff 其中,分别是矩阵,第 列的列向量， i b i u i fbufi1,2,in 则有： ()( )

21、( )ba vec uvec f 1112111 2122222 12 n n nnnnnn abababuf abababuf ufababab , 1 ,1,2, n ijji j ab uf in 12 , n aub aubaubf t aubf 从而原题得证。 2.2 chebyshev 多项式介绍及其性质 定义定义 2.3：在区间上的权函数以递归的形式定义的正交多项式 1,1 2 1 ( ) 1 x x 称为 chebyshev 多项式，它可写成：。( )cos( arccos ),0,1, n t xnx n chebyshev 多项式具有如下的性质： 性质(1) 31：正交性

22、1 1 2 2 1 0, ( ),( )( )( )(1),0 2 ,0 nmnm mn t x txt x txxdxmn mn 性质(2) 31：递推关系 11 01 ( )2( )( ) ( )1,( ) nnn txxt xtx t xt xx 11 01021 1 ( )2(1)( )( ) 1 ( )0,( )( ),( )4 ( ) nnn n txnt xtx n t xt xt x t xt x 性质(3) 31：是阶多项式，是阶多项式，是阶多项式，( ) n t xn ( ) n t x1n ( ) n tx2n 满足：，其中 1 0 1 ( )2( ) n nk k k

23、k nodd t xnt x c 1 22 0 1 ( )()( ) n nk k k k neven txn nk t x c 。 2,0 1,0 k k c k 定理定理 2.232：设，则: 2 ( )( )( ),( ),( ) ,( ),( ) kkkjkkjjkkj xt xtx axxbxx 2 (1)(2), 4 (1),2,4,6, 0,. jk jjkj ajkjjj kj or kjodd 1 , 2 ,2,2 2 0, k jk c kj bkjor kj otherwise 2.3 chebyshev 正交逼近的性质 我们讨论 chebyshev 逼近问题，需要借助带

24、权的 sobolev 空间，下面记以 为权的阶空间为，它 11 22 22 (1)(1)xy msobolev( 1,1) ( 1,1)( ) mm hh 的内积和范数定义分别为 (k)(k) , 0 ( , ), m m k u vu vdxdy 1 2 , , ( , )m m uu u 记。 2 11 22 0, ( , )( , ) l uu uu u 设区间是一个非空集，且是 lebesgue 可测的，记的范数为： n r( ) p l 1 () () , 1 p p l ffdxdyp 当时，。p () sup l x i fessf 接下来定义空间，设空间是有界的，且，有 p l

25、 n r1p () ( ): p p l lff 在空间和上的全体次连续可微的函数所构成的集合分别记为和k( ) k c 。记，( ) k c 12 (,) n n n 12n 12 12 n n dd dd 其中是广义导算子，接下来定义弱导数。 i i d x 定义定义 2.433：设，满足上面的式子，称 是的阶弱导数，记为 1 ,( ) loc u vlvu ，如果vd u 0 ( 1),( )uddxv dxc 有时又称在弱的意义下。vd u 下面定义空间，设区域是有界的，是非负整数， ,k p w n rk1p ，有：( ) p pl , ( ):( ),( ), k ppp wffl

26、d flk , () k p p w l k fd f 其中为空间上面的范数。 ,k p w f , ( ) k p w 在空间上的闭包记为，当时， 0 ( )c , ( ) k p w , 0, k p w 2p ,2 ( ) kk wh 。 ,2 0,0, ( ) kk wh 定理定理 2.3（gronwall 不等式）不等式）34：设和是上的非负的连续函数，并( ) t( )g t0, t 且在是可微的，如果存在常数满足，使得对任意的，都( ) t0, t0(0, )tt 有： ( ) ( )( )ttg t 或者等价的还有： 0 ( )(0)( )( ) t tgd 那么就有： ()

27、0 ( )(0)( ),0, t tt tegedtt 2.4 投影算子的性质 记是一个多项式空间，其最大自由度是，是到的正( )npn n p 2 ( )l( )np 交投影算子，是到的椭圆投 11 0, ( )( )()0hhv v 0 n p 1 0, ( )h ( )np 影算子，则有如下的定义和性质定理： 定义定义 2.535：空间 中从到的正交投影算子为： 2 ( )l 1 0, ( )h 2 n v n p 2 (, )0, nn p vvv 定义定义 2.635: 空间 中从到的椭圆投影算子为： 2 ( )l 1 0, ( )h ( )np 0 n p 02 ( (),()0,

28、 nwn p uuvvv 定理定理 2.436：对任意的非负整数 ，都有下面的不等式：s( ) s h 2 ()() s s n lh pcn 定理定理 2.537：对任意的非负整数 ，都有下面的不等式：s 1 0, ( )( ) s hh 12 00 () ()() s s nn h hl pnpcn 2.5 本章小结 本章给出了完成这篇论文所需要的一些必备的基础知识，首先介绍了克罗内 克积的定义以及性质，然后介绍了切比雪夫多项式性质以及一些重要的关系定理， 之后介绍了 chebyshev 正交逼近的性质，其中包括内积和范数的定义，空间的定 义和性质，最后还介绍了投影算子的定义和不等式性质，

29、为论文的进行做好准备 工作。 第第 3 章章 galerkin-chebyshev 谱谱方方法法和和边边界界值值法法 3.1 用 galerkin-chebyshev 谱方法求解椭圆型方程 考虑用 galerkin-chebyshev 谱方法来求解如下的椭圆型方程 (3-1) 22 22 ()( , )( , ) uu w x y uuf x y xy ( , ) 1,1 1,1x yi 边界条件是： ( 1, )(1, )0,11uyuyy ( , 1)( ,1)0,11u xu xx 由 chebyshev 多项式的定义和性质，设， 01 ( ),( ),( ) nn sspan t x

30、t xtx ，则方程(3-1)的, ( 1)0 nn vusu 2 ( )( ): ,0,1,2 nij vspanxyi jn galerkin-chebyshev 谱方法是求使得对任给的都满足 2 nn uv 2 n vv 22 22 ( (), )( ( , ), )(, )( ( , ), ) nn nn uu vw x y uvuvf x y v xy 其中，。( , )u vuv dxdy 11 22 22 (1)(1)xy 1,1 1,1 令,取 2 ,0 ( , )( )( ) n nkjkj k j ux yuxy ( )( ), ,0,1,2, lm vxy l mn 则有

31、 (, )(, )( ( , )( , ), ) nnn uvuvw x y uf x y v (,( )( )(,( )( )( ( , )( , ),( )( ) nlmnlmnlm uxyuxyw x y uf x yxy 将上式用矩阵表示即可写成 ()buaaubbubf 其中：，和满足定理 2.2 中的条件关系， ,0,1,2 () kjk jn aa ,0,1,2 () kjk jn bb kj a kj b ，且。 ,0,1,2 () kjk jn uu ,0,1,2 () kjk jn ff ( ( , )( , ),( )( ) kjnkj fw x y uf x yxy 由

32、定理 2.1 有 (3-2)()( )()( )( )abba vec ubb vec uvec f 对方程(3-2)进行求解，就可以求出其数值解，从而得到方程(3-1)的数值解 kj u 。( , ) n ux y 3.2 用边界值法求解常微分方程 边界值法是最近求解常微分方程数值解的常用方法，简称为 bvms，它是线 性多步法的一个推广，和其他常微分初值问题的数值解法相比较，bvms 具有高 精度和无条件稳定的特点38-41，是一个很好的方法。 考虑下面的初值问题 (3-3) . 0 ( )( , ( ),(0),0y tf t y tyyt 用步线性多步法离散上面的方程即可得到k (3-

33、4) 00 ,0,1,., kk rrjrrj rr ytfjnk 其中，为系数。(),( ,) rrrrr yy r t tr t ff ty r r 由泰勒展开有： 23 ()() 2!3! rjjjjj rtrt yyrt yyy 23 (4) ()() 2!3! rjrjjjjj rtrt fyyrt yyy 从而令： (1)( ) 01 00 kk qq rrjrrjjjqj rr lytfc yctyct y 则有： (3-5) 0 0 1 10 1 10 11 ,2,3, !(1)! k r r kk rr rr kk qq qrr rr c cr crrq qq 如果有次的连续

34、微商，那么就可以选取和使得( )y t2pk, jj ，即选取使其满足 01 0, p ccc 1 0 p c , rr (3-6) 0 11 1 11 0 0 11 0 !(1)! k r r kk rr rr kk pp rr rr r rr pp 此时就有 1(1)2 1 ( )() ppp p lctytot , 00 kk rrjrrjj k rr ytfr 其中为截断误差，略去，就得到了线性多步法(3-4)，该方法的精度是 , j k r , j k r 阶的。p 求解方程(3-4)需要个初始边界条件和个结尾边界条件，即我们需要k 和，初始边界条件可以由方程(3-3)得到。个初始

35、011 ,.,yyy 1,.,n kn yy 0 y1 边界条件和个结尾边界条件则来自于以下等式k (3-7) ( )( ) 00 ,1,.,1, kk jj rrrr rr ytfj 和 (3-8) ( )( ) 00 ,1,., kk jj rn k rrn k r rr ytfjnkn 其中系数和的选择，要满足使基于最初与最后的边界条件的方法的截断 ( ) j r ( ) j r 误差与基于公式(3-6)的方法具有相同的阶。 方程(3-4)(3-8)用矩阵形式表示可以写为 ( ,) eeeeee a ytb f ty 其中. 1(1) , nnn eeee tyra br (1)(1)(

36、1) 01 (1)(1)(1) 01 01 0 0 (1)(1) 0 ()() 0 (1) k k k e k n kn k k nn k nn a 用代替矩阵中的，即为矩阵，并且. e a e b 0101 (,.,) ,(,.,) tt enen yyyyffff 对进行划分, eee a bf 0000 , , , , tt eeee aaa yyybb bfff (1)(1)(1)(1) 00000000 ,0,0,0,0 tt ab 将第一列分离出来，可以得到的等价式( ,) eeeeee a ytb f ty (3-9) 0 ( , )aytbf t yg 其中是一个未知量，且有

37、n yr 11000000 ,( ,) tt nn yyyfffga ytb f ty 在这里我们用四阶 bvms 近似方程(3-3)，并取，由(3-6)得到3,2k 0123 1230123 123123 0 230 1227149 0 2438632 求解该方程组，得到其基础解系，取其中的三组解，一组代入到方程(3-4)，另两 组分别代入到(3-7)和(3-8)，即可以得到下面的关系式 (3-10) 32121 1 (99)() 122 jjjjjj t yyyyff 其中。 01230123 199111 ,0,0 1212121222 其对应的初始边界条件和结尾边界条件分别为 (3-1

38、1) 321010 1 (9917)(3) 244 t yyyyff 其中。 01230123 1799113 ,0,0 2424242444 (3-12) 3211 1 (9917)(3) 244 nnnnnn t yyyyff 其中。 01230123 1991731 ,0,0, 2424242444 把上面的三个式子化为等式(3-9)的形式，则其中的分别为 00 , ,a b a b 9/ 249/ 241/ 24 9/129/121/12 1/129/129/121/12 1/129/129/121/12 1/ 249/ 249/ 2417/ 24 a 3/ 4 1/ 21/ 2 1/

39、 21/ 2 3/ 41/ 4 b , 0 171 ,0,0 2412 t a 0 1 ,0,0 4 t b 如果我们考虑的是特殊的线性常微分方程组 (3-13) . 0 ( )( )( ), (0),0 x y tb y tg tyyt 其中是的矩阵,且 x bm m 1212 ( )( ),( ),.,( ) , ( )( ),( ),.,( ) tt mm y ty ty tytg tg tg tgt 那么我们可以将(3-13)化为矩阵形式如下 (3-14) 00000 ()()() mxmx aitbbyt bigt bb ygay 其中是的单位矩阵，且 m im m 00 ( )gg

40、 t 112111222212 ( ),( ),.,( ),( ),( ),.,( ),.,(),(),.,()t mmnnmn yy ty tyty ty tyty ty tyt 112111222212 ( ),( ),.,( ),( ),( ),.,( ),.,(),(),.,()t mmnnmn gg tg tgtg tg tgtg tg tgt 如果线性常微分方程组可以化为下面的形式 (3-15) . 0 ( )( )( ), (0),0 xx a y tb y tg tyy t 其中是的非奇异矩阵, 那么用四阶 bvm 法可以将方程(3-15)化为 x am m (3-16) 00

41、000 ()()() xxmxx aatbbyt bigt bb ygaa y 3.3 本章小结 本章给出了利用 galerkin-chebyshev 谱方法求解椭圆型方程的数值方法格式， 利用该方法将偏微分方程进行离散以后得到常微分方程组，再利用常微分方程的 解法对其进行求解即可达到将偏微分方程进行求解的目的，谱方法精度很高，稳 定性也很好，对于求解偏微分方程是一个非常好的方法。还给出了利用边界值法 求解微分方程组的过程，并给出了几种特殊形式的微分方程组的边界值法求解的 数值格式，边界值法求解微分方程具有很高的精度，对于求解常微分方程组也是 一个很好的方法。 第 4 章 求解二维薛定谔方程

42、4.1 区域和边界条件的处理 由于 galerkin-chebyshev 谱方法只能解决齐次边值条件的问题，故针对本文 的二维 schrdinger 方程问题，需要先进行区域的映射处理，对非齐次的边值问题 进行齐次化处理，将其转化成方程的一般形式，再进行求解。 4.1.1 区域的处理 原问题中，我们在这里对其进行一定的变换处理，使区域 ( , ),x ya ba b 变成。 1,11,1 令 ， 1 2 x xaba 1 2 y yaba 11 , , , 22 xy u x y tu abaabatu x y t , , ,u x y tu x y t tt 11 , , , 22 11 ,

43、 , , 22 22 xy u abaabat u x y tx xxx xy u abaabat babau x y t xx 2 2 22 2 222 1 , , , , ,() 2 44 x u abay t u x y tbau x y tba xxx 11 , , , 22 11 , , , 22 22 xy u abaabat u x y ty yyy xy u abaabat babau x y t yy 2 2 22 2 222 1 , , , , ,() 2 44 x u abay t u x y tbau x y tba yyy 于是方程(1-1)就转化为 (4-1) 22

44、 2222 , , ,44 (, , ) 11 (),() (, , ) 22 u x y tu x y tu ix y t txy baba xy w aba aba u x y t (, , )1,11,10,x y tt 初始条件为： 11 (, ,0)(),() 22 xy u x yaba aba 边界值条件为： 12 34 11 (, 1, )(), ),(,1, )(), ),0 22 11 ( 1, , )(), ),(1, , )(), ),0 22 xx u xtaba tu xtaba tt yy uy taba tuy taba tt 对方程(4-1)进行简化，可以将其

45、表示为 (4-2) 22 22 (, , ), , ,(, , ) uuu ix y tx y tx y tx y u x y t txy (, , )1,11,10,x y tt 初始条件为： (, ,0),u x yx y 边界值条件为： 12 34 (, 1, )(, ),(,1, )(, ),0 ( 1, , )( , ),(1, , )( , ),0 u xtx tu xtx tt uy ty tuy ty tt 其中 2 4 ba 11 ,(),() 22 xy x yw aba aba 11 ,(,() 22 xy x yaaba 1122 11 (, )(), ),(, )()

46、, ) 22 xx x taba tx taba t 3344 11 ( , )(), ),( , )(), ) 22 yy y taba ty taba t 4.1.2 边界条件的处理 在本文中对边界条件的处理过程，就是一个对边界条件进行齐次化的过程。 由方程(4-2)，我们很容易得到 ， 22 11 22 1 xx 22 22 22 1 xx 22 33 22 1 yy 22 44 22 1 yy ， 11 tt 32 tt 33 tt 44 tt 于是令，则有 1121 1 , ,(, )(, )(, ) 2 y ux y tx tx tx t a 111 (, )(, ), 1,0x

47、tx tuxt a 221 (, )(, ),1,0x tx tuxt a 331 ( , )( , )1, ,y ty tuy t a 441 ( , )( , )1, ,y ty tuy t 2222 1121 2222 (, , )1 2 ux y ty xxxx 1121 (, , )1 2 ux y ty tttt a 331 1, , u y t ttt a 441 1, , u y t ttt a 2222 3331 2222 1, , u y t yyyy a 2222 4414 2222 ! 1, , ! un y t yyyyrnr 令，即可得到 aaa 3432 1 ,

48、,( , )( , )( , ) 2 x ux y ty ty ty t aaa 2222 3432 2222 (, , )1 2 ux y tx yyyy aaa 3432( , , )1 2 ux y tx tttt aaa aa 3432 34 3141 3141 1 (, 1, )( 1, )( 1, )( 1, ) 2 11 ( 1, )( 1, ) 22 11 ( 1, )( 1, 1, )( 1, )(1, 1, ) 22 11 ( 1, )( 1, )( 1, )(1, )0 22 x uxtttt xx tt xx tuttut xx tttt 同理可得 2( ,1, )0u

49、xt 从而有 a 112 ( , )(, ), 1,0x tx tuxt a 222 (, )(, ),1,0x tx tuxt a 332 ( , )( , )1, ,0y ty tuy t a 442 ( , )( , )1, ,0y ty tuy t 再令，则有 12 uuuu a 12 41 141 4 (1, , ) (1, , )(1, , )(1, , ) (1, , )(1, , )( , ) (1, , )(1, , )( , )(1, , ) (1, , )( , )0 uy tuy tuy tuy t uy tuy ty t uy tuy ty tuy t uy ty t

50、 同理有 ， ( 1, , ) 0uy t ( ,1, ) 0u xt ( , 1, ) 0u xt 即 1, , 1,0uy tu xt 由可以得到 12 uuuu 12 uuuu tttt 2222 12 2222 uuuu xxxx 2222 12 2222 uuuu yyyy 于是方程(4-2)就转化为下面的方程 2222 1212 12 2222 , uuuuuuu ix yuuu tttxxyy 2222 12 2222 12 12 , , uuuuu ix y u txyxy uu x yuui tt 即可写成 (4-3) 22 22 , , uuu ix y uf x y t

51、txy (, , )1,11,10,x y tt 初始条件为： 0 (, ,0),u x yux y 边界值条件为： ( , 1, )0,0 ( 1, , ) 0,0 u xtt uy tt 其中 012 , ,0, ,0ux yx yux yux y 22 1212 12 22 , , uuuu f x y tx yuui xytt 于是原问题就转化成了标准问题。 接下来我们需要去复数域，令： , uriq 12 ffif 带入方程(4-3)得到 2222 12 2222 , rqrqrq iiiix yriqfif ttxxyy 根据复数的性质，即可以得到 (4-4) 22 2 22 22

52、 1 22 , , rqq x y qf txy qrr x y rf txy (, , )1,11,10,x y tt 初始条件为： 0000 (, )re(, ) ,(, )im(, )r x yux yq x yux y 边界值条件为： (, 1, )(, 1, )0,0 ( 1, , )( 1, , )0,0 r xtq xtt ry tqy tt 4.2 二维薛定谔方程的求解 对上面的方程组(4-4)进行移项变换，得到 (4-5) 22 2 22 22 1 22 , , qqr x y qf xyt rrq x y rf xyt 令 ， 2 ,0 n nkjkj k j rrtxy

53、2 ,0 n nkjkj k j qqtxy ， a 1 2 ,kj nkj fx y qfxy a 2 1 ,kj nkj fx y rfxy 0 0 , kjkj i rrx yxyd 0 0 , kjkj i qqx yxyd a a 00 0,1,20,1,2 0,1,2 , kj kjkj kjnkjn kjn rrqqrr aa aa 12 00 0,1,2 0,1,20,1,2 , kj kjkj kjn kjnkjn ffffqq 令，由第三章的方法即可以得到解的弱形式 lm vxy (4-6) 22 2 22 22 1 22 , , nnn n ii nnn n ii qqr

54、 x y qv dfv d xyt rrq x y rv dfv d xyt 00 , ,0,0, ,0,0 nn ii rx yrx yv dqx yqx yv d 将(4-6)用矩阵表示即可写成 (4-7) aaaa aaaa () () bqaaqbbr bf braarbbq bf ， a 0 (0)brbr a 0 (0)bqbq 其中 ， a a dr r dt a a dq q dt 0,1,2 kj kjn bb 利用定理 2.1，方程组(4-7)等价于 (4-8) aaa aaa ()()()( )() ()()()( )() bb vec rabba vec qvec f

55、bb vec qabba vec rvec f a 0 ()( (0)()bb vec rvec r a 0 ()( (0)()bb vec qvec q 然后利用第三章的边界值法对上面的常微分方程组进行求解，所得的数值解 即为原问题的数值解。 4.3 误差分析 由定理 2.4 可以得到 ， 2 h 000 s s n l rrcnr 2 h 000 s s n l qqcnq 那么问题(4-4)的弱形式是求使得 221 0, ,0, ;0, ;h,r qlt llt (4-9) 2 1 (), (), iii iii r qvdx y qv dfv d t q rvdx y rv dfv d

56、 t galerkin-chebyshev 谱方法是求,使得 2 ,v nnn rq 2 v,0, n tt (4-10) 2 1 0 0 (), (),. n nn iii n nn iii nn nn r qdx y qdfd t q rdx y rdfd t rp r qp q 令 ， 00 pp nnnn r rrrrr 00 pp nnnn q qqqqq ,，, 0 p nnn rr 0 pnrr 0 p n nn qq 0 pnqq 则有 ， nn r r n n q q 定理定理 4.14.1：假设，满足上面的关系，那么对任意的都有r n rq n q0, ,tt 22 nn

57、s ll rrqqcnm 其中 hh 1 22 2 ()()()() 0 2(0)(0) 2( )( ) ss s ss s t tt hhhh rqr tq tdm r tq t 证明：证明：由定理 2.5 可以得到 2 2 0 h ( )p s s n l l tr tr tcnr t 由定理 2.4 可以得到 22000 h 0,p s s nn ll rrrx yrcnr 同理可得 2 2 22 0 h 000 h p 0,p s s s n ll s nn ll tq tq tcnq t qqqx yqcnq 接下来估计和，令 2 n l 2 n l , ,(), ii a t g

58、vgvdx y gv d 由方程组(4-9)，(4-10)可以得到 2 (), (), iiii n nn iii r qvdx y qv dv df v d t r qvdx y q v dv d t 从而有 ()(),()0 n nn iii rr qqvdx yqqv dv d t 2 ( , , )( , )(, )(, )0,( ) n n n a tva tvvvvv tt 同理得到 2 ( , , )( , )(, )(, )0,( ) n nn a tva tvvvvv tt 由投影算子的性质，我们可以得到 0 0 0 0 0 0 2 0 0 ( , )(, )( , )(,

59、)(, ) ( , )(),(, ) ( , )(, ) ( , )(, )() n n n n n n n n nn n n n nn iii i p r r a tvva t p qqvvv ttt p r r a t p q vqvdx y q v dvv d tt p r a t p q vvf v d t p r a t p q vvqvd t 0 00 0 00 0 , (),(, )() , (, )(, )()(),() (, )(, )( n nn n nn n iii iii ii ii r x y qv dv d t p r p qvdx y p qv dvqvd t r

60、 x y qv dv d t p r r vvp qqvdx yp qq v d tt p r r vv tt , )v t 即 0 ( , )(, )(, )0,( ) n n n a tvvvvv tt 同理有 0 ( , )(, )(, )0,( ) n nn a tvvvvv tt 分别取，分别带入到上面两个式子中，就有 n v nv ( ,)(,)(,)0 n n nnn a t tt ( ,)(,)(,)0 n nnn n a t tt 由于，且 2 2 (,) nn l n 2 2 (,)(,)(,) 2(,)2(,) nnnnnnn nnnn l tttt tt 故得到 2 2