ある行列の固有ベクトル(numpy.ndarray
)を計算して転置したものがあります(各列が各固有値に対応しています)
その固有ベクトルの中から0
でない最も低い値をt
個抜き出してlist
に保存したいのですが、方法がわからずに困っています。
例として以下のような行列を用意しました。
import numpy as np
from numpy import linalg as LA
A = np.array([[3, 2, 0], [2, 9, 2], [0, 2, 7]])
w, v = LA.eig(L)
eigenvector = v.T
例の行列では固有ベクトルが0ベクトルになっていないですが、0ベクトルがある行列にも対応できように0をt
個分list
に加えるようにもしたいです。
よろしくお願いします。
np.sort(eigenvector[np.nonzero(eigenvector)])[:t]
が考えられます。np.apply_along_axis(lambda x: x[np.nonzero(x)][:t], 0, np.sort(eigenvector, axis=0)).T
でしょうか?(結果はリストではなくnumpy.ndarray
型になりますけれども)ValueError: could not broadcast input array from shape (2) into shape (1)
というエラーが出ました。明確な違いは非対称行列を扱ったくらいなのですが、対処法などはありますか?なお、非対称行列ですが固有値は実数です。n - t
よりも多くなっている事かと思います。0
をパディングするのであれば、np.apply_along_axis(lambda x: np.r_[(c:=x[np.nonzero(x)])[:t], np.zeros(max(0, t-c.shape[0]))], 0, np.sort(eigenvector, axis=0)).T
となりますが、:=
演算子(Walrus operator)を使用していますので Python 3.8 以降で実行する必要があります。