10.07.2015 Views

С.П. Шарый - Институт вычислительных технологий СО РАН

С.П. Шарый - Институт вычислительных технологий СО РАН

С.П. Шарый - Институт вычислительных технологий СО РАН

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

3.17. Численные методы для проблемы собственных значений 419не изменяется при умножениях на ортогональные матрицы слеваили справа.Доказательство. Напомним, что следом матрицы A = (a ij ), обозначаемымtrA, называется сумма всех её диагональных элементов:n∑trA = a ii .Нетрудно проверить, что привлечение понятия следа позволяет переписатьопределение фробениусовой нормы матрицы таким образом⎛ (n∑ n) ⎞ 1/2∑‖A‖ F = ⎝ a ij a ij⎠ = ( tr(A ⊤ A) ) 1/2.j=1i=1<strong>С</strong>ледовательно, для любой ортогональной матрицы Q справедливо‖QA‖ F =(tr ( (QA) ⊤ (QA) )) 1/2=i=1(tr ( A ⊤ Q ⊤ QA )) 1/2=(tr(A ⊤ A) ) 1/2= ‖A‖F .Для доказательства аналогичного соотношения с умножением на ортогональнуюматрицу справа заметим, что фробениусова норма не меняетсяпри транспонировании матрицы. <strong>С</strong>ледовательно,‖AQ‖ F = ∥ ∥ ( Q ⊤ A ⊤) ⊤∥∥F = ‖Q ⊤ A ⊤ ‖ F = ‖A ⊤ ‖ F = ‖A‖ F ,что завершает доказательство <strong>П</strong>редложения.<strong>С</strong>ледствие. Фробениусова норма матрицы не меняется при ортогональныхпреобразованиях подобия.Для более точного описания меры близости матриц A (k) , которыепорождаются конструируемым нами методом, к диагональной матрицевведём величину( ) 1/2 ∑ND(A) = a 2 ijj≠i— фробениусову норму внедиагональной части матрицы. Ясно, чтоматрица A диагональна тогда и только тогда, когда ND(A) = 0.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!