24小时热门版块排行榜    

查看: 748  |  回复: 4

husthc

木虫 (正式写手)

[求助] 求解本征态问题Au(x)=EBu(x)

首先薛定谔方程已经被矩阵化为如下形式:
Au(x)=EBu(x)


其中:A为一个对称矩阵,B为一个对角矩阵。E为能量本征值,u(x)为本证函数。

求教这个矩阵方程如何数值求解得到本征能量和本征函数?

具体算法? 谢谢

11.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhuisfp

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
husthc: 金币+20, 有帮助 2012-08-29 22:01:58
csgt0: 金币+1, 多谢应助 2012-08-29 23:11:55
MATLAB中
[e, v]=eigs(A)
就可以解得本征值本征矢
2楼2012-08-29 17:07:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

husthc

木虫 (正式写手)

引用回帖:
2楼: Originally posted by liuhuisfp at 2012-08-29 17:07:23
MATLAB中
=eigs(A)
就可以解得本征值本征矢

刚刚看了一下help,eigs(A,B)好像确实可以解这个广义本征值问题。另外,我想请问一下,matlab中eig和eigs两个函数有什么本质的区别吗?我看了help感觉没看出什么区别
3楼2012-08-29 20:27:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhuisfp

金虫 (小有名气)

引用回帖:
3楼: Originally posted by husthc at 2012-08-29 20:27:42
刚刚看了一下help,eigs(A,B)好像确实可以解这个广义本征值问题。另外,我想请问一下,matlab中eig和eigs两个函数有什么本质的区别吗?我看了help感觉没看出什么区别...

这个偶也不是很了解哦
4楼2012-08-29 20:36:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kngshuai

铁虫 (小有名气)

【答案】应助回帖


csgt0: 金币+1, 欢迎交流 2012-10-13 14:07:16
FORtran
LAPACK程序报  dsbgv
5楼2012-10-13 00:19:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 husthc 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 19:07 by 1js5n6upxm
[公派出国] +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 19:05 by 1js5n6upxm
[硕博家园] O(∩_∩)O +5 jk31rb6vyz 2026-06-12 7/350 2026-06-13 18:45 by lkmi35jdzt
[公派出国] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:33 by lkmi35jdzt
[考研] T^T +6 m584r2z4j3 2026-06-11 6/300 2026-06-13 18:29 by lkmi35jdzt
[文学芳草园] 咸菜 +3 myrtle 2026-06-10 3/150 2026-06-13 18:29 by ahsoarli
[考博] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:25 by lkmi35jdzt
[考研] o>_<o +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 18:22 by lkmi35jdzt
[考博] T^T +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 14:45 by 80j16ujcu8
[考研] O(∩_∩)O +3 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 14:42 by 80j16ujcu8
[考博] T^T +4 jk31rb6vyz 2026-06-12 4/200 2026-06-13 12:09 by 80j16ujcu8
[考博] T^T +3 m584r2z4j3 2026-06-11 3/150 2026-06-13 12:04 by 80j16ujcu8
[找工作] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 08:53 by t4hz8376wy
[考博] →_→ +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 08:49 by t4hz8376wy
[找工作] *囍 +3 m584r2z4j3 2026-06-12 3/150 2026-06-13 08:29 by t4hz8376wy
[基金申请] 面上项目没有好文章就没希望了吗? +19 1234567wang 2026-06-07 21/1050 2026-06-12 21:24 by littlezl
[考博] [急招]吉林大学招收26年9月入学考核制博士生1名,截止今晚24点前 +3 hehejin 2026-06-10 4/200 2026-06-12 18:24 by 冻雨天
[教师之家] 重庆理工大学副校长遇刺身亡 传涉案副教授疑因积怨行凶 +6 yexuqing 2026-06-07 7/350 2026-06-12 16:46 by wwwkkk83
[基金申请] 前几天时间戳更新了 +8 vito刘 2026-06-11 13/650 2026-06-12 14:17 by BlakeReary
[考博] 还有课题组有博士名额吗 +5 冻雨天 2026-06-07 6/300 2026-06-08 20:59 by RCS_tyro
信息提示
请填处理意见