matlab解薛定谔方程,定态薛定谔方程的MATLAB求解(一).doc
定态薛定谔方程的MATLAB求解(一)
利用矩阵法对定态薛定谔方程的MATLAB求解
摘要:本文首先对薛定谔方程的提出及发展做了一个简单介绍。然后,以在一维空间运动的粒子构成的谐振子的体系为例,详细介绍了矩阵法求解薛定谔方程的过程及公式推导。最后,通过MATLAB编程仿真实现了求解结果。
关键词:定态薛定谔方程求解 矩阵法 MATLAB仿真
薛定谔方程简介
1.1背景资料
薛定谔方程是由奥地利物理学家薛定谔提出的量子力学中的一个基本方程,是将物质波的概念和波动方程相结合建立的二阶偏微分方程,可描述微观粒子的运动,每个微观系统都有一个相应的薛定谔方程式,通过解方程可得到波函数的具体形式以及对应的能量,从而了解微观系统的性质。其仅适用于速度不太大的非相对论粒子,其中也没有包含关于粒子自旋的描述。当计及相对论效应时,薛定谔方程由相对论量子力学方程所取代,其中自然包含了粒子的自旋。
薛定谔方程建立于 1926年。它是一个非相对论的波动方程。它反映了描述微观粒子的状态随时间变化的规律,它在量子力学中的地位相当于牛顿定律对于经典力学一样,是量子力学的基本假设之一。设描述微观粒子状态的波函数为Ψ(r,t),质量为m的微观粒子在势场V(r,t)中运动的薛定谔方程为
在给定初始条件和边界条件以及波函数所满足的单值、有限、连续的条件下,可解出波函数Ψ(r,t)。由此可计算粒子的分布概率和任何可能实验的平均值(期望值)。当势函数V不依赖于时间t时,粒子具有确定的能量,粒子的状态称为定态。定态时的波函数可写成式中Ψ(r)称为定态波函数,满足定态薛定谔方程,这一方程在数学上称为本征方程,式中E为本征值,是定态能量,Ψ(r)又称为属于本征值E的本征函数。
量子力学中求解粒子问题常归结为解薛定谔方程或定态薛定谔方程。薛定谔方程揭示了微观物理世界物质运动的基本规律,被广泛地用于原子物理、核物理和固体物理,对于原子、分子、核、固体等一系列问题中求解的结果都与实际符合得很好。
定态薛定谔方程直角坐标系形式
定态薛定谔方程球坐标系形式
1.2定态薛定谔方程
条件
V(r,t)=V(r), 与t无关。
用分离变量法, 令Ψ=φ(r)f(t),代入薛定谔方程,得两个方程:
此称定态薛定谔方程
整个定态波函数形式:
特点:
波函数由空间部分函数与时间部分函数相乘;
B.时间部分函数是确定的。
定态波函数几率密度W与t无关,几率分布不随时间而变,因此称为定态。
1.3本征方程、本征函数与本征值
算符: 本征方程:
λ:本征值,有多个,甚至无穷多个
ψλ:本征值为λ的本征函数,也有多个,甚至无穷多个,有时一个本征值对应多个不同的本征函数,这称为简并。若一个本征值对应的不同本征函数数目为N,则称N重简并。
1.4 定态情况下的薛定谔方程一般解
1、定态薛定谔方程或不含时的薛定谔方程是能量本征方程,E就称为体系的能量本征值,而相应的解称为能量的本征函数。
2、当不显含时时,体系的能量是收恒量,可用分离变量。
3、解定态薛定谔方程,关键是写出哈密顿量算符。
2. 利用矩阵法求解薛定谔方程
以在一维空间运动的粒子构成的谐振子的体系为例。
该粒子的势能是,是谐振子的角频率,因此谐振子的哈密顿量为
。
当时,谐振子的势能变为无穷大,因此,粒子只能在有限的空间上运动,并且能量值谱是分立的。下面采用矩阵的方法,确定谐振子的能量分立值。
从运动方程出发 (1)
而势能 那么
又代入上式(1)得
即(2)
在矩阵形式下,该方程可以写为
含时坐标矩阵元 (3)
对它求导,我们得到
代入上式后,有
(4)
其中 (5)
所以,除了当或外,所有的坐标矩阵元都等于零
当时,由(5)式有
即 同理,
因此,只有变化时,才能得到频率即 所以不为零的坐标矩阵元为
根据定义[12-14]
对于存在的波函数,应为实数,所有的矩阵元也为实数,由厄密算符的性质得
为了计算坐标的矩阵元,由对易关系
又 代入上式易得
写为矩阵形式,有
根据矩阵的乘法规则,有
又,则有由前面的分析知,只有时,才存在矩阵元,代入上式,
从该方程我们可以得出
矩阵元不为零,但是当时,矩阵元则
即
又
依次类推,得出
最后,我们得到坐标矩阵元不为零的表达式
又谐振子的能量可以用来表示,且,计算该能量得
/amzdddc/ /ncylc/ /fbsylc/ /fbsylc/ /ywgjylc/ /ylgj/ /bhgjylc/ /nnyx/ /lzjlb/ /dzpkgw/
matlab解薛定谔方程,定态薛定谔方程的MATLAB求解(一).doc相关推荐
- matlab求薛定谔方程,定态薛定谔方程的MATLAB求解(一)
ψλ:本征值为λ的本征函数,也有多个,甚至无穷多个,有时一个本征值对应多个不同的本征函数,这称为简并.若一个本征值对应的不同本征函数数目为N,则称N重简并. 1.4 定态情况下的薛定谔方程一般解 1. ...
- matlab解方程最值点,MATLAB解方程与函数极值
1.线性方程数值求解 主要是用到了计算方法里的LU分解等不过是加快了求解速度而已相对于inv(A)*b或者A\b 2.非线性方程数值求解 1 单变量非线性方程求解 在MATLAB中提供了一个fzero ...
- matlab解方程x 2-x-2=0,matlab用三种方法求解二阶微分方程x''+0.2x'=0.4x=0.2u(t),u(t)是单位阶跃函数,初始状态为0...
问题描述: matlab用三种方法求解二阶微分方程x''+0.2x'=0.4x=0.2u(t),u(t)是单位阶跃函数,初始状态为0 1个回答 分类: 数学 2014-11-28 问题解答: 我来补答 ...
- matlab解微分方程组_MATLAB编程入门 求解常微分方程 通解 特解 数值解
一.求微分方程的通解 clear all; f1=dsolve('Dy-y=sin(x)') %默认情况下以t为自变量求解 求出的C2为积分常数 f2=dsolve('Dy-y=sin(x ...
- matlab解kmv方程,急:Matlab求解KMV模型,使用fsolve函数出现问题
最近也在做kmv模型 参考了你的和贴吧的改了一下代码,可是初始值设置总是有问题,参考了郑老师的初始值设置[1,1]可以搜索后,再改了以下代码,可是出的结果就只有1,1,但是用郑老师的代码改的话,出的结 ...
- matlab解二阶微分方程组ode,MATLAB解含参数方程、矩阵方程、二阶微分方程组
1.如下公式: 其中x.y.m.n为参数,a.b为未知数,利用MATLAB求解方程: syms x y m n a b; [a,b]=solve('x=m*cos(a)+n*cos(a+b)','y= ...
- 用matlab解出冲激响应的表达式,Matlab实验指导书
13 r2=linspace(5,0,length(t)); x2=r2.*cos(t); y2=r2.*sin(t); subplot(122),plot3(x2,y2,z) gridon 三.设计 ...
- 飞机飞行与碰撞问题matlab,飞行管理问题(建模与求解).doc
飞行管理问题(建模与求解).doc 飞行管理问题 建立了一个非线性规划模型,但由于模型求解过于复杂,设计了一个 算法,利用计算机求解,用模拟的方法建立了新的模型,并用它对给定的数据进行 计算,得出了合 ...
- matlab坡度计算公式,matlab解二阶微分方程怎么用matlab来解呢?x的定义域是(0,l),i是道路横坡坡度,w是降雨强度,k是路面横向...
共回答了16个问题采纳率:81.3% 你的方程即为: y'*y' + y''*y' +w/k =0 s = dsolve('Dy*D2y + Dy^2 + w/k ', 'y(l) = a', 'Dy ...
- matlab向量的模_基于MATLAB使用矩阵方法求解一维定态薛定谔方程
摘要:此文介绍了一种使用MATLAB求解一维定态薛定谔方程的方法.利用充分格式进行离散化,得出相应的矩阵方程,用MATLAB求解本征值和本征函数.此方法简单可靠,可以处理各种时间无关的束缚态问题.所用 ...
最新文章
- Powershell 音乐播放
- java jar log4j_java项目打包成可执行jar用log4j将日志写在jar所在目录操作
- 面试必会:HashMap 实现原理解读
- [iPhoneアプリ]iEscaper2攻略その6|龍の水晶
- c语言编程运动会分数统计系统报告,独家稿件:定稿运动会分数统计系统设计报告_完整版...
- linux学习wdlinux学堂
- 为何有些程序员总是想要“干掉”产品经理?
- ROW_NUMBER (Transact-SQL)
- SQL运行速度慢?查查中间件
- 【量产】波士顿动力机器狗,当警犬不错,上战场。。。
- java Socket实现简单在线聊天(二)
- 在VB中如何让背景图片铺满整个MDIForm
- ENSP实验五——三层交换机+二层交换机
- xynu 2139: 德莱联盟(判断线段是否相交 )
- 预留度数Rx表达式的推导
- 物流快递商家寄件运力查询接口API代码实例
- mysql创建日历表,可以按日或按月增加数据
- python手写数字识别实验报告_Python代码实现简单的MNIST手写数字识别(适合初学者看)...
- android studio 读取网络图片
- Android中拍照完就给图片加水印,Android调用相机拍照并添加水印