有限差分法模拟一维(二维)谐振子.doc

目目 录录 第第 1 1 章章 概述概述 4 第第 2 2 章章 有限差分方法有限差分方法 5 2 1 有限差分法基本思想 5 2 2 差分方程组的求解 5 2 2 1 高斯 赛德尔迭代法 5 2 2 2 逐次超松弛法 6 第第 3 3 章章 求解谐振子的微分方程求解谐振子的微分方程 7 3 1 一维谐振子 7 3 2 二维各向同性谐振子 9 第第 4 4 章章 总结总结 12 参考文献参考文献 13 附录附录 14 附 1 一维线性谐振子的程序设计 14 附 1 1 基态一维线性谐振子 14 附 1 2 第一激发态一维线性谐振子 15 附 1 3 第二激发态一维线性谐振子 15 附 2 二维线性谐振子的程序设计 16 2 第第 1 1 章章 概述概述 微分方程和积分微分方程数值解的方法 基本思想是把连续的定解区域 用有限个离散点构成的网格来代替 这些离散点称作网格的节点 把连续定 解区域上的连续变量的函数用在网格上定义的离散变量函数来近似 把原方 程和定解条件中的微商用差商来近似 积分用积分和来近似 于是原微分 方程和定解条件就近似地代之以代数方程组 即有限差分方程组 解此方 程组就可以得到原问题在离散点上的近似解 然后再利用插值方法便可以从 离散解得到定解问题在整个区域上的近似解 有限差分法可广泛用来求解偏微分方程的近似解 在电磁场中求解点位函 数的拉普拉斯方程时 可采用有限差分法的基本思想是 用网格将场域进行分 割 再把拉普拉斯方程用以各网格点处的点位作为未知数的差分方程式来进行 代换 将求解拉普拉斯方程解得问题变为求联立差分方程组的解得问题 在 1 差分网格非常多和情况下 利用并行计算方法对其进行区域分解 每个进程负 责运算一部分区域 区域边界之间进行必要地通信可有效提高计算速度 解决 更大规模的问题 往往只讨论它在静态场中的应用 即泊松方程或拉普拉斯方 程的有限差分形式 很少涉及到它在时谐场 即亥姆霍兹方程 中的应用 本 文重点讨论亥姆霍兹方程的有限差分形式以及它在时谐场中的应用 同时 有 限差分法 finite difference method 是基于差分原理的一种数值计算方法 在求解微分方程定解问题中广泛应用 有限差分法是以差分原理为基础的一种 数值计算法 它用离散的函数值构成的差商来近似逼近相应的偏导数 而所谓 的差商则是基于差分的应用的数值微分表达式 用离散的只含有有限个未知量 的差分方程组去近似代替连续变量的微分方程和定解条件 并把差分方程组的 姐作为威风方程定解问题的近似解 有限差分法可以处理几乎所有形式的势函数 且主程序不依赖于势函数的具体形式 对于多数两字体都可以进行相对准确的 计算 因此 将有限差分法应用于量子力学本征值问题的计算 有助于相对准 确地进行量子体系和形象直观地教学研究 32 量子力学教程中队一维无限深势阱 线性谐振子 氢原子等量子体系的薛 定谔方程进行了严格的求解 得到了描述体系状态的波函数和能量的精确解 多数量子体系的哈密顿算数比较复杂 薛定谔方程不能严格求解 因此 研究 和发展薛定谔方程的数值计算方法具有重要意义 4 3 第第 2 2 章章 有限差分方法有限差分方法 2 12 1 有限差分法基本思想有限差分法基本思想 有限差分法是解偏微分方程的主要数值方法之一 其基本思想是把连续的 问题离散化 即首先对求解区域作网格剖分 用有限个网格节点代替连续区域 其次将微分算子离散化 从而把微分方程的定解问题化为代数方程组的求解问 题 解方程组就可以得到原问题在离散点上的近似解 然后再利用插值方法便 可以从离散解得到定解问题在整个区域上的近似解 参照文献 给出有限差 5 分法数值计算的基本思想 1 区域的离散或子区域的划分 2 插值函数的选择 3 方程组的建立 4 方程组的求解 2 22 2 差分方程组的求解差分方程组的求解 利用有限差分要面临求解的问题 在实际应用中 可以采用逐次迭代的迭 代方法求解 这里介绍两种常用的迭代方法 高斯赛 德尔迭代法和逐次超松弛 迭代法 2 2 12 2 1 高斯高斯 赛德尔迭代法赛德尔迭代法 如图 1 1 采用正方形网格分割场域 这个方法是 先对节点选取 ji yx 初值 其中 0 表示 0 次近似值 下角标表示节点所在位置 即第 行 0 ij ji i 列的交点 再按下式j 其中 2 1 ij k ji k ji k ji k ji k ji h2 1 1 1 1 1 1 1 4 1 2 1 ji 反复迭代 一直进行到对所有节点满足下列条件为止 1 0 k 2 2 W k ji k ji 1 式中 是预定的最大允许误差 W 在高斯 赛德尔迭代中 网格节点一般按 自然顺序 排列 即先 从左到 右 再 从上到下 排列 如图 1 1 所示 4 y x 图 1 1 网格节点排列 1234 5678 2 2 22 2 2 逐次超松弛法逐次超松弛法 将式 2 1 在同一点上相邻两次迭代的差值计为 则可得 ji n R 2 3 2 1 1 1 1 1 1 1 4 1 k ji ij k ji k ji k ji k ji k ji k ji n hjiR 按式 2 1 迭代理想的收敛情况是所有内点的点位函数的余数为零 但这是不可 实现的 余数时正时负 时大时小 当网格很大时 收敛的速度很慢 ji n R 这就需要改小网格 但迭代的次数将随之增加 逐次超松弛法就是在高斯 赛德 尔迭代法中引入了加速收敛因子对其进行校正 即 2 4 2 1 1 1 1 1 1 1 4 4 k ji jik ji k ji k ji k ji n ji k ji k ji h 其中称为 加速收敛因子 是一个供选择的参数 其值在之间 该 21 方法的快慢与有着明显的关系 实践表明 如果选得好 可以较快的加速 迭代的速度 经验表明正方形场域由正方形网格划分 每边节点数为时 1 p 最佳的收敛因子为 0 2 5 p sin1 2 0 5 第第 3 3 章章 求解谐振子的微分方程求解谐振子的微分方程 3 13 1 一维谐振子一维谐振子 对一维线性谐振子 其能量本证方程为 3 1 Ex x h 22 2 22 2 1 d d 2 边界条件 3 2 x 0 x 去长度单位为 能量单位为 引入无量纲参量 2 ax 则式 3 1 化为无量纲形式 a 2 1 E 3 3 2 d d 2 令 1 1 式 2 3 可写为 i i i 2 m i1 mi 3 4 iiii 1 2 2 2 1 2 121 式 3 2 可写为 3 5 0 0 0 m 考虑式 3 5 式 3 4 可写为 其中 ii S 1 2 2 1 m m i 6 S 2 1 22 2 2 2 2 2 2 2 2 22 2 2 1 2 21 00 12 00 0 1 00 00 1 0 00 21 00 12 m m S 为三对角矩阵 是大型稀疏矩阵 用 Matlab 编程计算可同时得到矩阵 S 的本征值和本征矢 表 1 是取不同格点 步长值时一维谐振子能量 E 的数值计 算结果与精确解的比较 从表 1 可以看出有限差分法 m N 0N 1N 2 500 20 99752 98744 9673 1000 10 99942 99694 9919 2000 050 99982 99924 9980 计算结果比较精确 且计算精度随求解区域内格点数目增加而增高 主量 子数为 0 1 2 时 能量精确解分别为 简并度为 1 2 1 2 3 2 5 f 表 1 中 N 为量子数 取 m 200 0 05 时波函数图形如图 3 1 所示 表 1 中 N 为量子数 取 m 200 0 05 时波函数图形如图 3 1 所示 Matlab 程序见附 1 1 附 1 2 附 1 3 表 1 一维谐振子能量 以单 2 位 图 3 1 一维谐振子波函数 7 3 23 2 二维各向同性谐振子二维各向同性谐振子 对于二维各向同性谐振子 薛定谔方程为 3 6 Eyx yx 2 1 2 222 2 2 2 22 边界条件 3 7 0 yxx 0 yxy 取长度单位 能量单位为 引入无量纲参量 2 ayax 则式 3 6 可写为 a 2 1 E 3 8 22 2 2 2 2 令式 3 8 可写 11 11 2 2 njmi n j m i jijiji 为 jijijijijiji 1 1 2 22 2 1 2 1 1 1 2 1 3 9 式 3 7 可写为 3 10 0 0 0 0 0 0 nijmij 则在是 2 10 下 式 2 11 可写为 3 11 jjjj DCD 11 jm jm j j j D 1 2 2 1 2 2 1 0 0 1 8 22 1 222 2 22 1 222 2 22 1 222 2 22 1 22 11 2 1 000 111 2 1 00 00 111 2 1 00 111 2 j j j j C 则式 3 11 可进一步写成矩阵式 其中 jj S CD DC CD DC S 00 00 00 00 为三对角块矩阵 是大型稀疏矩阵 用Matlab编程计算可同时得到稀疏矩阵 的本征值和本征矢 表2给出了取不同格点 步长值时二维谐振子能量E的数 值计算结果与精确解的比较 从表2可以得到与表1相似的结论 主量子数为 0 1 2时 能量精确解分别为简并度分别为l 2 3 32 mn N 0N 1N 2 20200 50 51 96953 90465 7770 20200 50 5 3 91765 8054 20200 50 5 5 8810 40400 250 251 99593 98005 9525 40400 250 25 4 01455 9870 40400 250 25 6 1168 80800 1250 1252 00404 00005 9933 80800 1250 125 4 05076 0437 80800 1250 125 6 2166 表 2 二维谐振子的能量 以单位 2 取 m n 20 0 5 二维谐振子的波函数图形如图 3 3 所示 其中横 向的 2 个坐标分别为无量纲参量 纵向坐标为波函数 Matlab 程序 见附 2 9 图 3 3 二维谐振子波函数 10 第第 4 4 章章总结总结 计算物理学中 有限差分法 简称差分法 它以概念清晰 方法简单 直 观的特点应用于数值分析领域 因而应用广泛 无论是常微分方程还是偏微分 方程 各种类型的二阶线性方程 以至高阶或非线性方程 均可利用差分法转 换为代数方程组 然后利用计算机求其数值解 有限差分法是以差分原理为基础的一种数值解法 通过学习计算物理学的 有限差分法 我们对很多问题都可以得到足够高的计算精度 有限差分法是我 们较容易掌握的数值解法 我们在物理学中 可以利用有限差分法解任何偏微 分方程 利用有限差分法解边值问题时 首先将求解区域分为很多个网格和节 点 并用差商代替微商 然后 使区域中的偏微分方程转化为以节点的数值为 未知量的差分方程组 最后 解该方程组便可得到各离散点待求的数值解 该 数值解是近似解 但逼近区域的真实解 在我们的实验设计中 如果离散化的 点选择的足够密的话 我们得到的数值解与真实解的误差就能减小到可接受的 程度 11 参考文献参考文献 1 吴连坳 井孝功 丁慧明等 径向薛定谔方程的有限差分解法 J 吉林大学 自然科学报 1994 03 67 70 2 张志涌 精通 MATLAB6 5 北京 北京航空航天大学出版社 2003 3 刘建军 翟利学 有限差分法解能量本征方程 北京工业大学学报 2008 34 325 328 4 曾谨言 量子力学 M 北京 科学出版社 2000 12 附录附录 附附 1 1 一维线性谐振子的程序设计一维线性谐振子的程序设计 附附 1 11 1 基态一维线性谐振子基态一维线性谐振子 clc clear M 200 H 20 h H 2 M R 1 h 2 for i 1 2 M 1 for j 1 2 M 1 A i j 0 end end for p 1 2 M 1 A p p 2 R p M h 2 end for i 1 2 M 2 A i i 1 R end for i 2 2 M 1 A i i 1 R end v d eig A d eig A y v 1 z linspace H 2 H 2 2 M 1 y11 y 2 M2 trapz z y11 y12 1 M2 y11 plot z y12 y 1 M2 0 5 y plot z y hold on 13 附附 1 21 2 第一激发态一维线性谐振子第一激发态一维线性谐振子 clc clear M 200 H 20 h H 2 M R 1 h 2 for i 1 2 M 1 for j 1 2 M 1 A i j 0 end end for p 1 2 M 1 A p p 2 R p M h 2 end for i 1 2 M 2 A i i 1 R end for i 2 2 M 1 A i i 1 R end v d eig A d eig A y v 2 z linspace H 2 H 2 2 M 1 y11 y 2 M2 trapz z y11 y12 1 M2 y11 plot z y12 y 1 M2 0 5 y plot z y hold on 附附 1 31 3 第二激发态一维线性谐振子第二激发态一维线性谐振子 clc clear M 200 H 20 h H 2 M 14 R 1 h 2 for i 1 2 M 1 for j 1 2 M 1 A i j 0 end end for p 1 2 M 1 A p p 2 R p M h 2 end for i 1 2 M 2 A i i 1 R end for i 2 2 M 1 A i i 1 R end v d eig A d eig A y v 3 z linspace H 2 H 2 2 M 1 y11 y 2 M2 trapz z y11 y12 1 M2 y11 plot z y12 y 1 M2 0 5