一维谐振子Schrdinger定态波动方程的数值求解.pdfVIP

.. 一维谐振子 Schr o dinger 定态波动方程的数值求解 作者：丁长林， 指导教师：彭芳麟 单位：北京师范大学物理系 关键词：一维谐振子 数值解 摘要： 本文运用科学数学软件 MATLAB6.1 研究一维谐振子定态波动方程的数值求解。 一维谐振 子定态解是属于有边界条件的常微分方程问题， MATLAB 软件中的 bvp4c 指令正好可以用来求 解物理上的这一类问题，本文就用它求出了一维谐振子的本征值和本征函数的数值解，与我们 已经知道的解析解非常接近。全文分四部分：一、一维谐振子定态波动方程的处理；二、有关 的 MATLAB 指令；三、数值求解程序；四、用户界面设计。 Abstract In this paper, I investigate the numerical solution of the one-dimensional harmonic oscillator by using the mathematic software MATLAB6.1. The problem about the stationary state of the one-dimensional harmonic oscillator belongs to BVPs(boundary value problems) of ODE(ordinary differential equation), which can be solved by using the function BVP4C of MATLAB on the beam. I get the numerical solution of the one-dimensional harmonic oscillator from ground state to the eleventh excited state, and contrast it with the analytic solution known to us in eigenvalue and eigenfunction, they are most adjacent. The paper consists of four parts. The first part is the simplification of stationary state equation of the one-dimensional harmonic oscillator; The second part is about the usage of some MATLAB functions needed in programming; The third part is the ――――――――――――――――――――毕业论文―――――――――――――――――――――― programme; In the fourth part, I designed the graphical user interface. 一、 一维谐振子定态波动方程的简化处理 2 2 ( / 2 )h π d 一维谐振子定态波动方程为：[− 2 +V(x)]Ψ(x) = E Ψ(x) (1) 2m dx 1 2 2 势函数 V(x)= m ω x 2 边界条件是：|x| ∞ 时，Ψ(x) 0 .