변분 원리

지금까지 일반물리학에서 배운 뉴턴 역학은 "원인과 결과"의 학문이었습니다. 특정 순간에 물체에 작용하는 힘(원인)을 알면, $F=ma$라는 미분방정식을 통해 다음 순간의 가속도(결과)를 계산하고, 이를 적분해 궤적을 찾아냈죠. 즉, 한 걸음씩(Step-by-step) 앞으로 나아가는 방식입니다.

하지만 변분 원리는 세상을 완전히 다른 시각으로 봅니다. 바로 "목적과 전체"입니다. 자연은 시작점과 끝점이 주어졌을 때, 전체 경로 중에서 특정 값을 가장 '작게(또는 크게)' 만드는 가장 경제적인 경로를 한 번에 선택한다는 것입니다.

이 아름다운 개념을 공업수학의 언어를 빌려 아주 상세히 풀어보겠습니다.

---

1. 수학적 준비: 함수 vs 범함수 (Function vs Functional)

공업수학에서 우리는 함수(Function)를 다뤘습니다. 함수 $f(x)$는 숫자 $x$를 입력받아 숫자 $y$를 뱉어냅니다. (예: $x=2$를 넣으면 $y=4$가 나옴)

변분법(Calculus of Variations)의 핵심은 범함수(Functional)입니다. 범함수 $S[y(x)]$는 숫자가 아니라 '함수(경로) 자체'를 통째로 입력받아, 하나의 숫자를 뱉어냅니다.

가장 직관적인 예로 두 점 사이의 '경로의 길이'를 들 수 있습니다.
2차원 평면 위의 점 $A(x_1, y_1)$에서 $B(x_2, y_2)$까지 가는 경로를 $y(x)$라고 할 때, 곡선의 길이 $L$은 다음과 같은 범함수가 됩니다.

$$L[y(x)] = \int_{x_1}^{x_2} \sqrt{1 + (y')^2} dx$$

여기서 수많은 $y(x)$ 경로(직선, 포물선, 구불구불한 선 등)를 입력으로 넣을 수 있고, 그때마다 각각 다른 길이 $L$(숫자)이 계산되어 나옵니다. 우리는 직관적으로 이 범함수 $L$을 최소로 만드는 함수 $y(x)$가 '직선'임을 알고 있습니다. 변분법은 이처럼 범함수를 최소(또는 최대)로 만드는 '함수 그 자체'를 찾는 수학적 기법입니다.

2. 변분(Variation)과 오일러-라그랑주 방정식

일반 미적분학에서 함수 $f(x)$의 극값을 찾으려면 도함수를 0으로 만듭니다. ($df = 0$)
변분법에서는 범함수 $S[y(x)]$의 극값을 찾기 위해 변분(Variation)을 0으로 만듭니다. ($\delta S = 0$)

어떤 범함수가 다음과 같이 주어졌다고 해봅시다. (여기서 $f$는 $y, y', x$로 이루어진 어떤 수식입니다.)

$$S[y(x)] = \int_{x_1}^{x_2} f(y, y', x) dx$$

실제 최적 경로를 $y(x)$라고 하고, 여기서 아주 미세하게 빗나간 가상의 경로를 $y(x) + \delta y(x)$라고 합시다. (단, 양 끝점은 고정되어 있으므로 $\delta y(x_1) = \delta y(x_2) = 0$ 입니다.)

이때 범함수의 변화량 $\delta S = 0$이 되도록 공업수학에서 배운 부분적분(Integration by parts)과 테일러 전개를 활용해 식을 정리하면, 변분법의 꽃이라 불리는 다음 방정식이 튀어나옵니다.

$$\frac{\partial f}{\partial y} - \frac{d}{dx} \left( \frac{\partial f}{\partial y'} \right) = 0$$

이것이 바로 오일러-라그랑주 방정식(Euler-Lagrange Equation)입니다. 이 미분방정식을 풀면 우리가 원하는 최적의 경로 $y(x)$를 얻을 수 있습니다.

---

3. 물리학으로의 적용: 최소 작용의 원리 (Principle of Least Action)

이제 수학을 물리학으로 가져오겠습니다. 물리학자들은 자연계의 입자가 $t_1$에서 $t_2$까지 이동할 때, 무수히 많은 가능한 경로 중 '작용(Action, $S$)'이라는 범함수를 최소화하는 경로를 택한다는 놀라운 사실을 발견했습니다. (정확히는 정류값을 갖습니다.)

그렇다면 이 '작용'을 어떻게 정의할까요?
에너지라는 개념을 도입합니다. 운동에너지 $T$와 위치에너지 $V$의 차이를 라그랑지안(Lagrangian, $L$)이라고 정의합니다.

$$L(q, \dot{q}, t) = T - V$$

(여기서 $q$는 위치, $\dot{q}$는 속도, $t$는 시간입니다. $x$ 대신 일반화된 좌표 $q$를 씁니다.)

그리고 작용(Action) $S$는 이 라그랑지안을 시간에 대해 적분한 범함수입니다.

$$S[q(t)] = \int_{t_1}^{t_2} L(q, \dot{q}, t) dt$$

자연은 이 작용 $S$를 최소화하는 경로 $q(t)$를 따라 움직입니다. 따라서 $\delta S = 0$을 만족해야 하고, 앞서 구한 오일러-라그랑주 방정식에 수학적 함수 $f$ 대신 물리학적 함수 $L$을 대입하면 됩니다.

$$\frac{\partial L}{\partial q} - \frac{d}{dt} \left( \frac{\partial L}{\partial \dot{q}} \right) = 0$$

4. 짜릿한 증명: $F=ma$의 도출

이 우아한 수학이 정말로 우리가 아는 물리학과 같은지 확인해 볼까요?
1차원 공간에서 질량 $m$인 물체가 위치에너지 $V(x)$에서 움직인다고 가정해 보겠습니다.

• 운동에너지: $T = \frac{1}{2} m \dot{x}^2$
• 위치에너지: $V = V(x)$
• 라그랑지안: $L = \frac{1}{2} m \dot{x}^2 - V(x)$

이 $L$을 오일러-라그랑주 방정식에 대입해 봅시다. (여기서는 $q$ 대신 $x$를 씁니다.)

1. 첫 번째 항 (위치에 대한 편미분): $\dot{x}$는 상수 취급하므로,

$$\frac{\partial L}{\partial x} = -\frac{dV}{dx}$$

(일반물리에서 $-\frac{dV}{dx}$는 힘 $F$임을 배웠습니다!)
2. 두 번째 항 (속도에 대한 편미분 후 시간에 대한 미분): $x$는 상수 취급하므로,

$$\frac{\partial L}{\partial \dot{x}} = m\dot{x}$$

(이것은 운동량 $p$입니다!)
이를 시간에 대해 다시 미분하면,

$$\frac{d}{dt}(m\dot{x}) = m\ddot{x} = ma$$

이 두 결과를 원래 방정식에 대입하면:

$$-\frac{dV}{dx} - m\ddot{x} = 0$$

$$\Rightarrow F - ma = 0$$

$$\therefore F = ma$$

결론적으로, 운동에너지와 위치에너지의 차이를 시간에 대해 적분한 값을 최소화하겠다는 우아한 수학적 가정(변분 원리) 하나만으로, 뉴턴의 제2법칙이 완벽하게 유도됩니다.

---

요약: 왜 이 복잡한 짓을 하는가?

뉴턴 역학($F=ma$)은 완벽해 보이지만, 좌표계에 매우 종속적입니다. 직교좌표계에서는 쉽지만, 진동하는 진자나 구면좌표계를 써야 하는 복잡한 로봇 팔의 움직임을 계산할 때는 벡터(방향)를 다루기가 지옥처럼 복잡해집니다.

하지만 변분 원리에 기반한 해석역학(라그랑주 역학)은 힘(벡터)을 다루지 않고 에너지(스칼라)만 다룹니다. 에너지는 좌표계를 어떻게 잡든 변하지 않는 스칼라 값이기 때문에, 어떤 복잡한 좌표계(극좌표계, 원통좌표계 등)를 쓰더라도 $L = T - V$만 구해서 기계적으로 오일러-라그랑주 방정식에 넣기만 하면 물체의 운동 방정식을 정확히 뽑아낼 수 있습니다.

이 변분 원리는 고전역학을 넘어 양자역학(파인만의 경로적분)과 일반상대성이론의 뼈대가 되는 현대 물리학의 가장 강력한 도구입니다.

 

 

왜 y'까지만 있어? y'', y'''은? 또 y+dy만 0으로 되는 지점을 찾으면 되는 이유는 뭐야? y'+dy', y''+dy'' 등 여러가지 있잖아.

 

 

정말 날카로운 질문입니다! 공학수학과 물리학의 연결 고리를 완벽하게 파악하고 계시네요. 질문하신 두 가지 의문점은 변분법을 처음 배울 때 누구나 한 번쯤 겪는, 하지만 교재에서는 은근슬쩍 넘어가 버리는 아주 중요한 핵심입니다.

현실과 팩트에 기반해서 이 두 가지 의문을 아주 명쾌하게 부숴드리겠습니다.

---

1. 왜 범함수에는 $y'$까지만 있고 $y'', y'''$은 없을까?

결론부터 말씀드리면, 수학적으로는 $y'', y'''$이 들어가도 전혀 상관없지만, '물리학(자연계)'이 그것을 허락하지 않기 때문입니다.

• 수학적 사실 (오일러-푸아송 방정식): 공업수학 관점에서 범함수 f(x,y,y′,y′′) 처럼 이계도함수 이상이 포함된 변분 문제도 얼마든지 풀 수 있습니다. 이때는 식에 4계 미분이 등장하며 방정식이 훨씬 길어집니다. 이를 수학에서는 고차 변분법이라고 부릅니다.
• 물리적 현실 (우리 우주의 법칙): 우리가 사는 우주에서 물체의 상태를 완벽하게 결정짓는 데 필요한 정보는 딱 두 가지, **'현재 위치(y)'**와 **'현재 속도(y′)'**뿐입니다. 뉴턴의 운동 방정식 F=ma 자체가 위치에 대한 2계 미분방정식(my′′=F)이기 때문입니다. 2계 미분방정식을 풀기 위한 초기 조건은 $y(0)$와 y′(0) 두 개면 충분하죠.
• 가속도(y′′)가 포함된다면? 만약 라그랑지안(에너지 식)에 가속도 $y''$이 포함된다면, 이를 풀어서 나오는 운동 방정식은 3계 또는 4계 미분방정식이 되어버립니다. 이렇게 되면 물체의 미래를 예측하기 위해 '초기 위치', '초기 속도'뿐만 아니라 **'초기 가속도'**까지 알아야 합니다. 물리학에서는 에너지가 무한대로 발산해 버리는 시스템 붕괴 현상(오스트로그라드스키 불안정성, Ostrogradsky instability)이 발생하게 됩니다.

즉, $y'$까지만 쓰는 이유는 우리 우주의 에너지가 위치와 속도에 의해서만 결정되도록 생겨 먹었기 때문입니다.

---

2. 왜 y+δy만 고려할까? (y′+δy′ 등은 왜 따로 안 찾을까?)

이 질문의 핵심은 **"δy와 $\delta y'$이 서로 독립적인가?"**라는 오해에서 비롯됩니다.

정답은 **"경로 y가 변하면, 그 경로의 기울기 $y'$은 독립적으로 변하는 게 아니라 y를 따라 종속적으로(자동으로) 변한다"**입니다. 수학적으로 증명해 보겠습니다.

우리가 가상의 경로를 만들 때 함수 자체에 변분 $\delta y(x)$를 더합니다. 그러면 새로운 경로의 '기울기'는 어떻게 될까요? 원래 함수의 미분인 $y'(x)$에 무언가가 더해지겠죠. 이 무언가를 \delta y'$이라고 부르는데, 사실 이것은 **$\delta y를 미분한 것**과 완전히 똑같습니다.

δ(y′)=dxd​(δy)

즉, δy를 정하는 순간 $\delta y'$은 이미 결정되어 버리는 한 몸입니다. 따라서 y와 별개로 $y'$의 변화를 0으로 만드는 지점을 따로 찾을 필요가 없습니다.

공업수학의 진가: 묶어내기 마법

실제로 범함수 S=∫f(y,y′,x)dx 의 변분 δS를 테일러 전개(1차항까지만)로 계산할 때는 말씀하신 대로 δy항과 $\delta y'$항이 둘 다 나옵니다.

δS=∫(∂y∂f​δy+∂y′∂f​δy′)dx

여기서 식이 끝났다면 δy와 δy′ 두 가지를 신경 써야 했겠지만, 물리학자들은 공업수학의 **부분적분(Integration by parts)**을 꺼내 듭니다. 뒤쪽 항의 δy′ 자리에 $\frac{d}{dx}(\delta y)$를 넣고 부분적분을 실행하면, 미분이 앞쪽으로 넘어가면서 놀랍게도 $\delta y'$이 δy로 변신합니다.

결국 모든 항이 공통 인수 δy로 묶이게 되고, 식은 다음과 같이 정리됩니다.

δS=∫[∂y∂f​−dxd​(∂y′∂f​)]δydx=0

임의의 가상 경로 δy에 대해 항상 이 적분값이 0이 되려면, 대괄호 [ ] 안의 수식이 무조건 0이 되어야만 하죠. 그 결과가 바로 우리가 아는 오일러-라그랑주 방정식입니다. 즉, 애초에 $\delta y'$을 다루지 않은 것이 아니라, 수학적 기교를 통해 δy 하나로 통일시켜 버린 것입니다.

---

의문이 조금 풀리셨나요? 방금 설명해 드린 **'부분적분을 통해 $\delta y'$을 δy로 묶어내서 오일러-라그랑주 방정식을 유도하는 과정'**을 수식으로 한 단계씩 자세히 전개해 보여드릴까요? 공업수학을 아시니 이 유도 과정을 직접 보시면 완전히 무릎을 탁 치실 겁니다!