用MATLAB语言解氢原子与类氢离子的定态薛定谔方程

氢原子薛定谔方程引言薛定谔方程是量子力学的基石之一，描述了微观粒子的行为。

而氢原子是最简单的原子系统，因此研究其薛定谔方程有助于我们理解量子力学的基本原理。

本文将深入探讨氢原子薛定谔方程，从基本概念到具体计算，全面分析该方程的背景、推导和解析。

薛定谔方程简介薛定谔方程是描述量子系统的一维时间无关定态的方程。

对于一个粒子的波函数ψ(x)、能量E和势能V(x)，薛定谔方程可以写作：Ĥψ(x)=Eψ(x)其中，Ĥ是哈密顿算符，定义为Ĥ=−ℏ22md2dx2+V(x)，ℏ是约化普朗克常数，m是粒子的质量，x是粒子的位置。

对于氢原子，势能V(x)由于原子核和电子之间的相互作用而产生。

氢原子的薛定谔方程氢原子是由一个质子和一个电子构成的，因此氢原子的薛定谔方程是描述电子在氢原子中的运动。

使用球坐标系，薛定谔方程可以重写为：[−ℏ22m(1r2ddr(r2ddr)−L̂22mr2)+V(r)]ψ(r,θ,ϕ)=Eψ(r,θ,ϕ)其中，L̂2是角动量算符的平方，定义为L̂2=−ℏ2(1sinθddθ(sinθddθ)+1sin2θd2dϕ2)。

氢原子的径向方程为了简化氢原子的薛定谔方程，我们考虑分离变量，假设波函数可以表示为一个径向部分和一个角向部分的乘积：ψ(r,θ,ϕ)=R(r)Y(θ,ϕ)。

代入薛定谔方程并分离变量，可以得到径向方程和角向方程。

径向方程的推导通过分离变量，我们将薛定谔方程转化为径向方程和角向方程。

径向方程可以通过将薛定谔方程乘以r2并对角度积分得到。

经过一系列数学推导，可以得到氢原子的径向方程为：[−ℏ22md2dr2+ℏ22ml(l+1)r2+V(r)−E]R(r)=0其中，l是角量子数，通过求解该方程可以得到径向波函数R(r)和能量E。

解析解与数值解氢原子的薛定谔方程可以通过解析方法求解，得到精确的解析解。

然而，尽管存在解析解，推导和计算过程非常复杂，通常需要使用数值方法来近似求解。

氢原子的薛定谔方程在量子力学中，薛定谔方程是描述微观粒子运动的基本方程之一。

对于氢原子来说，薛定谔方程起着至关重要的作用，它能够描述氢原子中电子的运动状态和能级分布，为我们理解氢原子的结构和性质提供了重要依据。

氢原子由一个质子和一个围绕质子运动的电子组成。

在薛定谔方程中，波函数描述了电子的运动状态，包括位置和动量等信息。

通过求解薛定谔方程，我们可以得到氢原子中电子的能级和波函数，从而揭示出氢原子的量子性质。

薛定谔方程的解可以分解为径向部分和角向部分，分别描述了电子在氢原子中径向和角向的运动。

径向部分的解决定了氢原子中电子的轨道半径和能级，而角向部分则描述了电子在轨道上的运动方式。

通过这两部分的解，我们可以全面了解氢原子中电子的运动规律。

薛定谔方程的一个重要应用是计算氢原子的能级结构。

通过求解薛定谔方程，我们可以得到氢原子中不同能级的能量和波函数。

这些能级决定了氢原子的光谱线，可以用来解释氢原子在不同波长下的吸收和发射现象。

因此，薛定谔方程的解不仅可以帮助我们理解氢原子的内部结构，还可以解释氢原子的光谱特性。

除了氢原子外，薛定谔方程还可以应用于其他原子和分子系统的研究。

通过对薛定谔方程的求解，我们可以得到不同原子和分子系统的波函数和能级，从而揭示它们的量子性质和相互作用规律。

这为我们研究原子和分子的结构、性质和反应机制提供了重要的理论基础。

总的来说，薛定谔方程是量子力学中的重要方程之一，对于理解氢原子和其他微观粒子系统的性质和行为具有重要意义。

通过求解薛定谔方程，我们可以揭示微观世界的奥秘，探索物质世界的微观规律，为科学技术的发展提供重要支持。

希望未来能有更多科学家通过对薛定谔方程的研究，揭示出更多微观世界的奥秘，推动人类对自然界的认识和探索。

matlab数值薛定谔方程摘要：I.引言- 介绍薛定谔方程- 介绍matlab 数值求解方法II.薛定谔方程的数值求解方法- 有限差分法- 有限元法- 谱方法III.matlab 数值求解薛定谔方程的步骤- 准备薛定谔方程的数值模型- 选择数值求解方法- 编写matlab 代码- 运行代码，分析结果IV.结果与讨论- 结果展示- 结果分析- 结果验证V.结论- 总结matlab 数值求解薛定谔方程的方法- 展望未来的研究方向正文：I.引言薛定谔方程是量子力学中的一个基本方程，用于描述一个微观粒子在给定势能场中的运动状态。

然而，由于薛定谔方程本身是一个偏微分方程，它的求解在大多数情况下是非常困难的。

matlab 作为一种强大的科学计算软件，可以用于数值求解薛定谔方程。

本文将介绍薛定谔方程的数值求解方法，以及如何使用matlab 进行数值求解。

II.薛定谔方程的数值求解方法薛定谔方程的数值求解方法主要有以下几种：1.有限差分法：将薛定谔方程的解表示为离散的点，通过差分代替微分，将方程转化为一个线性代数方程组，从而求解薛定谔方程。

2.有限元法：将薛定谔方程的解表示为有限个基函数的线性组合，通过插值或逼近基函数，将方程转化为一个线性代数方程组，从而求解薛定谔方程。

3.谱方法：通过在一组基函数上将薛定谔方程进行投影，将方程转化为一个线性代数方程组，从而求解薛定谔方程。

III.matlab 数值求解薛定谔方程的步骤使用matlab 进行数值求解薛定谔方程的步骤如下：1.准备薛定谔方程的数值模型：首先需要根据实际问题建立薛定谔方程的数值模型，包括势能场、边界条件等。

2.选择数值求解方法：根据问题的特点和求解需求，选择合适的数值求解方法，如有限差分法、有限元法或谱方法。

3.编写matlab 代码：根据所选方法，编写matlab 代码，实现薛定谔方程的数值求解。

4.运行代码，分析结果：运行编写的matlab 代码，得到薛定谔方程的数值解。

氢原子的薛定谔方程
氢原子是最简单的原子之一，由一个质子和一个电子组成。

在量子力学中，描述氢原子的运动状态的数学模型就是薛定谔方程。

薛定谔方程是量子力学的基础方程之一，它描述了微观粒子在势场中的运动规律。

薛定谔方程是一个偏微分方程，描述了波函数随时间和空间的演化。

波函数包含了粒子的所有信息，包括位置、动量等。

在氢原子的情况下，薛定谔方程可以被简化为一个径向部分和一个角向部分的乘积。

径向部分描述了电子在原子核周围运动的距离，角向部分描述了电子在不同方向上的概率分布。

通过求解薛定谔方程，我们可以得到氢原子的能级和波函数，从而进一步研究原子的性质和行为。

薛定谔方程的求解需要考虑原子核和电子之间的相互作用，以及外加的势场对电子的影响。

通过引入适当的近似和数值方法，可以求解薛定谔方程并得到氢原子的能级和波函数。

氢原子的能级是量子化的，即只能取离散的数值。

能级越高，电子离原子核越远，能量也越大。

每个能级对应一个波函数，描述了电子在原子周围的分布情况。

薛定谔方程的求解不仅可以用于氢原子，还可以推广到其他原子和
分子系统。

通过求解薛定谔方程，我们可以理解原子和分子的结构、性质和反应规律，为化学和物理学的发展提供重要的理论基础。

薛定谔方程是描述氢原子和其他微观粒子运动的重要方程，它揭示了量子力学世界的奥秘。

通过求解薛定谔方程，我们可以深入理解原子和分子的微观世界，为科学研究和技术应用提供重要支持。

希望未来能够进一步探索量子力学的奥秘，推动科学的发展和进步。

氢原子的薛定谔方程
薛定谔方程是一个著名的电子结构理论，可以用来描述一个原子的电子状态。

它是一个带有四个变量的复合实现方程，被称为薛定谔方程。

它由20世纪伟大的物理学家Ernst Schrdinger发明，他是量子力学的创始人。

当谈到氢原子时，薛定谔方程还可以用来解释它的电子状态。

氢原子只有一个电子，因此为了解释它的电子状态，只需要一个薛定谔方程。

薛定谔方程可以如下表达：
iψ/t = ^2/2m·^2ψ + Vψ
其中，ψ表示波函数；i是虚数单位；表示普朗克常数，ψ/t表示时间导数；m是电子的质量；^2表示laplace算符；V表示电子的势能。

薛定谔方程简写为：
Hψ = εψ
其中，H表示哈密顿量，ε表示电子的能量。

对氢原子的薛定谔方程可以写为：
[^2/2m·^2+ V(r)E]ψ(r) = 0
其中，V(r)表示电子势能，E表示电子能量，r表示电子的位置半径。

解决氢原子的薛定谔方程需要一些技巧——定义一个适应性正交基函数组，利用拉普拉斯算符变换到正交空间，然后使用矩阵方法解决。

有时，哈密顿量可以被简化为一个对角矩阵，这一点取决于电
子势能的类型。

任何时候，电子能量的计算都是从在某个特定的位置的电子的能量开始的。

氢原子可以通过薛定谔方程来解释，并且可以计算出它的电子能量，解释的结果可以用来解释它的原子结构。

薛定谔方程对氢原子的电子状态起着至关重要的作用。

薛定谔方程求解氢原子
氢原子的薛定谔方程为：（−h¯22m∇2+V）ψ=Eψ（−h28π2m∇2−Ze24πε0r）ψ=Eψ。

薛定谔方程（Schrödinger equation），又称薛定谔波动方程（Schrodinger wave equation），是由奥地利物理学家薛定谔提出的量子力学中的一个基本方程，也是量子力学的一个基本假定。

它是将物质波的概念和波动方程相结合建立的二阶偏微分方程，可描述微观粒子的运动，每个微观系统都有一个相应的薛定谔方程式，通过解方程可得到波函数的具体形式以及对应的能量，从而了解微观系统的性质。

在量子力学中，粒子以概率的方式出现，具有不确定性，宏观尺度下失效可忽略不计。

薛定谔方程（Schrodinger equation）在量子力学中，体系的状态不能用力学量（例如x）的值来确定，而是要用力学量的函数Ψ（x，t），即波函数（又称概率幅，态函数）来确定。

力学量取值的概率分布如何，这个分布随时间如何变化，这些问题都可以通过求解波函数的薛定谔方程得到解答。

这个方程是奥地利物理学家薛定谔于1926年提出的，它是量子力学最基本的方程之一。

氢原子或类氢离子的薛定谔方程一般解的探讨根据关于薛定谔方程的定义，以及它的应用，本文讨论的是氢原子或类氢离子的薛定谔方程一般解。

首先，我们来谈谈薛定谔方程的定义。

薛定谔方程是一个常微分方程，它的解描述了量子系统的运动状态。

它既可以用于描述物理系统，也可以用于描述物理状态，比如电子，氢原子等。

根据这个方程，得出的结果是描述物理系统的运动状态在不同时间点上的解析式。

这个方程也可以用于描述复杂的流体流动状态，以及微观世界中物质间的相互作用。

其次，让我们来看看薛定谔方程在氢原子或类氢离子中的应用。

这些氢原子或类氢离子被称为单子系统，因为它们都是由一个复杂的几何结构组成的，而这个结构可以用一个常微分方程来描述。

薛定谔方程可以用来描述氢原子或类氢离子中电子质量和其他参数的变化，从而可以得出氢原子或类氢离子中电子能量层次的变化，以及电子轨道结构的变化。

通过薛定谔方程，可以得出在不同的电子能量层次上电子的位置以及它们的概率分布。

这一结果可以为量子力学提供新的见解。

最后，让我们来看看薛定谔方程一般解的探讨。

通常情况下，可以使用标准的数学方法解决薛定谔方程，包括方程组求解、数值解以及符号运算。

有时，也可以采用克服难题的改进方法，比如用一组有穷基组求解，或者使用李萨克法求解，以及扩展薛定谔方程等方法，以解决薛定谔方程的一般解。

总之，通过本文，我们可以了解到薛定谔方程在氢原子或类氢离子中的应用，以及薛定谔方程一般解的探讨。

本文就此主题做了深入浅出的探讨，以期能够为读者更好地认识氢原子或类氢离子的薛定谔方程的基础及其一般解的解析式的获得提供便利。

matlab数值薛定谔方程薛定谔方程是描述量子力学中粒子的行为的基本方程。

在数值计算中，我们可以使用数值方法来求解薛定谔方程。

下面我将从多个角度来回答关于在MATLAB中数值求解薛定谔方程的问题。

1. 数值方法的选择：在MATLAB中，我们可以采用多种数值方法来求解薛定谔方程，其中常用的方法包括有限差分法、有限元法和谱方法等。

选择合适的数值方法取决于问题的特点和计算资源的可用性。

2. 离散化：在数值计算中，我们需要将薛定谔方程离散化为有限个点上的代数方程。

通常，我们会将空间离散化为网格，并在每个网格点上计算波函数的值。

时间离散化则是通过迭代的方式逐步求解时间演化。

3. 有限差分法：有限差分法是一种常见的数值方法，它将导数近似为有限差分。

在薛定谔方程中，我们可以将二阶导数近似为中心差分，然后使用差分方程来求解离散化的薛定谔方程。

4. 有限元法：有限元法是一种广泛应用于偏微分方程求解的数值方法。

在薛定谔方程中，我们可以使用有限元法将波函数表示为一组基函数的线性组合，并通过求解线性方程组来确定系数。

5. 谱方法：谱方法是一种基于函数展开的数值方法，它使用一组特定的基函数来表示波函数。

在薛定谔方程中，我们可以使用傅里叶级数或其他正交多项式作为基函数，并通过求解线性方程组来确定系数。

6. 边界条件：在数值求解薛定谔方程时，我们需要指定合适的边界条件。

常见的边界条件包括固定边界条件和周期性边界条件，具体取决于问题的物理背景。

7. 算法实现：在MATLAB中，我们可以使用内置的数值计算函数和工具箱来实现数值求解薛定谔方程。

例如，可以使用MATLAB的PDE Toolbox来求解偏微分方程，或者使用MATLAB的FFT函数来进行傅里叶变换。

总结起来，数值求解薛定谔方程是一个复杂而重要的问题，需要根据具体情况选择合适的数值方法并进行适当的离散化和边界条件处理。

MATLAB提供了丰富的数值计算工具和函数，可以帮助我们实现数值求解薛定谔方程的算法。

MATLAB数值薛定谔方程介绍薛定谔方程是量子力学中描述微观粒子行为的基本方程之一。

它描述了粒子的波函数随时间的演化。

在实际研究中，常常需要通过数值方法来求解薛定谔方程，特别是对于复杂的体系或无法通过解析方法求解的情况。

MATLAB作为一种强大的数值计算工具，提供了丰富的函数和工具箱，可以帮助我们求解薛定谔方程。

本文将介绍如何使用MATLAB进行数值求解，并给出一些示例代码和注意事项。

数值方法求解薛定谔方程通常需要使用数值方法，其中最常用的方法之一是有限差分法。

有限差分法将波函数离散化为网格点上的数值，通过近似微分来代替薛定谔方程中的导数项，从而转化为一个矩阵方程。

具体来说，我们可以将一维薛定谔方程表示为：iℏ∂Ψ(x,t)∂t=−ℏ22m∂2Ψ(x,t)∂x2+V(x)Ψ(x,t)其中，Ψ(x,t)是波函数，m是粒子的质量，V(x)是势能函数。

为了使用有限差分法求解，我们将空间坐标x离散化为网格点x i，时间t离散化为时间步长Δt，波函数Ψ(x,t)在网格点上的值用Ψi n表示，其中i表示网格点的索引，n表示时间步的索引。

将导数项用中心差分近似表示，我们可以得到：iℏΨi n+1−Ψi nΔt=−ℏ22mΨi+1n−2Ψi n+Ψi−1nΔx2+V i nΨi n其中，Δx是空间步长，V i n表示势能函数在网格点x i上的值。

通过这个差分方程，我们可以逐步更新波函数的值，从而得到波函数随时间的演化。

MATLAB代码示例下面是一个简单的MATLAB代码示例，演示如何使用有限差分法求解一维薛定谔方程。

% 定义参数hbar = 1; % 约化普朗克常数m = 1; % 粒子质量L = 10; % 空间范围N = 1000; % 网格点数dx = L/N; % 空间步长dt = 0.01; % 时间步长% 初始化波函数x = linspace(-L/2, L/2, N); % 空间坐标psi = exp(-x.^2); % 初始波函数% 求解薛定谔方程for n = 1:1000% 计算势能函数V = 0.5*x.^2;% 更新波函数psi = psi - 1i*dt*(hbar/(2*m))*(circshift(psi,-1,2)-2*psi+circshift(psi,1, 2))/(dx^2) - 1i*dt*V.*psi;% 绘制波函数随时间的演化plot(x, abs(psi).^2);xlim([-L/2, L/2]);ylim([0, 1]);xlabel('x');ylabel('|\psi|^2');title(['Time step ', num2str(n)]);drawnow;end在这个示例中，我们假设粒子质量m=1，空间范围L=10，网格点数N=1000。

写出一份matlab对氢原子波函数仿真模型的报告氢原子波函数仿真模型报告1. 引言在量子物理领域中，研究氢原子波函数是非常重要的。

波函数用于描述氢原子的能级和态的特性，对理解氢原子的电子结构和行为具有重要意义。

本报告旨在介绍使用MATLAB编写的氢原子波函数仿真模型。

2. 模型原理氢原子波函数基本理论氢原子的波函数可以使用Schrodinger方程进行描述，波函数的形式可通过求解Schrodinger方程得出。

在模型中，我们使用量子数n 和l来描述氢原子的态，其中n表示主量子数，l表示角量子数。

模型工具和算法本模型使用MATLAB进行编写和实现。

主要使用的算法包括数值解法和绘图函数。

具体而言，我们使用数值方法求解氢原子的Schrodinger方程，并利用MATLAB的绘图函数可视化波函数的形状。

3. 模型实现数据处理首先，我们需要定义氢原子的参数，包括质量、电荷、Planck常数等。

然后，使用数值方法求解氢原子的Schrodinger方程，得到波函数的数值解。

绘图根据数值解得到的结果，我们使用MATLAB的绘图函数绘制波函数在三维空间中的形状。

可以根据不同的量子数n和l，绘制不同的波函数图像。

4. 模型应用与结果分析模型应用本模型可以应用于研究氢原子的能级和态的特性。

可以通过对不同量子数的波函数进行分析，得到氢原子的电子分布和结构。

结果分析通过对氢原子波函数的仿真，我们可以观察到不同量子数对应的波函数形状和能级分布。

这些结果与已知的氢原子理论相吻合，验证了模型的正确性和可靠性。

5. 总结与展望本报告介绍了一种基于MATLAB编写的氢原子波函数仿真模型。

通过数值求解Schrodinger方程和绘图函数，实现了对氢原子波函数的模拟。

模型在氢原子波函数的研究中具有重要应用价值。

未来可以进一步优化模型的求解算法，提高模型的计算效率和精度。

参考文献[1] Griffiths, D. J. (2005). Introduction to quantum mechanics. Pearson Education.[2] Matlab Documentation. (2021). Retrieved from。

在MATLAB中模拟定态氢原子波函数通常涉及解决定态薛定谔方程。

对于氢原子，波函数可以表示为径向部分和角部分（通常使用球谐函数来描述角部分）的乘积。

以下是使用MATLAB解决定态氢原子问题的基本步骤：1. 定义氢原子的势能函数V(r) = -1/r，其中r是电子与核之间的距离。

2. 写出定态薛定谔方程并将其简化为径向薛定谔方程，因为波函数是径向对称的。

3. 将简化后的薛定谔方程转换为一个一维二阶微分方程，然后求解该方程。

4. 计算波函数的概率密度函数，即径向部分的模平方|R(r)|^2。

5. （可选）绘制波函数的图形，展示其在不同量子数n、l（角动量量子数）和m（磁量子数）下的变化。

以下是一个简单的MATLAB代码示例，用于计算并绘制氢原子的第一定态（n=1）的径向波函数：氢原子第一定态的径向波函数势能常数（氢原子核电荷）Z = 1; 对于氢原子，核电荷数Z=1势能函数（以Z为单位）V = -1/sqrt(pi); V(r) = -Z/r，取Z=1能量（对于n=1，l=0，氢原子第一定态）E = -13.6eV; E = -13.6 eV * n^2，取n=1计算径向波函数的参数k = sqrt(E);定义计算范围r = linspace(0.01, 10, 1000); 从0.01到10个单位长度，共1000个点计算径向波函数R(r)Rr = exp(-V*r) * sin(k*r);计算概率密度函数|R(r)|^2prob_density = Rr * conj(Rr);绘图figure;plot(r, prob_density);xlabel('Radius (bohr)');ylabel('Probability density');title('Probability density of the ground state hydrogen atom');grid on;请注意，这个代码示例仅计算和绘制了第一定态（n=1，l=0）的径向波函数的绝对平方，它代表了在给定半径下找到电子的概率密度。

定态薛定谔方程的MATLAB求解(一)利用矩阵法对定态薛定谔方程的MATLAB求解摘要：本文首先对薛定谔方程的提出及发展做了一个简单介绍。

然后，以在一维空间运动的粒子构成的谐振子的体系为例，详细介绍了矩阵法求解薛定谔方程的过程及公式推导。

最后，通过MATLAB编程仿真实现了求解结果。

关键词：定态薛定谔方程求解矩阵法MATLAB仿真薛定谔方程简介1.1背景资料薛定谔方程是由奥地利物理学家薛定谔提出的量子力学中的一个基本方程，是将物质波的概念和波动方程相结合建立的二阶偏微分方程，可描述微观粒子的运动，每个微观系统都有一个相应的薛定谔方程式，通过解方程可得到波函数的具体形式以及对应的能量，从而了解微观系统的性质。

其仅适用于速度不太大的非相对论粒子，其中也没有包含关于粒子自旋的描述。

当计及相对论效应时，薛定谔方程由相对论量子力学方程所取代，其中自然包含了粒子的自旋。

薛定谔方程建立于1926年。

它是一个非相对论的波动方程。

它反映了描述微观粒子的状态随时间变化的规律，它在量子力学中的地位相当于牛顿定律对于经典力学一样，是量子力学的基本假设之一。

设描述微观粒子状态的波函数为Ψ（r，t），质量为m的微观粒子在势场V （r，t）中运动的薛定谔方程为在给定初始条件和边界条件以及波函数所满足的单值、有限、连续的条件下，可解出波函数Ψ（r，t）。

由此可计算粒子的分布概率和任何可能实验的平均值（期望值）。

当势函数V不依赖于时间t时，粒子具有确定的能量，粒子的状态称为定态。

定态时的波函数可写成式中Ψ（r）称为定态波函数，满足定态薛定谔方程，这一方程在数学上称为本征方程，式中E为本征值，是定态能量，Ψ（r）又称为属于本征值E的本征函数。

量子力学中求解粒子问题常归结为解薛定谔方程或定态薛定谔方程。

薛定谔方程揭示了微观物理世界物质运动的基本规律，被广泛地用于原子物理、核物理和固体物理，对于原子、分子、核、固体等一系列问题中求解的结果都与实际符合得很好。

⼀维定态薛定谔⽅程求解的两种⽅法（matlab）量⼦⼒学中，薛定谔⽅程是核⼼。

薛定谔的猫描述了态的概念，但实际研究中，要想细致地研究⼀个原⼦，分⼦，甚⾄⼀块物质，都需要从薛定谔⽅程的求解开始。

下⾯将会以我的⼀次作业的题⽬为例，向⼤家展⽰整个求解过程。

薛定谔⽅程的完整形式为：以上⽅程有对时间的微分，还有对空间的微分。

⽽对于定态的薛定谔⽅程，我们只需考虑某⼀时刻的波函数，所以直接可将能量算符替代为E（⼀个常数）。

(1)分段势能法对于空间的梯度，如果只是⼀维情况的话，可以直接将梯度算符改为微分。

所以⼀维定态薛定谔⽅程就显得很简单：就是⼀个简单的⼆阶微分⽅程。

此⽅程的解想必⼀眼就可以看出来。

就是这个解是假设U(x)与x⽆关，是⼀个常数才得出这个⾃由波的解。

类似与微积分中的⽅法，对于⼀个任意势场函数，我们可以假设在某⼀个极⼩的dt范围内，势函数是不变的，因此可以将任意⼀个势函数⽤有限个⼀定宽度的恒定势场来代替。

如下图所⽰：其中的各个⼩段的波函数就可以表⽰为这样就会有2N个⽅程，然后利⽤内部的n-1个边界条件（界⾯处波函数连续，波函数的倒数连续），和两端的衔接（假设⼊射为1，则A1=1, B1=r；且最终透射端没有反射波，AN=t, BN=0. ），就可以写出2N个线性⽆关的⽅程，从⽽可以将系数都求解出来。

注意，这种情况下，我们⽆从得知基态的能量值，以及能量的分⽴的特性。

但是从这种⾓度出发，我们可以很容易计算出波在这样的势函数中传输特性，可以计算出⼊射端的反射系数R，以及不同能量所对应的⼊射波的透射系数T。

下⾯将以⼀个例⼦应⽤上述关系。

根据上图中所⽰的势函数求解薛定谔⽅程，得到透射系数和反射系数随温度的变化关系为(2)差分法现在我们从另外⼀个⾓度出发，⼀维定态薛定谔⽅程如下在这⾥，我们要求的是，可以将分为N份，采⽤数值计算⽅法，将微分⽅程变成差分⽅程。

参考相应书籍可知可以化为对于上述波函数也可以转化为类似的形式，所以可以由矩阵T的特征值对应能量，特征向量对应于波函数在每⼀个节点的解。

如何用matlab解薛定谔方程？数值求解的无量纲化技术我前面讲了 matlab解二次微分方程的方法。

薛定谔方程当然是个二次微分方程. 所以，上一讲的matlab 的ode函数是可以解薛定谔方程的。

不过，在求解之前，我们还有个工作必须先做。

薛定谔方程中有个hbar，它的数值是如此之小，而且还要平方。

还有电子电荷e，光速c, 电子质量m 这样的数值也是如此。

这样的数值是不适合在计算程序中出现的。

凡是天文数值都不适合在计算程序中出现。

有个很优美的技术来消除它们，就是无量纲化。

这个技术是我们做计算的时候必须做的，所以，我在这里讲讲这个事情。

无量纲化，就是用一些特征的长度做长度单位，用一些特征的能量做能量单位。

假设我们研究的问题是原子，我们就可以用玻尔半径a = 0.529埃为长度单位，以氢原子基态能量的绝对值 13.6eV 为能量单位。

为了用它们做无量纲化，我们需要它们的公式形式：a = hbar^2 /me^2, |E0| = e^2/2a。

氢原子的径向波函数满足的薛定谔方程是[（-hbar^2/2m) （d^2/dr^2 + (2/r) d/dr) - e^2/r] R(r) = E R(r).把这方程两边除上述|E0|，得到 [a^2 （d^2/dr^2 + (2/r) d/dr) - 2a/r] R(r) = E R(r)。

这里的E 是以|E0| = 13.6eV 为单位的。

然后，把 a 除到导数下面的r上，方程就变成[（d^2/dr^2 + (2/r) d/dr) - 2/r] R(r) = E R(r)，这里的r 是以a 为单位的。

这个方程里面的每个量都仅仅是一些无量纲的数了，方程大大简化了。

我们最后需要求解的方程是这个无量纲化的薛定谔方程：[（d^2/dr^2 + (2/r) d/dr) - 2/r] R(r) = E R(r)。

这方程怎么解，上一讲已经讲过了。

好不好懂，请给个意见。

摘要：此文介绍了一种使用MATLAB求解一维定态薛定谔方程的方法。

利用充分格式进行离散化，得出相应的矩阵方程，用MATLAB求解本征值和本征函数。

此方法简单可靠，可以处理各种时间无关的束缚态问题。

宏观物体遵循的是牛顿运动方程，给定初始条件以及受力条件，我们就可以求出任意时刻粒子的状态。

而在原子尺度上，所有粒子都表现出波的行为，粒子的状态用波函数Ψ(x,y,z,t)来描述，描述微观粒子运动的方程是由薛定谔于1925年提出薛定谔方程。

微观粒子的运动与其所处的势场相关，当势场不随时间变化时，称为定态，一维情况下的定态薛定谔方程为−ℏ22md2ψ(x)dx2+U(x)ψ(x)=Eψ(x) (1)其中，U(x)表示粒子所处的势场。

由定态波函数可以得出总波函数Ψ(x,t)=ψn(x)e−Ent/ℏ1926年，Born提出，波函数模的平方|Ψ(x,t)|2代表在位置x，时刻t寻找到电子的概率，这就是波函数的条统计解释。

可以由遵循微分方程的波函数表示，薛定谔波方程涉及空间坐标和时间。

对于一维情况，在时刻t时，在x1和x2之间的某处找到电子的概率由下式给出P(t)=∫x1x2Ψ(x,t)Ψ∗(x,t)dx时间无关薛定谔方程(1)是系统的能量本征方程。

该特征值方程通常由一组特定的函数ψ1(x),ψ2(x),ψ3(x),...和相应的一组常数E1,E2,E3,...满足解的条件，被称为是哈密顿算符的本征函数和相应的本征值。

当测量处于ψn(x)状态的系统的能量时，结果将始终为En。

对于束缚态情况，必须有ψn(x)=0,ifx→±∞一维定态薛定谔方程是二阶微分方程，但是，能够解析求解的情况屈指可数，如氢原子，谐振子，无限深势阱等。

随着计算机技术的发展，我们可以数值上进行求解。

本文利用MATLAB软件，使用矩阵方法求解束缚态的本征值问题。

对于原子系统，以nm和eV的能量测量长度更方便。

我们可以使用缩放因子Length:Lse=1×10−9,convert m to nm;Energy:Ese=1.6×10−19,convert J to eV.因此我们可以将等式(1)写成[Csed2dx2+U(x)]ψ(x)=Eψ(x)whereCse=−ℏ22m(1Lse2Ese)为了求解上述方程式，首先进行离散化处理。