Альтернативный способ задать дифференциал в геометрической алгебре / forpes.ru

Главная
Альтернативный способ задать дифференциал в геометрической алгебре

Альтернативный способ задать дифференциал в геометрической алгебре +2

29.03.2026 12:08

Exlt8 10 7400 Источник

Привет, Хабр!

В геометрической алгебре достаточно абстрактно введен дифференциал, здесь предлагается наглядный численный метод — выразить дифференциал через ориентированный объём и геометрическое произведение.

Это даёт возможность интерполировать значения функции вне сетки и отдельно учитывать параллельную и ортогональную составляющие приращения.

Написана статья с целью собрать мнения специалистов о достоинствах и недостатках такого подхода. В общем буду рад комментариям.

1. Обозначения

Δ — конечные приращения на сетке (шаги Δx,Δy,Δz и разности функции fΔxf,Δyf,Δzf).
δ — изменяемые приращения для произвольной точки в окрестности узла.
σ1,σ2,σ3 — орты осей в матричном представлении (матрицы Паули):
$σ _1 =\begin{bmatrix} 0 & 1 \\ 1 & 0 \end{bmatrix}, \quad σ _2 =\begin{bmatrix} 0 & -i \\ i & 0 \end{bmatrix}, \quad σ _3 =\begin{bmatrix} 1 & 0 \\ 0 & -1 \end{bmatrix}.$

Умножение с точкой ⋅ — матричное (геометрическое).
Умножение с крестиком × — векторное.
I=σ1⋅σ2⋅σ3 — псевдоскаляр (аналог мнимой единицы в ГА).

Вектор в таком базисе записывается как

Обратите внимание: σ3 соответствует оси — x, σ2 — y, σ1 — z. Так сделано по мотивам моей предыдущей статьи Давайте объединим линейную и геометрическую алгебры. На простом примере. Часть 1. В разделе "6. Базис Клиффорда для матриц 2х2" можно видеть, что расположение осей несколько иное, чем в декартовых координатах. Так же такое расположение еще делает привычнее некоторые представления, напишу как нибудь.

2. Ориентированный объём и вектор конечных разностей

Ориентированный объём, образованный площадями граней ячейки и смещениями δx,δy,δz:

$ΔS_x =Δy⋅Δz,\quad ΔS_y =Δx⋅Δz, \quad ΔS_z =Δx⋅Δy$

Разделив на объём ячейки Δv=ΔxΔyΔz, получим относительный ориентированный объём:

$δV ^′ = \frac {δV} {Δv} = \frac {Δx}{δx} ⋅ σ _3 + \frac {Δy}{δy} ⋅ σ _2 + \frac {Δz}{δz} ⋅ σ _1 .$

Как то так как на обложке статьи будет это выглядеть

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

$\Delta(f) = \Delta_x f\,\sigma_3 + \Delta_y f\,\sigma_2 + \Delta_z f\,\sigma_1$ $\Delta_x f = f(x+\Delta x,y,z)-f(x,y,z),\quad$ $\Delta_y f = f(x,y+\Delta y,z)-f(x,y,z)$ $\Delta_z f = f(x,y,z+\Delta z)-f(x,y,z)$

3. Скалярное и обобщённое приращения

Скалярное приращение (классический дифференциал, аппроксимированный конечными разностями) — это скалярная часть геометрического произведения δV′ и Δ(f):

$\delta f_{\|} = \frac{\Delta_x f}{\Delta x}\,\delta x + \frac{\Delta_y f}{\Delta y}\,\delta y + \frac{\Delta_z f}{\Delta z}\,\delta z$

При

$\Delta x,\Delta y,\Delta z\to0 \quad и \quad \delta x,\delta y,\delta z\to0$

это выражение стремится к обычному дифференциалу df.

Полное геометрическое произведение даёт обобщённый дифференциал — сумму скалярной и бивекторной частей:

$\delta_g f = \delta f_{\|} + I\;\delta f_{\perp}$ $\delta f_{\perp} = \left( \Delta_y(f) \; \frac{\delta z}{\Delta z} - \Delta_z(f) \; \frac{\delta y}{\Delta y} \right) \; \sigma_3 + \left( \Delta_x(f) \; \frac{\delta z}{\Delta z}-\Delta_z(f) \; \frac{\delta x}{\Delta x} \right) \; \sigma_2 + \left( \Delta_x(f) \; \frac{\delta y}{\Delta y} - \Delta_y(f) \; \frac{\delta x}{\Delta x} \right) \cdot \sigma_1$

Бивекторная часть δf⊥ равна нулю тогда и только тогда, когда вектор относительного объёма δV′ параллелен вектору конечных разностей Δ(f) — то есть когда смещение (δx,δy,δz) коллинеарно градиенту. Когда бивекторная часть не равна нулю происходит следующее:

Бивекторная часть кодирует нормаль к плоскости перемноженных векторов. Плоскость в которой расположены компоненты диффенциала. Так же она кодирует объект в этой плоскости, ортогональный вектору на который происходит проецирование.
Скалярная часть ортогональна бивекторной части и нормали к плоскости.
Вместе они образуют аналог комплексного числа в этой плоскости.

4. Зачем это нужно

Интерполяция вне сетки. Зная в узле сетки значения Δxf,Δyf,Δzf (например, из численного решения), можно вычислить приращение δgfдля любого смещения (δx,δy,δz) в окрестности узла. Никакой дополнительной сетки не требуется.
Разделение параллельной и ортогональной составляющих. Скалярная часть δf∥ — классическое изменение вдоль направления смещения. Бивекторная часть δf⊥ — мера того, насколько смещение отклоняется от градиента.
Выделение траекторий, где обобщённый дифференциал совпадает с классическим. Если потребовать δf⊥=0, получим условия коллинеарности, которые задают прямые, параллельные Δ(f). Вдоль таких траекторий δgf=δf∥.

5. Пример применения

Пусть в узле сетки (0,0,0)(0,0,0) известны:

$\Delta x = 0.1,\; \Delta y = 0.1,\; \Delta z = 0.1$ $\Delta_x f = 0.2,\; \Delta_y f = 0.3,\; \Delta_z f = 0.4$

Это соответствует частным производным: 2,3,4.

Требуется найти значение функции в точке (δx,δy,δz)=(0.02,0.03,0.04) — внутри ячейки.

Вычисляем скалярное приращение:

$\delta f_{\|} = \frac{0.04}{0.1}\cdot0.3 + \frac{0.3}{0.1}\cdot0.03 + \frac{0.4}{0.1}\cdot0.04 = 0.29$

Вычисляем бивекторную часть (коэффициенты при σ3,σ2,σ1):

$c_3 = 0.3\cdot\frac{0.04}{0.1} - 0.4\cdot\frac{0.03}{0.1} =0$ $c_2 = 0.4\cdot\frac{0.02}{0.1} - 0.2\cdot\frac{0.04}{0.1} =0$ $c_1 = 0.2\cdot\frac{0.03}{0.1} - 0.3\cdot\frac{0.02}{0.1} = 0$

Здесь смещение оказалось точно параллельным вектору конечных разностей, поэтому бивекторная часть обнулилась. Обобщённый дифференциал совпал с классическим: δgf=0.29

Если взять смещение (δx,δy,δz)=(0.03,0.02,0.04), то:

Скалярное приращение:

$\delta(f) = 2\cdot0.03 + 3\cdot0.02 + 4\cdot0.04 = 0.28$

Бивекторные коэффициенты:

$c_3 = 0.3\cdot0.4 - 0.4\cdot0.2 = 0.04$ $c_2 = 0.4\cdot0.3 - 0.2\cdot0.4 =0.04$ $c_1 = 0.2\cdot0.2 - 0.3\cdot0.3 = -0.05$

Бивекторная часть не равна нулю — смещение не параллельно градиенту. Полное приращение δgf будет содержать как скалярную, так и бивекторную компоненту, что позволяет оценить степень «неградиентности» направления.

6. Итог

Предложен способ выразить обобщённый дифференциал через конечные разности и геометрическое произведение.
Формулы работают непосредственно с сеточными данными и позволяют интерполировать значения функции в произвольную точку.
Бивекторная часть даёт дополнительную информацию — меру отклонения направления смещения от градиента.
Для вычислений достаточно скалярных арифметических операций, матрицы не требуются (если не считать исходного определения σ).

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

Пишите ваши комментарии!

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

runaway
29.03.2026 13:09
#29741032
Когда я учился в универе, меня сильно били по рукам линейкой за формализм:

dx -- это когда мы говорим о бесконечно малых, которые могут расти или убывать (например, приращение координаты);

δх -- это когда мы говорим о бесконечно малых, которые не могут убывать (например, приращение совершённой работы в процессе).

У вас в статье есть некая фривольность в этом плане.
1. Exlt8 Автор
  29.03.2026 13:09
  #29741074
  Так мне показалось читать проще. Как пледложите назвать фиксированное и изменяемое, но не бесконечно малое приращение, чтобы не путаться в многообразии дельт?
  1. runaway
    29.03.2026 13:09
    #29741362
    Стерильный формализм, конечно -- dx. Это подходит ко всему. δх -- когда хотят подчеркнуть, что дифференциал не убывает.
    
    Exlt8 Автор
    29.03.2026 13:09
    #29741388
    Так это у меня конечное приращение, а не дифференциал. В статье нет дифференциалов

MasterMentor
29.03.2026 13:09
#29741434
В статье развиты идеи нестандартного анализа.

https://ru.wikipedia.org/wiki/Гипервещественное_число

Термин «гипервещественное число» предложен американским математиком Эдвином Хьюиттом в 1948 году. Теорию поля гипервещественных чисел как расширения поля вещественных чисел опубликовал в 1960-е годы Абрахам Робинсон, который назвал её «нестандартным анализом».

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

Ключевые особенности НА:

-- Принцип переноса: Любое математическое утверждение верное для обычных чисел, остается верным и для гипервещественных. Это позволяет переносить правила алгебры на бесконечно малые.

-- Стандартная часть - st: Операция st() «округляет» гипервещественное число до ближайшего обычного. Что заменяет процедуру вычисления предела.

Выгоды НА:

– Функция непрерывна, если бесконечно малое изменение аргумента вызывает бесконечно малое изменение значения

– Дифференциалы dx и dy становятся реальными числами, а не символами. Их можно делить, умножать и сокращать, как в обычной школе.

– Геометрическая наглядность: Криволинейная трапеция действительно состоит из бесконечно тонких прямоугольников. Это больше не метафора.

– Связь с производной: Основная теорема анализа (Ньютона-Лейбница) доказывается почти автоматически: мы просто суммируем бесконечно малые приращения функции df, и их сумма дает общее приращение f(b)-f(a)

В статье используется инвариантность формы первого дифференциала. В нестандартном анализе это следствие того, что dy = f’(x)dx — это прямое равенство между числами. Это позволяет строить геометрические выводы (через ориентированные объемы и векторы), оперируя приращениями как физическими векторами, а не переменными “стремящимися к нулю”.
1. MasterMentor
  29.03.2026 13:09
  #29741584
  И немного что значит «инвариантность формы» первого дифференциала.
  
  Первый дифференциал df это результат скалярного произведения градиента ∇f на вектор перемещения dr. Имеем df = (∇f, dr). Поскольку скалярное произведение двух векторов - число не зависящее от выбора системы координат (это числовой инвариант), значит df будет неизменным.
  
  В обычном матанализе инвариантность формы — это теорема о правиле дифференцирования сложной функции. В нестандартном анализе акцент смещается на геометрию:
  
  – Мы смотрим на функцию в бесконечно малой окрестности («под микроскопом»).
  
  – Там поверхность функции выглядит как плоскость.
  
  – У этой плоскости есть наклон (градиент).
  
  – Формула df = (∇f, dr) — это просто уравнение этой плоскости.
  
  – Так как геометрия плоскости не зависит от системы координат, то и формула её описания (первый дифференциал) сохраняет свой вид.
  
  PS Это же, вероятно (проверить!), будет справедливо и для криволинейных координат: «внутри» них это скалярное произведение будет давать одно и то же числовое значение приращения df.
  
  Для второго дифференциала это уже не работает т.к. он “меряет” кривизну самой плоскости (или самой системы координат). В функции d(dy) - появится не нулевое (второе) слагаемое, которое нельзя отбросить.
  1. Exlt8 Автор
    29.03.2026 13:09
    #29741948
    Саму идею взял не от сюда) Взял из англоязычной статьи в Вики о теореме о дивергенции, там дано второе определение дивергенции через предел проекции на относительный объем. В русскоязычной статье этого нет, так что хорошо, что догадался заглянуть
    
    MasterMentor
    29.03.2026 13:09
    #29742716
    В "русскоязычной Википедии" - много чего нет. И многое намеренно написано так, чтобы делать из читателя - идиота. Рекомендую всегда пользоваться английской, тем более переводчики - позволяют.
    https://translate.google.com/?hl=ru&sl=auto&tl=ru&op=websites

alexshvets
29.03.2026 13:09
#29743606
По существу здесь интересна не нотация и не нестандартный анализ, а сама дискретная геометрическая конструкция. Вы берёте конечные разности, собираете из них локальный «градиентный» вектор и умножаете его на произвольное смещение внутри ячейки. Скалярная часть даёт обычное направленное приращение первого порядка, бивекторная — естественную меру непараллельности смещения и дискретного градиента. Это живая идея, и развивать стоит именно её.
1. Exlt8 Автор
  29.03.2026 13:09
  #29743612
  Спасибо, ровно это и имел ввиду

Альтернативный способ задать дифференциал в геометрической алгебре +2

1. Обозначения

2. Ориентированный объём и вектор конечных разностей

3. Скалярное и обобщённое приращения

4. Зачем это нужно

5. Пример применения

6. Итог

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

runaway

Exlt8 Автор

runaway

Exlt8 Автор

MasterMentor

MasterMentor

Exlt8 Автор

MasterMentor

alexshvets

Exlt8 Автор