Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Chislennye_metody.doc
Скачиваний:
27
Добавлен:
09.05.2015
Размер:
1 Mб
Скачать

3.Главное собственное число

Пусть А=(aij) – симметричная матрица. Метод максимизации столбцов дает следующее:

Перемножим равенства:

максимальное собственное число матрицы A2(из метода максимизации).

λ=±b– главное собственное число симметричной матрицы А.

§8. Метод вращения

Рассмотрим симметричную матрицу А=(аij). Обозначим ортогональную матрицу вращенияUpq(α):

uqq=upp=c=cosα, -upq =uqp= -s= -sinα, остальные диагональные элементы равны 1, а недиагональные – нулю.

Рассморим вращение матрицы в плоскости Opq:

Матрица В=AUpqотличается от А столбцамиbpиbq:

bp=c·ap+s·aq

bq= -s·ap+c·aq

bj=aj, j≠p,q

Матрица отличается от В р-ой иq-ой строками:

dpj=c·bpj+s·bqj

dqj= -s·bpj+c·bqj

Остальные строки не изменяют.

Тогда,

Разобьем Sна диагональную и недиагональную части:

Заметим:

При элементарном вращении Dнедиагональные элементы аpi,aqiи аip,aiq(i≠p,q) меняются так, что попарные суммы квадратов их модулей сохраняются.

Кроме этих элементов вне диагонали меняется элемент аpq.

Т.е. величина S2меняется при элементарном вращении настолько, насколько изменится |apq|2.

Будем подбирать вращения так, чтобы S2максимально уменьшалась.

Положим dpq=0:

dpq=c·bpq+s·bqq=c(-s·app+c·apq)+s(-s·aq p+c·aqq)=

αвыберем следующим образом:

1) при аppcos2α=0 =>α=π/4

2) иначе

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]