正の半定行列の平方根はユニークな結果ですか?


8

私は観測の 時系列を分散共分散構造とランダム系列に分解しようとしています。nvcn×nv

したがって、の自己相関関数から分散共分散行列を導出できます。これは、半正定値のテプリッツ行列になります。したがって、適切な行列を計算して、相関系列をランダムな信号に変換できます。 vc12

v=12vc

これは、MATLABのsqrt(m)関数を使用して実行できますが、分散共分散行列のコレスキー分解を検出し、これを使用して相関を誘導することもできます。ただし、sqrtm法とコレスキー法を使用したランダムシリーズの結果は異なります(ただし多少似ています)。

さまざまな行列の平方根を確認する方法を決定するためにいくつかのテキストを読み、固有値分解法などを調べてきました。特定の所定の条件下では一意のソリューションしかないことがわかりますが、これらの一意のソリューションはまだ多くのルーツの1つにすぎないと思いますか?

私の質問はこれです。ある特定の平方根が別の平方根よりも好ましいと主張する方法はありますか。そうでない場合、すべての可能な解を抽出して、すべての可能なランダム関数を取得できる方法はありますか?

回答:


7

行列に「平方根」とます。あれは、VAB

V=AA=BB.

簡単にするために、元の行列が可逆であると仮定します (これは、仮定の下で正定であることに相当します)。次に、、、およびそれらの転置も反転可能でなければなりません。VAB

I=V1V=V1AA=(V1A)A

左逆を示し、も反転可能であることを意味します。もちろん、同じ議論がにも当てはまります。これらの逆を利用して書き込みますAAB

(B1A)(B1A)=B1(AA)B1=B1(V)B1=B1(BB)B1=I I=I,

ことを示すである直交行列:であり、。このようなマトリックスの形態2つの平滑実次元のマニホールドの集合あるによって。幾何学的に、直交行列は、行列式の符号に応じて、回転または回転に続く反射に対応します。O=B1AOO=In(n1)/2Vnn

逆に、が平方根である場合、同様の(しかしより簡単な)計算により、も直交行列平方根であることがわかりますが可逆かどうかはここでは関係ありません。AVAOOA

直交行列(と等しくない)による乗算が、可逆行列の平方根を実際に変更することも簡単にわかります。結局のところ、すぐに意味します。これは、正定行列の平方根が直交行列と1対1で対応できることを示しています。IAO=AO=A1A=I

これは、正定行列の平方根が直交行列による乗算までしか決定されないことを示しています。 半定値の場合、状況はより複雑になりますが、少なくとも、直交行列を乗算すると平方根であるという特性が維持されます。

平方根に追加の基準を適用する場合は、一意の基準を特定するか、少なくともあいまいさを絞り込むことができる場合があります。それは特定の設定に依存します。


5
(+1)@hydrologist:whuberの回答の補足として:一意性につながる可能性のある1つの基準は、平方根自体が半正定であることを主張することです。一意性は、が正の半正定値であるという弱い条件の下で成り立ちます。「うまくいかない」ことができるか見るために有益な例は、の可能平方根を見ることである、でもちょうど対角のもの!:)AVI
枢機卿の

@cardinal:ご回答いただきありがとうございます。これは非常に役に立ち、大変感謝しています!
水文学者、2013年

@whuber:ご協力ありがとうございます。これは非常に便利です。
水文学者

1
これは平方根の単一の自由
kjetil b halvorsen '24年
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.