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


Схема установки для исследования методом сопротивления: 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.

Формула (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)$

Это общая формула скалярного произведения векторов (взаимных норм).

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

Выражения вида (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)$.



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

Для измерения такого («кажущегося») сопротивления источник тока (напряжение) прикладывается к одним узлам (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 вершины различны:



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

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

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



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

$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-х вершинах, показали его основные свойства. В следующей статье рассмотрим пространство графа, определим на нем точечный базис и разберемся с его возможностями.

Комментарии (6)


  1. R9A_019
    13.10.2017 14:40
    #10471160

    Возможно я чего то не понял, но грубо говоря скалярное произведение — это машина где на входе два вектора а на выходе число. Так вот, то что вы написали в формуле 3.5 называется "спариванием", а не скалярным произведением. Если выражаться выражаться через тензоры -скалярное произведение это тензорное произведение ковариантного метрического тензора ранга (0,2) на первый контравариантный вектор и на второй контравариантный вектор, вот здесь и наблюдается свертка тензора и двух векторов. А у вас произведение вектора на ковектор, повторюсь эта процедура называется спариванием и не требует существование метрики того пространства где они существуют.


    1. dmagin Автор
      14.10.2017 11:44
      #10472526

      Да, вы все верно написали. Только сколько новых терминов вам пришлось использовать, чтобы донести мысль об отличии «спаривания» векторов от их скалярного произведения? Вряд ли в контексте серии эта разница существенна, но в любом случае замечание полезное.

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


  1. IvanKor2017
    14.10.2017 09:51
    #10472436

    Помню, взял китайский тестер и два шампура. Воткнул их в грунт и получил что сопротивление больше 1 МОм, удивился. Разнес шампура и воткнул в грунт на расстоянии 40 метров,
    вуаля, между ними только 600 Ом. Вот что значит тензор.


    1. dmagin Автор
      14.10.2017 12:06
      #10472548

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


      1. IvanKor2017
        14.10.2017 12:50
        #10472584

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


        1. dmagin Автор
          14.10.2017 20:28
          #10473008

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