問題タブ [qr-decomposition]

LAPACK のDGEQRF および SGEQRFは、QR 分解の Q 部分をパック形式で返します。開梱には手順が必要なようでO(k^3)(k 個の低ランク製品)、それほど簡単ではないようです。kさらに、逐次乗算の数値安定性は私にはわかりません。

LAPACK には Q をアンパックするためのサブルーチンが含まれていますか? そうでない場合は、どのようにすればよいですか?

問題: この問題では、As と呼ばれる行列のリストが与えられます。あなたの仕事は、それぞれの QR 因数分解を見つけることです。

qr_by_gram_schmidt の実装: この関数は行列 A を入力として取り、QR 分解を計算し、2 つの変数 Q と R を返します。ここで、A=QR であり、Q は直交し、R は対角線より下にゼロです。

A は、n≧m の n×m 行列です (つまり、列より行が多い)。

修正された Gram-Schmidt 手順を使用して、この関数を実装する必要があります。

入力：

As: 配列のリスト

出力：

Qs: qr_by_gram_schmidt によって出力された Q 行列のリスト (As と同じ順序)。形状が n×m の行列 A の場合、Q は形状が n×m でなければなりません。
Rs: qr_by_gram_schmidt が As と同じ順序で出力した R 行列のリスト。形状が n×m の行列 A の場合、R は形状が m×m でなければなりません

私は正しいと信じているQR因数分解のコードを書きました：

As の各行列の QR 分解を計算し、それらをその順序で格納するループを作成するにはどうすればよいですか?

編集: コードにもエラーがあります。デバッグにご協力いただければ幸いです。

ありがとう

QR 分解を学習しようとしていますが、従来の行列計算に頼らずに beta_hat の分散を取得する方法がわかりません。私はirisデータセットを使って練習しています.これが私がこれまでに持っているものです:

ご協力いただきありがとうございます！

コレスキー分解を行う方法と、TensorFlow で QR 法を使用して線形システムを解く方法があることを確認しましたが、TensorFlow で QR 分解を実装する方法が見つかりません。

TensorFlow で QR 分解を実行するにはどうすればよいですか?

Rで行列のQR因数分解を見つけるためのコードに取り組んでいます.

しかし、ここではH、世帯主の反射を表す行列をすべてのステップで計算していません。またA、すべてのステップで行列を計算していません。

としてH = I - 2 v v'、 を掛けるQと得られます

これで、この操作はすべてのステップで機能するはずです。ただし、最初のマトリックスHと2番目のマトリックスを考慮すると、H1これらのマトリックスは最初のマトリックスよりも小さくなります。それを避けるために、次のコード行を使用しました。

しかし、すべてのステップでゼロa_rの最初のエントリを持つ新しいベクトルを生成すると、コードがうまく機能する理由がわかりません。k

QR 分解を使用して重回帰を解くための関数を作成しようとしています。入力: y ベクトルと X 行列。出力: b、e、R^2。これまでのところ、私はこれを手に入れましたが、ひどく立ち往生しています。私はすべてを複雑にしすぎたと思います：

Eigenに行列を再直交化する組み込みの方法はありますか? 多くの回転を乗算すると、最終的に行列を再直交化する必要があります。SVD を使用するなどの標準的な手法があり、手動で記述して最適化するのに 1 日を費やすことができます。ただし、Eigen には組み込みの機能がほとんどないことを願っています。以前に関連する質問がありましたが、Eigen のコンテキストで明確な回答がありませんでした。