方程的时域差分数值解法-安徽师范大学学报.PDF

第 卷 期 安 徽 师 范 大 学 学 报 （自然科学版） ４０ ２ 犞狅犾．４０犖狅．２ 年 月 （ ） ２０１７ ３ 犑狅狌狉狀犪犾狅犳犃狀犺狌犻犖狅狉犿犪犾犝狀犻狏犲狉狊犻狋 犖犪狋狌狉犪犾犛犮犻犲狀犮犲 犕犪狉．２０１７ 狔 ： ／ 犇犗犐１０．１４１８２犑．犮狀犽犻．１００１－２４４３．２０１７．０２．００６ 犛犮犺狉犱犻狀犲狉方程的时域差分数值解法 犵 王广保 （滁州职业技术学院 传媒工程系，安徽 滁州 ２３９０００） 摘 要：将时域差分方法应用于求解量子体系的 方程 详细推导了含 方程 犛犮犺狉犱犻狀犲狉 ． 犛犮犺狉犱犻狀犲狉 犵 犵 的离散化公式，得到了耦合微分方程的差分格式，编写了相应的数值计算程序 以谐振子势为例，对 ． 程序代码作了检验，结果令人满意． 关键词： 方程；时域差分；数值解； 变换 犛犮犺狉犱犻狀犲狉 犉狅狌狉犻犲狉 犵 中图分类号： 文献标志码： 文章编号： （ ） 犗５６ 犃 １００１－２４４３２０１７０２－０１３５－０４ 引 言 在束缚态能级计算中，需要求解 方程 由于具有精确解的势函数只有无限深阱势、谐振子 犛犮犺狉犱犻狀犲狉 ． 犵 势、中心力场势、 势、 势等少数几种情形，因此实际应用中近似的数值解法就十分重要 通 犕狅狉狊犲 犘狅狊犮犺犾犜犲犾犾犲狉 ． ［，］ 常的算法是 犆狅狅犾犾犲等提出的有限差分方法 １２ ，即先把波函数作 犜犪犾狅狉展开，取有限项（多为二次项）近似， 狔 再采用 两步格式，通过前向、后向迭代波函数在匹配点的光滑连接，而得到能量本征值 在量子数 犖狌犿犲狉狅狏 ． 较低时，该方法数值稳定性和计算结果都相当好，但对高量子数却不能令人满意 代表性的改进方案，主要有 ． ［，］ ［，］ ３４ ５６ 分立变量表象方法（ ） 和时间演化方法（ ） 等 等人则将有限差分方法和时间演化方法 犇犞犚 犜犘犉犜 ．犢犲犲 ［］