Вычислительные методы линейной алгебры. Горбаченко В.И - 84 стр.

UptoLike

ВУЗ: 

ПГУ | Пенза

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

Рубрика: 

84
-149 -50 -154
537 180 546
-27 -9 -25
>> [U,T]=schur(A)
U =
-0.3162 0.6529 0.6882
0.9487 0.2176 0.2294
-0.0000 -0.7255 0.6882
T =
1.0000 -7.1119 815.8706
0 2.0000 55.0236
0 0 3.0000
Собственные значения, равные 1, 2 и 3, расположены на диагонали матрицы
T. Вектор собственных значений может быть получен следующим образом:
>> V = ordeig (T)
V =
1.0000
2.0000
3.0000
В случае комплексных собственных значений получаем:
>> A=[0 2 4; 1 0 0; 0 1 0]
A =
0 2 4
1 0 0
0 1 0
>> [U,T]=schur(A)
U =
0.8729 0.2611 -0.4122
0.4364 -0.0399 0.8988
0.2182 -0.9645 -0.1488
T =
2.0000 -3.3320 1.0659
0 -1.0000 -0.5365
0 1.8639 -1.0000
Матрицу Шура в комплексной форме можно получить следующим образом:
>> [U,T]=schur(A,'complex')
U =
0.8729 -0.3633 + 0.1234i -0.2301 + 0.1949i
0.4364 0.7921 - 0.0189i 0.0352 - 0.4249i
0.2182 -0.1311 - 0.4560i 0.8499 + 0.0703i
T =
2.0000 0.9393 - 1.5752i 2.9361 - 0.5039i
0 -1.0000 + 1.0000i -1.3274
0 0 -1.0000 - 1.0000i
Аналогичный результат можно получить, применив функцию rsf2csf:
>> [U,T]=schur(A)
U =
0.8729 0.2611 -0.4122 
  -149   -50 -154
   537   180   546
   -27    -9   -25
>> [U,T]=schur(A)
U =
   -0.3162    0.6529      0.6882
    0.9487    0.2176      0.2294
   -0.0000   -0.7255      0.6882
T =
    1.0000   -7.1119    815.8706
         0    2.0000     55.0236
         0         0      3.0000
Собственные значения, равные 1, 2 и 3, расположены на диагонали матрицы
T. Вектор собственных значений может быть получен следующим образом:
>> V = ordeig (T)
V =
    1.0000
    2.0000
    3.0000

    В случае комплексных собственных значений получаем:
>> A=[0 2 4; 1 0 0; 0 1 0]
A =
     0     2     4
     1     0     0
     0     1     0
>> [U,T]=schur(A)
U =
    0.8729    0.2611   -0.4122
    0.4364   -0.0399    0.8988
    0.2182   -0.9645   -0.1488
T =
    2.0000   -3.3320    1.0659
         0   -1.0000   -0.5365
         0    1.8639   -1.0000
Матрицу Шура в комплексной форме можно получить следующим образом:
>> [U,T]=schur(A,'complex')
U =
   0.8729            -0.3633   + 0.1234i   -0.2301 + 0.1949i
   0.4364             0.7921   - 0.0189i    0.0352 - 0.4249i
   0.2182            -0.1311   - 0.4560i    0.8499 + 0.0703i
T =
   2.0000             0.9393   - 1.5752i    2.9361 - 0.5039i
        0            -1.0000   + 1.0000i   -1.3274
        0                  0               -1.0000 - 1.0000i
Аналогичный результат можно получить, применив функцию rsf2csf:
>> [U,T]=schur(A)
U =
    0.8729    0.2611     -0.4122

                                                                     84

Страницы