7.42 비동차 상미분 방정식과 특수해

1. 비동차 상미분 방정식의 정의

비동차(nonhomogeneous) 상미분 방정식은 우변에 독립변수의 함수가 포함된 미분 방정식이다. $n$ 차 선형 비동차 상미분 방정식의 일반적 형태는 다음과 같다.

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

여기서 $g(t) \neq 0$ 인 함수를 비동차 항(nonhomogeneous term) 또는 강제 함수(forcing function)라 한다. 로봇 동역학에서 이 강제 함수는 외부 토크, 외란(disturbance), 또는 제어 입력에 해당한다.

상수 계수를 갖는 2차 선형 비동차 상미분 방정식은 로봇공학에서 가장 빈번하게 등장하며, 다음과 같은 형태를 취한다.

$\ddot{y} + a\dot{y} + by = g(t)$

이 방정식에 대응하는 동차 방정식(homogeneous equation)은 $g(t) = 0$ 으로 놓은 것이다.

$\ddot{y} + a\dot{y} + by = 0$

2. 일반해의 구조

비동차 선형 상미분 방정식의 일반해(general solution) $y(t)$ 는 다음 두 성분의 합으로 구성된다.

$y(t) = y_h(t) + y_p(t)$

여기서 $y_h(t)$ 는 대응하는 동차 방정식의 일반해(여해, complementary solution)이고, $y_p(t)$ 는 비동차 방정식을 만족하는 하나의 특수해(particular solution)이다.

이 구조가 성립하는 이유는 선형성에 기인한다. $y_1(t)$ 과 $y_2(t)$ 가 모두 비동차 방정식의 해이면, 그 차 $y_1(t) - y_2(t)$ 는 동차 방정식의 해가 된다. 따라서 비동차 방정식의 모든 해는 하나의 특수해에 동차 방정식의 일반해를 더한 형태로 표현할 수 있다.

$n$ 차 동차 방정식의 일반해가 $n$ 개의 선형 독립인 해 $\{y_1, y_2, \ldots, y_n\}$ 의 선형 결합으로 표현되므로, 비동차 방정식의 일반해는 다음과 같다.

$y(t) = c_1 y_1(t) + c_2 y_2(t) + \cdots + c_n y_n(t) + y_p(t)$

여기서 $c_1, c_2, \ldots, c_n$ 은 초기 조건에 의해 결정되는 임의 상수이다.

3. 미정 계수법

3.1 기본 원리

미정 계수법(method of undetermined coefficients)은 상수 계수 선형 비동차 상미분 방정식에서 강제 함수 $g(t)$ 가 다항식, 지수함수, 삼각함수 또는 이들의 곱으로 구성될 때 적용 가능한 방법이다. 핵심 아이디어는 $g(t)$ 의 형태로부터 특수해 $y_p(t)$ 의 형태를 추정하고, 미지의 계수를 방정식에 대입하여 결정하는 것이다.

3.2 강제 함수에 따른 특수해의 추정 형태

다음 표는 상수 계수 선형 비동차 방정식에서 강제 함수의 유형에 따른 특수해의 추정 형태를 정리한 것이다.

강제 함수 $g(t)$	추정 특수해 $y_p(t)$
$P_n(t)$ ( $n$ 차 다항식)	$A_n t^n + A_{n-1} t^{n-1} + \cdots + A_1 t + A_0$
$e^{\alpha t}$	$A e^{\alpha t}$
$\sin(\beta t)$ 또는 $\cos(\beta t)$	$A \cos(\beta t) + B \sin(\beta t)$
$e^{\alpha t} P_n(t)$	$e^{\alpha t}(A_n t^n + \cdots + A_0)$
$e^{\alpha t} \sin(\beta t)$ 또는 $e^{\alpha t} \cos(\beta t)$	$e^{\alpha t}(A \cos(\beta t) + B \sin(\beta t))$

3.3 수정 규칙

추정한 특수해의 항이 동차 방정식의 해와 중복되는 경우, 특수해 추정 형태에 $t^s$ 를 곱하여 수정하여야 한다. 여기서 $s$ 는 해당 항이 동차 해의 근과 일치하는 중복도(multiplicity)이다.

예를 들어, 다음 방정식을 고려하라.

$\ddot{y} - 4\dot{y} + 4y = e^{2t}$

특성 방정식 $r^2 - 4r + 4 = 0$ 의 근은 $r = 2$ (이중근)이다. 동차 해가 $y_h = (c_1 + c_2 t)e^{2t}$ 이므로, $e^{2t}$ 와 $te^{2t}$ 는 이미 동차 해에 포함된다. 따라서 특수해는 다음과 같이 추정한다.

$y_p(t) = At^2 e^{2t}$

이를 원래 방정식에 대입하면 $A = \frac{1}{2}$ 을 얻는다.

3.4 중첩 원리의 활용

강제 함수가 여러 항의 합으로 구성된 경우, 중첩 원리(superposition principle)를 적용할 수 있다. $g(t) = g_1(t) + g_2(t) + \cdots + g_m(t)$ 이면, 각 $g_i(t)$ 에 대한 특수해 $y_{p_i}(t)$ 를 구한 후 이들을 합산하면 전체 특수해가 된다.

$y_p(t) = y_{p_1}(t) + y_{p_2}(t) + \cdots + y_{p_m}(t)$

4. 매개변수 변환법

4.1 기본 원리

매개변수 변환법(variation of parameters)은 미정 계수법보다 일반적인 방법으로, 강제 함수 $g(t)$ 의 형태에 제한이 없다. 이 방법은 동차 방정식의 해를 이미 알고 있다는 전제 하에, 동차 해의 상수를 함수로 대체하여 특수해를 구한다.

2차 방정식 $\ddot{y} + p(t)\dot{y} + q(t)y = g(t)$ 에 대하여, 동차 해가 $y_1(t)$ 과 $y_2(t)$ 일 때 특수해를 다음과 같이 가정한다.

$y_p(t) = u_1(t)y_1(t) + u_2(t)y_2(t)$

4.2 론스키안과 해의 공식

$u_1(t)$ 과 $u_2(t)$ 를 결정하기 위하여 다음 연립 방정식을 풀어야 한다.

$u_1'y_1 + u_2'y_2 = 0$

$u_1'y_1' + u_2'y_2' = g(t)$

이 연립 방정식의 해는 론스키안(Wronskian) $W(t) = y_1 y_2' - y_2 y_1'$ 을 이용하여 다음과 같이 구해진다.

$u_1'(t) = -\frac{y_2(t)g(t)}{W(t)}, \quad u_2'(t) = \frac{y_1(t)g(t)}{W(t)}$

따라서 특수해는 다음과 같다.

$y_p(t) = -y_1(t)\int \frac{y_2(t)g(t)}{W(t)}dt + y_2(t)\int \frac{y_1(t)g(t)}{W(t)}dt$

4.3 $n$ 차 방정식으로의 확장

$n$ 차 선형 비동차 방정식에 대한 매개변수 변환법의 일반화된 형태에서, $n$ 개의 동차 해 $\{y_1, y_2, \ldots, y_n\}$ 이 알려진 경우, 특수해는 다음과 같다.

$y_p(t) = \sum_{k=1}^{n} y_k(t) \int \frac{W_k(t)}{W(t)} g(t) \, dt$

여기서 $W(t)$ 는 $\{y_1, y_2, \ldots, y_n\}$ 의 론스키안이며, $W_k(t)$ 는 $W(t)$ 의 $k$ 번째 열을 $(0, 0, \ldots, 0, 1)^T$ 로 대체한 행렬식이다.

5. 그린 함수를 이용한 특수해 표현

5.1 그린 함수의 정의

그린 함수(Green’s function)는 선형 미분 연산자의 역연산자를 나타내는 핵(kernel) 함수이다. 선형 미분 연산자 $L$ 에 대하여 $Ly = g(t)$ 의 해를 다음과 같이 적분으로 표현할 수 있다.

$y_p(t) = \int_{t_0}^{t} G(t, \tau) g(\tau) \, d\tau$

여기서 $G(t, \tau)$ 가 그린 함수이다. 상수 계수 2차 방정식의 경우, 그린 함수는 충격 응답(impulse response)과 동일하며, 동차 해를 이용하여 다음과 같이 구할 수 있다.

$G(t, \tau) = \frac{y_1(\tau)y_2(t) - y_1(t)y_2(\tau)}{W(\tau)}$

5.2 로봇 제어 시스템에서의 해석

로봇 제어 시스템에서 그린 함수는 물리적으로 단위 충격 입력(Dirac delta function $\delta(t - \tau)$ )에 대한 시스템 응답을 의미한다. 선형 시불변(LTI) 시스템의 경우, 그린 함수는 $G(t, \tau) = h(t - \tau)$ 의 형태를 취하며, 여기서 $h(t)$ 는 임펄스 응답 함수이다. 이를 통해 임의의 입력 $g(t)$ 에 대한 출력은 합성곱(convolution)으로 표현된다.

$y_p(t) = \int_0^t h(t - \tau) g(\tau) \, d\tau = (h * g)(t)$

이 관계는 라플라스 변환 영역에서 전달 함수와 직접적으로 연결되며, 로봇 관절 구동기의 시간 영역 응답 해석에 핵심적으로 활용된다.

6. 로봇 동역학에서의 비동차 상미분 방정식

6.1 단일 관절의 구동 방정식

로봇 매니퓰레이터의 단일 관절에 대한 운동 방정식은 전형적인 2차 비동차 상미분 방정식이다. 관절 각도 $\theta(t)$ 에 대하여 다음과 같이 표현된다.

$J\ddot{\theta} + b\dot{\theta} + k\theta = \tau(t)$

여기서 $J$ 는 관성 모멘트, $b$ 는 점성 마찰 계수, $k$ 는 강성 계수, $\tau(t)$ 는 인가된 토크이다. 이 방정식에서 $\tau(t)$ 가 강제 함수이며, 이 강제 함수의 형태에 따라 특수해의 물리적 의미가 달라진다.

6.2 계단 입력에 대한 응답

제어 토크가 계단 함수(step function) $\tau(t) = \tau_0 \cdot u(t)$ 로 주어지는 경우, 특수해는 정상 상태(steady-state) 응답에 해당한다. $k \neq 0$ 일 때 미정 계수법을 적용하면 특수해는 다음과 같다.

$y_p = \frac{\tau_0}{k}$

이는 관절이 최종적으로 도달하는 평형 위치를 나타낸다. 동차 해는 과도 응답(transient response)으로, 시간이 지남에 따라 감쇠하여 소멸한다.

6.3 정현파 입력에 대한 응답

강제 함수가 정현파 $\tau(t) = F_0 \sin(\omega t)$ 인 경우, 특수해는 강제 진동(forced vibration)의 정상 상태 응답을 나타낸다. 미정 계수법을 적용하면 특수해는 다음과 같다.

$y_p(t) = A\cos(\omega t) + B\sin(\omega t)$

미지 계수를 원래 방정식에 대입하여 정리하면 다음을 얻는다.

$A = \frac{-b\omega F_0}{(k - J\omega^2)^2 + (b\omega)^2}, \quad B = \frac{(k - J\omega^2)F_0}{(k - J\omega^2)^2 + (b\omega)^2}$

진폭과 위상의 관점에서 특수해를 다시 쓰면 다음과 같다.

$y_p(t) = M \sin(\omega t - \phi)$

여기서 진폭 $M$ 과 위상 지연 $\phi$ 는 각각 다음과 같다.

$M = \frac{F_0}{\sqrt{(k - J\omega^2)^2 + (b\omega)^2}}, \quad \phi = \arctan\left(\frac{b\omega}{k - J\omega^2}\right)$

강제 주파수 $\omega$ 가 고유 진동수 $\omega_n = \sqrt{k/J}$ 에 접근하면 진폭이 급격히 증가하는 공진(resonance) 현상이 발생한다. 이는 로봇 구조물 설계 및 제어기 설계 시 반드시 고려하여야 하는 현상이다.

6.4 비감쇠 공진의 특수해

감쇠가 없는 경우( $b = 0$ )에 강제 주파수가 고유 진동수와 정확히 일치하면( $\omega = \omega_n$ ), 통상적인 추정 형태의 특수해가 동차 해와 중복된다. 수정 규칙을 적용하여 $t$ 를 곱하면 특수해는 다음과 같다.

$y_p(t) = \frac{F_0}{2J\omega_n} t \sin(\omega_n t)$

이 해는 시간에 비례하여 진폭이 무한히 증가하는 순수 공진 현상을 나타내며, 실제 로봇 시스템에서는 비선형 효과와 구조적 한계에 의해 제한된다.

7. 특수해의 존재성과 유일성

비동차 선형 상미분 방정식의 해의 존재성과 유일성은 다음 정리에 의해 보장된다.

정리 (피카르-린델뢰프, Picard-Lindelof). 구간 $I$ 에서 $a_i(t)$ ( $i = 0, 1, \ldots, n$ )와 $g(t)$ 가 연속이고 $a_n(t) \neq 0$ 이면, 초기 조건 $y(t_0) = y_0, y'(t_0) = y_0', \ldots, y^{(n-1)}(t_0) = y_0^{(n-1)}$ 을 만족하는 해가 구간 $I$ 에서 유일하게 존재한다.

이 정리는 로봇 동역학 시뮬레이션에서 수치 해법의 이론적 근거를 제공한다. 관성 행렬이 양정치(positive definite)이면 최고차 계수가 영이 아님이 보장되므로, 로봇 동역학 방정식의 해는 언제나 존재하고 유일하다.

8. 요약

비동차 상미분 방정식의 일반해는 동차 해(과도 응답)와 특수해(정상 상태 응답)의 합으로 구성된다. 미정 계수법은 상수 계수 방정식에서 강제 함수가 특정 형태일 때 효율적인 방법이며, 매개변수 변환법은 보다 일반적인 강제 함수에 적용 가능하다. 그린 함수 접근법은 임펄스 응답을 통한 체계적인 해 표현을 제공한다. 로봇 동역학에서 비동차 항은 외부 토크와 제어 입력에 대응하며, 특수해의 분석은 관절 응답 특성, 공진 회피, 정상 상태 정밀도 등 제어 시스템 설계의 핵심 정보를 제공한다.

참고 문헌

Boyce, W. E. & DiPrima, R. C. (2012). Elementary Differential Equations and Boundary Value Problems (10th ed.). Wiley.
Zill, D. G. (2017). Advanced Engineering Mathematics (6th ed.). Jones & Bartlett Learning.
Craig, J. J. (2005). Introduction to Robotics: Mechanics and Control (3rd ed.). Pearson.
Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.
Ogata, K. (2010). Modern Control Engineering (5th ed.). Prentice Hall.

version: 0.1.0