用numpy对角化酉矩阵不会产生正交本征向量

dgiusagp  于 2022-11-29  发布在  其他
关注(0)|答案(1)|浏览(119)

我试图用numpy对角化一个酉矩阵,特别是numpy.linalg.eig函数。由于矩阵是酉矩阵,它的特征向量应该形成一个标准正交基。然而,情况似乎并非如此:

import numpy as np
from qiskit.circuit.library import QFT
from qiskit.quantum_info import Operator

op = Operator(QFT(num_qubits=4, do_swaps=True)).data
eigvals, eigvecs = np.linalg.eig(op) # Compute the eigenvalues
eigvecs = eigvecs.T # Since the eigenvectors are arranged in columns
lambda_i = eigvals[2]
lambda_j = eigvals[-1]
v_i = eigvecs[2].reshape((-1, 1))
v_j = eigvecs[-1].reshape((-1, 1))
print(np.linalg.norm(op @ v_i - lambda_i * v_i)) # Should be close to 0 by definition, actually yields 6.706985734026871e-16, which is fine
print(np.linalg.norm(op @ v_j - lambda_j * v_j)) # Should be close to 0 by definition, actually yields 8.151878100248519e-16, which is fine
print(v_j.T.conj() @ v_i) # Should be close to zero since the basis is supposed to be orthonormal but actually yields array([[-0.15147621-0.06735767j]]), which is not fine
print(np.linalg.norm(op.T.conj() @ op - np.eye(16))) # Should be around zero if and only if op.data is unitary, actually yields 2.3337334181537826e-15, which is ine

我已经用qiskit库测试过了,但是我的问题纯粹是与数值相关的。如果需要,可以使用下面的. npy文件加载相同的矩阵:

b"\x93NUMPY\x01\x00v\x00{'descr': '<c16', 'fortran_order': False, 'shape': (16, 16), }                                                       \n\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00C\x8d2\xcfk\x90\xcd?`\xa9\xae\xa6\xe2}\xb8?\xcb;\x7ff\x9e\xa0\xc6?\xca;\x7ff\x9e\xa0\xc6?c\xa9\xae\xa6\xe2}\xb8?C\x8d2\xcfk\x90\xcd?\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?^\xa9\xae\xa6\xe2}\xb8\xbfC\x8d2\xcfk\x90\xcd?\xc9;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?A\x8d2\xcfk\x90\xcd\xbfc\xa9\xae\xa6\xe2}\xb8?\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00C\x8d2\xcfk\x90\xcd\xbf`\xa9\xae\xa6\xe2}\xb8\xbf\xcb;\x7ff\x9e\xa0\xc6\xbf\xca;\x7ff\x9e\xa0\xc6\xbfc\xa9\xae\xa6\xe2}\xb8\xbfC\x8d2\xcfk\x90\xcd\xbf\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf^\xa9\xae\xa6\xe2}\xb8?C\x8d2\xcfk\x90\xcd\xbf\xc9;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbfA\x8d2\xcfk\x90\xcd?c\xa9\xae\xa6\xe2}\xb8\xbf\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xcb;\x7ff\x9e\xa0\xc6?\xca;\x7ff\x9e\xa0\xc6?\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xca;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xcb;\x7ff\x9e\xa0\xc6\xbf\xca;\x7ff\x9e\xa0\xc6\xbf\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xca;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbf\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xcb;\x7ff\x9e\xa0\xc6?\xca;\x7ff\x9e\xa0\xc6?\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xca;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xcb;\x7ff\x9e\xa0\xc6\xbf\xca;\x7ff\x9e\xa0\xc6\xbf\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xca;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbf\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00a\xa9\xae\xa6\xe2}\xb8?C\x8d2\xcfk\x90\xcd?\xca;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?D\x8d2\xcfk\x90\xcd\xbf[\xa9\xae\xa6\xe2}\xb8\xbf\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbfC\x8d2\xcfk\x90\xcd?c\xa9\xae\xa6\xe2}\xb8\xbf\xcb;\x7ff\x9e\xa0\xc6?\xc9;\x7ff\x9e\xa0\xc6?Z\xa9\xae\xa6\xe2}\xb8\xbfC\x8d2\xcfk\x90\xcd?\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00a\xa9\xae\xa6\xe2}\xb8\xbfC\x8d2\xcfk\x90\xcd\xbf\xca;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbfD\x8d2\xcfk\x90\xcd?[\xa9\xae\xa6\xe2}\xb8?\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?C\x8d2\xcfk\x90\xcd\xbfc\xa9\xae\xa6\xe2}\xb8?\xcb;\x7ff\x9e\xa0\xc6\xbf\xc9;\x7ff\x9e\xa0\xc6\xbfZ\xa9\xae\xa6\xe2}\xb8?C\x8d2\xcfk\x90\xcd\xbf\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00^\xa9\xae\xa6\xe2}\xb8\xbfC\x8d2\xcfk\x90\xcd?\xcb;\x7ff\x9e\xa0\xc6\xbf\xca;\x7ff\x9e\xa0\xc6\xbfC\x8d2\xcfk\x90\xcd?d\xa9\xae\xa6\xe2}\xb8\xbf\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?C\x8d2\xcfk\x90\xcd\xbf]\xa9\xae\xa6\xe2}\xb8\xbf\xc9;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbfd\xa9\xae\xa6\xe2}\xb8?A\x8d2\xcfk\x90\xcd?\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00^\xa9\xae\xa6\xe2}\xb8?C\x8d2\xcfk\x90\xcd\xbf\xcb;\x7ff\x9e\xa0\xc6?\xca;\x7ff\x9e\xa0\xc6?C\x8d2\xcfk\x90\xcd\xbfd\xa9\xae\xa6\xe2}\xb8?\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbfC\x8d2\xcfk\x90\xcd?]\xa9\xae\xa6\xe2}\xb8?\xc9;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?d\xa9\xae\xa6\xe2}\xb8\xbfA\x8d2\xcfk\x90\xcd\xbf\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xc9;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xcb;\x7ff\x9e\xa0\xc6?\xc9;\x7ff\x9e\xa0\xc6?\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xc9;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbf\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xcb;\x7ff\x9e\xa0\xc6\xbf\xc9;\x7ff\x9e\xa0\xc6\xbf\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xc9;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xcb;\x7ff\x9e\xa0\xc6?\xc9;\x7ff\x9e\xa0\xc6?\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xc9;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbf\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xcb;\x7ff\x9e\xa0\xc6\xbf\xc9;\x7ff\x9e\xa0\xc6\xbf\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00C\x8d2\xcfk\x90\xcd\xbfc\xa9\xae\xa6\xe2}\xb8?\xca;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbfZ\xa9\xae\xa6\xe2}\xb8\xbfD\x8d2\xcfk\x90\xcd?\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbfd\xa9\xae\xa6\xe2}\xb8?C\x8d2\xcfk\x90\xcd?\xcb;\x7ff\x9e\xa0\xc6\xbf\xc9;\x7ff\x9e\xa0\xc6\xbfC\x8d2\xcfk\x90\xcd?Y\xa9\xae\xa6\xe2}\xb8?\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00C\x8d2\xcfk\x90\xcd?c\xa9\xae\xa6\xe2}\xb8\xbf\xca;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?Z\xa9\xae\xa6\xe2}\xb8?D\x8d2\xcfk\x90\xcd\xbf\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?d\xa9\xae\xa6\xe2}\xb8\xbfC\x8d2\xcfk\x90\xcd\xbf\xcb;\x7ff\x9e\xa0\xc6?\xc9;\x7ff\x9e\xa0\xc6?C\x8d2\xcfk\x90\xcd\xbfY\xa9\xae\xa6\xe2}\xb8\xbf\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00C\x8d2\xcfk\x90\xcd\xbf`\xa9\xae\xa6\xe2}\xb8\xbf\xcb;\x7ff\x9e\xa0\xc6?\xca;\x7ff\x9e\xa0\xc6?c\xa9\xae\xa6\xe2}\xb8\xbfC\x8d2\xcfk\x90\xcd\xbf\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?^\xa9\xae\xa6\xe2}\xb8?C\x8d2\xcfk\x90\xcd\xbf\xc9;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?A\x8d2\xcfk\x90\xcd?c\xa9\xae\xa6\xe2}\xb8\xbf\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00C\x8d2\xcfk\x90\xcd?`\xa9\xae\xa6\xe2}\xb8?\xcb;\x7ff\x9e\xa0\xc6\xbf\xca;\x7ff\x9e\xa0\xc6\xbfc\xa9\xae\xa6\xe2}\xb8?C\x8d2\xcfk\x90\xcd?\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf^\xa9\xae\xa6\xe2}\xb8\xbfC\x8d2\xcfk\x90\xcd?\xc9;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbfA\x8d2\xcfk\x90\xcd\xbfc\xa9\xae\xa6\xe2}\xb8?\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xcb;\x7ff\x9e\xa0\xc6\xbf\xca;\x7ff\x9e\xa0\xc6\xbf\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xca;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbf\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xcb;\x7ff\x9e\xa0\xc6?\xca;\x7ff\x9e\xa0\xc6?\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xca;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xcb;\x7ff\x9e\xa0\xc6\xbf\xca;\x7ff\x9e\xa0\xc6\xbf\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xca;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbf\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xcb;\x7ff\x9e\xa0\xc6?\xca;\x7ff\x9e\xa0\xc6?\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xca;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00a\xa9\xae\xa6\xe2}\xb8\xbfC\x8d2\xcfk\x90\xcd\xbf\xca;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?D\x8d2\xcfk\x90\xcd?[\xa9\xae\xa6\xe2}\xb8?\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbfC\x8d2\xcfk\x90\xcd\xbfc\xa9\xae\xa6\xe2}\xb8?\xcb;\x7ff\x9e\xa0\xc6?\xc9;\x7ff\x9e\xa0\xc6?Z\xa9\xae\xa6\xe2}\xb8?C\x8d2\xcfk\x90\xcd\xbf\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00a\xa9\xae\xa6\xe2}\xb8?C\x8d2\xcfk\x90\xcd?\xca;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbfD\x8d2\xcfk\x90\xcd\xbf[\xa9\xae\xa6\xe2}\xb8\xbf\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?C\x8d2\xcfk\x90\xcd?c\xa9\xae\xa6\xe2}\xb8\xbf\xcb;\x7ff\x9e\xa0\xc6\xbf\xc9;\x7ff\x9e\xa0\xc6\xbfZ\xa9\xae\xa6\xe2}\xb8\xbfC\x8d2\xcfk\x90\xcd?\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00^\xa9\xae\xa6\xe2}\xb8?C\x8d2\xcfk\x90\xcd\xbf\xcb;\x7ff\x9e\xa0\xc6\xbf\xca;\x7ff\x9e\xa0\xc6\xbfC\x8d2\xcfk\x90\xcd\xbfd\xa9\xae\xa6\xe2}\xb8?\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?C\x8d2\xcfk\x90\xcd?]\xa9\xae\xa6\xe2}\xb8?\xc9;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbfd\xa9\xae\xa6\xe2}\xb8\xbfA\x8d2\xcfk\x90\xcd\xbf\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00^\xa9\xae\xa6\xe2}\xb8\xbfC\x8d2\xcfk\x90\xcd?\xcb;\x7ff\x9e\xa0\xc6?\xca;\x7ff\x9e\xa0\xc6?C\x8d2\xcfk\x90\xcd?d\xa9\xae\xa6\xe2}\xb8\xbf\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbfC\x8d2\xcfk\x90\xcd\xbf]\xa9\xae\xa6\xe2}\xb8\xbf\xc9;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?d\xa9\xae\xa6\xe2}\xb8?A\x8d2\xcfk\x90\xcd?\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xc9;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbf\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xcb;\x7ff\x9e\xa0\xc6\xbf\xc9;\x7ff\x9e\xa0\xc6\xbf\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xc9;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xcb;\x7ff\x9e\xa0\xc6?\xc9;\x7ff\x9e\xa0\xc6?\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00\xc9;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbf\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbf\xcb;\x7ff\x9e\xa0\xc6\xbf\xc9;\x7ff\x9e\xa0\xc6\xbf\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00\xc9;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?\xcb;\x7ff\x9e\xa0\xc6?\xc9;\x7ff\x9e\xa0\xc6?\xfd\xff\xff\xff\xff\xff\xcf?\x00\x00\x00\x00\x00\x00\x00\x00C\x8d2\xcfk\x90\xcd?c\xa9\xae\xa6\xe2}\xb8\xbf\xca;\x7ff\x9e\xa0\xc6?\xcb;\x7ff\x9e\xa0\xc6\xbfZ\xa9\xae\xa6\xe2}\xb8?D\x8d2\xcfk\x90\xcd\xbf\x05\\\x143&\xa6q\xbc\xfd\xff\xff\xff\xff\xff\xcf\xbfd\xa9\xae\xa6\xe2}\xb8\xbfC\x8d2\xcfk\x90\xcd\xbf\xcb;\x7ff\x9e\xa0\xc6\xbf\xc9;\x7ff\x9e\xa0\xc6\xbfC\x8d2\xcfk\x90\xcd\xbfY\xa9\xae\xa6\xe2}\xb8\xbf\xfd\xff\xff\xff\xff\xff\xcf\xbf\x00\x00\x00\x00\x00\x00\x00\x00C\x8d2\xcfk\x90\xcd\xbfc\xa9\xae\xa6\xe2}\xb8?\xca;\x7ff\x9e\xa0\xc6\xbf\xcb;\x7ff\x9e\xa0\xc6?Z\xa9\xae\xa6\xe2}\xb8\xbfD\x8d2\xcfk\x90\xcd?\x05\\\x143&\xa6q<\xfd\xff\xff\xff\xff\xff\xcf?d\xa9\xae\xa6\xe2}\xb8?C\x8d2\xcfk\x90\xcd?\xcb;\x7ff\x9e\xa0\xc6?\xc9;\x7ff\x9e\xa0\xc6?C\x8d2\xcfk\x90\xcd?Y\xa9\xae\xa6\xe2}\xb8?"

这种行为似乎只发生在纯复数(即实部为零)的特征值上,而对于其他大多数特征向量,标积的结果确实为0,但我可以在其他特征向量上重复这种行为。
我不认为这是由于某种数值近似,因为最终的标量积相当大。
我不认为这是因为我忘记了一些共轭,作为一个健全的检查,我已经确保v_iv_j确实是op的特征向量,这意味着它们是正确的,也意味着它们应该是正交的。

dm7nw8vv

dm7nw8vv1#

推理中的错误来自这一部分:
作为一个健全的检查,我已经确保v_i和v_j确实是op的特征向量,这意味着它们是正确的,也意味着它们应该是正交的。
如果v_i和v_j与相同的特征值lambda相关联,则(v_i+v_j)/sqrt(2)是与lambda相关联的特征向量,并且其不与v_i或v_j正交。
为了确保特征向量正交,只需执行:

eigvecs = np.linalg.qr(eigvecs)[0]

eigvecs现在将包含op的标准正交本征向量。

相关问题