Геометрия данных 3. Скалярное произведение векторов

Любые системы координат — лишь подпорки для описания «положения вещей». Поэтому важно отличать понятия, связанные с системой координат и независимые от нее — инварианты. В данной статье рассмотрим один из таких инвариантов — скалярное произведение векторов. Одна из целей серии (начало здесь) — обновить сложившиеся стереотипы, поэтому постараемся выйти за рамки привычного геометрического смысла.

59db677295fa5901506865.jpeg
Схема установки для исследования методом сопротивления: A и B — питающие заземления; M и N — измерительные заземления; 1 — измерительный прибор (из книги «Электроразведка», Якубовский Ю. B., M., 1980).

На рисунке показана схема измерения скалярного произведения векторов $\vec{AB}$ и $\vec{MN}$ на грунте.

Тензоры координат


Напомним, что пространства точек и пространство векторов — это два разных пространства. В точечном базисе это различие выражено в структуре координат — ди-координаты векторов имеют нулевую скалярную компоненту, а би-координаты векторов являются балансовыми (сумма их компонент равна нулю).

Если координаты точек пространства собрать в список (массив), то получим матрицу. Строки данной матрицы будут соответствовать точкам пространства, а столбцы — компонентам координат, то есть базовым точкам. Обозначим массив ди-координат точек как $X_{ia}$. Индекс $i$ означает точку пространства, а индекс $a$ — точку базиса. Метрический тензор определяется на вершинах базиса, поэтому при «жонглировании индексами» (переходе от ди-координат к би-координатам $X_i^a$ и обратно) индекс точки положение не меняет:

$X_i^a=X_{ib}g^{ab}$ и $X_{ia}=X_i^b g_{ab} \quad(3.1)$

Взаимные нормы точек пространства будут матрицей:

$N_{ij}=X_{ia} X_j^a \quad(3.2)$

Как показано в предыдущей части (см. 2.17.2, 2.17.3), нормы точек связаны с дистанцией до пространства базиса, и по этой причине не являются инвариантами.

Нормы векторов


Вектор представляет собой разность координат двух точек, соответственно его тензор имеет два индекса точек и один индекс базиса:

$V_{ij,a}=X_{ja} - X_{ia}, \space \space V^{a}_{ij}=X^a_j - X^a_i \quad(3.3)$

Здесь определены массивы координат ди- $V_{ij,a}$ и би- $V^{a}_{ij}$ координат векторов соответственно. Тензоры координат антисимметричны относительно перестановок индексов точек: $V_{ij,a}=-V_{ji,a}$.

Дистанция между точками $i$ и $j$ будет двумерной матрицей $D_{ij}$ норм векторов:

$D_{ij} = V_{ij,a} V^{a}_{ij} \quad(3.4)$

В (3.4) учтено, что правило свертки индексов применяется только к индексам разного уровня.

Дистанция между точками является инвариантом — не зависит от выбора базиса.

Скалярное произведение векторов


Под скалярным произведением понимается свертка координат векторов по индексу базиса. Результат свертки двух векторов — скаляр, не зависит от выбора системы координат (инвариант). Если $\vec{u}$ и $\vec{v}$ — это два неких вектора, то их скалярное произведение в тензорной форме можно записать как:

$\vec{u} \cdot \vec{v}=u_a v^a \quad(3.5)$

Выражение аналогично скалярному произведению координат точек, но тут речь о перемножении векторов. Далее будем называть также результат скалярного произведения векторов взаимной нормой векторов, чтобы отличать понятие операции скалярного произведения (свертки любых координат) от результата ее применения к векторам.

Поскольку умножать можно любые векторы пространства, то в общем случае выражение для взаимной нормы векторов представляет собой тензор:

$D_{ij,kl} = V_{ij,a} V^{a}_{kl} \quad(3.6)$

Индексы $ij$ задают координаты одного вектора, а индексы $kl$ — другого.

Связанные пары — три вершины, теорема косинусов


Согласно определению вектор можно представить как сумму других векторов, начало и конец которых совпадают:

$ V_{ij,a}=V_{ik,a}+V_{kj,a}=-V_{ki,a}+V_{kj,a} $

$ V^a_{ij}=V^a_{ik}+V^a_{kj}=-V^a_{ki}+V^a_{kj} \quad(3.7) $

Если теперь раскрыть норму вектора $\vec{V_{ij}}$, подставив (3.7) в (3.4), то получим выражение для взаимной нормы (скалярного произведения) векторов $\vec{V_{ki}}$ и $\vec{V_{ki}}$:

$ D_{ki,kj}=V_{ki,a} V^{a}_{kj}=(D_{ik} + D_{jk} - D_{ij}) / 2 \quad(3.8)$

Здесь $D_{ik}$ — норма вектора $\vec{V_{ik}}$, совпадает с дистанцией между точками i и k.
59dc80c3f1816969388285.png
Формула (3.8) — это обычная теорема косинусов для треугольника. Здесь векторы, между которыми определена взаимная норма, имеют общую точку k, (смежные векторы).

Квадрат скалярного произведения на 3-х вершинах связан с площадью образуемого ими треугольника $S_{ijk}$ (формула Герона):

$ 4S_{ijk}^2 = D_{ik} D_{jk} - (D_{ik,jk})^2 \quad(3.9)$

Независимые пары — четыре вершины


В общем случае взаимная норма может быть определена для векторов, не имеющих общих точек. Вывод формулы аналогичен (3.8), только вместо одной промежуточной точки используем две. Получаем:

$D_{ij,kl}=(D_{il} + D_{jk} - D_{ik} - D_{jl}) / 2 \quad(3.10)$

Это общая формула скалярного произведения векторов (взаимных норм).
59dc87b4b9af3103252327.png
Геометрически значение взаимной нормы отражает произведение косинуса угла между векторами на их длину. Порядок индексов важен — задает направление векторов.

Выражения вида (3.10) (четыре слагаемых — два положительных, два отрицательных) встречаются в разных разделах математики. Его присутствие обычно означает, что скорее всего есть и пара векторов, скалярное произведение которых задает структуру данного выражения.

Симметрия тензора норм векторов


Перечислим свойства тензора взаимных норм $D_{ij,kl}$.

1) Очевидно, что он антисимметричен относительно перестановки индексов $ij$ или $kl$:

$D_{ij,kl}=-D_{ji,kl}=D_{ij,lk} \quad(3.11.1)$

2) Не зависит от перестановки векторов:

$D_{ij,kl}=D_{kl,ij} \quad(3.11.2)$

3) Существует только две независимых нормы на заданных 4-х точках ввиду тождества:

$D_{ij,kl} + D_{ik,lj} + D_{il,jk} = 0 \quad(3.11.3)$

Знатоки математики должны увидеть в формуле (3.11.3) первое (алгебраическое) тождество Бьянки. Из чего можно сделать вывод, что структуры тензора кривизны (Римана) и тензора взаимных норм векторов — подобны.

4) Векторная норма на 4-х вершинах может быть выражена через разность норм на 3-х вершинах:

$D_{ij,kl} = D_{ij,il} - D_{ij,ik} \quad(3.11.4)$

Это полезное тождество, которое нам еще понадобится.

Скалярное произведение векторов на графах


На графах скалярное произведение векторов — это сопоставление двух пар вершин. Допустим, задан граф с известной матрицей смежности, а, значит, и лапласианом $L$. Тогда скалярным произведением между вектором разности узлов $i, j$ с одной стороны и вектором на узлах $k, l$ с другой является значение кофактора 2-го порядка, деленное на скалярный потенциал лапласиана (алгебраическое определение):

$D_{ij,kl}=U_{ij,kl}(L) / u(L) \quad(3.12.1)$

Напомним, что кофактором $U_{ij,kl}$ называется определитель минора квадратной матрицы (с учетом знака). Скалярный потенциал $u(L)$ — это кофактор 1-го порядка от лапласиана (см. (1.12) из первой статьи). То есть для получения значения взаимной номы векторов на графе необходимо из матрицы лапласиана удалить столбцы, соответствующие одному вектору (в формуле это i-й и j-й столбцы), и строки, соответствующие другому (k и l), после чего разделить определитель получившегося минора на скалярный потенциал $u(L)=det'(L)$.

59df66765c650961401182.png

Если удаляемые столбцы и строки лапласиана — одни и те же, то получим значение эффективного сопротивления между узлами графа (резистивную дистанцию). Таким образом скалярное произведение векторов связано с понятием обобщенного сопротивления в электрической цепи.

Для измерения такого («кажущегося») сопротивления источник тока (напряжение) прикладывается к одним узлам (A и B), а разность потенциалов измеряется между другими (M и N). (Это основной метод измерений в электрометрии — см. КДПВ). Взаимная норма (скалярное произведение) векторов равна нулю, если внешняя разность потенциалов не приводит к разности потенциалов на измеряемом векторе.

В следующей статье разберем подробнее, почему отношение заданной и измеряемой разностей потенциалов узлов связано со скалярным произведением векторов.

Обращение минора лапласиана


Подматрицу значений взаимных норм можно получить обращением минора лапласиана. Обозначим лапласиан, из которого удалены a-я строка и b-й столбец как $L^{i(a),j(b)}$. Тогда имеет место тождество

$D_{ia,jb}=(L^{i(a),j(b)})^{-1} \quad(3.12.2)$

Формула (3.12.2) — один из способов получить дистанционную матрицу по заданному лапласиану. Удаляем из лапласиана графа какой-либо из узлов (пусть будет a) и обращаем. Получаем матрицу взаимных норм векторов $D_{ia,ja}$. Значение индекса a фиксировано — базовый узел. Матрица $D_{ia,ja}$ представляет собой скалярные произведения векторов, начало которых находится в базовом узле a, а конец пробегает по точкам пространства.

Для восстановления дистанционной матрицы $D_{ij}$ добавим к множеству узлов базовый (удаленный из лапласиана), то есть просто расширим матрицу $D_{ia,ja}$ нулевой строкой и столбцом. Далее применяем к матрице преобразование дистанции:

$F_{ij}(X) = X_{ii} + X_{jj} - 2X_{ij}$.

В результате получим дистанционную матрицу $D_{ij}=F_{ij}(D_{ia,ja})$.

Квадрат скалярного произведения, матрица Якоби


Если в графе изменить значение проводимости ребра (элемент лапласиана), то очевидно, что изменятся и все дистанции между вершинами (нормы векторов) $D_{ij}$. При увеличении проводимости дистанции должны сократиться (уменьшиться). Подарок богов в том, что можно оценить изменение дистанций не только качественно, но и количественно. Обозначим производную матрицы дистанций по лапласиану как

$Q_{ij,kl}=\partial D_{ij} / \partial L^{kl} \quad(3.13)$

Тензор $Q_{ij,kl}$ — это матрица Якоби, то есть выражение изменений значений дистанционной матрицы $D_{ij}$ при изменении значений лапласиана $L^{kl}$. Оказывается, что данный тензор выражается через квадрат скалярного произведения векторов (квадрат норм):

$ Q_{ij,kl}=(D_{ij,kl})^2=(D_{il} + D_{jk} - D_{ik} - D_{jl})^2/4 \quad(3.14)$

Круто, что сущности не множатся. Вариация лапласиана связана с вариацией дистанционной матрицы через известное понятие скалярного произведения векторов. Этим подчеркивается роль данного тензора. Кроме того, тензор квадратов скалярного произведения $(D_{ij,kl})^2$ можно обращать.

Тензор квадросвязности вершин


Выражение (3.13) можно представить в таком эквивалентном виде:

$\partial D_{ij} = Q_{ij,kl} \space \partial L^{kl} \quad(3.13')$

Данную формулу можно трактовать как отклик $\partial D_{ij}$ на воздействие $\partial L^{kl}$. Тензор $Q_{ij,kl}$ играет роль передаточной функции (реакции на воздействие).

Возможна и обратная ситуация, при которой воздействие и отклик меняются местами. Прямая и обратная $Q^{ij,kl}$ передаточные функции связаны соотношением:

$Q_{ij,kl}Q^{kl,mn}=\delta_{ij}^{mn} \quad(3.15)$

Снова удача — тензор $Q^{ij,kl}$ можно выразить через лапласиан:

$Q^{ij,kl}=\partial L^{ij} / \partial D_{kl}=( L^{il} L^{jk} + L^{ik} L^{jl})/2 \quad(3.16)$

Если важен тензор взаимных норм, то и тензор $Q^{ij,kl}$ должен быть не менее значим. Назовем его тензором квадросвязности и кратко рассмотрим его свойства.

Свойства квадросвязности вершин


Значения тензора определены через значения лапласиана для 4-х вершин. Будем считать данные вершины вершинами графа. Допустим, что все 4 вершины различны:

59df7763079d8201832515.png

Здесь векторами обозначены пары вершин, между которыми считается квадросвязность. Пары имеют отличную от нуля квадросвязность только тогда, когда их вершины попарно связаны (необходимые связи показаны на рисунке одинаковым цветом). Если все связи в графе положительны, то квадросвязность между разными вершинами также всегда больше или равна нулю.

Если пары имеют общую вершину, то смысл квадросвязности меняется. Это связано с тем, что диагональные элементы лапласиана не равны нулю (как в дистанционной матрице), а отражают общую связность (проводимость) узла.

$Q^{ki,kj}=( L^{ki} L^{kj} + L^{kk} L^{ij})/2 \quad(3.16.1)$

59df7f683c5d9759958700.png

Пары вершин могут совпадать — диагональные элементы тензора квадросвязности.

$Q^{ij}=Q^{ij,ij}=(L^{ii} L^{jj} + (L^{ij})^2)/2 \quad(3.16.2)$

Здесь $Q^{ij}$ характеризует связь двух узлов i и j. Считается как сумма произведения их суммарной проводимости (степени вершины) $L^{ii} L^{jj}$ и квадрата связи между узлами $(L^{ij})^2$.

Несмотря на то, что формально для тензора $Q^{ij,kl}$ могут вычисляться элементы вида $Q^{ij,ll}$ (одна из пар вырождена), данные (вырожденные) элементы являются линейно-зависимыми от остальных. Могут быть вычислены через сумму тензора по одному из индексов пары:

$Q^{ij,ll} = -Q^{ij,kl} 1_k \quad(3.16.3)$
___

На этом пока все. Мы определили инвариант скалярного произведения векторов на 4-х вершинах, показали его основные свойства. В следующей статье рассмотрим пространство графа, определим на нем точечный базис и разберемся с его возможностями.

© Habrahabr.ru