Apakah determinan kecil menyiratkan pengondisian matriks?

29

det(A)0

Apakah kebalikannya juga benar? Apakah matriks yang dikondisikan buruk memiliki determinan yang hampir nol?

Ini adalah sesuatu yang saya coba di 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
Pemeriksaan resmi
sumber
1
Penentu menunjukkan apakah sebuah matriks adalah reguler atau tunggal. Itu tidak menunjukkan apakah itu terkondisi dengan baik atau buruk.
Allan P. Engsig-Karup
5
Besarnya penentu tidak dapat mencerminkan kondisi buruk: tetapi . κ(A)=κ(A1)det(A1)=(detA)1
faleichik
Haruskah ada atau suatu tempat?
Pemeriksaan
3
Jika Anda tertarik untuk belajar lebih banyak tentang efek floating-point matematika pada matriks spektrum, Anda harus memeriksa buku Nick Trefethen ini: Spectra dan Pseudospectra: The Perilaku tidak normal Matriks dan Operator dan Pseudospectra Gateway .
Aron Ahmadia

Jawaban:

38

Ini adalah besarnya jumlah kondisi yang mengukur kedekatan dengan singularitas, bukan kelenturan penentu.κ(A)

Sebagai contoh, matriks diagonal memiliki determinan kecil, tetapi dikondisikan dengan baik.1050I

Di sisi lain, perhatikan keluarga berikut dari matriks segitiga atas, karena Alexander Ostrowski (dan juga dipelajari oleh Jim Wilkinson):

U=(122121)

Determinan dari matriks U selalu 1 , namun rasio terbesar dengan nilai singular terkecil (yaitu jumlah kondisi 2-norma κ 2 ( U ) = σ 1n×nU1 ) ditunjukkan oleh Ostrowski sama dengancot2πκ2(U)=σ1σn , yang terlihat meningkat untuk meningkatkann.cot2π4nn

JM
sumber
1
@Nunoxic: pasti tidak; sebelum saya memulai ke detail, apakah Anda sudah terbiasa dengan dekomposisi nilai singular?
JM
2
ΣΣκ
2
logbκbκ1013
1
Ya, tapi itu bukan metode yang disarankan untuk menentukan nomor kondisi (penjelasan untuk pertanyaan lainnya). Saya kira Anda tahu bagaimana membalikkan matriks diagonal, bukan?
JM
2
"Regd. Kehilangan digit, bisakah kamu memberi saya referensi untuk ini?" - Saya bisa, tetapi ini benar-benar salah satu hal yang harus Anda coba sendiri di lingkungan komputasi untuk penguatan.
JM
17

det(kA)=kndetA

Dengan demikian tidak pernah menggunakan determinan untuk menilai kondisi atau kedekatan dengan singularitas.

Di sisi lain, untuk hampir semua masalah numerik yang diposisikan dengan baik, kondisi ini terkait erat dengan jarak ke singularitas, dalam arti gangguan relatif terkecil yang diperlukan untuk membuat masalah tersebut menjadi masalah. Secara khusus, ini berlaku untuk sistem linear.

Arnold Neumaier
sumber