Анализ функций нескольких переменных · Множитель Лагранжа

Множитель Лагранжа

Чтобы найти наибольшее или наименьшее значение f вдоль кривой ограничения g = c, ищите точки, где при ∇g ≠ 0 векторы ∇f и ∇g параллельны: ∇f = λ∇g. Здесь λ — скорость изменения оптимального значения при изменении c. Перемещайте P по кривой и следите за составляющей ∇f вдоль неё.

Множитель Лагранжа в этой визуализации

Нужно найти наибольшее или наименьшее значение функции f(x,y)f(x, y), когда точка не свободна, а обязана лежать на кривой ограничения g(x,y)=cg(x, y) = c. Метод множителей Лагранжа указывает, где искать. Пусть у ff и gg непрерывные первые частные производные в окрестности точки P на кривой, ограничение ff на кривую имеет в P локальный максимум или минимум и ∇g(P)≠0\nabla g(P) \neq 0. Тогда существует число λ\lambda — множитель Лагранжа, для которого

Уравнения для точек-кандидатов
∇f(x,y)=λ ∇g(x,y),g(x,y)=c\nabla f(x, y) = \lambda\, \nabla g(x, y), \qquad g(x, y) = c

Это три уравнения для x, y и λ\lambda. Их решения — кандидаты, а не готовые ответы: каждого ещё нужно проверить, а точки с ∇g=0\nabla g = 0 исследовать отдельно.

Карта показывает линии уровня ff, кривую ограничения жирной линией и точку P на ней. В P синяя стрелка — ∇f\nabla f, полая стрелка показывает направление ∇g\nabla g, а ∇f\nabla f разложен на жёлтую составляющую вдоль кривой и пунктирную поперёк неё. Под картой профиль показывает ff в зависимости от пройденного вдоль кривой расстояния.

Почему это работает? Вдоль кривой gg постоянно, поэтому единичная касательная TT перпендикулярна ∇g\nabla g. При движении с единичной скоростью ff меняется со скоростью ∇f⋅T\nabla f\cdot T — это наклон профиля в P, а жёлтая стрелка равна (∇f⋅T) T(\nabla f\cdot T)\,T. В наивысшей или наинизшей точке профиля, где ∇g≠0\nabla g \neq 0, этот наклон нулевой, поэтому у ∇f\nabla f нет составляющей вдоль кривой. На плоскости такой вектор кратен ∇g\nabla g, значит, ∇f=λ∇g\nabla f = \lambda\nabla g; при ∇f=0\nabla f = 0 имеем λ=0\lambda = 0. Там, где жёлтая стрелка ненулевая, движение по ней ещё увеличивает ff, поэтому P не может быть ни максимумом, ни минимумом.

Если также ∇f(P)≠0\nabla f(P) \neq 0, то линия уровня ff через P перпендикулярна ∇f\nabla f и имеет общую касательную с ограничением. В условном максимуме или минимуме кривая ограничения касается линии уровня, не переходя через неё.

У λ\lambda есть собственный смысл. Обозначим через f∗(c)f^*(c) оптимальное значение на кривой g=cg = c. Если точка, в которой оно достигается, плавно движется при изменении cc, то λ\lambda — скорость изменения этого значения:

Множитель как скорость изменения
λ=df∗dc\lambda = \frac{d f^*}{d c}

На окружности x2+y2=cx^2 + y^2 = c наибольшее значение xyxy равно c/2c/2, а в обоих максимумах λ=1/2\lambda = 1/2. Окно λ показывает f∗f^* в зависимости от cc с этим наклоном. Скорость зависит от записи ограничения. Если записать 2g=C2g = C при C=2cC = 2c, то для той же кривой множитель уменьшится вдвое: теперь он измеряет изменение на единицу C. Кнопка Записать как 2g показывает это.

Две кривые показывают границы метода. На y−x3=cy - x^3 = c у функции f=yf = y в точке (0,c)(0, c) градиенты параллельны, но вдоль кривой f=x3+cf = x^3 + c лишь выравнивается и продолжает расти. Параллельности недостаточно. На y2−x3=0y^2 - x^3 = 0 наименьшее значение f=xf = x достигается в острие (0,0)(0, 0), где ∇g=(0,0)\nabla g = (0, 0), а ∇f=λ∇g\nabla f = \lambda\nabla g не имеет решения. Метод его пропускает, поэтому такие точки нужно проверять отдельно.

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

На что обратить внимание

  • Знак λ\lambda показывает, направлен ли ∇f\nabla f в одну сторону с ∇g\nabla g или в противоположную. На окружностях он положителен в максимумах и отрицателен в минимумах.
  • Касание линии уровня не означает, что кривая остаётся с одной стороны от неё. В точке-кандидате кубическая кривая касается линии уровня и всё же пересекает её.

Связанные понятия