Процесс Грама ― Шмидта

Процесс Грама ― Шмидта преобразует последовательность линейно независимых векторов ${𝐚}_1,\dots ,{𝐚}_n$ в ортонормированную систему векторов ${𝐞}_1,\dots ,{𝐞}_n$, причём так, что каждый вектор ${𝐞}_j$ есть линейная комбинация ${𝐚}_1,\dots ,{𝐚}_j$.

Классический процесс Грама — Шмидта

Алгоритм

Пусть имеются линейно независимые векторы ${𝐚}_1,\dots ,{𝐚}_n$ и пусть ${𝐩𝐫𝐨𝐣}_{𝐛}𝐚$ — оператор проекции вектора $𝐚$ на вектор $𝐛$, определённый как

${𝐩𝐫𝐨𝐣}_{𝐛}𝐚=\frac{⟨𝐚,𝐛⟩}{⟨𝐛,𝐛⟩}𝐛,$

где $⟨𝐚,𝐛⟩$ — скалярное произведение векторов $𝐚$ и $𝐛$.

Классический процесс Грама — Шмидта выполняется следующим образом:

$\begin{array}{cccc}{𝐛}_1& =& {𝐚}_1& (1)\\ {𝐛}_2& =& {𝐚}_2-{𝐩𝐫𝐨𝐣}_{{𝐛}_1}{𝐚}_2& (2)\\ {𝐛}_3& =& {𝐚}_3-{𝐩𝐫𝐨𝐣}_{{𝐛}_1}{𝐚}_3-{𝐩𝐫𝐨𝐣}_{{𝐛}_2}{𝐚}_3& (3)\\ & \dots & & \\ {𝐛}_n& =& {𝐚}_n-{𝐩𝐫𝐨𝐣}_{{𝐛}_1}{𝐚}_n-{𝐩𝐫𝐨𝐣}_{{𝐛}_2}{𝐚}_n-\dots -{𝐩𝐫𝐨𝐣}_{{𝐛}_{n-1}}{𝐚}_n& (n)\end{array}$

На основе каждого вектора ${𝐛}_j(j=1\dots n)$ может быть получен нормированный вектор ${𝐞}_j$ единичной длины, определённый как

${𝐞}_j=\frac{{𝐛}_j}{\Vert {𝐛}_j\Vert }$

Результаты процесса Грама — Шмидта:

${𝐛}_1,\dots ,{𝐛}_n$ — система ортогональных векторов либо

${𝐞}_1,\dots ,{𝐞}_n$ — система ортонормированных векторов.

Вычисление ${𝐛}_1,\dots ,{𝐛}_n$ носит название ортогонализации Грама — Шмидта, а ${𝐞}_1,\dots ,{𝐞}_n$ — ортонормализации Грама — Шмидта.

Геометрическая интерпретация

Файл:Gram-Schmidt-step2.svg

Рис. 1. Второй шаг процесса Грама — Шмидта

Рассмотрим формулу (2) — второй шаг алгоритма. Её геометрическое представление изображено на рис. 1:

1. получение проекции вектора ${𝐚}_2$ на ${𝐛}_1$;
2. вычисление ${𝐚}_2-{𝐩𝐫𝐨𝐣}_{{𝐛}_1}{𝐚}_2$, то есть перпендикуляра, которым выполняется проецирование ${𝐚}_2$ на ${𝐛}_1$. Этот перпендикуляр — вычисляемый в формуле (2) вектор ${𝐛}_2$;
3. перемещение полученного на шаге 2 вектора ${𝐛}_2$ в начало координат. Это перемещение сделано на рисунке лишь для наглядности;

На рисунке видно, что вектор ${𝐛}_2$ ортогонален вектору ${𝐛}_1$, так как ${𝐛}_2$ является перпендикуляром, по которому ${𝐚}_2$ проецируется на ${𝐛}_1$.

Рассмотрим формулу (3) — третий шаг алгоритма — в следующем варианте:

$\begin{array}{ccc}{𝐛}_3={𝐚}_3-({𝐩𝐫𝐨𝐣}_{{𝐛}_1}{𝐚}_3+{𝐩𝐫𝐨𝐣}_{{𝐛}_2}{𝐚}_3)& & (6)\end{array}$

Её геометрическое представление изображено на рис. 2:

Файл:Gram-schmidt-step3base.svg

Рис. 2. Третий шаг процесса Грама — Шмидта

1. получение проекции вектора ${𝐚}_3$ на ${𝐛}_1$;
2. получение проекции вектора ${𝐚}_3$ на ${𝐛}_2$;
3. вычисление суммы ${𝐩𝐫𝐨𝐣}_{{𝐛}_1}{𝐚}_3+{𝐩𝐫𝐨𝐣}_{{𝐛}_2}{𝐚}_3$, то есть проекции вектора ${𝐚}_3$ на плоскость, образуемую векторами ${𝐛}_1$ и ${𝐛}_2$. Эта плоскость закрашена на рисунке серым цветом;
4. вычисление ${𝐚}_3-({𝐩𝐫𝐨𝐣}_{{𝐛}_1}{𝐚}_3+{𝐩𝐫𝐨𝐣}_{{𝐛}_2}{𝐚}_3)$, то есть перпендикуляра, которым выполняется проецирование ${𝐚}_3$ на плоскость, образуемую векторами ${𝐛}_1$ и ${𝐛}_2$. Этот перпендикуляр — вычисляемый в формуле (6) вектор ${𝐛}_3$;
5. перемещение полученного ${𝐛}_3$ в начало координат. Это перемещение сделано на рисунке лишь для наглядности. Оно не является математическим действием и поэтому не отражается в формуле (6).

На рисунке видно, что вектор ${𝐛}_3$ ортогонален векторам ${𝐛}_1$ и ${𝐛}_2$, так как ${𝐛}_3$ является перпендикуляром, по которому ${𝐚}_3$ проецируется на плоскость, образуемую векторами ${𝐛}_1$ и ${𝐛}_2$.

Таким образом, в процессе Грама — Шмидта для вычисления ${𝐛}_j$ выполняется проецирование ${𝐚}_j$ ортогонально на гиперплоскость, натянутую на векторы ${𝐛}_1,\dots ,{𝐛}_{j-1}$. Вектор ${𝐛}_j$ затем вычисляется как разность между ${𝐚}_j$ и его проекцией. То есть ${𝐛}_j$ — это перпендикуляр от ${𝐚}_j$ к гиперплоскости, натянутой на векторы ${𝐛}_1,\dots ,{𝐛}_{j-1}$. Поэтому ${𝐛}_j$ ортогонален векторам, образующим эту гиперплоскость.

Особые случаи

Процесс Грама — Шмидта может применяться также к бесконечной последовательности линейно независимых векторов.

Кроме того, процесс Грама — Шмидта может применяться к линейно зависимым векторам. В этом случае он выдаёт $\mathrm{𝟎}$ (нулевой вектор) на шаге $j$, если ${𝐚}_j$ является линейной комбинацией векторов ${𝐚}_1,\dots ,{𝐚}_{j-1}$. Для сохранения ортогональности выходных векторов и для предотвращения деления на ноль при ортогонализации алгоритм должен отбрасывать нулевые векторы. Количество векторов, выдаваемых алгоритмом, будет равно размерности подпространства, порождённого векторами (то есть количеству линейно независимых векторов, которые можно выделить среди исходных векторов).

Свойства

• Матрица перехода от ${𝐚}_1,\dots ,{𝐚}_j$ к ${𝐛}_1,\dots ,{𝐛}_j$ — верхнетреугольная матрица.

• Произведение длин ${𝐛}_1,\dots ,{𝐛}_j$ равно объёму параллелепипеда, построенного на векторах системы ${𝐚}_1,\dots ,{𝐚}_j$ как на рёбрах.

• Процесс Грама ― Шмидта для векторов-столбцов матрицы из $GL(n)$ определяет гомотопическую эквивалентность $GL(n)\to O(n)$.

Дополнительные толкования

Процесс Грама ― Шмидта может быть истолкован как разложение невырожденной квадратной матрицы в произведение ортогональной (или унитарной в случае эрмитова пространства) и верхнетреугольной матрицы с положительными диагональными элементами ― QR-разложение, что является частным случаем разложения ИвасавыШаблон:Нет АИ.

Литература

• Шаблон:Книга

Ссылки

• https://www.youtube.com/watch?v=GhhpG1tCQLU Анимация процесса на плоскости

Шаблон:Rq