Кривизна Гаусса

Материал из Википедии — свободной энциклопедии
Перейти к: навигация, поиск
Слева направо: поверхность с отрицательной гауссовой кривизной (гиперболоид), поверхность с нулевой гауссовой кривизной (цилиндр), и поверхность с положительной гауссовой кривизной (сфера).

Гауссова кривизна — мера искривления поверхности в окрестности какой-либо её точки. Гауссова кривизна является объектом внутренней геометрии поверхностей, в частности, не изменяется при изометрических изгибаниях.

Определение[править | править вики-текст]

Кривизна Гаусса для двумерной поверхности[править | править вики-текст]

Обозначим нормальные кривизны в главных направлениях (главные кривизны) в рассматриваемой точке поверхности \kappa_1 и \kappa_2. Величина:

K=\kappa_1\kappa_2

называется гауссовой кривизной, полной кривизной или просто кривизной поверхности. Встречается также термин скаляр кривизны, который подразумевает результат свёртки тензора кривизны; при этом скаляр кривизны вдвое больше, чем гауссова кривизна.

Гауссова кривизна может быть вычислена через метрику поверхности, и поэтому она является объектом внутренней геометрии (отметим, что главные кривизны к внутренней геометрии не относятся). По знаку кривизны можно классифицировать точки поверхности (см. рисунок). Кривизна плоскости равна нулю. Кривизна сферы радиуса R всюду равна \frac{1}{R^2}. Существует и поверхность постоянной отрицательной кривизны — псевдосфера.

Кривизна Гаусса для гиперповерхности[править | править вики-текст]

Кривизна n-мерной гиперповерхности в точке полностью описывается её главными кривизнами k^{(1)}, k^{(2)}, \dots, k^{(n)} и соответствующими главными направлениями.

Рассмотрим (с точностью до знака) симметрические многочлены, составленные из чисел

(1) \qquad k^{(1)}, k^{(2)}, \dots, k^{(n)} :



(2) \qquad 
\left\{\begin{array}{rcl}
K^{[1]} &=& - (k^{(1)} + k^{(2)} + \dots + k^{(n)}) = - \sum\limits_{i} k^{(i)} \\
K^{[2]} &=& k^{(1)} k^{(2)} + k^{(1)} k^{(3)} + \dots + k^{(n-1)} k^{(n)} = \sum\limits_{i < j} k^{(i)} k^{(j)} \\
\cdots &\cdots& \cdots \cdots \cdots \cdots \cdots  \cdots  \cdots  \cdots  \\
K^{[n]} &=& (-1)^n k^{(1)} k^{(2)} \cdots k^{(n)} \\
\end{array}\right.

Назовем вышеприведенные величины кривизнами Гаусса соответствующей степени. Общая формула кривизны Гаусса степени m запишется так:

(3) \qquad K^{[m]} = \sum_{i_1 < i_2 < \dots < i_m} k^{(i_1)} k^{(i_2)} \cdots k^{(i_m)}

Кривизны Гаусса являются коэффициентами характеристического многочлена для матрицы тензора полной кривизны гиперповерхности:

(4) \qquad \det(\lambda \delta^i_j - b^i_j) = \lambda^n + K^{[1]} \lambda^{n-1} + \dots + K^{[n-1]} \lambda + K^{[n]}

Тензорная формула для кривизны Гаусса[править | править вики-текст]

Формула (3) определяет кривизну Гаусса через собственные числа тензора полной кривизны гиперповерхности b_{ij}. Попробуем выразить эти величины через компоненты самого тензора b_{ij} в любой системе координат. Для вычисления определителя произвольного тензора второго ранга мы имеем такую формулу с использованием тензора метрической матрешки (см. Абсолютно антисимметричный единичный тензор):

(5) \qquad \det(a^i_j) = {1 \over n!} g^{i_1 i_2 \dots i_n}_{j_1 j_2 \dots j_n} a^{j_1}_{i_1} a^{j_2}_{i_2} \cdots a^{j_n}_{i_n}

Подставим в эту формулу a^i_j = \lambda \delta^i_j - b^i_j, чтобы вычислить левое выражение формулы (4), тогда имеем:

(6) \qquad n! \det (\lambda \delta^i_j - b^i_j) = g^{i_1 i_2 \dots i_n}_{j_1 j_2 \dots j_n} (\lambda \delta^{j_1}_{i_1} - b^{j_1}_{i_1}) \cdots (\lambda \delta^{j_n}_{i_n} - b^{j_n}_{i_n})

Раскроем скобки в формуле (6). Поскольку тензор метрической матрешки g^{i_1 i_2 \dots i_n}_{j_1 j_2 \dots j_n} не меняется при синхронной перестановке верхних и нижних индексов, то все слагаемые при одинаковой степени \lambda^m будут одинаковыми (их количество равно биномиальному коэффициенту C^m_n), и мы получаем:

(7) \qquad n! \det (\lambda \delta^i_j - b^i_j) = \lambda^n g^{s_1 s_2 \dots s_n}_{s_1 s_2 \dots s_n} - C^1_n \lambda^{n-1} g^{i s_2 \dots s_n}_{j s_2 \dots s_n} b^j_i + C^2_n \lambda^{n-2} g^{i_1 i_2 \dots s_n}_{j_1 j_2 \dots s_n} b^{j_1}_{i_1} b^{j_2}_{i_2} - \dots

Поскольку последовательные свертки тензора метрической матрешки равны:

(8) \qquad g^{i_1 i_2 \dots i_m s_{m+1} s_{m+2} \dots s_n}_{j_1 j_2 \dots j_m s_{m+1} s_{m_2} \dots s_n} = 
(n-m)! \, g^{i_1 \dots i_m}_{j_1 \dots j_m}

То из формулы (7) и формулы для биномиальных коэффициентов C^m_n = {n! \over m! (n-m)!} находим такую формулу для характеристического многочлена (разделив обе стороны уравнения (7) на n!):

(9)  \qquad  \det (\lambda \delta^i_j - b^i_j) = \lambda^n - {\lambda^{n-1} \over 1!} g^i_j b^j_i + 
{\lambda^{n-2} \over 2!} g^{i_1 i_2}_{j_1 j_2} b^{j_1}_{i_1} b^{j_2}_{i_2} - \dots

Сравнивая формулы (9) и (4), находим такую формулу для кривизны Гаусса:

(10) \qquad K^{[m]} = {(-1)^m \over m!} g^{i_1 i_2 \dots i_m}_{j_1 j_2 \dots j_m} b^{j_1}_{i_1} b^{j_2}_{i_2} \dots b^{j_m}_{i_m}

Выражение через тензор Римана[править | править вики-текст]

Для скалярной кривизны гиперповерхности мы имеем такую ​​формулу

(11) \qquad R = g^{ik} g^{jl} R_{ijkl} = 2 \sum_{i < j} k^{(i)} k^{(j)} = 2 K^{[2]}

Чтобы обобщить эту формулу для более высоких степеней, попробуем заменить произведение двух метрических тензоров в формуле (11) на тензор метрической матрешки четвертого ранга:

(12) \qquad g^{ijkl} R_{ijkl} = \begin{vmatrix} g^{ik} & g^{il} \\ g^{jk} & g^{jl} \end{vmatrix} R_{ijkl} = 
(g^{ik} g^{jl} - g^{il} g^{jk}) R_{ijkl} = 2 R = 4 K^{[2]}

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

(13) \qquad g_{ij} = \delta_{ij} = \begin{cases} 1, & i=j \\ 0, & i \ne j \end{cases}

а потому мы можем численно не различать ковариантные и соответствующие контравариантные компоненты тензоров (верхние и нижние индексы). Тензор Римана в точке P будет в некотором смысле диагональным, а именно, его ненулевые компоненты будут равны:

(14) \qquad R_{ijij} = -R_{ijji} = k^{(i)} k^{(j)} \qquad (i \ne j)

и равны нулю все те компоненты R_{ijkl}, где вторая пара индексов (kl) не совпадает с (ij) с точностью до перестановки в паре.

Левая часть формулы (12) является линейной формой от тензора Римана, а коэффициентами этой формы служат компоненты тензора метрической матрешки. Очевидным обобщением является рассмотрение билинейной формы и форм высших степеней от компонента тензора Римана. Проведем вычисления формулы (12) еще раз и таким образом, чтобы эти вычисления можно было легко обобщить. Имеем, учитывая диагональность тензора Римана:

(15) \qquad g^{ij}_{kl} R^{kl}_{ij} = \sum_{i,j} \left ( \sum_{k,l} g^{ij}_{kl} R^{kl}_{ij} \right ) = \sum_{i, j} \left ( g^{ij}_{ij} R^{ij}_{ij} + g^{ij}_{ji} R^{ji}_{ij} \right )

Далее, два слагаемых в правой части формулы (15) одинаковы вследствие антисимметрии по индексам внутри пары как тензора метрической матрешки, так тензора Римана. Кроме того, диагональная компонента метрической матрешки равна единице, поскольку (в следующей формуле сложения по одинаковым индексам не производится, а индексы i, j разные):

(16) \qquad g^{ij}_{ij} = \begin{vmatrix} \delta^i_i & \delta^i_j \\ \delta^j_i & \delta^j_j \end{vmatrix} = \begin{vmatrix} 1 & 0 \\ 0 & 1 \end{vmatrix} = 1

Учитывая вышесказанное и формулу (14), превращаем формулу (15) далее:

(17) \qquad g^{ij}_{kl} R^{kl}_{ij} = 2 \sum_{i \ne j} 1 \cdot k^{(i)} k^{(j)} = 2 \cdot 2! \sum_{i < j} k^{(i)} k^{(j)} = 2 \cdot 2! K^{[2]}

Теперь перейдем к вычислению следующей квадратичной формы:

(18) \qquad \Phi_2(R) = g^{i_1 j_1 i_2 j_2}_{k_1 l_1 k_2 l_2} R^{k_1 l_1}_{i_1 j_1} R^{k_2 l_2}_{i_2 j_2}

Коэффициентами этой формы служат компоненты тензора метрической матрешки восьмого ранга. Этот тензор имеет две группы индексов, и является антисимметричным по перестановке индексов внутри этих групп. Вычисляем аналогично формуле (15).

(19) \qquad \Phi_2(R) = \sum_{i_1, j_1, i_2, j_2} \left ( \sum_{k_1, l_1, k_2, l_2} g^{i_1 j_1 i_2 j_2}_{k_1 l_1 k_2 l_2} R^{k_1 l_1}_{i_1 j_1} R^{k_2 l_2}_{i_2 j_2} \right ) = 2^2 \sum_{i_1, j_1, i_2, j_2} g^{i_1 j_1 i_2 j_2}_{i_1 j_1 i_2 j_2} R^{i_1 j_1}_{i_1 j_1} R^{i_2 j_2}_{i_2 j_2}

Обозначим индексы i_1, j_1, i_2, j_2 как i,j,k,l для упрощения записи:

(19a) \qquad \Phi_2(R) = 2^2 \sum_{i, j, k, l} g^{ijkl}_{ijkl} R^{ij}_{ij} R^{kl}_{kl} = 2^2 4! \sum_{i, j, k, l \over all\; different} g^{ijkl}_{ijkl} k^{(i)} k^{(j)} k^{(k)} k^{(l)}

Все четыре индекса i,j,k,l должны быть попарно различными, поскольку компоненты тензора метрической матрешки равны нулю при наличии двух одинаковых индексов в одной группе. В правой сумме формулы (19a) стоят диагональные компоненты тензора метрической матрешки, которые равны единице (аналогично формуле 16).

(19b) \qquad \Phi_2(R) = 2^2 \sum_{i, j, k, l \over all\; different} k^{(i)} k^{(j)} k^{(k)} k^{(l)} = 2^2 \cdot 4! \sum_{i < j < k < l} k^{(i)} k^{(j)} k^{(k)} k^{(l)} = 2^2 \cdot 4! K^{[4]}

Множитель 4! при переходе ко второй сумме в формуле (19a) возник вследствие того, что для одного слагаемого в правой сумме, характеризующегося фиксированным набором четырех различных чисел i < j < k < l, соответствует 4! = 24 одинаковых по величине слагаемого в левой сумме, характеризующихся перестановками этих четырех чисел.

Формулы (19), (19a), (19b) легко обобщаются на формы высших степеней. Таким образом получаем общую формулу для нахождения кривизны Гаусса парной степени 2 m:

(20) \qquad K^{[2 m]} = {1 \over 2^m (2 m)!} g^{i_1 j_1 \dots i_m j_m}_{k_1 l_1 \dots k_m l_m} R^{k_1 l_1}_{i_1 j_1} \cdots R^{k_m l_m}_{i_m j_m}

Альтернативный вывод формулы кривизны Гаусса для парной степени[править | править вики-текст]

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

(21) \qquad R^{kl}_{ij} = b^k_i b^l_j - b^k_j b^l_i

и начнем в формуле (10) группировать сомножители по два, например начиная с первых двух (здесь мы считаем, что степень 2 m кривизны Гаусса не меньше двух (m \ge 1), и для упрощения записи опустим обозначения m):

(22) \qquad (2 m)! K = g^{i j \dots}_{k l \dots} b^k_i b^l_j \cdots = - g^{j i \dots}_{k l \dots} b^k_i b^l_j \cdots

Последнее преобразование справедливо вследствие антисимметрии тензора метрической матрешки относительно индексов в верхней группе. Далее, в последнем выражении поменяем местами индексы i, j:

(23) \qquad (2 m)! K = - g^{i j \dots}_{k l \dots} b^k_j b^l_i \cdots

Теперь добавим уравнение (22) и (23), при этом учтя (21). Получаем, опять изменив обозначение индексов:

(24) \qquad 2 (2 m)! K^{[2 m]} = g^{i_1 j_1 i_2 j_2 \dots i_m j_m}_{k_1 l_1 k_2 l_2 \dots k_m l_m} R^{k_1 l_1}_{i_1 j_1} b^{k_2}_{i_2} \cdots b^{k_m}_{i_m} b^{l_m}_{j_m}

Множитель 2 в левой части уравнения (24) появился в результате группировки двух множителей b^{k_1}_{i_1} b^{l_1}_{j_1}. Очевидно, мы можем аналогичным образом сгруппировать попарно и остальные сомножители, тогда в левой части мы получим множитель 2^m, а в правой - выражение, в котором участвует только тензор Римана и тензор метрической матрешки, т.е. мы получим формулу (20).

Кривизна Гаусса нечетной степени[править | править вики-текст]

Кривизна Гаусса нечетной степени также связана с тензором Римана, но более сложными формулами, чем (20). К тому же из этих формул кривизна Гаусса выражается неоднозначно.

Значение кривизны Гаусса[править | править вики-текст]

В начале было дано определение кривизны Гаусса только для гиперповерхности (формулы 2, 3). Но формула (20), как и формулы для нахождения кривизны Гаусса нечетной степени, позволяют распространить это понятие на произвольные (абстрактные) многообразия. Таким образом мы можем рассматривать кривизны Гаусса как скалярные инварианты тензора Римана.

Внутренняя кривизна многообразия полностью описывается тензором Римана.

Кривизну Гаусса как скаляр можно интегрировать по объему всего многообразия (смотрите статью Интегралы Гаусса). Интеграл от K [n] является топологическим инвариантом n-мерного многообразия (не меняется при непрерывной деформации многообразиях).

Формула Бриоски для двумерной поверхности[править | править вики-текст]

ds^2 = Edu^2+2Fdudv+Gdv^2 \,
и их производные первого и второго порядков по так называемой формуле Бриоски[1]:
 K =\frac{\begin{vmatrix} -\frac{1}{2}E_{vv} + F_{uv} - \frac{1}{2}G_{uu} & \frac{1}{2}E_u & F_u-\frac{1}{2}E_v\\F_v-\frac{1}{2}G_u & E & F\\\frac{1}{2}G_v & F & G \end{vmatrix}- \begin{vmatrix} 0 & \frac{1}{2}E_v & \frac{1}{2}G_u\\\frac{1}{2}E_v & E & F\\\frac{1}{2}G_u & F & G \end{vmatrix}}{(EG-F^2)^2}
  • В случае F=0 (так называемая ортогональная параметризация поверхности) эту формулу можно переписать в виде
K = -\frac{1}{2\sqrt{EG}}\left(\frac{\partial}{\partial u}\frac{G_u}{\sqrt{EG}} + \frac{\partial}{\partial v}\frac{E_v}{\sqrt{EG}}\right).
  • В случае, если поверхность задана как график дифференцируемой функции z = f(x,y) в трёхмерном евклидовом пространстве с метрикой ds^2 = dx^2+dy^2+dz^2, эта формула принимает вид
K = \frac{f_{xx}\cdot f_{yy}- f_{xy}^2}{(1+f_x^2+ f_y^2)^2}
  • Если поверхность в трёхмерном евклидовом пространстве с метрикой ds^2 = dx^2+dy^2+dz^2 задана уравнением f(x,y,z) = 0, формула принимает вид[2]
K=\frac{[f_z(f_{xx}f_z-2f_xf_{xz})+f_x^2f_{zz}][f_z(f_{yy}f_z-2f_yf_{yz})+f_y^2f_{zz}]-[f_z(-f_xf_{yz}+f_{xy}f_z-f_{xz}f_y)+f_xf_yf_{zz}]^2}{f_z^2(f_x^2+f_y^2+f_z^2)^2}
  • В случае, если первая квадратичная форма конформно эквивалентна евклидовой, то естьE = G = e^{u} с некоторой функцией u=u(x,y) и F= 0, формула для кривизны Гаусса принимает вид
 K = -\frac{1}{2e^u}\Delta u,
где \Deltaоператор Лапласа.

См. также[править | править вики-текст]

Примечания[править | править вики-текст]

Литература[править | править вики-текст]