- •Глава1. Проблема аппроксимации
- •§1. Полиномиальная апппроксимация
- •§2. Интерполяционный полином в форме Лагранжа
- •§3. Интерполяционный полином в форме Ньютона
- •§4. Аппроксимация сплайнами
- •§5. Метод наименьших квадратов
- •§6. Полиномиальная интерполяция с кратными узлами
- •§7. Свойства разделенных разностей
- •§8. Задача Чебышева. Разрешимость системы
- •§9. Теорема Чебышева
- •§10. Многочлены Чебышева
- •Глава2. Численное дифференцирование
- •Глава3. Численное интегрирование
- •§1. Интерполяционные квадратурные формулы
- •1.Интерполяционные квадратурные формулы
- •2.Интерполяционные квадратурные формулы наилучшей алгебраической точности
- •3.Ортогональные многочлены и их свойства
- •§2. Применение квадратурных формул
- •§3. Метод Монте-Карло (метод статистических испытаний)
- •§4. Правило Рунге практической оценки погрешности
- •Глава4. Алгебраическая проблема собсвенных значений
- •§1. Ортогональные матрицы
- •1.Ортогональные матрицы
- •2.Матрица элементарного поворота
- •§2. Вариационное свойство собственных значений
- •§3. Приведение симметричной матрицы к диагональному виду
- •§4. Сингулярное разложение матрицы
- •§5. Сопряженная матрица
- •§6. Частная спектральная задача
- •1.Вариационный метод
- •2.Степенной метод
- •§7. Метод максимизации столбцов
- •1.Максимизация первого столбца
- •2.Алгоритм сингулярного разложения
- •3.Главное собственное число
- •§8. Метод вращения
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) иначе