Методы вычислений. Часть I. Численные методы алгебры. Курцева К.П - 45 стр.

UptoLike

Составители: 

45
1
2
0 ... 0
0 ... 0
... ... ... ...
000
n
λ
λ
λ
⎛⎞
⎜⎟
⎜⎟
Λ=
⎜⎟
⎜⎟
⎜⎟
⎝⎠
,
00
00
1
00
iii
i
ee
i
λ
λ
⎛⎞
⎜⎟
⎜⎟
⎜⎟
Λ=Λ = =
⎜⎟
⎜⎟
⎜⎟
⎝⎠
MM
Собственные числа матрицы
A
i
λ
, но мы получили приближение к
диагональной матрице отсюда следует, что собственные числа матрицы
A
()k
ii
a
, что
собственные векторы матрицы
A
()k
iiii
Ue a
λ
()k
ii
a
11 12 1 1
21 22 2 2
12
...
0
...
0
... ... ... ... ...
1
...
0
ni
ni
i
n n nn ni
uu u u
uu u u
Ue
i
uu u u
⎛⎞
⎛⎞
⎜⎟
⎜⎟
⎜⎟
⎜⎟
==
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎝⎠
⎝⎠
То есть
i-тому собственному числу будет соответствовать i-тый столбец
матрицы
U (определяющий собственный вектор).
Пример: Найти собственные числа и собственные векторы матрицы А методом
вращения.
11
1
321 cos sin 0
231 sin cos 0
113 0 0 1
ii
AU
ϕ
ϕ
ϕϕ
==
⎛⎞
⎜⎟
==
⎜⎟
⎜⎟
⎝⎠
матрица вращения осей i j
(1, 2)ij==
2
22
tg 2
33 4
ij
ii jj
a
aa
π
ϕϕ
====
−−
11 11
00
22 22
11 11
00
11
22 22
001 001
T
UU
⎛⎞⎛⎞
⎜⎟⎜⎟
⎜⎟⎜⎟
⎜⎟⎜⎟
⎜⎟⎜⎟
⎜⎟⎜⎟
⎜⎟⎜⎟
⎜⎟⎜⎟
⎝⎠⎝⎠
==
1 1 32 32 5 1
011
22 22 22
321
1 1 23 23 5 1
231 0 1 1
1
22 2 2 22
113
001 2 2
03 03
22
AU
⎛⎞
+−+
⎛⎞
⎜⎟
⎜⎟
⎜⎟
⎛⎞
⎜⎟
⎜⎟
⎜⎟
+−+
⎜⎟
==
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎝⎠
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎝⎠
⎝⎠
=
                                                                    ⎛0⎞ ⎛ 0 ⎞
        ⎛ λ1 0 ... 0 ⎞                                              ⎜0⎟ ⎜ 0 ⎟
        ⎜ 0 λ ... 0 ⎟                                               ⎜ ⎟ ⎜ ⎟
      Λ=⎜       2         ⎟,                                Λei = Λ ⎜ M ⎟ = ⎜ M ⎟ = λ i ei
        ⎜ ... ... ... ... ⎟                                         ⎜ ⎟ ⎜ ⎟
        ⎜⎜                ⎟⎟                                        ⎜1⎟ i ⎜λ i ⎟
         ⎝ 0 0 0 λn ⎠                                               ⎜0⎟ ⎜ 0 ⎟
                                                                    ⎝ ⎠ ⎝ ⎠
      Собственные числа матрицы A –                                λi , но мы получили приближение к
                                                                                                     (k )
диагональной матрице отсюда следует, что собственные числа матрицы A – aii , что
                                                              (k )
собственные векторы матрицы A – Uei ⇔ aii                            ≈ λi aii( k )
            ⎛ u11 u12                 ... u1n ⎞ ⎛ 0 ⎞ ⎛ u1i ⎞
            ⎜u     u                  ... u2 n ⎟⎟ ⎜⎜ 0 ⎟⎟ ⎜⎜ u2i ⎟⎟
      Uei = ⎜ 21 22                                       =
            ⎜ ... ...                 ... ... ⎟ ⎜ 1 ⎟ i ⎜ ... ⎟
            ⎜⎜                                  ⎟⎜ ⎟ ⎜ ⎟
             ⎝ un1 un 2               ... unn ⎟⎠ ⎝ 0 ⎠ ⎜⎝ uni ⎟⎠
      То есть i-тому собственному числу будет соответствовать i-тый столбец
матрицы U (определяющий собственный вектор).

      Пример: Найти собственные числа и собственные векторы матрицы А методом
вращения.
          ⎛ 3 2 1 ⎞ i =1                         ⎛ cos ϕ     − sin ϕ        0 ⎞ i =1
     A = ⎜⎜ 2 3 1 ⎟⎟                       U1 = ⎜⎜ sin ϕ     cos ϕ          0 ⎟⎟ –матрица вращения осей i j
          ⎜ 1 1 3⎟                               ⎜ 0            0           1 ⎟⎠
          ⎝        ⎠                             ⎝
                                                                                              (i = 1, j = 2)

                      2aij            2⋅2          π
      tg 2ϕ =                     =       = ∞ ⇒ϕ =
                 aii − a jj           3−3          4
            ⎛    1           1         ⎞                ⎛ 1            1      ⎞
            ⎜           −             0⎟                ⎜ 2                  0⎟
            ⎜     2           2        ⎟                ⎜               2     ⎟
      U1 = ⎜⎜    1          1          ⎟
                                      0⎟          U1 = ⎜⎜ − 1
                                                   T                   1      ⎟
                                                                             0⎟
                  2          2                               2          2
            ⎜                          ⎟                ⎜                     ⎟
            ⎜⎜   0          0         1⎟                ⎜⎜ 0           0     1⎟
                                       ⎟                                      ⎟
             ⎝                         ⎠                 ⎝                    ⎠

                         ⎛            1           1     ⎞ ⎛ 3+ 2     −3+ 2   ⎞ ⎛ 5      −1    ⎞
                         ⎜                   −         0⎟ ⎜                 1⎟ ⎜             1⎟
             ⎛ 3 2 1⎞ ⎜                2           2    ⎟ ⎜ 2           2    ⎟ ⎜ 2       2    ⎟
                                                        ⎟ ⎜ 2+3              ⎟ ⎜ 5            ⎟
       AU1 = ⎜⎜ 2 3 1 ⎟⎟ ⎜⎜           1          1
                                                       0⎟ = ⎜
                                                                     −2 + 3
                                                                            1⎟ = ⎜
                                                                                        1
                                                                                             1⎟
              ⎜        ⎟               2          2           2         2    ⎟ ⎜ 2       2
             ⎝ 1 1 3⎠ ⎜                                 ⎟ ⎜                                   ⎟
                         ⎜            0          0     1⎟ ⎜ 2                ⎟ ⎜ 2            ⎟
                         ⎜                              ⎟ ⎜            0    3⎟ ⎜        0    3⎟
                         ⎝                              ⎠ ⎝ 2                ⎠ ⎝ 2            ⎠




                                                                                                            45