4.10.16Advanced Topics (Elite Level)

Isoperimetric problems — constraints (Lagrange multipliers in variational sense)

1,868 words8 min readdifficulty · medium

WHAT is an isoperimetric problem?

The name comes from the original problem: among all closed curves of fixed perimeter (K=K= length ==\ell), find the one enclosing maximum area (J=J= area). Answer: a circle. "Iso-perimetric" = "same perimeter".

Recall WHY can't we just use the plain Euler–Lagrange equation?

Because plain E–L finds the stationary point over all admissible functions. Here we are only allowed to wander on the surface K[y]=K[y]=\ell. Free E–L would walk off that surface. We need to restrict variations to those that keep KK fixed — that's exactly what the multiplier handles.


HOW to derive the rule (from first principles)

We want δJ=0\delta J=0 for variations δy\delta y that also keep KK fixed, i.e. δK=0\delta K=0.

Step 1 — Use a TWO-parameter family. A single ϵη(x)\epsilon\eta(x) won't work: one knob can't both vary JJ and respect the constraint. So take

y(x)=yˉ(x)+ϵ1η1(x)+ϵ2η2(x),y(x)=\bar y(x)+\epsilon_1\eta_1(x)+\epsilon_2\eta_2(x),

with η1(a)=η1(b)=η2(a)=η2(b)=0\eta_1(a)=\eta_1(b)=\eta_2(a)=\eta_2(b)=0 (boundary values are fixed, so the bumps vanish at the ends).

Why this step? Two parameters give us enough freedom: we'll spend one degree of freedom satisfying the constraint and still have one left to demand stationarity.

Step 2 — Reduce to ordinary calculus. Now JJ and KK become ordinary functions of two numbers: Φ(ϵ1,ϵ2)=J[y],Ψ(ϵ1,ϵ2)=K[y].\Phi(\epsilon_1,\epsilon_2)=J[y],\qquad \Psi(\epsilon_1,\epsilon_2)=K[y]. At yˉ\bar y we sit at ϵ1=ϵ2=0\epsilon_1=\epsilon_2=0. We want Φ\Phi stationary subject to Ψ=\Psi=\ell. This is now a finite-dimensional constrained extremum — solved by ordinary Lagrange multipliers!

Why this step? We've converted a hard calculus-of-variations constraint into the familiar 2-variable multiplier problem we already know.

Step 3 — Apply ordinary multiplier rule. There exists λ\lambda such that (ΦλΨ)=0at (0,0),\nabla(\Phi-\lambda\Psi)=0 \quad\text{at }(0,0), i.e. ϵi(ΦλΨ)0=0,i=1,2.\frac{\partial}{\partial\epsilon_i}\big(\Phi-\lambda\Psi\big)\Big|_{0}=0,\quad i=1,2.

Step 4 — Compute the derivatives. Using fϵi=fyηi+fyηi\dfrac{\partial f}{\partial\epsilon_i}=f_y\eta_i+f_{y'}\eta_i' and integrating fyηif_{y'}\eta_i' by parts (boundary terms die because ηi\eta_i vanish at ends):

Φϵi0=ab(fyddxfy)ηidx,\frac{\partial\Phi}{\partial\epsilon_i}\Big|_0=\int_a^b\Big(f_y-\frac{d}{dx}f_{y'}\Big)\eta_i\,dx, and similarly for Ψ\Psi with gg. So the condition ϵi(ΦλΨ)=0\partial_{\epsilon_i}(\Phi-\lambda\Psi)=0 reads ab[(fyddxfy)λ(gyddxgy)]ηidx=0.\int_a^b\Big[(f_y-\tfrac{d}{dx}f_{y'})-\lambda(g_y-\tfrac{d}{dx}g_{y'})\Big]\eta_i\,dx=0.

Step 5 — Fundamental Lemma. Define F=fλgF=f-\lambda g. The bracket is (FyddxFy)\big(F_y-\frac{d}{dx}F_{y'}\big). Since η2\eta_2 is essentially arbitrary (after η1\eta_1 is used to fix the constraint), the fundamental lemma of CoV forces the bracket to vanish everywhere.

Figure — Isoperimetric problems — constraints (Lagrange multipliers in variational sense)

Worked Example 1 — Catenary: hanging chain of fixed length

A chain of fixed length \ell hangs between two posts. It settles to minimise potential energy \propto height of centre of mass, i.e. minimise J=aby1+y2dxsubject to lengthK=ab1+y2dx=.J=\int_a^b y\sqrt{1+y'^2}\,dx \quad\text{subject to length}\quad K=\int_a^b\sqrt{1+y'^2}\,dx=\ell.

Step 1 — Augment. f=y1+y2f=y\sqrt{1+y'^2}, g=1+y2g=\sqrt{1+y'^2}, so F=fλg=(yλ)1+y2.F=f-\lambda g=(y-\lambda)\sqrt{1+y'^2}. Why? FF folds the constraint in; now we solve unconstrained.

Step 2 — Use Beltrami (since FF has no explicit xx). FyFy=CF-y'F_{y'}=C. Compute Fy=(yλ)y1+y2F_{y'}=(y-\lambda)\frac{y'}{\sqrt{1+y'^2}}, so (yλ)1+y2(yλ)y21+y2=yλ1+y2=C.(y-\lambda)\sqrt{1+y'^2}-(y-\lambda)\frac{y'^2}{\sqrt{1+y'^2}}=\frac{y-\lambda}{\sqrt{1+y'^2}}=C. Why Beltrami? It's a first integral that saves us one integration whenever xx is absent.

Step 3 — Solve. yλ=C1+y2y-\lambda=C\sqrt{1+y'^2}. Separate and integrate to get y(x)=λ+Ccosh ⁣(xx0C).y(x)=\lambda+C\cosh\!\Big(\frac{x-x_0}{C}\Big). The shape is a catenary (a cosh\cosh). Constants λ\lambda (vertical shift), CC, x0x_0 are fixed by the two endpoints and the length condition 1+y2=\int\sqrt{1+y'^2}=\ell.


Worked Example 2 — Dido's problem (maximum area, fixed perimeter)

Among curves y(x)y(x) from (a,0)(-a,0) to (a,0)(a,0) of fixed arc-length \ell, maximise enclosed area J=aaydxsubject toK=aa1+y2dx=.J=\int_{-a}^{a}y\,dx\quad\text{subject to}\quad K=\int_{-a}^{a}\sqrt{1+y'^2}\,dx=\ell.

Step 1 — Augment. F=yλ1+y2F=y-\lambda\sqrt{1+y'^2}. Why? f=yf=y, g=1+y2g=\sqrt{1+y'^2}.

Step 2 — Beltrami (xx absent): FyFy=CF-y'F_{y'}=C. Fy=λy1+y2F_{y'}=-\lambda\frac{y'}{\sqrt{1+y'^2}}, so yλ1+y2+λy21+y2=yλ1+y2=C.y-\lambda\sqrt{1+y'^2}+\lambda\frac{y'^2}{\sqrt{1+y'^2}}=y-\frac{\lambda}{\sqrt{1+y'^2}}=C.

Step 3 — Solve. λ1+y2=yC\dfrac{\lambda}{\sqrt{1+y'^2}}=y-C. Let u=yCu=y-C. Then 1+y2=λ2/u2y2=(λ2u2)/u21+y'^2=\lambda^2/u^2\Rightarrow y'^2=(\lambda^2-u^2)/u^2, giving uduλ2u2=±dxλ2u2=±(xx0).\frac{u\,du}{\sqrt{\lambda^2-u^2}}=\pm dx \Rightarrow -\sqrt{\lambda^2-u^2}=\pm(x-x_0). Square: (xx0)2+(yC)2=λ2(x-x_0)^2+(y-C)^2=\lambda^2. A circular arc! Why this confirms intuition: fixed perimeter → maximum area → circle, exactly the classic isoperimetric result. The multiplier λ\lambda turned out to be the radius.


Common Mistakes (Steel-manned)


Recall Feynman: explain to a 12-year-old

Imagine you have a fixed length of string and you want to enclose the biggest possible playground. You try squares, triangles, blobs… and it turns out a perfect circle wins every time. Now, how does math find that? There's a clever trick: instead of carefully keeping your string length exactly right while you search (hard!), you pretend there's a "fine" λ\lambda for every extra bit of string used. You add this fine to the thing you're maximising, then you're allowed to search freely. After finding the best shape, you tune the fine λ\lambda until the string length comes out exactly right. That fine is the Lagrange multiplier, and the free search is just the ordinary Euler–Lagrange equation.


Flashcards

What is an isoperimetric problem?
Make a functional J[y]=fdxJ[y]=\int f\,dx stationary subject to another integral functional K[y]=gdx=K[y]=\int g\,dx=\ell being fixed.
The augmented integrand for an isoperimetric problem is?
F=fλgF=f-\lambda g, and you solve FyddxFy=0F_y-\frac{d}{dx}F_{y'}=0.
Why can't plain Euler–Lagrange be used directly with a constraint?
It searches over all functions, including those violating K[y]=K[y]=\ell; we must restrict to constraint-preserving variations.
Why do we need a TWO-parameter variation yˉ+ϵ1η1+ϵ2η2\bar y+\epsilon_1\eta_1+\epsilon_2\eta_2?
One knob can't both keep KK fixed and probe stationarity; two parameters reduce it to ordinary 2-variable Lagrange multipliers.
How many unknowns and how many conditions in the final ODE?
Three unknowns (C1,C2,λC_1,C_2,\lambda); three conditions (two boundary conditions + the constraint integral).
Solution shape of: minimise y1+y2\int y\sqrt{1+y'^2} with fixed length?
A catenary, y=λ+Ccoshxx0Cy=\lambda+C\cosh\frac{x-x_0}{C}.
Solution shape of: maximise enclosed area with fixed perimeter?
A circular arc — the isoperimetric inequality's equality case is the circle.
Physical meaning of the multiplier λ\lambda?
Shadow price / marginal value: λ=dJ/d\lambda=dJ^*/d\ell, the rate of change of the optimum with the constraint budget.
Integral constraint vs pointwise constraint multiplier?
Integral (isoperimetric) → single constant λ\lambda; pointwise g(x,y)=0g(x,y)=0 → multiplier function λ(x)\lambda(x).
Which first integral simplifies these problems when xx is absent from FF?
Beltrami identity FyFy=CF-y'F_{y'}=C.

Connections

Concept Map

optimise

subject to

classic case

handled by

forms

derive via

reduces to

introduces

yields

fixes

solve for

Isoperimetric problem

Maximise J of y

Constraint K of y equals l

Add lambda times constraint

Auxiliary h equals f minus lambda g

Two-parameter family

Ordinary Lagrange multiplier

Euler-Lagrange on h

Optimal y and lambda

Fixed perimeter gives circle

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, isoperimetric problem ka matlab simple hai: tumhe ek functional J[y]=fdxJ[y]=\int f\,dx ko maximise ya minimise karna hai, lekin ek constraint ke saath — doosra integral K[y]=gdxK[y]=\int g\,dx fixed rehna chahiye (jaise rope ki length fixed, ya perimeter fixed). Classic example: utni hi fence se sabse bada khet banao — answer hamesha circle aata hai.

Ab trick yeh hai — constraint ko zabardasti maintain karte hue search karna mushkil hai. Toh hum bribe karte hain: ek number λ\lambda (Lagrange multiplier) leke, naya integrand banate hain F=fλgF=f-\lambda g. Iske baad bilkul normal Euler–Lagrange equation laga do: FyddxFy=0F_y-\frac{d}{dx}F_{y'}=0. Bas! Constraint apne aap handle ho jaata hai. Yeh same Lagrange multiplier idea hai jo tumne functions ke max-min me dekha tha, bas yahan points ki jagah poore function pe kaam kar raha hai.

Counting bhi yaad rakho: solution me teen unknown aayenge — do constants (C1,C2C_1,C_2) aur ek λ\lambda. Inhe fix karne ke liye teen condition: do boundary condition aur ek constraint gdx=\int g\,dx=\ell. λ\lambda ka physical matlab "shadow price" hai — agar tum budget \ell thoda badhao, toh optimum kitna improve hoga, woh rate λ\lambda batata hai.

Do mast examples: (1) chain hang karti hai fixed length ke saath, energy minimise karke shape catenary (cosh\cosh) banta hai. (2) Fixed perimeter me area maximise karo, toh circle (arc) milta hai — multiplier λ\lambda to seedha radius nikalta hai! Jab bhi FF me explicit xx na ho, Beltrami identity FyFy=CF-y'F_{y'}=C use karke ek integration bacha lo. Mnemonic: Bribe, Apply EL, Boundary, Constraint = BAB-C.

Go deeper — visual, from zero

Test yourself — Advanced Topics (Elite Level)

Connections