行列の 2 つの最小の固有ベクトルを取得しようとしています。
[v,c]=eigs(lap,2,'sm');
結果vは、約 66% の確率で「正しい」です。もちろん、私が正しいと言うとき、私が解決しようとしている問題に関して「正しく見える」ことを意味します。他の部分では、さまざまなベクトルを取得します。
私は eigs が数値ソルバーを使用していることを知っており、その最初の推測はランダムであるため、それが説明されています。私を悩ませているのは、matlab のドキュメントによると、停止する基準として使用される許容範囲がeps最初に設定されていることがわかり、増加しようとopts.maxit=10000000;しましたが、結果や実行時間には影響しないようです。最大反復回数に達する前に。
一貫した結果を得るにはどうすればよいですか? 計算時間的には問題ありません。
eigマトリックスは非常に大きくてまばらなので、 では作業できず、 でのみ作業できることに注意してくださいeigs