Un determinante minuscolo implica mal condizionamento di una matrice?


29

det(A)0

È vero anche il contrario? Una matrice mal condizionata ha un determinante quasi zero?

Ecco qualcosa che ho provato in 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
Il determinante mostra se una matrice è regolare o singolare. Non mostra se è ben condizionato o mal condizionato.
Allan P. Engsig-Karup,

5
L'entità del determinante non può riflettere il mal condizionamento: ma . κ(A)=κ(A1)det(A1)=(detA)1
Faleichik,

Dovrebbe esserci un o da qualche parte?
Inchiesta il

3
Se sei interessato a saperne di più sugli effetti della matematica in virgola mobile sugli spettri a matrice, dai un'occhiata al libro di Nick Trefethen: Spectra e Pseudospectra: il comportamento delle matrici e degli operatori non anormali e il Pseudospectra Gateway .
Aron Ahmadia,

Risposte:


38

κ(UN)

10-50io

D' altro canto , considera la seguente famiglia di matrici triangolari superiori quadrate, dovuta ad Alexander Ostrowski (e anche studiato da Jim Wilkinson):

U=(122121)

n×nU1κ2(U)=σ1σnculla2π4nn


1
@Nunoxic: sicuramente no; prima di lanciarmi nei dettagli, hai già familiarità con la scomposizione del valore singolare?
JM,

2
ΣΣκ

2
logBκBκ1013

1
Sì, ma questo non è il metodo raccomandato per determinare il numero della condizione (una spiegazione per un'altra domanda). Presumo tu sappia come invertire una matrice diagonale, no?
JM,

2
"Regd. La perdita di cifre, potresti darmi un riferimento per questo?" - Potrei, ma questa è davvero una di quelle cose che dovresti sperimentare da solo in un ambiente informatico per il rinforzo.
JM,

17

det(KUN)=KndetUN

Pertanto, non utilizzare mai il determinante per valutare la condizione o la vicinanza alla singolarità.

D'altra parte, per quasi tutti i problemi numerici ben posati, la condizione è strettamente correlata alla distanza dalla singolarità, nel senso della più piccola perturbazione relativa necessaria per rendere il problema mal posto. In particolare, questo vale per i sistemi lineari.

Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.