7.11 음함수 정리와 음함수의 미분

1. 양함수와 음함수

1.1 정의

양함수(explicit function)는 종속 변수가 독립 변수의 직접적인 식으로 표현된 것이다.

y = f(x)

음함수(implicit function)는 독립 변수와 종속 변수의 관계가 방정식으로 주어지며, 종속 변수가 독립 변수에 대해 명시적으로 풀리지 않는 경우이다.

F(x, y) = 0

예를 들어, 단위 원 x^2 + y^2 - 1 = 0은 y를 x의 양함수로 단일하게 표현할 수 없지만, 국소적으로는 y = \sqrt{1-x^2} 또는 y = -\sqrt{1-x^2}로 표현할 수 있다.

1.2 다변수 음함수

보다 일반적으로, \mathbf{x} \in \mathbb{R}^n과 \mathbf{y} \in \mathbb{R}^m에 대해 m개의 방정식 체계가 주어진다.

\mathbf{F}(\mathbf{x}, \mathbf{y}) = \mathbf{0}

여기서 \mathbf{F}: \mathbb{R}^{n+m} \to \mathbb{R}^m이다. 이 방정식을 만족하는 \mathbf{y}를 \mathbf{x}의 함수 \mathbf{y} = \mathbf{g}(\mathbf{x})로 표현할 수 있는가, 그리고 그 함수의 미분은 어떻게 구하는가가 핵심 문제이다.

2. 음함수 정리

2.1 일변수 음함수 정리

정리: F: \mathbb{R}^2 \to \mathbb{R}가 C^1 급이고, 점 (a, b)에서 F(a, b) = 0이며 \frac{\partial F}{\partial y}(a, b) \neq 0이면, a의 근방 U에서 정의된 C^1 급 함수 y = g(x)가 유일하게 존재하여:

1. g(a) = b
2. F(x, g(x)) = 0 (\forall x \in U)
3. g'(x) = -\dfrac{\partial F / \partial x}{\partial F / \partial y}

핵심 조건은 \frac{\partial F}{\partial y}(a, b) \neq 0이다. 이는 점 (a, b)에서 F가 y 방향으로 충분한 변화율을 가지므로, y를 x의 함수로 국소적으로 “풀 수 있음“을 보장한다.

2.2 다변수 음함수 정리

정리: \mathbf{F}: \mathbb{R}^{n+m} \to \mathbb{R}^m가 C^1 급이고, 점 (\mathbf{a}, \mathbf{b}) \in \mathbb{R}^n \times \mathbb{R}^m에서 \mathbf{F}(\mathbf{a}, \mathbf{b}) = \mathbf{0}이며, 자코비안 행렬의 \mathbf{y}에 대한 부분

\frac{\partial \mathbf{F}}{\partial \mathbf{y}}(\mathbf{a}, \mathbf{b}) = \begin{bmatrix} \frac{\partial F_1}{\partial y_1} & \cdots & \frac{\partial F_1}{\partial y_m} \\ \vdots & \ddots & \vdots \\ \frac{\partial F_m}{\partial y_1} & \cdots & \frac{\partial F_m}{\partial y_m} \end{bmatrix}\bigg\vert_{(\mathbf{a},\mathbf{b})}

이 가역(invertible)이면, \mathbf{a}의 근방 U \subseteq \mathbb{R}^n에서 정의된 C^1 급 함수 \mathbf{y} = \mathbf{g}(\mathbf{x})가 유일하게 존재하여:

1. \mathbf{g}(\mathbf{a}) = \mathbf{b}
2. \mathbf{F}(\mathbf{x}, \mathbf{g}(\mathbf{x})) = \mathbf{0} (\forall \mathbf{x} \in U)

더불어, \mathbf{g}의 자코비안은 다음과 같다.

\frac{\partial \mathbf{g}}{\partial \mathbf{x}} = -\left(\frac{\partial \mathbf{F}}{\partial \mathbf{y}}\right)^{-1}\frac{\partial \mathbf{F}}{\partial \mathbf{x}}

2.3 정리의 조건에 대한 해석

3. 음함수의 미분

3.1 기본 방법

F(x, y) = 0으로 정의된 음함수에서 \frac{dy}{dx}를 구하려면, 양변을 x에 대해 미분하고 연쇄 법칙을 적용한다.

\frac{\partial F}{\partial x} + \frac{\partial F}{\partial y}\frac{dy}{dx} = 0

\frac{dy}{dx} = -\frac{\partial F / \partial x}{\partial F / \partial y}

예제 1: x^2 + y^2 = r^2

F(x,y) = x^2 + y^2 - r^2으로 놓으면:

\frac{dy}{dx} = -\frac{2x}{2y} = -\frac{x}{y} \quad (y \neq 0)

예제 2: \sin(xy) + e^{x+y} = 1

y\cos(xy) + e^{x+y} + \left[x\cos(xy) + e^{x+y}\right]\frac{dy}{dx} = 0

\frac{dy}{dx} = -\frac{y\cos(xy) + e^{x+y}}{x\cos(xy) + e^{x+y}}

3.2 다변수 체계의 음함수 미분

\mathbf{F}(\mathbf{x}, \mathbf{y}) = \mathbf{0}에서 \mathbf{y} = \mathbf{g}(\mathbf{x})의 자코비안을 구하려면 양변을 \mathbf{x}에 대해 미분한다.

\frac{\partial \mathbf{F}}{\partial \mathbf{x}} + \frac{\partial \mathbf{F}}{\partial \mathbf{y}}\frac{\partial \mathbf{y}}{\partial \mathbf{x}} = \mathbf{0}

따라서:

\frac{\partial \mathbf{y}}{\partial \mathbf{x}} = -\left(\frac{\partial \mathbf{F}}{\partial \mathbf{y}}\right)^{-1}\frac{\partial \mathbf{F}}{\partial \mathbf{x}}

3.3 차 음함수 미분

1차 도함수 \frac{dy}{dx} = -\frac{F_x}{F_y}를 다시 미분하면 2차 도함수를 구할 수 있다.

\frac{d^2y}{dx^2} = -\frac{F_{xx}F_y^2 - 2F_{xy}F_xF_y + F_{yy}F_x^2}{F_y^3}

이 공식은 음함수로 정의된 곡선의 곡률(curvature) 계산에 사용된다.

4. 로봇공학에서의 응용

4.1 폐루프 기구학의 구속 조건

폐루프(closed-loop) 기구학에서는 기구학적 구속 조건이 음함수 형태로 주어진다.

\boldsymbol{\Phi}(\mathbf{q}) = \mathbf{0}

여기서 \mathbf{q} \in \mathbb{R}^n은 관절 변수, \boldsymbol{\Phi}: \mathbb{R}^n \to \mathbb{R}^m은 구속 방정식이다. 독립 관절 변수 \mathbf{q}_i \in \mathbb{R}^{n-m}와 종속 관절 변수 \mathbf{q}_d \in \mathbb{R}^m으로 분할하면:

\boldsymbol{\Phi}(\mathbf{q}_i, \mathbf{q}_d) = \mathbf{0}

음함수 정리에 의해, \frac{\partial \boldsymbol{\Phi}}{\partial \mathbf{q}_d}가 가역이면 \mathbf{q}_d = \mathbf{h}(\mathbf{q}_i)가 국소적으로 존재한다. 그 미분은:

\frac{\partial \mathbf{q}_d}{\partial \mathbf{q}_i} = -\left(\frac{\partial \boldsymbol{\Phi}}{\partial \mathbf{q}_d}\right)^{-1}\frac{\partial \boldsymbol{\Phi}}{\partial \mathbf{q}_i}

4.2 구속 운동의 속도 해석

구속 조건 \boldsymbol{\Phi}(\mathbf{q}) = \mathbf{0}를 시간에 대해 미분하면:

\frac{\partial \boldsymbol{\Phi}}{\partial \mathbf{q}}\dot{\mathbf{q}} = \boldsymbol{\Phi}_\mathbf{q}\dot{\mathbf{q}} = \mathbf{0}

이는 속도 수준의 구속 조건으로, 허용 가능한 속도 \dot{\mathbf{q}}는 자코비안 \boldsymbol{\Phi}_\mathbf{q}의 영공간(null space)에 놓여야 한다.

\dot{\mathbf{q}} \in \text{null}(\boldsymbol{\Phi}_\mathbf{q})

한 번 더 미분하면 가속도 수준의 구속 조건을 얻는다.

\boldsymbol{\Phi}_\mathbf{q}\ddot{\mathbf{q}} + \dot{\boldsymbol{\Phi}}_\mathbf{q}\dot{\mathbf{q}} = \mathbf{0}

\boldsymbol{\Phi}_\mathbf{q}\ddot{\mathbf{q}} = -\dot{\boldsymbol{\Phi}}_\mathbf{q}\dot{\mathbf{q}}

4.3 역기구학의 존재성과 유일성

순기구학 \mathbf{x} = f(\mathbf{q})에서 역기구학 해의 존재를 F(\mathbf{q}, \mathbf{x}) = f(\mathbf{q}) - \mathbf{x} = \mathbf{0}으로 정식화할 수 있다. 음함수 정리에 의해 자코비안 \mathbf{J} = \frac{\partial f}{\partial \mathbf{q}}가 정방이고 가역인 경우 (n = m이고 \det(\mathbf{J}) \neq 0), 현재 해 (\mathbf{q}_0, \mathbf{x}_0) 근방에서 역기구학 해 \mathbf{q} = g(\mathbf{x})가 유일하게 존재한다.

\frac{\partial \mathbf{q}}{\partial \mathbf{x}} = \mathbf{J}^{-1}(\mathbf{q})

이는 자코비안의 역행렬이 역기구학의 국소적 미분임을 나타낸다. 특이점 \det(\mathbf{J}) = 0에서는 음함수 정리의 조건이 위배되므로, 역기구학 해의 유일성이나 존재성이 보장되지 않는다.

4.4 접촉 구속과 마찰

로봇 손(gripper)이 물체와 접촉할 때, 접촉 구속 조건은 음함수 형태로 표현된다.

\Phi(\mathbf{q}_r, \mathbf{q}_o) = 0

여기서 \mathbf{q}_r은 로봇의 형상, \mathbf{q}_o는 물체의 형상이다. 접촉점에서의 법선 방향과 접선 방향의 구속을 음함수 미분으로 해석할 수 있다.

4.5 라그랑주 승수법과의 관계

구속 조건 \boldsymbol{\Phi}(\mathbf{q}) = \mathbf{0} 하에서의 최적화 문제:

\min_{\mathbf{q}} f(\mathbf{q}) \quad \text{subject to} \quad \boldsymbol{\Phi}(\mathbf{q}) = \mathbf{0}

의 필요 조건은 라그랑주 조건:

\nabla f(\mathbf{q}) + \boldsymbol{\lambda}^\top \frac{\partial \boldsymbol{\Phi}}{\partial \mathbf{q}} = \mathbf{0}, \qquad \boldsymbol{\Phi}(\mathbf{q}) = \mathbf{0}

이 연립 방정식 체계도 음함수 정리의 틀 안에서 해석할 수 있으며, 정칙 조건(regularity condition) \text{rank}\left(\frac{\partial \boldsymbol{\Phi}}{\partial \mathbf{q}}\right) = m은 음함수 정리의 가역 조건에 대응한다.

5. 역함수 정리와의 관계

음함수 정리와 밀접하게 관련된 것이 역함수 정리(inverse function theorem)이다.

역함수 정리: \mathbf{f}: \mathbb{R}^n \to \mathbb{R}^n이 C^1 급이고, 점 \mathbf{a}에서 자코비안 행렬 \mathbf{J}_{\mathbf{f}}(\mathbf{a})가 가역이면, \mathbf{a}의 근방에서 \mathbf{f}는 가역이며 역함수 \mathbf{f}^{-1}도 C^1 급이다.

(\mathbf{f}^{-1})'(\mathbf{b}) = [\mathbf{J}_{\mathbf{f}}(\mathbf{a})]^{-1}, \quad \mathbf{b} = \mathbf{f}(\mathbf{a})

음함수 정리는 역함수 정리로부터 유도할 수 있으며, 그 역도 가능하다. 로봇공학에서 역함수 정리는 순기구학의 국소적 역변환(역기구학)의 존재를 보장하는 이론적 근거를 제공한다.

참고 문헌

• Rudin, W. (1976). Principles of Mathematical Analysis. 3rd ed. McGraw-Hill.
• Apostol, T. M. (1974). Mathematical Analysis. 2nd ed. Addison-Wesley.
• Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.
• Murray, R. M., Li, Z., & Sastry, S. S. (1994). A Mathematical Introduction to Robotic Manipulation. CRC Press.

v 0.1