Jacobin menetelmä ominaisarvoille

Jacobin ominaisarvojen menetelmä  on iteratiivinen algoritmi todellisen symmetrisen matriisin ominaisarvojen ja ominaisvektorien laskemiseen . Nimetty Carl Gustav Jacob Jacobin mukaan, joka ehdotti tätä menetelmää vuonna 1846 [1] , vaikka menetelmä otettiin käyttöön vasta 1950-luvulla tietokoneiden myötä [2] .

Kuvaus

Olkoon  symmetrinen matriisi ja anna  olla rotaatiomatriisi . Sitten

on symmetrinen ja matriisimainen .

Lisäksi se sisältää seuraavat komponentit:

missä ja .

Koska  on ortogonaalinen matriisi, matriiseilla ja on samat Frobenius-normit (kaikkien komponenttien neliösumman juuret), ja voimme valita niin, että , ja tässä tapauksessa on suurempi diagonaalielementtien neliöiden summa:

Kun tämä lasketaan nollaan, saamme

Jos , niin

Optimaalisen vaikutuksen saavuttamiseksi on vaadittava, että se on suurin itseisarvoltaan poikkeava alkio, ns. peruselementti .

Jacobi-menetelmä ominaisarvoille pyörii, kunnes matriisi on lähes diagonaalinen. Sitten diagonaalin elementit approksimoivat matriisin ominaisarvoja .

Muistiinpanot

  1. Jacobi, CGJ Über ein leichtes Verfahren, die in der Theorie der Säkularstörungen vorkommenden Gleichungen numerisch aufzulösen  (saksa)  // Crelle's Journal . - 1846. - T. 30 . - S. 51-94 .
  2. Golub, G.H.; van der Vorst, HA Omaarvolaskenta 1900-luvulla  //  Journal of Computational and Applied Mathematics : päiväkirja. - 2000. - Voi. 123 , nro. 1-2 . - s. 35-65 . - doi : 10.1016/S0377-0427(00)00413-1 .