Antaa vuoronsa

Kokeneet kirjoittajat eivät ole vielä tarkistaneet sivun nykyistä versiota, ja se voi poiketa merkittävästi 7.11.2021 tarkistetusta versiosta . tarkastukset vaativat 2 muokkausta .

Annettu rotaatio - lineaarisessa algebrassa , lineaarinen operaattori vektorin kiertämiseksi tietyllä kulmalla .

Annettu matriisi [1] [2] [3]

Givens-matriisilla on seuraava muoto: ${\displaystyle G_{kl))$

$G_{kl}={\begin{bmatrix}1&\cdots &0&\cdots &0&\cdots &0\\\vdots &\ddots &\vdots &&\vdots &&\vdots \\0&\cdots &\cos {\ phi }&\cdots &-\sin {\phi }&\cdots &0\\\vdots &&\vdots &\ddots &\vdots &&\vdots \\0&\cdots &\sin {\phi }&\cdots &\ cos {\phi }&\cdots &0\\\vdots &&\vdots &&\vdots &\ddots &\vdots \\0&\cdots &0&\cdots &0&\cdots &1\end{bmatrix))$

Tämä matriisi eroaa identiteettimatriisista vain alimatriisin suhteen

$M(\phi )={\begin{bmatrix}\cos {\phi }&-\sin {\phi }\\\sin {\phi }&\cos {\phi }\end{bmatrix))$

sijaitsevat riveillä ja sarakkeilla, joissa on numerot ja . On ortogonaalinen. $k$ $l$

Jos vektori , annetaan , niin valitsee $a={\begin{bmatrix}a_{1}&\ldots &a_{n}\end{bmatrix}}^{T}\in \mathbb {R} ^{n}$ $s={\sqrt {a_{k}^{2}+a_{l}^{2))}\neq 0$

cos ⁡ ϕ = a k a k 2 + a l 2 {\displaystyle \cos {\phi }={\frac {a_{k}}{\sqrt {a_{k}^{2}+a_{l}^{2))))} $\cos {\phi }={\frac {a_{k}}{\sqrt {a_{k}^{2}+a_{l}^{2))))$ synti ⁡ ϕ = − a l a k 2 + a l 2 {\displaystyle \sin {\phi }={\frac {-a_{l}}{\sqrt {a_{k}^{2}+a_{l}^{2))))} $\sin {\phi }={\frac {-a_{l}}{\sqrt {a_{k}^{2}+a_{l}^{2))))$

voit asettaa vektorin :nnen komponentin nollaksi : $l$ $a$

[ cos ⁡ ϕ − synti ⁡ ϕ synti ⁡ ϕ cos ⁡ ϕ ] [ a k a l ] = [ cos ⁡ ϕ ⋅ a k − synti ⁡ ϕ ⋅ a l synti ⁡ ϕ ⋅ a k + cos ⁡ ϕ ⋅ a l ] = [ a k 2 + a l 2 a k 2 + a l 2 − a l ⋅ a k + a k ⋅ a l a k 2 + a l 2 ] = [ a k 2 + a l 2 0 ] {\displaystyle {\begin{bmatrix}\cos {\phi }&-\sin {\phi }\\\sin {\phi }&\cos {\phi }\end{bmatrix)){\begin{bmatrix} a_{k}\\a_{l}\end{bmatrix}}={\begin{bmatrix}\cos {\phi }\cdot a_{k}-\sin {\phi }\cdot a_{l}\\ \sin {\phi }\cdot a_{k}+\cos {\phi }\cdot a_{l}\end{bmatrix))={\begin{bmatrix}{\frac {a_{k}^{2} +a_{l}^{2}}{\sqrt {a_{k}^{2}+a_{l}^{2}}}}\\{\frac {-a_{l}\cdot a_{k }+a_{k}\cdot a_{l}}{\sqrt {a_{k}^{2}+a_{l}^{2}}}}\end{bmatrix}}={\begin{bmatrix} {\sqrt {a_{k}^{2}+a_{l}^{2}}}\\0\end{bmatrix}}}

{\begin{bmatrix}\cos {\phi }&-\sin {\phi }\\\sin {\phi }&\cos {\phi }\end{bmatrix)){\begin{bmatrix} a_{k}\\a_{l}\end{bmatrix}}={\begin{bmatrix}\cos {\phi }\cdot a_{k}-\sin {\phi }\cdot a_{l}\\ \sin {\phi }\cdot a_{k}+\cos {\phi }\cdot a_{l}\end{bmatrix))={\begin{bmatrix}{\frac {a_{k}^{2} +a_{l}^{2}}{\sqrt {a_{k}^{2}+a_{l}^{2}}}}\\{\frac {-a_{l}\cdot a_{k }+a_{k}\cdot a_{l}}{\sqrt {a_{k}^{2}+a_{l}^{2}}}}\end{bmatrix}}={\begin{bmatrix} {\sqrt {a_{k}^{2}+a_{l}^{2}}}\\0\end{bmatrix}}

Givensin rotaatioiden avulla voidaan laskea matriisien QR-hajoaminen ja piirtää hermiittiset matriisit kolmikulmaiseen muotoon .

Givens-matriisien käyttö tridiagonalisointiin

Haluamme pelkistää symmetrisen matriisin kolmikulmaiseen muotoon: $A={\begin{bmatrix}a_{11}&\cdots &a_{1p}&\cdots &a_{1q}&\cdots &a_{1n}\\\vdots &&\vdots &&\vdots &&\vdots \ \a_{p1}&\cdots &a_{pp}&\cdots &a_{pq}&\cdots &a_{pn}\\\vdots &&\vdots &&\vdots &&\vdots \\a_{q1}&\cdots &a_{ qp}&\cpisteet &a_{qq}&\cpisteet &a_{qn}\\\vpisteet &&\vpisteet &&\vpisteet &&\vpisteet \\a_{n1}&\cpisteet &a_{np}&\cpisteet &a_{nq}& \cdots &a_{nn}\end{bmatrix}}$

Missä . Sitten kerromme sen Givensin rotaatiomatriisilla: . on transponoitu matriisi. Tämä muuttaa vain elementit , ja $a_{pq}=a_{qp}$ $G'_{pq}(\theta )AG_{pq}(\theta )$ $G$ $a_{pp}$ ${\displaystyle a_{pq))$ ${\displaystyle a_{qq))$

${\displaystyle a'_{pp}=c^{2}a_{pp}+2csa_{pq}+s^{2}a_{qq))$

${\displaystyle a'_{pq}=sc(a_{qq}-a_{pp})+(c^{2}-s^{2})a_{pq))$

${\displaystyle a'_{qq}=s^{2}a_{pp}-2csa_{pq}+c^{2}a_{qq))$

Tässä alkuluku tarkoittaa elementtiä, joka tulee näkyviin kierron jälkeen. Valitaan kertoimet ja niin, että diagonaalista poikkeava elementti asetetaan nollaan ja suhde ja kanssa ja $c$ $s$ $c$ $s$ $\cos \phi$ $\sin \phi$

Sitten: $\phi =1/2\tan ^{-1}(2a_{pq}/(a_{pp}-a_{qq}))$

$c=\cos \phi$

$s=\sin\phi$

Tällaista kiertoa sovelletaan peräkkäin kaikkien ensimmäisen rivin elementtien nollaamiseksi kahta ensimmäistä lukuun ottamatta. Eli (1,2), (1,3), (1,4)...(1,n) Sitten rinnakkain toinen rivi (2,3),(2, 4)...(2 ,n )

C++ koodi:

for ( etumerkitön int i = 0 ; i < N -1 ; ++ i ) { for ( etumerkitön int j = i + 2 ; j < N ; ++ j ) { t = 2 * matr [ i ][ j ] / ( matr [ i ][ i ] - matr [ j ][ j ]); phi = 0,5 * atan ( t ); c = cos ( phi ); s = sin ( phi ); bii = c * c * matr [ i ][ i ] + 2 * c * s * matr [ i ][ j ] + s * s * matr [ j ][ j ]; bij = s * c * ( matr [ j ][ j ] - matr [ i ][ i ]) + matr [ i ][ j ] * ( c * c - s * s ); bjj = s * s * matr [ i ][ i ] + c * c * matr [ j ][ j ] - 2 * c * s * matr [ i ][ j ]; bji = bij ; matr [ i ][ i ] = bii ; matr [ i ][ j ] = bij ; matr [ j ][ i ] = bji ; matr [ j ][ j ] = bjj ; } }

Muistiinpanot

↑ Tyrtyshnikov E. E. Numeerisen analyysin menetelmät. - M. , 2006. - S. 73-74.
↑ Björck, Åke, 1934-. Numeeriset menetelmät vähintään neliösumma - tehtäville . - Philadelphia: SIAM, 1996. - S. 121-123. — xvii, 408 sivua, s. - ISBN 0-89871-360-9 , 978-0-89871-360-2.
↑ Demmel, James W. Sovellettu numeerinen lineaarinen algebra . - Philadelphia: Society for Industrial and Applied Mathematics, 1997. - S. 53-56. — xi, 419 sivua s. - ISBN 0-89871-389-7 , 978-0-89871-389-3, 0-89871-361-7, 978-0-89871-361-9.