应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求解本征态问题Au(x)=EBu(x)
51/1返回列表
查看: 683  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

husthc

木虫 (正式写手)

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

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


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

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

具体算法？ 谢谢

11.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-08-29 10:58:44
已阅   回复此楼   关注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 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭