7.39 상미분 방정식의 분류와 차수

1. 상미분 방정식의 정의

상미분 방정식(ordinary differential equation, ODE)은 하나의 독립 변수(일반적으로 시간 $t$ )에 대한 미지 함수와 그 도함수들 사이의 관계를 나타내는 방정식이다. 일반적인 형태는

$F\bigl(t, y, y', y'', \dots, y^{(n)}\bigr) = 0$

이며, 여기서 $y = y(t)$ 는 미지 함수, $y' = dy/dt$ , $y'' = d^2y/dt^2$ , $y^{(n)} = d^n y/dt^n$ 은 각각 1차, 2차, $n$ 차 도함수이다.

상미분 방정식은 독립 변수가 하나인 점에서 편미분 방정식(partial differential equation, PDE)과 구별된다. 편미분 방정식은 둘 이상의 독립 변수에 대한 편도함수를 포함한다.

2. 차수와 계수

2.1 차수

상미분 방정식의 차수(order)는 방정식에 나타나는 최고차 도함수의 계수(階數)이다.

1차 ODE: $y' = f(t, y)$
2차 ODE: $y'' = f(t, y, y')$
$n$ 차 ODE: $y^{(n)} = f(t, y, y', \dots, y^{(n-1)})$

로봇공학에서 뉴턴-오일러 운동 방정식이나 라그랑주 동역학 방정식은 일반화 좌표의 2차 도함수(가속도)를 포함하므로 2차 상미분 방정식에 해당한다.

2.2 계수

상미분 방정식의 계수(degree)는 최고차 도함수의 거듭제곱 차수이다. 예를 들어

$\left(\frac{d^2 y}{dt^2}\right)^3 + \frac{dy}{dt} = 0$

은 2차, 3계의 상미분 방정식이다. 계수가 정의되려면 방정식이 도함수에 대한 다항식 형태로 표현 가능해야 한다.

3. 선형성에 의한 분류

3.1 선형 상미분 방정식

$n$ 차 선형 상미분 방정식(linear ODE)은 미지 함수 $y$ 와 그 도함수가 모두 1차의 형태로만 나타나며, 서로의 곱이 존재하지 않는 방정식이다.

$a_n(t) y^{(n)} + a_{n-1}(t) y^{(n-1)} + \cdots + a_1(t) y' + a_0(t) y = g(t)$

여기서 계수 $a_0(t), a_1(t), \dots, a_n(t)$ 와 비동차 항(inhomogeneous term) $g(t)$ 는 $t$ 의 함수이다. 선형 ODE의 핵심 성질은 중첩 원리(superposition principle)가 성립한다는 것이다. 즉, $y_1$ 과 $y_2$ 가 동차 방정식( $g = 0$ )의 해이면, 임의의 상수 $c_1, c_2$ 에 대해 $c_1 y_1 + c_2 y_2$ 도 해이다.

3.2 비선형 상미분 방정식

미지 함수나 그 도함수의 비선형 항(예: $y^2$ , $y y'$ , $\sin y$ )을 포함하는 방정식을 비선형 상미분 방정식(nonlinear ODE)이라 한다.

$y'' + \sin y = 0 \quad \text{(단진자 방정식)}$

$y'' + y y' + y^3 = 0$

비선형 ODE에서는 중첩 원리가 성립하지 않으며, 해의 존재성, 유일성, 안정성 분석이 선형의 경우보다 본질적으로 복잡하다. 로봇 동역학 방정식은 원심력, 코리올리력, 중력 항의 비선형성으로 인해 비선형 ODE에 해당한다.

4. 동차성에 의한 분류

4.1 동차 방정식

선형 ODE에서 비동차 항 $g(t) = 0$ 이면 동차 방정식(homogeneous equation)이라 한다.

$a_n(t) y^{(n)} + a_{n-1}(t) y^{(n-1)} + \cdots + a_0(t) y = 0$

동차 방정식의 해 공간은 벡터 공간을 형성하며, $n$ 차 동차 선형 ODE의 해 공간의 차원은 $n$ 이다.

4.2 비동차 방정식

$g(t) \neq 0$ 인 경우를 비동차 방정식(inhomogeneous equation 또는 nonhomogeneous equation)이라 한다.

$a_n(t) y^{(n)} + \cdots + a_0(t) y = g(t)$

비동차 방정식의 일반해(general solution)는 대응하는 동차 방정식의 일반해 $y_h$ 와 비동차 방정식의 하나의 특수해(particular solution) $y_p$ 의 합이다.

$y = y_h + y_p$

5. 계수의 성질에 의한 분류

5.1 상수 계수 방정식

계수 $a_0, a_1, \dots, a_n$ 이 모두 상수인 선형 ODE를 상수 계수 선형 ODE(constant coefficient linear ODE)라 한다.

$a_n y^{(n)} + a_{n-1} y^{(n-1)} + \cdots + a_1 y' + a_0 y = g(t)$

상수 계수 방정식은 특성 방정식(characteristic equation)을 이용하여 체계적으로 풀 수 있다. 로봇공학에서 선형화된 동역학 모형이나 제어 시스템의 전달 함수 모형은 상수 계수 선형 ODE로 표현된다.

5.2 변수 계수 방정식

계수가 독립 변수 $t$ 의 함수인 경우를 변수 계수 선형 ODE(variable coefficient linear ODE)라 한다. 오일러-코시 방정식(Euler-Cauchy equation)

$t^2 y'' + t y' + y = 0$

이 대표적인 예이다. 변수 계수 방정식은 일반적으로 닫힌 형태(closed-form)의 해가 존재하지 않으며, 멱급수 해법(power series method)이나 수치 해법을 필요로 한다.

6. 양함수 형태와 음함수 형태

6.1 양함수 형태

최고차 도함수에 대해 풀린 형태를 양함수 형태(explicit form) 또는 정규 형태(normal form)라 한다.

$y^{(n)} = f(t, y, y', \dots, y^{(n-1)})$

양함수 형태에서는 피카르-린들뢰프 정리(Picard-Lindelöf theorem)에 의한 해의 존재와 유일성 조건을 직접 확인할 수 있다.

6.2 음함수 형태

최고차 도함수에 대해 풀릴 수 없는 형태를 음함수 형태(implicit form)라 한다.

$F(t, y, y', \dots, y^{(n)}) = 0$

음함수 형태의 방정식은 해의 다가성(multivaluedness)이나 특이점(singular point)을 수반할 수 있다.

7. 초기값 문제와 경계값 문제

7.1 초기값 문제

$n$ 차 ODE에 $n$ 개의 초기 조건이 동일 시점 $t_0$ 에서 주어지는 문제를 초기값 문제(initial value problem, IVP)라 한다.

$y^{(n)} = f(t, y, y', \dots, y^{(n-1)}), \quad y(t_0) = y_0, \; y'(t_0) = y_0', \; \dots, \; y^{(n-1)}(t_0) = y_0^{(n-1)}$

로봇 동역학 시뮬레이션은 초기 관절 위치와 속도가 주어진 상태에서 시간에 따른 운동을 추적하는 초기값 문제이다.

7.2 경계값 문제

조건이 서로 다른 두 시점(또는 그 이상)에서 주어지는 문제를 경계값 문제(boundary value problem, BVP)라 한다.

$y'' = f(t, y, y'), \quad y(t_0) = \alpha, \; y(t_f) = \beta$

로봇 궤적 계획에서 초기 위치와 최종 위치를 동시에 만족하는 궤적을 구하는 문제는 경계값 문제에 해당한다.

8. 연립 상미분 방정식

$n$ 차 스칼라 ODE는 $n$ 개의 1차 ODE로 이루어진 연립 방정식으로 변환할 수 있다. 상태 변수 $x_1 = y$ , $x_2 = y'$ , $\dots$ , $x_n = y^{(n-1)}$ 을 도입하면

$\dot{x}_1 = x_2, \quad \dot{x}_2 = x_3, \quad \dots, \quad \dot{x}_n = f(t, x_1, x_2, \dots, x_n)$

이 되어, 벡터 표기로

$\dot{\mathbf{x}} = \mathbf{f}(t, \mathbf{x})$

의 1차 연립 ODE 형태를 얻는다. 이 변환은 수치 해법과 상태 공간 제어 이론의 기초가 된다.

9. 로봇공학에서 나타나는 ODE의 분류

로봇공학에서 주요하게 다루는 상미분 방정식의 유형을 정리하면 다음과 같다.

로봇공학 모형	ODE 분류	차수	선형성
관절 동역학	연립 비선형 ODE	2차	비선형
선형화 동역학	연립 선형 ODE (상수 계수)	2차	선형
제어 시스템 상태 방정식	연립 1차 ODE	1차	선형/비선형
DC 모터 전기 동역학	1차 선형 ODE	1차	선형
감쇠 진동 모형	2차 선형 ODE (상수 계수)	2차	선형
단진자	2차 비선형 ODE	2차	비선형

로봇 매니퓰레이터의 운동 방정식

$\mathbf{M}(\mathbf{q}) \ddot{\mathbf{q}} + \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \dot{\mathbf{q}} + \mathbf{g}(\mathbf{q}) = \boldsymbol{\tau}$

은 $n$ 개의 2차 비선형 연립 ODE이다. 관성 행렬 $\mathbf{M}(\mathbf{q})$ 가 위치에 의존하고, 코리올리/원심력 항 $\mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \dot{\mathbf{q}}$ 가 속도의 이차 형식이며, 중력 항 $\mathbf{g}(\mathbf{q})$ 가 삼각 함수를 포함하므로 비선형성이 발생한다.

10. 해의 존재와 유일성

10.1 피카르-린들뢰프 정리

1차 ODE의 초기값 문제

$y' = f(t, y), \quad y(t_0) = y_0$

에서 $f$ 가 직사각형 영역 $R = \{(t, y) : \lvert t - t_0 \rvert \leq a, \lvert y - y_0 \rvert \leq b\}$ 에서 연속이고, $y$ 에 대해 립시츠 조건(Lipschitz condition)

$\lvert f(t, y_1) - f(t, y_2) \rvert \leq L \lvert y_1 - y_2 \rvert, \quad \forall\, (t, y_1), (t, y_2) \in R$

을 만족하면, $t_0$ 를 포함하는 구간에서 초기값 문제의 해가 존재하며 유일하다(Picard, 1890; Lindelöf, 1894).

립시츠 조건은 $f$ 가 $y$ 에 대해 편미분 가능하고 $\lvert \partial f / \partial y \rvert$ 가 유계이면 자동으로 만족된다. 로봇 동역학 방정식에서 관성 행렬 $\mathbf{M}(\mathbf{q})$ 가 양정치이면 $\ddot{\mathbf{q}}$ 에 대해 유일하게 풀리므로, 양함수 형태의 연립 ODE로 변환 가능하며 해의 존재와 유일성이 보장된다.

11. 참고 문헌

Boyce, W. E., & DiPrima, R. C. (2012). Elementary Differential Equations and Boundary Value Problems. 10th ed. Wiley.
Coddington, E. A., & Levinson, N. (1955). Theory of Ordinary Differential Equations. McGraw-Hill.
Khalil, H. K. (2002). Nonlinear Systems. 3rd ed. Prentice Hall.
Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.
Ascher, U. M., & Petzold, L. R. (1998). Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations. SIAM.

v 0.2