質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.47%
Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Q&A

解決済

1回答

1501閲覧

固有ベクトルが正しく求まらない

Hirot

総合スコア4

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

0グッド

0クリップ

投稿2021/06/06 17:04

linalg.eig()で固有ベクトルを求めたいのですが正しく求まりません。なにか問題があるのでしょうか。
H=numpy.zeros((4,4),complex)
H[0,2]=1
H[1,3]=1j
H[2,0]=-1
H[3,1]=-1j
eigval, eigvec= linalg.eig(H)
print(eigvec[0])
このようなコードでHの固有ベクトル(例えばゼロ番目)を表示すると下のベクトルが返ってきました。
[0.70710678+0.j 0.70710678+0.j 0. +0.j 0. +0.j]
このベクトルはHの固有ベクトルになっていないと思うのですがなにがおかしいんでしょうか。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

取り出す次元が違っています

python

1import numpy as np 2H=np.zeros((4,4),complex) 3H[0,2]=1 4H[1,3]=1j 5H[2,0]=-1 6H[3,1]=-1j 7eigval, eigvec= np.linalg.eig(H) 8 9eigval_0 = eigval[0] 10eigvec_0 = eigvec[:, 0] 11print(H @ eigvec_0) 12print(eigval_0 * eigvec_0)

実行結果

python

1>>> print(H @ eigvec_0) 2[ 0. +0.70710678j 0. +0.j -0.70710678+0.j 3 0. +0.j ] 4>>> print(eigval_0 * eigvec_0) 5[ 0. +0.70710678j 0. +0.j -0.70710678+0.j 6 0. +0.j ]

投稿2021/06/06 22:57

ppaul

総合スコア24666

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

Hirot

2021/06/07 07:25

助かりました。ありがとうございます。
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.47%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問