私の目的は、行列の固有ベクトルを見つけることです。Matlabには、[V,D] = eig(M)次を使用して行列の固有ベクトルを取得するためのがあります[V,D] = eig(M)。または、WebサイトWolframAlphaを使用して結果を再確認しました。

と10X10呼ばれる行列がありMます：

D：

V：

次の結果が得られました。

1. 元のマトリックス：

1. WolframAlphaの結果：

2. Matlab Eigの結果：

D（固有値）

V（固有ベクトル）

固有ベクトルのさまざまな解を得ることが可能ですか、それとも独自の答えである必要があります。この概念を明確にすることに興味があります。

こんにちは、みんな！！！

次の形式の一般化された固有分解を計算したい：

Lf=ラムダAf

scipy.sparse.linalg.eigs関数を使用しますが、次のエラーが発生します：

/usr/local/lib/python2.7/dist-packages/scipy/linalg/decomp_lu.py:61：RuntimeWarning：対角線の数65は正確にゼロです。特異行列。RuntimeWarning）**DLASCLパラメーター番号4への入力時に不正な値がありました

対角行列、正半定値（PSD）行列、数値K（最初のK個の固有値）の3つの引数を渡します。Matlabのeigs関数は、同じ入力パラメーターを使用してうまく機能しますが、私が理解しているように、SciPyでは、PSDで計算するために、シグマパラメーターも指定する必要があります。

だから、私の質問は、MatLabのようにシグマパラメータの設定を回避する方法はありますか、そうでない場合は、シグマ値を取得する方法はありますか？

アドバイスやヒントをお待ちしております...よろしくお願いします！

次の API 関数は、OpenCV ライブラリの PCA のソース コードです。

この関数は Eigenobjects.cpp からのものです

そして、私は次のコードを理解していません:

PCA は共分散行列とその固有値と固有ベクトルを計算することであることがわかっていますが、関数の一部が内部関数 icvCalcCovarMatrixEx_8u32fR() を呼び出して固有値と固有ベクトルを計算すると、「上のコードの目的は何ですか?」 [その他の部分]"? [その他の部分] は固有ベクトルを計算することだと思いますが、関数 icvCalcCovarMatrixEx_8u32fR() がそれを行ったので、このコードの目的は何かを答えるのを手伝ってくれる人はいますか? ありがとう！！！

こんにちは、私はJamaライブラリを持っていますが、このライブラリはDouble数値でしか動作しません..そして非常に遅いです。Androidアプリの場合..そして最後に、eig decompの高精度は必要ありません..フロートnumを持つJAva libarayはありますか.jamaとの構文が似ていますか? なぜなら...私は440行のコードを書き直したくありません。ありがとう。eig。転置、逆変換などの基本的な線形代数演算..

または、同じ Java ライブラリ固有値が存在します。スレッドで？

Matlab を使用してサンプル共分散行列の固有値の統計的分散を調査しようとしています。明確にするために、各サンプル共分散行列は有限数のベクトル スナップショット (ランダムなホワイト ガウス ノイズに悩まされている) から構築されます。次に、多数の試行にわたって、多数のそのような行列が生成され、固有値の理論的統計を推定するために固有分解されます。

いくつかの情報源 ([1, Eq.3] および [2, Eq.11] などを参照) によると、各サンプルの固有値の分散は、理論上の固有値の 2 乗を使用したベクトル スナップショットの数で割った値に等しいはずです。各共分散行列について。ただし、Matlab から取得した結果は、それほど近いものではありません。

これは私のコードの問題ですか? マトラブで？（同様の問題に取り組むのにこれほど苦労したことはありません）。

非常に簡単な例を次に示します。

参考文献:

[1]フリードランダー、B。ワイス、AJ。、「標本共分散行列の固有ベクトルの二次統計について」、Signal Processing、IEEE Transactions on 、vol.46、no.11、pp.3136-3139、1998 年 11 月 [2] Kaveh, M.; Barabell、A.; 、「雑音中の平面波を解決する際の MUSIC および最小ノルム アルゴリズムの統計的性能」、Acoustics, Speech and Signal Processing、IEEE Transactions on 、vol.34、no.2、pp.331-341、1986 年 4 月

私はC++を初めて使用し、LAPACKを使用して無限帯行列の固有値を見つける方法を理解しようとしています（非調和振動子の問題）。値を確認したところ、すべて一致しているので、マトリックスを正しく計算していることがわかります。ただし、値をサブルーチンに正しく渡しているのか、または返される固有値が期待どおりではないために何かが混同されているのかはわかりません。これを計算するためにdsbtrdサブルーチンを使用しています。そのためのマニュアルは次のとおりです。http：//www.netlib.org/lapack/explore-html/d0/d62/dsbtrd_8f.html

私がどこで間違っているのかについてのアイデアはありますか？

以下のような未知の定数がたくさんある行列があります。

お気づきかもしれませんが、対角線に関して対称であるため、対角線の値はすべて正です。すべての定数が0より大きい。

matlabの固有値についてこれを解決したいと思います。どうすればこれを行うことができますか？a、b、c、d、eの値がわかりません。私はこのようなことをしたいと思います：

ただし、eig関数は関数ハンドルを受け入れません。

Matlabで関数を使用します：

VとDは常に昇順でソートされていることがわかります。それはいつもそれが好きですか、それとも私が値を取得した後にそれらをソートする必要がVありDますか？

これは、IronPythonユーザーにとって頻繁に発生する問題である可能性があります。しかし、私はPythonを初めて使用します。

IronPython2.7を使用して約50x50行列の固有値と固有ベクトルを計算したいNumpyとScipyを使用する可能性を検討しましたが、IronPythonではサポートされていません。

私が必要とするものを達成するためのより良い方法はありますか？

行列が与えられた場合、固有値と対応する固有ベクトルに興味があります。

Jamaを使用して、固有値と固有ベクトルを取得できますが、2つの間の相関関係は定義されていません。各固有ベクトルを対応する固有値にマップします。

その方法を教えていただけますか？自分で実装しようとしたのですが、やっかいになりました。

ありがとう ：）

私は許可された答えを探していますが、今のところ、私が行った実験と観察によれば、固有ベクトルとevigenValuesは対応しているようです。