logo

편미분방정식에서 라그랑지안과 오일러-라그랑주 방정식 📂편미분방정식

편미분방정식에서 라그랑지안과 오일러-라그랑주 방정식

정의1

  • 라그랑지안Lagrangian

    스무스 함수 L:Rn×RnRL : \mathbb{R}^n \times \mathbb{R}^n \to \mathbb{R}이 주어졌다고 하자. 이를 라그랑지안 이라 부르고 다음과 같이 표기한다.

    L=L(v,x)=L(v1,,vn,x1,,xn)v,xRnDvL=(Lv1,,Lvn),DxL=(Lx1,,Lxn) L = L(v,x)=L(v_{1}, \dots, v_{n}, x_{1}, \dots, x_{n}) \quad v,x\in \mathbb{R}^{n} \\ D_{v}L = (L_{v_{1}}, \dots, L_{v_{n}}), \quad D_{x}L = (L_{x_{1}}, \dots, L_{x_{n}})

    여기서 변수를 v,xv, x로 쓰는 이유는 실제로 물리학에서 각 변수가 속도와 위치를 의미하기 때문이다.

  • 액션action, 어드미서블 클래스admissible class

    고정된 두 점 x,yRnx,y \in \mathbb{R}^{n}와 시간 t>0t>0에 대해서 다음과 같이 정의된 범함수 II액션 이라 한다.

    I[w()]:=0tL(w˙(s),w(s))ds( ˙=dds) I[ \mathbf{w}(\cdot)] := \int_{0}^tL(\dot{\mathbf{w}}(s), \mathbf{w}(s) ) ds \quad \left( \dot{ }=\dfrac{d}{ds}\right)

    이 때 함수 w()=(w1(),,wn())\mathbf{w}(\cdot)=\big( w^1(\cdot), \cdots, w^n(\cdot) \big)는 다음과 같이 정의되는 어드미서블 클래스 A\mathcal{A}의 원소이다.

    A:={w()=C2([0,t];Rn)  w(0)=y,w(t)=x} \mathcal{A} := \left\{ \mathbf{w}(\cdot) = \in C^2 \big([0,t];\mathbb{R}^n \big) \ \big| \ \mathbf{w}(0)=y, \mathbf{w}(t)=x\right\}

    다시 말해 A\mathcal{A}는 시간이 00부터 tt까지 흐를 때, 위치가 yy에서 시작해 xx로 끝나는, 2번 연속적으로 미분가능한, 모든 경로를 모아놓은 집합을 의미힌다.

설명

변분법calculus of variations의 목적은 액션 II의 적분값이 최소가 되게 하는 xA\mathbf{x} \in \mathcal{A}를 찾는 것이다. 이때의 x\mathbf{x}II미니마이저minimiaer라고 한다.

I[x()]=infw()AI[w()] I[ \mathbf{x} (\cdot) ] = \inf_{\mathbf{w}(\cdot)\in \mathcal{A}} I[\mathbf{w}(\cdot)]

이러한 x\mathbf{x}를 왜 찾느냐하면 라그랑지안의 액션이 최소가 되는 경로가 실제로 물체가 운동하는 경로이기 때문이다. 다시 말해 물체의 운동에 대해서 알고 싶기 때문이고, 본질적으로 이는 F=maF=ma를 푸는 것과 같다. 고젼역학에서 라그랑지안은 구체적으로 운동에너지와 퍼텐셜에너지의 차이로 주어진다.

미니마이저의 판별에 대해서 다음의 정리가 있다.

정리

x()A\mathbf{x}(\cdot) \in \mathcal{A}를 액션 II의 미니마이저라고 가정하자. 그러면 x()\mathbf{x}(\cdot)는 아래의 식을 만족한다.

dds[DvL(x˙(s),x(s))]+DxL(x˙(s),x(s))=0(0st) -\dfrac{d}{ds} \Big[ D_{v}L\big( \dot{\mathbf{x}}(s), \mathbf{x}(s) \big) \Big] + D_{x}L\big( \dot{\mathbf{x}}(s), \mathbf{x}(s)\big)=0 \quad (0 \le s \le t)

위 식을 오일러-라그랑주 방정식Euler-Lagrange equations이라 한다.


주의해야할 점은 미니마이저는 오일러-라그랑주 방정식을 만족하지만, 오일러-라그랑주 방정식을 만족한다고해서 미니마이저인 것은 아니라는 것이다. 최솟값을 가지는 점에서 미분하면 00이지만 미분해서 00인 점이라고해서 최솟값을 가지는 점은 아닌 것과 같다. 이같은 센스로 오일러-라그랑주 방정식을 만족하는 x()A\mathbf{x}(\cdot) \in \mathcal{A}II극점ciritical point 이라 부른다. 따라서 미니마이저는 극점이지만 극점이라고 해서 미니마이저인 것은 아니다.

증명

xA\mathbf{x} \in \mathcal{A}를 액션 II의 미니마이저라고 가정하자.

  • Step 1.

    함수 y:[0,t]Rn,y()=(y1(),,yn())\mathbf{y} : [0,t] \to \mathbb{R}^{n}, \mathbf{y}(\cdot) = (y^1(\cdot), \cdots, y^n(\cdot) )가 아래의 식을 만족하는 스무스 함수라고 하자.

    y(0)=y(t)=0 \begin{equation} \mathbf{y}(0)=\mathbf{y}(t)=\mathbf{0} \label{eq1} \end{equation}

    그리고 임의의 τR\tau \in \mathbb{R}에 대해서 w()\mathbf{w}(\cdot)를 다음과 같이 정의하자.

    w():=x()+τy()A \mathbf{w}(\cdot) : = \mathbf{x}(\cdot) + \tau \mathbf{y}(\cdot) \in \mathcal{A}

    그러면 w\mathbf{w}x\mathbf{x}와 시점, 종점의 값은 같지만 그 사이에서는 τy()\tau \mathbf{y}(\cdot)만큼 차이가 나는 경로를 갖는다. 또한 x()\mathbf{x}(\cdot)II의 미니마이저이므로 아래의 식이 성립한다.

    I[x()]I[w()]=I[x()+τy()]=:i(τ) I[\mathbf{x}(\cdot)] \le I[\mathbf{w}(\cdot)]=I[\mathbf{x}(\cdot) + \tau \mathbf{y}(\cdot)] =: i(\tau)

    또한 함수 ii는 미니마이저의 정의에 의해 τ=0\tau=0에서 최솟값을 가진다. 따라서 ii의 미분이 존재한다면, i(0)=0i^{\prime}(0)=0이다.

  • Step 2.

    위에서 정의한바에 따라 ii는 다음과 같다.

    i(τ)=0tL(x˙(s)+τy˙(s),x(s)+τy(s))ds i(\tau) = \int_{0} ^t L\big( \dot{\mathbf{x}}(s) + \tau \dot{\mathbf{y}}(s), \mathbf{x}(s)+ \tau \mathbf{y}(s) \big)ds

    LL, y\mathbf{y}는 스무스 함수이므로 ii의 미분을 구하면 다음과 같다.

    i(τ)=0ti=1n[Lvi(x˙+τy˙,x+τy)y˙i+Lxi(x˙+τy˙,x+τy)yi]ds i^{\prime}(\tau) = \int_{0}^t \sum_{i=1}^{n} \left[ L_{v_{i}} ( \dot{\mathbf{x}} + \tau \dot{\mathbf{y}}, \mathbf{x}+ \tau \mathbf{y} )\dot{y}^i + L_{x_{i}} ( \dot{\mathbf{x}} + \tau \dot{\mathbf{y}}, \mathbf{x}+ \tau \mathbf{y} ) y^i \right] ds

    τ=0\tau=0을 대입하면 Step 1. 의 결과에 의해 다음을 얻는다.

    0=i(0)=0ti=1n[Lvi(x˙,x)y˙i+Lxi(x˙,x)yi]ds 0=i^{\prime}(0) = \int_{0}^t \sum_{i=1}^{n} \left[ L_{v_{i}} ( \dot{\mathbf{x}} , \mathbf{x})\dot{y}^i + L_{x_{i}} ( \dot{\mathbf{x}} , \mathbf{x} )y^i \right] ds

    Lviy˙iL_{v_{i}}\dot{y}^i의 각 항에 부분적분을 취하면, 가정 y(0)=y(t)=0\mathbf{y}(0)=\mathbf{y}(t)=\mathbf{0}에 의해서 다음을 얻는다.

    0tLviy˙ids=Lviyi]0t0tddsLviyi=0tddsLviyi \int_{0}^t L_{v_{i}}\dot{y}^i ds = \left. L_{v_{i}}y^i \right]_{0}^t- \int_{0}^t \dfrac{d}{ds}L_{v_{i}}y^i=\int_{0}^t-\dfrac{d}{ds} L_{v_{i}}y^i

    따라서 다음의 식이 성립한다.

    0=i(0)=i=1n0t[ddsLvi(x˙,x) +Lxi(x˙,x)]yids 0=i^{\prime}(0) = \sum_{i=1}^{n} \int_{0}^t \left[ -\dfrac{d}{ds} L_{v_{i}} ( \dot{\mathbf{x}} , \mathbf{x})\ + L_{x_{i}} ( \dot{\mathbf{x}} , \mathbf{x} ) \right]y^i ds

Step 2. 의 결과는 (eq1)\eqref{eq1}을 만족하는 모든 스무스 함수 y:[0,t]Rn\mathbf{y} : [0,t] \to \mathbb{R}^n에 대해서 성립한다. 그러므로 각괄호 안의 값이 00이어야만 한다. 그러므로 다음이 성립한다.

ddsLvi(x˙,x)+Lxi(x˙,x)=0 -\dfrac{d}{ds} L_{v_{i}}( \dot{\mathbf{x}}, \mathbf{x} ) +L_{x_{i}}( \dot{\mathbf{x}}, \mathbf{x}) =0


이 결과로부터 해밀턴 방정식을 이끌어낼 수 있다.

같이보기


  1. Lawrence C. Evans, Partial Differential Equations (2nd Edition, 2010), p115-117 ↩︎