首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >问答首页 >相同hermitian矩阵的numpy's "linalg.eig()“和"linalg.eigh()”

相同hermitian矩阵的numpy's "linalg.eig()“和"linalg.eigh()”
EN

Stack Overflow用户
提问于 2022-11-23 23:18:06
回答 1查看 26关注 0票数 -1

这个问题是由于误解造成的。见下面的答案.

numpy.linalg方法eig()和eigh()似乎为相同的hermitian矩阵返回不同的特征向量。在这里,代码:

代码语言:javascript
运行
复制
import numpy as np

H = [[0.6 , -1j, 0], [1j, 0.4, 0], [0, 0, -1]]

evals, evects = np.linalg.eig(H)
print('\nOutput of the eig function')
for i in range(0,3):
    print('evect for eval=',evals[i],'\n',evects[i,0],'\n',evects[i][1],'\n',evects[i][2])
    
evals, evects = np.linalg.eigh(H)
print('\nOutput of the eigh function')
for i in range(0,3):
    print('evect for eval=',evals[i],'\n',evects[i,0],'\n',evects[i][1],'\n',evects[i][2])
EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2022-11-24 13:42:34

张贴这篇文章是为了帮助那些可能和我有过同样误解的人:

特征向量是两个函数的生成矩阵的。故障发生在原始代码中,从特征向量矩阵中提取行而不是列。正确的代码如下。

代码语言:javascript
运行
复制
H = [[0.6 , -1j, 0], [1j, 0.4, 0], [0, 0, -1]]

evals, evects = np.linalg.eig(H)
print('\nOutput of the eig function')
for i in range(0,3):
    print('evect for eval=',evals[i],'\n',evects.T[i,0],'\n',evects.T[i][1],'\n',evects.T[i][2])
    
evals, evects = np.linalg.eigh(H)
print('\nOutput of the eigh function')
for i in range(0,3):
    print('evect for eval=',evals[i],'\n',evects.T[i,0],'\n',evects.T[i][1],'\n',evects.T[i][2])
票数 0
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/74553962

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档