7.41 2차 선형 상미분 방정식과 특성 방정식

1. 차 선형 상미분 방정식의 일반 형태

2차 선형 상미분 방정식(second-order linear ODE)의 일반적인 형태는

$a(t) y'' + b(t) y' + c(t) y = g(t)$

이며, $a(t) \neq 0$ 인 영역에서 양변을 $a(t)$ 로 나누어 표준 형태(standard form)

$y'' + p(t) y' + q(t) y = r(t)$

로 정리할 수 있다. $r(t) = 0$ 이면 동차 방정식(homogeneous equation), $r(t) \neq 0$ 이면 비동차 방정식(nonhomogeneous equation)이다.

2. 동차 방정식의 해 구조

2.1 중첩 원리와 해 공간

2차 동차 선형 ODE

$y'' + p(t) y' + q(t) y = 0$

의 해 집합은 2차원 벡터 공간을 형성한다. $y_1(t)$ 과 $y_2(t)$ 가 선형 독립인 두 해이면, 일반해(general solution)는

$y(t) = c_1 y_1(t) + c_2 y_2(t)$

이며, 상수 $c_1, c_2$ 는 초기 조건 $y(t_0) = y_0$ , $y'(t_0) = y_0'$ 에 의해 유일하게 결정된다.

2.2 론스키안과 선형 독립성

두 해 $y_1, y_2$ 의 선형 독립성은 론스키안(Wronskian)

$W(y_1, y_2)(t) = \begin{vmatrix} y_1(t) & y_2(t) \\ y_1'(t) & y_2'(t) \end{vmatrix} = y_1(t) y_2'(t) - y_2(t) y_1'(t)$

으로 판정된다. $p(t)$ 와 $q(t)$ 가 연속인 구간에서, $W \neq 0$ 이면 $y_1$ 과 $y_2$ 는 선형 독립이다. 아벨 정리(Abel’s theorem)에 의해 론스키안은

$W(t) = W(t_0) \exp\left(-\int_{t_0}^{t} p(\tau) \, d\tau\right)$

를 만족하므로, $W(t_0) \neq 0$ 이면 론스키안은 해당 구간 전체에서 0이 아니다.

3. 상수 계수 동차 방정식

3.1 특성 방정식의 유도

상수 계수 2차 동차 선형 ODE

$ay'' + by' + cy = 0, \quad a, b, c \in \mathbb{R}, \; a \neq 0$

에서 $y = e^{\lambda t}$ 형태의 해를 시행하면

$a\lambda^2 e^{\lambda t} + b\lambda e^{\lambda t} + c e^{\lambda t} = 0$

$e^{\lambda t} \neq 0$ 이므로

$a\lambda^2 + b\lambda + c = 0$

을 얻는다. 이 2차 대수 방정식을 특성 방정식(characteristic equation)이라 하며, 그 근을 특성근(characteristic root)이라 한다.

3.2 판별식에 의한 분류

특성 방정식의 판별식(discriminant) $\Delta = b^2 - 4ac$ 의 부호에 따라 해의 형태가 세 가지로 분류된다.

3.2.1 경우 1: 서로 다른 두 실수근 ( $\Delta > 0$ )

두 특성근 $\lambda_1 \neq \lambda_2$ 가 모두 실수이면, 일반해는

$y(t) = c_1 e^{\lambda_1 t} + c_2 e^{\lambda_2 t}$

이다. 두 해 $e^{\lambda_1 t}$ 와 $e^{\lambda_2 t}$ 의 론스키안은 $W = (\lambda_2 - \lambda_1) e^{(\lambda_1 + \lambda_2)t} \neq 0$ 이므로 선형 독립이 보장된다.

3.2.2 경우 2: 중복 실수근 ( $\Delta = 0$ )

$\lambda_1 = \lambda_2 = \lambda = -b/(2a)$ 인 경우, $e^{\lambda t}$ 만으로는 두 개의 선형 독립 해를 구성할 수 없다. 계수 변환법(reduction of order)에 의해 두 번째 선형 독립 해가 $te^{\lambda t}$ 임을 보일 수 있다. 일반해는

$y(t) = (c_1 + c_2 t) e^{\lambda t}$

이다.

3.2.3 경우 3: 복소 켤레근 ( $\Delta < 0$ )

두 특성근이 복소 켤레쌍 $\lambda = \alpha \pm i\beta$ 이면, 여기서 $\alpha = -b/(2a)$ , $\beta = \sqrt{4ac - b^2}/(2a)$ 이다. 오일러 공식(Euler’s formula) $e^{i\beta t} = \cos(\beta t) + i \sin(\beta t)$ 를 적용하면, 실수값 일반해는

$y(t) = e^{\alpha t} \bigl(c_1 \cos(\beta t) + c_2 \sin(\beta t)\bigr)$

이다. 이를 진폭-위상(amplitude-phase) 형태로 다시 쓰면

$y(t) = A e^{\alpha t} \cos(\beta t - \phi)$

이며, $A = \sqrt{c_1^2 + c_2^2}$ , $\phi = \arctan(c_2/c_1)$ 이다.

4. 해의 물리적 해석

4.1 특성근과 시스템 응답

특성근의 실수부 $\alpha$ 와 허수부 $\beta$ 는 시스템 응답의 성격을 결정한다.

특성근의 형태	$\alpha$	$\beta$	응답 특성
서로 다른 음의 실수근	$\alpha_1, \alpha_2 < 0$	$0$	비진동 감쇠 (overdamped)
중복 음의 실수근	$\alpha < 0$	$0$	임계 감쇠 (critically damped)
음의 실수부를 가진 복소근	$\alpha < 0$	$\beta \neq 0$	감쇠 진동 (underdamped)
순허수근	$0$	$\beta \neq 0$	비감쇠 진동 (undamped)
양의 실수부를 가진 근	$\alpha > 0$	임의	불안정 (발산)

4.2 감쇠비와 고유 진동수

표준 2차 시스템

$\ddot{y} + 2\zeta \omega_n \dot{y} + \omega_n^2 y = 0$

에서 $\omega_n > 0$ 은 비감쇠 고유 진동수(undamped natural frequency), $\zeta \geq 0$ 은 감쇠비(damping ratio)이다. 특성 방정식은

$\lambda^2 + 2\zeta \omega_n \lambda + \omega_n^2 = 0$

이며, 특성근은

$\lambda = -\zeta \omega_n \pm \omega_n \sqrt{\zeta^2 - 1}$

이다.

$\zeta > 1$ : 과감쇠 (overdamped), 두 음의 실수근
$\zeta = 1$ : 임계 감쇠 (critically damped), 중복 음의 실수근 $\lambda = -\omega_n$
$0 < \zeta < 1$ : 부족 감쇠 (underdamped), 복소근 $\lambda = -\zeta \omega_n \pm i \omega_d$
$\zeta = 0$ : 비감쇠, 순허수근 $\lambda = \pm i \omega_n$

부족 감쇠의 경우, 감쇠 진동수(damped natural frequency)는

$\omega_d = \omega_n \sqrt{1 - \zeta^2}$

이다.

5. 고차 상수 계수 방정식으로의 확장

$n$ 차 상수 계수 동차 선형 ODE

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

의 특성 방정식은

$a_n \lambda^n + a_{n-1} \lambda^{n-1} + \cdots + a_1 \lambda + a_0 = 0$

이다. $n$ 차 다항식의 근의 개수(중복도 포함)가 $n$ 개이므로, 다음의 규칙에 따라 일반해를 구성한다.

단순 실수근 $\lambda_k$ : $e^{\lambda_k t}$
중복도 $m$ 의 실수근 $\lambda_k$ : $e^{\lambda_k t}, te^{\lambda_k t}, \dots, t^{m-1}e^{\lambda_k t}$
복소 켤레근 $\alpha \pm i\beta$ (단순): $e^{\alpha t}\cos(\beta t)$ , $e^{\alpha t}\sin(\beta t)$
복소 켤레근 $\alpha \pm i\beta$ (중복도 $m$ ): $t^k e^{\alpha t}\cos(\beta t)$ , $t^k e^{\alpha t}\sin(\beta t)$ , $k = 0, 1, \dots, m-1$

6. 로봇공학에서의 응용

6.1 질량-스프링-댐퍼 모형

로봇 관절의 기계적 모형은 회전 관성 $J$ , 점성 감쇠 계수 $B$ , 강성 계수 $K$ 에 의한 2차 선형 ODE

$J\ddot{\theta} + B\dot{\theta} + K\theta = 0$

로 표현된다. 특성 방정식은

$J\lambda^2 + B\lambda + K = 0$

이며, 비감쇠 고유 진동수 $\omega_n = \sqrt{K/J}$ , 감쇠비 $\zeta = B/(2\sqrt{JK})$ 이다. 감쇠비의 값에 따라 관절의 자유 응답이 진동적인지 비진동적인지가 결정된다.

6.2 PD 제어 시스템의 폐루프 동역학

단일 관절에 PD 제어기 $\tau = -K_p (\theta - \theta_d) - K_d \dot{\theta}$ 를 적용하면, 폐루프 동역학은

$J\ddot{\theta} + (B + K_d)\dot{\theta} + K_p \theta = K_p \theta_d$

가 된다. 동차 부분의 특성 방정식

$J\lambda^2 + (B + K_d)\lambda + K_p = 0$

의 근에 의해 폐루프 시스템의 과도 응답 특성이 결정된다. 제어 이득 $K_p$ 와 $K_d$ 의 선택은 원하는 고유 진동수와 감쇠비를 달성하도록 설계된다.

$\omega_n = \sqrt{\frac{K_p}{J}}, \quad \zeta = \frac{B + K_d}{2\sqrt{J K_p}}$

6.3 안정성 판정

특성 방정식의 모든 근의 실수부가 음수이면, 시스템은 점근적으로 안정(asymptotically stable)하다. 2차 시스템 $a\lambda^2 + b\lambda + c = 0$ 에서 모든 계수가 양수, 즉 $a > 0$ , $b > 0$ , $c > 0$ 이면 두 근의 실수부가 모두 음수임이 보장된다. 이 조건은 로봇 제어 시스템의 안정성을 확보하기 위한 이득 선택의 기초가 된다.

7. 참고 문헌

Boyce, W. E., & DiPrima, R. C. (2012). Elementary Differential Equations and Boundary Value Problems. 10th ed. Wiley.
Ogata, K. (2010). Modern Control Engineering. 5th ed. Prentice Hall.
Nise, N. S. (2015). Control Systems Engineering. 7th ed. Wiley.
Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.
Strang, G. (2014). Differential Equations and Linear Algebra. Wellesley-Cambridge Press.

v 0.2