Ортогональная матрица. Ортогонализация Грамма-Шмидта.

Квадратная матрица ||aij||  i=1,2,...n,   j=1,2,...n называется ортогональной, если выполняется условие

AAT=E,

где E - единичная матрица, AT- транспонированная матрица.

Из условия AAT=E следует, что AT является обратной к матрице A:

A-1=AT.

Для любой невырожденной квадратной матрицы можно построить ортогональную матрицу. Для построения ортогональной матрицы применяют метод ортогонализации Грамма-Шмидта. Затем нормируют полученные векторы строки. Эти две процедуры вместе называют методом ортонормализации Грамма-Шмидта.

Ортогонализация Грамма-Шмидта

Пусть задана некоторая квадратная матрица, строки которой являются векторы

a_1, a_2, ...a_n

и пусть эти векторы линейно независимы (т.е. матрица построенная этими векторами строками невырождена). Требуется получить взаимно ортогональные n векторы

b_1,b_2,b_3

Суть метода заключается в следующем:

1. Выбирается некоторая строка (пусть это будет a1). b1 выбирается равным a1.

2. Вектор b2 получается пректированием a2 на нуль-пространство b1.

3. Вектор b3 получается проектированием a3 на нуль-пространство матрицы

b_1,b_2

И т.д.

Рассмотрим подробнее процесс ортогонализации.

На первом шаге выбираем b1=a1.

На втором шаге вычисляем нуль-пространство b1:

N(b1)=E-b1^+b1

где E- единичная матрица порядка nxn,-псевдообратная к b1. Так как b1 является вектором-строкой (матрицей-строкой), то

b1^+=b1^T*1/b1*b1^T

Для пректирования a2 на нуль-пространство b1 вычисляем

a2-a2b1^Tb1/b1b1^T

На третьем шаге вычисляем b3:

Так как векторы b1 и b2 ортогональны, то

Следовательно

Для n -го вектора получим:

Таким образом, процедура ортогонализации Грамма-Шмидта имеет следующий вид:

Ортонормализация Грамма-Шмидта

Суть метода ортонормализации Грамма-Шмидта заключается в ортогонализации методом Грамма-Шмидта а затем нормализации полученных векторов строк:

b1/||b1||,...

Полученная матрица

(c1, c2, ...,cn)^T

является ортогональной матрицей т.к.