Sylvesterovo kritérium

Z testwiki
Skočit na navigaci Skočit na vyhledávání

Šablona:Podrobně

Symetrická matice je pozitivně definitní, právě když ona i všechny její vedoucí hlavní podmatice mají kladný determinant

Sylvesterovo kritérium pozitivní definitnosti, pojmenované po Jamesi Sylvesterovi, uvádí nutnou a postačující podmínku pro rozhodnutí, zda komplexní hermitovská matice je pozitivně definitní.

Kritérium je možné zobecnit pro určení definitnosti matic, tj. pro rozpoznání pozitivně semidefinitních, negativně definitních, negativně semidefinitních i indefinitních matic.

Znění kritéria

Hermitovská matice je pozitivně definitní, právě když má všechny vedoucí hlavní subdeterminanty kladné.

Konkrétně, pro matici

𝑨=(a11a12a1na21a22a2nan1an2ann)

mají být kladné determinanty:

|a11|,|a11a12a21a22|,,|a11a12a1na21a22a2nan1an2ann|=det𝑨.

Kritérium platí v nezměněné formě pro reálné symetrické matice, protože tyto matice jsou speciálním případem hermitovských matic.

Použitím vhodné permutace na řádky i sloupce matice 𝑨 lze ukázat, 𝑨 je pozitivně definitní, právě když subdeterminanty v libovolné posloupnosti n do sebe vnořených hlavních podmatic matice 𝑨 jsou všechny kladné.

Ukázky

Reálná symetrická matice

𝑨=(210132024)

je pozitivně definitní, protože

|a11|=|2|=2>0,|a11a12a21a22|=|2113|=5>0 a det𝑨=|210132024|=12>0.

Namísto všech vedoucích hlavních podmatic by bylo možné použít i jinou posloupnost tří do sebe vnořených hlavních podmatic, např.

|a33|=|4|=4>0,|a11a13a31a33|=|2004|=8>0 a det𝑨=|210132024|=12>0.

Naopak matice

𝑩=(210132024)

není pozitivně definitní, protože alespoň jeden z determinantů vedoucích hlavních podmatic není kladný:

|b11|=|2|=2>0,|b11b12b21b22|=|2113|=7<0 a det𝑩=|210132024|=20>0

Matice 𝑩 zároveň ilustruje nutnost vzít do sebe vnořené hlavní podmatice. Kdyby namísto uvedené matice řádu 2 byla použita jiná hlavní podmatice, jejíž determinant je:

|b22b23b32b33|=|3224|=16>0,

byly by všechny determinanty tří hlavních podmatic různých řádů kladné, ale přesto matice 𝑩 není pozitivně definitní.

Zobecnění

Pozitivně semidefinitní matice

Pozitivně semidefinitní matice lze charakterizovat analogicky:

Hermitovská matice je pozitivně semidefinitní, právě když má všechny hlavní subdeterminanty nezáporné.

Oproti pozitivně definitním maticím nestačí uvážit pouze n hlavních vedoucích subdeterminantů, ale je třeba otestovat všech 2n1 vedoucích subdeterminantů.

Ukázky

Reálná symetrická matice

𝑨=(100042021)

je pozitivně semidefinitní, protože

|a11|=|a33|=|1|=10, |a22|=|4|=40, |a11a12a21a22|=|1004|=40, |a11a13a31a33|=|1001|=10, |a22a23a32a33|=|4221|=0, a det𝑨=|100042021|=0.

Reálná symetrická matice.

𝑩=(001010100)

má hodnoty vedoucích hlavních subdeterminantů nezáporné:

|0|=0,|0001|=0 a det𝑩=|001010100|=1>0

Tento výpočet ovšem není dostatečný pro test pozitivní semidefinitnosti matice 𝑩. Pro vektor 𝒙=(0,2,0)T platí: 𝒙T𝑩𝒙=4<0, čili matice 𝑩 není pozitivně semidefinitní.

Negativně (semi)definitní matice

Hermitovská matice je negativně definitní právě když znaménko každého vedoucího hlavního subdeterminantu řádu k je (1)k.

U negativně semidefinitních matic je třeba otestovat všech 2n1 vedoucích subdeterminantů, jejichž znaménko musí být buď 0 nebo (1)k, kde k je řád subdeterminantu.

V ostatních případech je matice indefinitní.

Idea důkazu

Dopřednou implikaci lze dokázat například následujícím způsobem: Je-li hermitovská 𝑨 řádu n pozitivně definitní a 𝑩 je její vedoucí hlavní podmatice řádu kn, potom je 𝑩 také pozitivně definitní, protože libovolný nenulový vektor 𝒗k lze doplnit nk nulami na nenulový vektor 𝒖n, pro nějž pak platí:

𝒗H𝑩𝒗=𝒖H𝑨𝒖>0

Pozitivně definitní matice mají kladný determinant, kupříkladu protože mají všechna vlastní čísla kladná anebo protože mají své odmocniny (viz. též Choleského rozklad).

Pro obrácenou implikaci je možné využít algoritmus pro výpočet Choleského rozkladu 𝑨=𝑳𝑳H. Ten by selhal jen v případě, kdy by některý prvek lkk na diagonále byl dán odmocninou z nekladného čísla s. V tuto chvíli vypočtený rozklad by měl splňovat 𝑨k=𝑳k𝑳kH i pro vedoucí hlavní podmatice matic 𝑨 a 𝑳 řádu k. Pokud by s0, nastal by v takovém případě spor s předpokladem, že matice 𝑨k má kladný determinant, protože jej lze podle věty o determinantu součinu matic vyjádřit výrazem:

det𝑨k=si=1k1lii20

Kromě uvedeného argumentu lze podat i přímý, ale pracný důkaz matematickou indukcí, jenž je založen na výpočtu signatury kvadratické formy Gaussovou eliminací upravenou pro řádky i sloupce.

Numerické záležitosti

Při výpočtu všech determinantů zároveň Gaussovou eliminací, při níž je zakázáno měnit pořadí řádků, lze spočítat všechny determinanty v čase O(n3), což je asymptoticky stejně rychlé jako test pozitivní definitnosti pomocí Choleského rozkladu.

Pokud by byl každý determinant počítán zvlášť, časová složitost vychází O(n4) pro test pozitivně definitních matic.

Složitost testu pozitivně semidefinitních matic prostřednictvím Sylvestreova kritéria je O(2nn3), a proto je vhodnější použít jiné postupy.

Odkazy

Reference

Šablona:Překlad

Literatura

Související články

Šablona:Autoritní data Šablona:Portály