matlab解薛定谔方程,定态薛定谔方程的MATLAB求解(一)
定态薛定谔方程的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]
对于存在的波函数,应为实数,所有的矩阵元也为实数,由厄密算符的性质得
为了计算坐标的矩阵元,由对易关系
又 代入上式易得
写为矩阵形式,有
根据矩阵的乘法规则,有
又,则有由前面的分析知,只有时,才存在矩阵元,代入上式,
从该方程我们可以得出
矩阵元不为零,但是当时,矩阵元则
即
又
依次类推,得出
最后,我们得到坐标矩阵元不为零的表达式
又谐振子的能量可以用来表示,且,计算该能量得
matlab解薛定谔方程,定态薛定谔方程的MATLAB求解(一)相关推荐
- 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坡度计算公式,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求解本征值和本征函数.此方法简单可靠,可以处理各种时间无关的束缚态问题.所用 ...
- matlab解薛定谔方程,有限差分法解薛定谔方程与MATLAB实现
第30卷 第3期 高 师 理 科 学 刊 Vol. 30 No.3 2010年 5月 Journal of Science of Teachers′ College and University Ma ...
最新文章
- 打印机打印一次出三张_一次性帮你搞定3D打印机堵塞问题
- java httpclient put_[工具类-HttpClientUtils]HttpClient之GET PUT DELETE POST
- ITK:分段线性曲线的数据结构
- 05模块-barcode【条码扫描】
- 感染性的木马病毒分析之样本KWSUpreport.exe
- python使用指南(part1)--去重
- Qt for Android环境配置
- ASP+COM技术在嵌入式Webserver中的应用
- BZOJ 2527 Meteors | 整体二分
- python 目标检测 训练_YOLOv3目标检测有了TensorFlow实现,可用自己的数据来训练
- 判断某个字符串是否是数字(简单方法)
- VMWaer克隆centos后网络的问题解决
- 人工智能能为旅游业带来多大的想象空间?
- Python代码实现验证码识别
- 赤峰中考计算机考试软件,2017年内蒙古赤峰中考信息技术操作考试实施细则
- 硕思闪客精灵v2021如何用Flash制作动画短片流程分享
- Excel函数——小数点计算
- 【深度之眼cs231n第七期】笔记(二十七)
- 统计给定的n个数中,负数、零和正数的个数。
- Linux(xUbuntu18.04)安装VMware workstation15.1 解决报错:Before you can runVMware