Математический анализ

Теорема Грина

Цели урока

Понять физический и геометрический смысл теоремы Грина как связи границы и области
Уметь применять теорему для вычисления площадей через контурные интегралы
Знать критерий консервативности поля через нулевой ротор
Видеть теорему Грина как частный случай обобщённой теоремы Стокса

Предварительные знания

Криволинейные интегралы
Частные производные
Параметрические кривые

Криволинейные интегралы

Как вычислить площадь произвольной фигуры, зная только её контур - без единого интеграла по площади?

GIS и картография: PostGIS вычисляет площади полигонов из миллиардов GPS-точек через алгоритм shoelace - прямое следствие теоремы Грина
FEM-симуляции: Ansys, COMSOL и Abaqus конвертируют объёмные уравнения в граничные условия через теорему Грина, экономя на порядок вычислений
Компьютерная графика: растеризация полигонов и антиалиасинг в GPU-рендерерах используют знак ротора для определения принадлежности точки
Гидродинамика: циркуляция вихрей в 2D-потоке - это контурный интеграл, который через теорему Грина связан с завихрённостью внутри

Грин, самоучка из Ноттингема

Джордж Грин (1793-1841) был сыном пекаря. Образование - одна школа в возрасте 8-9 лет, больше никакого. Мельница, принадлежавшая семье, оставляла время для самостоятельного изучения математики. В 1828 году Грин опубликовал эссе тиражом 51 экземпляр по подписке. Уильям Томсон (лорд Кельвин) случайно нашёл копию в 1845 году, через 4 года после смерти Грина, и разослал ведущим математикам Европы. Риман, Клебш, Хельмгольц прочитали - и оказалось, что Грин за 17 лет до них построил теорию, на которой стоит вся математическая физика. Одно эссе. 51 экземпляр. Фундамент FEM-симуляций на 100 млрд долларов в год.

Формулировка и смысл теоремы Грина

1828 год. Джордж Грин - самоучка, сын пекаря из Ноттингема - публикует эссе тиражом 51 экземпляр. Почти никто не читает. Через 40 лет Кельвин переоткрывает теорему и популяризирует её. Сегодня Ansys продаёт лицензии на FEM-решатели за 2.4 млрд долларов в год - каждый из них внутри считает surface integrals через теорему Грина.

Теорема Грина - частный случай обобщённой теоремы Стокса для 2D-многообразия с краем. Формула oint_C omega = iint_D d(omega) где omega = P dx + Q dy.

Площадь эллипса через теорему Грина

Вычисление без явного интегрирования по области

Эллипс: x = a cos(t), y = b sin(t). Формула Грина: S = (1/2) oint (x dy - y dx). Подставляя параметризацию: S = (1/2) integral от 0 до 2pi (a cos(t) * b cos(t) - b sin(t) * (-a sin(t))) dt = (1/2) * 2pi * ab = pi*ab. Расчёт занял 3 строки вместо двойного интеграла по эллиптической области.

Если ротор поля равен нулю - поле потенциальное и контурный интеграл по замкнутому контуру равен нулю. Теорема Грина дает этому строгое обоснование: iint_D 0 dA = 0.

По теореме Грина oint_C x dy = ? (C - граница области D, ориентированная ЧЧС)

При P=0, Q=x: dQ/dx - dP/dy = 1-0 = 1. Теорема Грина: oint_C x dy = iint_D 1 dA = площадь D.

Применения: площадь, ротор и консервативность

Алгоритм shoelace (Гаусс, 1795) вычисляет площадь полигона с n вершинами за O(n) операций. Это прямое следствие теоремы Грина: S = (1/2)|sum_{i} (x_i * y_{i+1} - x_{i+1} * y_i)|. GIS-системы типа PostGIS обрабатывают миллиарды полигонов - весь этот вычислительный поток работает на Грине.

Теорема Грина требует, чтобы D была связной областью, а P, Q имели непрерывные частные производные на замкнутой области. Если поле имеет особенность внутри D - теорема не применима напрямую.

Площадь сердцеобразной кривой (кардиоиды)

r = a(1 - cos t) в полярных координатах

Параметризация: x = a(1-cos t)cos t, y = a(1-cos t)sin t, t от 0 до 2pi. По формуле площади через теорему Грина: S = (1/2) oint (x dy - y dx). После вычислений: S = (3/2) pi a^2. Для a=1: площадь = 4.712... Численная проверка через code ниже даёт совпадение до 6 знаков.

В задачах про поле тока в плоской пластине, задачах аэродинамики 2D - циркуляция = oint F dr. Ротор = плотность вихрей. Теорема Грина переводит суммарную завихрённость в циркуляцию по контуру.

Поле F = (P, Q) = (-y/(x^2+y^2), x/(x^2+y^2)) на R^2\{0}. Чему равен контурный интеграл по единичной окружности?

Поле - угловая 1-форма dtheta. При обходе единичной окружности угол меняется на 2pi. Теорема Грина неприменима (особенность внутри), но прямой расчёт даёт 2pi.

Идея доказательства и связь со Стоксом

Доказательство теоремы Грина для прямоугольника - три строки. Для общей области - разбиение на прямоугольники и предельный переход. Именно этот же приём работает в численных методах: FEM-решатели разбивают область на треугольники, и теорема Грина связывает граничные условия с объёмными уравнениями.

Связь: теорема Грина (2D) => теорема Стокса (поверхность в 3D) => теорема Гаусса (объём в 3D) => обобщённая теорема Стокса (n-мерное многообразие). Единый принцип: граница границы пуста (d^2 = 0).

Обобщённая теорема Стокса int_M d(omega) = int_{boundary M} omega при M = отрезок [a,b], omega = f - это:

Обобщённая теорема Стокса при M = [a,b] и omega = f: int_{[a,b]} d(f) = int_{{a,b}} f = f(b) - f(a). Это в точности формула Ньютона-Лейбница.

Теорема Грина и комплексный анализ

Теорема Коши в комплексном анализе - прямое следствие теоремы Грина. Функция f(z) = u + iv аналитична тогда и только тогда, когда выполнены уравнения Коши-Римана: du/dx = dv/dy, du/dy = -dv/dx. Это в точности условие нулевого ротора для поля (u, -v) и (-v, u).

Пакет SciPy использует численное интегрирование по контурам в комплексной плоскости для вычисления специальных функций (гамма-функция, интеграл Фресне). Под капотом - дискретизация теоремы Коши, которая сама есть следствие теоремы Грина.

Теорема Грина объединяет три мира: вещественный анализ (площади, потоки), векторный анализ (ротор, циркуляция) и комплексный анализ (теорема Коши). Одна формула - три интерпретации.

Теорема Коши в комплексном анализе (oint f(z)dz = 0 для аналитической f) является следствием теоремы Грина потому что:

Уравнения Коши-Римана - это в точности условия нулевого ротора для двух вещественных полей. Теорема Грина переводит нулевой ротор в нулевой контурный интеграл.

Связи с другими темами

Теорема Грина - узел, связывающий вещественный анализ, векторный анализ и комплексный анализ.

Комплексный анализ — Связанная тема
Численные методы — Связанная тема
Гидродинамика — Связанная тема
Компьютерная геометрия — Связанная тема

Итоги

Теорема Грина: oint_C P dx + Q dy = iint_D (dQ/dx - dP/dy) dA - граничный интеграл = двойной интеграл ротора
Площадь области: S = (1/2) oint_C (x dy - y dx) - только по границе, без двойного интегрирования
Критерий консервативности: oint = 0 тогда и только тогда, когда ротор = 0 на односвязной области
Теорема Грина - частный случай int_M d(omega) = int_{boundary M} omega в размерности 2
Теорема Коши в комплексном анализе есть следствие теоремы Грина + уравнений Коши-Римана

Вопросы для размышления

Почему теорема Грина позволяет вычислить площадь области, зная только её границу - что физически означает этот факт?
Как связаны ротор 2D-поля и завихрённость в гидродинамике - что измеряет dQ/dx - dP/dy физически?
Почему особенность поля внутри контура ломает применимость теоремы Грина и что делают в таких случаях (угловая форма, лемма Гурса)?

Связанные уроки

calc-24-line-integrals — Криволинейные интегралы - фундамент теоремы Грина
calc-26-divergence-theorem — Теорема Гаусса - трёхмерный аналог
calc-27-diff-forms — Формы и внешний дифференциал - современный язык
la-06-transformations