24小时热门版块排行榜    

查看: 749  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 函评 +3 布布和一二 2026-06-13 3/150 2026-06-13 23:35 by 芋圆香飘飘
[硕博家园] T^T +6 m584r2z4j3 2026-06-11 7/350 2026-06-13 23:10 by 3eaav2pea1
[找工作] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 19:07 by 1js5n6upxm
[公派出国] T^T +5 q1zg6x4lf3 2026-06-12 5/250 2026-06-13 18:53 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
[考博] 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
[考研] ╭(╯ε╰)╮ +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 14:22 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
[论文投稿] →_→ +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 08:13 by t4hz8376wy
[基金申请] 面上项目没有好文章就没希望了吗? +19 1234567wang 2026-06-07 21/1050 2026-06-12 21:24 by littlezl
[教师之家] 重庆理工大学副校长遇刺身亡 传涉案副教授疑因积怨行凶 +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
[教师之家] AASSDDFF +3 m584r2z4j3 2026-06-11 3/150 2026-06-12 11:28 by zs48m4lz8g
[考博] 博士申请 +3 张城斐 2026-06-07 3/150 2026-06-08 09:15 by A大朋同学
信息提示
请填处理意见