- •Глава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. Метод вращения
§1. Ортогональные матрицы
1.Ортогональные матрицы
Рассмотрим два ортонормированных базиса в пространстве Rn: . Т.е..
Разложим вектора второго базиса по векторам первого:
Обозначим как U=(uik).
Рассмотрим скалярное произведение:
Таким образом, сумма Cikтрактуется как произведениеiой строки матрицыUнаkый столбец матрицыUT.Cдругой стороны, из свойств ортонормированного базиса матрицаC=(cik)=E.
Отсюда получим: UUT=EилиUTU=E.
Из выведенных формул следуют равенства:
1)U-1 =UT(из определения обратной матрицы)
2) (Ux,Uy) = (x,y) для любых векторовx,yизRn(т.е. скалярные произведения векторов и их образов относительноUсовпадают).
Действительно, (Ux,Uy) = (UUTx,y) = (x,y)
3)Отсюда в частности следует: ||x|| = ||Ux||
Действительно, .
Матрица V, удовлетворяющая любому из этих трех равенств,ортогональная. Она обладает следующими свойствами:
1)переводит ортогональный базис в ортонормиррованный
2)сохраняет углы между векторами, а также их нормы
Представление матрицы-оператора в другом базисе
Из теорем линейных операторов:
1)Всякой матрице А можно сопоставить некоторый линейный оператор
Ω: x→y=Ax
2)Всякий линейный ограниченный оператор вRnпредставим в виде матрицы.
Рассмотрим некоторый ортонормированный базис и найдем представление в нем матрицы А.
Пусть x=(x1..xn)Tв базисе .
Т.е. .
Рассмотрим оператор Ω: x→yв . Т.е. оператор(β – оператор Ω в новом базисе).
- представление матрицы А в базисе .
2.Матрица элементарного поворота
Рассмотрим матрицу элементарного поворота в двух пространствах:
1)в пространствеR2:
Рассмотрим матрицу, где с=cosα,s=sinα.
Тогда .
Докажем, что U– ортогональная матрица:
;
А следовательно, U– ортогональная матрица.
Задача: доказать, что линейный оператор φ, соответствующий матрице U, есть оператор поворота на угол α.
2)в пространствеRn:
Рассмотрим матрицу, где с=cosα,s=sinα.
Т.е. upp =uqq =c,upq = -uqp =s,uii= |i≠p,q| =1,uij= |i,j≠p,q| = 0.
Задача: доказать, что линейный оператор Φ, соответствующий матрице U, есть оператор поворота на угол α вn-мерном пространстве.
§2. Вариационное свойство собственных значений
Матрица А неотрицательно определена, если.
Лемма Собственные числа неотрицательно определенной матрицы неотрицательны.
Доказательство
Пусть А – неотрицательно определенная матрица, λ1– её собственное число.
Тогда .
Т.к. .
Что и требовалось доказать.
Матрица А симметричная , если.
Лемма 1 Если А симметрична, неотрицательно определена и (Ay,y)=0, тоAy=0. Геометрически: А не может сделать вектор у ортогональным самому себе, а может лишь обратить его в ноль.
Доказательство
Рассмотрим выражение (A(y+tz),y+tz),
где z– произвольный вектор;
t– вещественное, малое число (t<<1).
Если бы (Ay,z)≠0, то знак всего выражения можно было бы сделать отрицательным, выбирая знакt, что противоречило бы условию неотрицательной определенности А.
Значит, (Ay,z)=0.
Т.к. элеент zвыбирался произвольно, тоAy=0.
Что и требовалось доказать.
Рассмотрим выражение: .
Функция nпеременныхF(y)=(Ay,y) непрерывна на единичной сфере. Следовательно, по первой теореме ВейерштрассаF(y) достигает наS1своих точных граней. Это означает, что.
Лемма 2 (вариационное свойство) Если А – симметричная матрица, то- собственное число матрицы А.
Доказательство
Из приведенных выше рассуждений и достигается в некоторой точке.
Т.к. , то ||y||=1, а значит.
Получим: .(*)
Из выбора точку y’ верно: (Ay,y)≤λ1для любого вектораyизS1.
Т.е. Последнее утверждение верно и для любого вектораxиз пространстваRn. Т.е.- неотрицательно определенная матрица.
Следовательно, из (*), используя Лемму 1, получим:
, т.е. λ1– собственное число матрицы А.
Что и требовалось доказать.