7.62 최단 강하선 문제와 등주 문제

1. 최단 강하선 문제

1.1 문제의 정식화

최단 강하선 문제(brachistochrone problem)는 중력장에서 마찰 없는 곡선을 따라 질점이 점 A = (0, 0)에서 점 B = (x_f, y_f)까지 활강할 때, 이동 시간을 최소화하는 곡선을 찾는 문제이다. y축을 아래 방향 양으로 설정한다.

에너지 보존에 의하여 높이 y에서의 속도는 v = \sqrt{2gy}이다. 곡선 y = y(x)를 따른 미소 호의 길이는 ds = \sqrt{1 + y'^2}\,dx (y' = dy/dx)이므로, 총 이동 시간은 다음의 범함수이다.

T[y] = \int_0^{x_f} \frac{\sqrt{1 + y'^2}}{\sqrt{2gy}}\,dx

이 범함수를 최소화하는 함수 y(x)를 찾아야 한다.

1.2 오일러-라그랑주 방정식의 적용

라그랑지안 L(y, y') = \sqrt{(1 + y'^2)/(2gy)}x를 명시적으로 포함하지 않으므로, 벨트라미 항등식을 적용한다.

L - y'\frac{\partial L}{\partial y'} = C

이를 계산하여 정리하면 다음을 얻는다.

\frac{1}{\sqrt{y(1 + y'^2)}} = C_1 \quad \text{(상수)}

여기서 C_1 = 1/\sqrt{2R}으로 놓으면 y(1 + y'^2) = 2R이 된다.

1.3 사이클로이드 해

위의 미분 방정식을 매개변수 \theta를 도입하여 풀면, 해는 사이클로이드(cycloid)이다.

x(\theta) = R(\theta - \sin\theta), \quad y(\theta) = R(1 - \cos\theta)

이는 반지름 R인 원이 x축을 따라 구를 때, 원 위의 한 점이 그리는 궤적이다. 상수 R은 경계 조건 y(x_f) = y_f에 의해 결정된다.

1.4 물리적 의미

사이클로이드는 직선보다 긴 경로이지만, 초반에 급격히 하강하여 빠르게 속도를 얻으므로 전체 이동 시간이 최소가 된다. 이 문제는 로봇 궤적 최적화에서 “최적 경로가 반드시 최단 경로는 아니다“라는 원리를 보여주는 고전적 예시이다.

2. 등주 문제

2.1 문제의 정식화

등주 문제(isoperimetric problem)는 적분 제약 조건(integral constraint) 하에서 범함수를 최적화하는 문제이다. 기본 형태는 다음과 같다.

\min_{y} J[y] = \int_{t_0}^{t_f} L(t, y, \dot{y})\,dt

제약 조건:

K[y] = \int_{t_0}^{t_f} M(t, y, \dot{y})\,dt = \ell \quad \text{(주어진 상수)}

경계 조건: y(t_0) = y_0, y(t_f) = y_f

2.2 고전적 등주 문제

이 문제 유형의 원형은 다음과 같다: 주어진 둘레 \ell을 갖는 평면 폐곡선 중에서 내부 넓이를 최대화하는 곡선을 찾아라. 답은 원이며, 이로부터 “등주 문제(isoperimetric problem)“라는 명칭이 유래하였다.

2.3 라그랑주 승수법에 의한 풀이

등주 문제는 라그랑주 승수(Lagrange multiplier) \lambda를 도입하여 비제약 문제로 변환한다. 증강된 범함수(augmented functional)를 다음과 같이 정의한다.

\bar{J}[y] = J[y] + \lambda(K[y] - \ell) = \int_{t_0}^{t_f} [L + \lambda M]\,dt - \lambda\ell

\lambda는 상수이므로, 증강된 라그랑지안 \bar{L} = L + \lambda M에 대하여 오일러-라그랑주 방정식을 적용한다.

\frac{\partial \bar{L}}{\partial y} - \frac{d}{dt}\frac{\partial \bar{L}}{\partial \dot{y}} = 0

이 방정식과 제약 조건 K[y] = \ell을 연립하여 극값 함수 y^*(t)와 라그랑주 승수 \lambda를 동시에 결정한다.

2.4 다중 적분 제약

m개의 적분 제약 조건이 존재하는 경우, m개의 라그랑주 승수 \lambda_1, \lambda_2, \ldots, \lambda_m을 도입한다.

\bar{L} = L + \sum_{k=1}^{m} \lambda_k M_k

3. 로봇공학에서의 등주 문제 응용

3.1 에너지 제약 하의 최소 시간 궤적

로봇 궤적 최적화에서 전형적인 등주 문제는 다음과 같다: 총 에너지 소비가 주어진 상한 E_{max}를 초과하지 않으면서 이동 시간을 최소화하라.

\min J = \int_{t_0}^{t_f} 1\,dt = t_f - t_0

\text{subject to} \quad \int_{t_0}^{t_f} \boldsymbol{\tau}^T\boldsymbol{\tau}\,dt \leq E_{max}

3.2 궤적 평활도 제약

경로 길이를 제한하면서 궤적의 평활도를 최대화하는 문제도 등주 문제의 형태를 갖는다.

\min J = \int_{t_0}^{t_f} \lVert\ddot{\mathbf{q}}\rVert^2\,dt

\text{subject to} \quad \int_{t_0}^{t_f} \lVert\dot{\mathbf{q}}\rVert\,dt = L \quad \text{(경로 길이)}

4. 최단 거리 문제

4.1 문제와 풀이

두 점을 잇는 최단 거리 곡선을 찾는 문제는 가장 단순한 변분 문제이다. 유클리드 평면에서 두 점 (x_0, y_0)(x_f, y_f) 사이의 곡선 길이는 다음과 같다.

J[y] = \int_{x_0}^{x_f} \sqrt{1 + y'^2}\,dx

라그랑지안 L = \sqrt{1 + y'^2}xy를 포함하지 않으므로, \partial L/\partial y = 0에서 다음을 얻는다.

\frac{\partial L}{\partial y'} = \frac{y'}{\sqrt{1 + y'^2}} = C

이로부터 y' = \text{const}, 즉 해는 직선 y = ax + b이다. 유클리드 공간에서 두 점 사이의 최단 경로가 직선이라는 자명한 결과가 변분법에 의해 엄밀하게 증명된다.

4.2 측지선으로의 확장

곡면 위에서의 최단 거리 문제는 측지선(geodesic) 문제로 확장된다. 로봇 매니퓰레이터의 형상 공간(configuration space)이 비유클리드 다양체인 경우, 관절 공간에서의 최단 경로는 리만 측지선이 된다. 관성 행렬이 리만 계량(Riemannian metric)의 역할을 하며, 이때의 측지선은 자유 운동(무중력, 무마찰)에서의 자연스러운 궤적에 해당한다.

5. 요약

최단 강하선 문제와 등주 문제는 변분법의 고전적이고 핵심적인 응용이다. 최단 강하선 문제의 사이클로이드 해는 벨트라미 항등식을 통해 유도되며, 최적 경로가 최단 경로와 다를 수 있음을 보여준다. 등주 문제는 적분 제약 조건 하의 범함수 최적화 문제로, 라그랑주 승수를 도입하여 비제약 문제로 환원하여 풀 수 있다. 이러한 고전적 문제의 수학적 구조는 로봇 궤적 최적화에서 에너지 제약, 경로 길이 제약, 평활도 제약 등 다양한 제약 조건을 갖는 최적화 문제의 이론적 기초를 형성한다.


참고 문헌

  • Gelfand, I. M. & Fomin, S. V. (2000). Calculus of Variations. Dover Publications.
  • Liberzon, D. (2012). Calculus of Variations and Optimal Control Theory: A Concise Introduction. Princeton University Press.
  • Goldstein, H., Poole, C., & Safko, J. (2002). Classical Mechanics (3rd ed.). Addison-Wesley.
  • Kirk, D. E. (2004). Optimal Control Theory: An Introduction. Dover Publications.
  • Sussmann, H. J. & Willems, J. C. (1997). “300 years of optimal control: From the brachystochrone to the maximum principle.” IEEE Control Systems Magazine, 17(3), 32–44.

version: 0.1.0