Изопериметрические задачи и множители Лагранжа

Классическая загадка Дидоны

По преданию, финикийская царевна Дидона бежала в Северную Африку и попросила у местного вождя столько земли, «сколько можно охватить шкурой быка». Она разрезала шкуру на тонкие полосы и огородила ими участок — взяв прямую линию берега как одну сторону. Какую форму выбрала Дидона? Конечно, полуокружность. Именно такие задачи — максимизировать площадь при фиксированном периметре — называются изопериметрическими. Строгое математическое доказательство того, что окружность является ответом, потребовало изобретения вариационного исчисления.

Постановка изопериметрической задачи

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

В аналитической форме: max J₁[y] = ∫₀ᵃ y dx при J₂[y] = ∫₀ᵃ √(1+y'²) dx = L.

Это функционал с ограничением-функционалом — аналог задачи условной оптимизации в конечномерном пространстве.

Метод множителей Лагранжа для функционалов

Теорема: для задачи min J[y] = ∫F(x,y,y')dx при J₂[y] = ∫G(x,y,y')dx = C экстремаль является экстремалью вспомогательного функционала:

H[y] = ∫(F(x,y,y') + λ G(x,y,y')) dx

где λ — числовой множитель Лагранжа, то есть экстремаль удовлетворяет уравнению ЭЛ для F + λG:

(F + λG)y − d/dx (F + λG){y'} = 0

Как находить λ? Из условия J₂[y] = C — подставляем найденную экстремаль в ограничение и решаем уравнение для λ.

Задача о катенарии (цепной линии)

Постановка: гибкая цепь фиксированной длины L висит между точками A и B. Под действием силы тяжести она принимает форму, минимизирующую потенциальную энергию. Найти форму цепи.

Потенциальная энергия: J₁ = ∫y√(1+y'²) dx (высота центра масс элемента дуги × длина элемента). Длина: J₂ = ∫√(1+y'²) dx = L.

Вспомогательный функционал с множителем λ: H = ∫(y+λ)√(1+y'²) dx.

Уравнение ЭЛ: функционал H не зависит от x явно → используем первый интеграл:

H − y'H_{y'} = C₁ → (y+λ)/√(1+y'²) = C₁.

Разделяя переменные: dy/√((y+λ)²/C₁² − 1) = dx. Интегрируя: y + λ = C₁ cosh((x−C₂)/C₁).

Ответ: y = a cosh(x/a + b) − λ, где a, b — константы из граничных условий. Это цепная линия (катенарий).

Физический смысл: форма цепи не является параболой (как думали до Гюйгенса и Лейбница) — это гиперболический косинус! Разница хорошо видна при больших провесах.

Задача Дидоны: полуокружность

Задача: огородить максимальную территорию с прямолинейным берегом, используя верёвку длиной L.

Граничные условия: y(0) = y(a) = 0. max ∫₀ᵃ y dx при ∫₀ᵃ √(1+y'²) dx = L.

Вспомогательный функционал: H = ∫(y + λ√(1+y'²)) dx.

Уравнение ЭЛ: 1 − λ d/dx(y'/√(1+y'²)) = 0.

d/dx(y'/√(1+y'²)) = 1/λ → y'/√(1+y'²) = x/λ + C.

Из граничного условия y'(0) нечётный: C = 0. Решение: y = √(λ² − x²) — полуокружность радиуса λ.

Из условия длины: πλ = L → λ = L/π, площадь = πλ²/2 = L²/(2π).

Это теорема об изопериметрическом неравенстве: для кривой длиной L площадь ≤ L²/(4π), причём равенство достигается для окружности.

Общая постановка с несколькими ограничениями

min J[y₁,...,yₙ] при K₁[y] = C₁,...,Kₘ[y] = Cₘ

Вспомогательный функционал: H = F + λ₁G₁ + ... + λₘGₘ. Система ЭЛ по каждой yᵢ + m уравнений для нахождения m множителей λⱼ из условий Kⱼ[y] = Cⱼ.

Полный разбор: задача Эйлера о колонне

Задача: найти форму упругой колонны, несущей максимальную осевую нагрузку P до потери устойчивости. Длина колонны фиксирована.

Функционал энергии: J = ∫₀ᴸ (EI/2) κ²(s) ds, где κ — кривизна, EI — жёсткость. Ограничение: ∫₀ᴸ 1 ds = L (длина колонны).

Задача Штурма-Лиувилля: EI y'' + P y = 0. Критические нагрузки: Pₙ = n²π²EI/L² (формула Эйлера).

Наименьшая критическая нагрузка (потеря устойчивости): P₁ = π²EI/L² — критическая нагрузка Эйлера. Это минимальная сила, при которой прямолинейная конфигурация перестаёт быть минимумом энергии.

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

Принцип Лагранжа: общая формулировка

Изопериметрическое правило обобщается на широкий класс задач условной оптимизации функционалов. Принцип Лагранжа в вариационном исчислении: для задачи минимизации J[y] при ограничениях K_j[y] = c_j (j = 1,...,m) необходимым условием оптимальности является существование чисел λ_j таких, что y* — экстремаль вспомогательного функционала J[y] + Σ λ_j K_j[y]. Множители Лагранжа имеют смысл «теневых цен»: ∂J*/∂c_j = −λ_j (производная оптимума по правой части ограничения).

Условия второго порядка с ограничениями

Чтобы экстремаль действительно была минимумом, нужно проверять знакоопределённость второй вариации на касательном пространстве к ограничениям — наборе вариаций η, удовлетворяющих линеаризованным ограничениям. Это аналог приведённой матрицы Гессе в конечномерной оптимизации.

Современные изопериметрические задачи

• Спектральная оптимизация: какая форма мембраны имеет наименьшее основное собственное значение лапласиана при заданной площади? Ответ — круг (Faber-Krahn inequality)
• Изопериметрические неравенства в высоких размерностях: в ℝⁿ объём шара V(r) даёт минимальную поверхность Σ при заданном объёме (теорема Брунна-Минковского)
• Изопериметрия на многообразиях: на сфере геодезические шары минимизируют площадь границы (Леви, Громов)
• Изопериметрия в теории графов (expansion): минимизация числа рёбер на границе подмножества вершин при заданном размере — основа теории случайных блужданий и спектральных методов

Применения

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