小さな行列式は、行列の悪条件を意味しますか?


29

正方の可逆行列があり、その行列式を取得し、そのを見つけた場合、det(A)0これは行列の条件が悪いことを意味しますか?

その逆も本当ですか?悪条件の行列の行列式はほぼゼロですか?

Octaveで試したものを次に示します。

a = rand(4,4);
det(a) %0.008
cond(a)%125
a(:,4) = 1*a(:,1) + 2*a(:,2) = 0.000000001*ones(4,1);
det(a)%1.8E-11
cond(a)%3.46E10

1
行列式は、行列が正則か特異かを示します。条件が良いか悪いかは表示されません。
アランP. Engsig-Karup

5
行列式の大きさは、悪条件を反映することはできません:κ(A)=κ(A1)が、det(A1)=(detA)1
-faleichik

どこかにまたは\ neが必要ですか?
調査

3
マトリックススペクトルに対する浮動小数点演算の影響について詳しく知りたい場合は、ニックトレフェテンの著書:Spectra and Pseudospectra:The Behavior of Nonnormal Matrices and Operators and the Pseudospectra Gatewayを参照してください
アロンアーマディア

回答:


38

行列式の小ささではなく、特異点への近さを測定するのは条件数 \ kappa(\ mathbf A)の大きさですκ(A)

たとえば、対角行列1050Iは小さな行列式がありますが、条件が整っています。

反対に、Alexander Ostrowski(およびJim Wilkinsonによる研究)により、次の正方形の上三角行列のファミリーを検討します。

U=(122121)

行列式行列常にが、最小特異値に対する最大の比(すなわち、2ノルムの条件数)に等しくなるようにOstrowskiによって示された増大させるために増加することがわかる、。n×nU1κ2(U)=σ1σncot2π4nn


1
@Nunoxic:確かにそうではありません。詳細を説明する前に、特異値分解にすでに精通していますか?
JM

2
とても良い。あなたが知る必要があるのはそれだけです。考え方は、コンディショニングに関する非常に重要な情報が集中しているということです。特に、その行列の対角線で最大値と最小値(対角要素が非負になるように分解が定義されていることを思い出してください)を探したいでしょう。最大対最小の対角要素の比率は、条件番号です。考慮すべき条件番号のサイズは、作業中のマシンによって異なりますΣΣκ
JM

2
...しかし、一般に、その行列で線形方程式を解くと、解で base桁を失うことになります。これは、条件番号の大まかな目安です。したがって、16桁のみで作業している場合、が問題の原因になります。logbκbκ1013
JM

1
はい。ただし、条件番号を決定するための推奨される方法ではありません(その説明は別の質問に対するものです)。対角行列を反転する方法を知っていると思いますか?
JM

2
「それでは、桁落ちです。これについての参考にしてください」-私はできましたが、これは本当に強化のためにコンピューティング環境で自分で実験すべきものの1つです。
JM

17

、行列式は、(条件数を変化させない)、単純な再スケーリングによって任意に大きく又は小さくすることができます。特に高次元では、2の無邪気な係数によるスケーリングでさえ、決定要因を膨大に変化させます。det(kA)=kndetA

したがって、条件または特異点への近さを評価するために決定要因を使用しないでください。

一方、ほとんどすべての適切な数値問題では、問題を不適切にするために必要な最小の相対的摂動という意味で、条件は特異点までの距離に密接に関連しています。特に、これは線形システムに当てはまります。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.