부분분수 분해법을 통한 역변환

이론적 개요

라플라스 변환에서 유리함수 형태의 변환식을 역변환하기 위해 가장 널리 활용되는 기법 중 하나가 부분분수 분해이다. 실함수 영역에서의 부분분수 분해와 유사하나, 라플라스 변환에서는 특정 형태의 기저 함수(지수 함수, 사인/코사인 등)를 기반으로 식을 재구성하므로, 보다 정밀하게 분해 과정을 다루어야 한다. 부분분수 분해를 통해 얻어진 각 항은 표준 라플라스 변환 표에 대응하는 간단한 형태가 되어, 역변환 과정이 명료해진다.

분해를 적용하기 위해서는 먼저 라플라스 변수인 $s$에 대한 분모 다항식을 완전히 인수분해해야 한다. 인수분해가 가능하다는 것은 실수 영역에서 일차 혹은 이차 인수로 나타낼 수 있음을 의미한다. 일차 인수로 완전히 인수분해되는 경우, 서로 다른 단순근(distinct real roots)이거나 중근(repeated roots)일 수 있다. 만약 허근(complex roots)이 존재하면 공액 쌍을 이루는 이차 다항식 인수로 나타낼 수 있다. 이 각각의 경우에 대해 부분분수 분해 형태는 조금씩 달라진다.

단순근을 갖는 분자·분모 다항식의 경우

단순근을 갖는 분모의 경우를 먼저 살펴보자. 예를 들어, 다음과 같은 $s$ 영역의 함수 $F(s)$가 있다고 가정한다.

$$F(s) = \frac{2s + 3}{s(s+1)}.$$

분모 $s(s+1)$는 서로 다른 일차 인수 $s$와 $s+1$로 이미 인수분해된 상태이다. 이때 분자 차수가 분모 차수보다 낮으므로(분자 차수가 1, 분모 차수가 2) 별도의 다항식 나눗셈 없이 바로 분해를 시도할 수 있다. 부분분수 분해의 표준 형태에 따라,

$$\frac{2s + 3}{s(s+1)} = \frac{A}{s} + \frac{B}{s+1}.$$

이 식을 양변에서 통분하여 얻는 항등식을 활용해 $A$와 $B$를 구한다.

$$2s + 3 = A(s+1) + B s.$$

양변이 $s$에 대한 다항식으로 동일하므로, 계수를 비교하면,

$$2s + 3 = (A + B)s + A.$$

이로부터

$$A + B = 2, \quad A = 3.$$

을 얻는다. 따라서

$$B = 2 - A = 2 - 3 = -1.$$

결과적으로

$$F(s) = \frac{3}{s} - \frac{1}{s+1}.$$

이제 표준 라플라스 변환 식을 이용하면,

$$\mathcal{L}^{-1}\{ \frac{1}{s} \} = 1, \quad \mathcal{L}^{-1}\{ \frac{1}{s+1} \} = e^{-t}.$$

를 알고 있으므로,

$$\mathcal{L}^{-1}\{F(s)\} = \mathcal{L}^{-1}\Big\{\frac{3}{s} - \frac{1}{s+1}\Big\} = 3 \mathcal{L}^{-1}\Big\{\frac{1}{s}\Big\} - 1 \cdot \mathcal{L}^{-1}\Big\{\frac{1}{s+1}\Big\} = 3 - e^{-t}.$$

이와 같이 서로 다른 단순근으로 이루어진 분모에 대해서는 적절한 부분분수 항을 구성하고, 항등식 계수 비교를 통해 분할 상수를 구한 뒤, 그 각각에 대한 역변환을 합산함으로써 최종적인 시간영역 함수를 얻는다.

중근을 갖는 경우

분모가 중근(repeated root)을 포함하면, 부분분수 분해는 단순 일차 항 외에 높은 차수까지 고려한다. 예를 들어, 분모에 $(s-a)^2$ 형태의 인수가 등장하면,

$$\frac{P(s)}{(s-a)^2} \quad\text{대신}\quad \frac{A}{s-a} + \frac{B}{(s-a)^2}$$

와 같이 분해한다. 분자 $P(s)$의 차수가 충분히 낮다는 전제하에서 이와 같이 각각의 거듭제곱 항에 대해 모든 차수를 포함하는 항을 따로 두는 방식으로 전개한다. 이는 서로 다른 단순근의 경우보다 조금 더 복잡한 식 전개를 야기하나, 기본 원리는 동일하다. 궁극적으로, 분해된 각각의 항은 $\frac{1}{(s-a)^n}$ 형태(또는 계수 배)로 나타나며, 표준 라플라스 역변환 식

$$\mathcal{L}^{-1}\big\{\frac{1}{(s-a)^n}\big\} = \frac{t^{n-1}}{(n-1)!} e^{at}$$

을 사용하여 시간 영역으로 옮길 수 있다. 따라서 중근 이상이 존재하더라도, 올바른 형태로 항들을 나누기만 하면 동일하게 역변환을 처리할 수 있다.

이차 인수를 갖는 경우

허근을 포함하는 이차 인수도 중요하다. 일반적으로 $s^2 + \omega^2$나 $s^2 + 2\zeta\omega s + \omega^2$ 형태의 인수는 흔히 나타난다. 이 경우에는

$$\frac{P(s)}{(s^2 + \omega^2)} = \frac{As + B}{s^2 + \omega^2}$$

와 같은 식으로 표현되어야 하며, $A$, $B$를 적절히 찾는다. 이차 인수가 공액 복소근을 포함한다면, 실제역(real domain)에서의 역변환 공식을 활용하기 위해 사인, 코사인 형태가 등장한다는 점에 유의한다. 예컨대,

$$\mathcal{L}^{-1}\big\{\frac{s}{s^2 + \omega^2}\big\} = \cos(\omega t), \quad \mathcal{L}^{-1}\big\{\frac{\omega}{s^2 + \omega^2}\big\} = \sin(\omega t).$$

과 같이 사인∙코사인 성분으로 분해되므로, 결과적으로 시간 영역에서의 운동 방정식이나 회로 방정식이 가진 주기 특성을 간단히 해석할 수 있다.

예시를 통한 분해 절차 요약

라플라스 역변환에서 부분분수 분해 절차를 체계적으로 수행하기 위해서는 분자 및 분모 차수 관계를 먼저 확인하고, 필요하다면 다항식 나눗셈을 실시하여 분자 차수를 낮춘다. 이후 분모를 가능한 모든 일차 혹은 이차 인수로 인수분해한 뒤, 각 인수에 대응하는 항들을 설정한다. 서로 다른 단순근이든 중근이든, 또는 이차 인수이든, 모든 경우에 대해 계수를 비교하면서 항들에 대한 상수(혹은 일차식)를 찾는다. 최종적으로 얻어진 각각의 항은 라플라스 변환 표에 수록된 표준 역변환 식을 통해 간단하게 시간 영역으로 전환할 수 있다.

Heaviside의 Cover-up 방법

부분분수 분해를 수행할 때, 항등식 계수 비교 외에도 Heaviside가 제안한 간편한 “cover-up” 방법을 활용할 수 있다. 이 방법은 라플라스 변환식의 분모에 포함된 일차 인수들에 대한 분해 계수를 빠르게 찾아내는 기법으로 잘 알려져 있다. 예를 들어,

$$F(s) = \frac{5s + 3}{(s+1)(s+2)}$$

와 같이 분모가 $(s+1)(s+2)$ 형태로 인수분해되어 있을 경우, 부분분수 분해는

$$\frac{5s + 3}{(s+1)(s+2)} = \frac{A}{s+1} + \frac{B}{s+2}$$

와 같은 형태가 된다. Heaviside의 cover-up 방법에 따르면 $A$는 분모에서 $(s+1)$ 인수를 제외하고 나머지 부분에 $s=-1$을 대입하여 구할 수 있다. 즉,

$$A = \left.\frac{5s+3}{s+2}\right|_{s=-1} = \frac{5(-1)+3}{-1+2} = \frac{-5+3}{1} = -2.$$

마찬가지로 $B$는 $(s+2)$ 인수를 제외한 나머지 부분에 $s=-2$를 대입하면,

$$B = \left.\frac{5s+3}{s+1}\right|_{s=-2} = \frac{5(-2)+3}{-2+1} = \frac{-10+3}{-1} = 7.$$

이로부터,

$$F(s) = -2\frac{1}{s+1} + 7\frac{1}{s+2}$$

임을 얻는다. 이 방식을 사용하면 항등식 전개 및 계수 비교를 거치지 않고도 간단하게 분해 계수를 구할 수 있다. 다만, 인수가 일차형태로 분명히 구분되어 있어야 하고, 중근 이상으로 반복되는 인수나 이차 인수에 대해서는 별도의 처리 과정을 거쳐야 한다.

중근 이상에서의 cover-up 방법

중근(repeated root)을 포함하는 경우 Heaviside cover-up 방법을 직접적으로 적용하기는 까다롭다. 예를 들어,

$$F(s) = \frac{1}{s(s+1)^2}$$

와 같은 함수를 생각해 볼 수 있다. 분모가 $s$와 $(s+1)^2$를 인수로 갖는데, $(s+1)$이 중근이므로 다음과 같은 형태로 분해해야 한다.

$$\frac{1}{s(s+1)^2} = \frac{A}{s} + \frac{B}{s+1} + \frac{C}{(s+1)^2}.$$

이때 $A$를 구하기 위해 cover-up 방식을 단순 적용하려 하면, $s=0$을 대입하는 순간 분모에 $(s+1)^2$가 남아 있어 그 결과가 $B$, $C$를 통합한 형태가 되므로 $B$와 $C$를 분리해내기 어렵다. 따라서 중근 이상의 경우에는 일반 항등식 계수 비교를 사용하거나, 부분적으로 cover-up을 적용하되 남은 항에 대해 다시 부분분수 분해를 수행하는 방식을 택하기도 한다.

중근 이상은 표준형이

$$\frac{A_1}{(s-a)} + \frac{A_2}{(s-a)^2} + \cdots + \frac{A_n}{(s-a)^n}$$

처럼 전개된다는 사실을 기억해야 하며, 각 항마다 개별 계수를 찾아낸 뒤 역변환 공식을 활용하여 시간 영역 함수를 도출한다.

다항식 나눗셈의 중요성

분자와 분모의 차수가 같거나 분자가 더 높을 때는, 먼저 다항식 나눗셈(polynomial long division)을 통해 분자를 ‘낮은 차수 + 나머지’ 형태로 만들고, 그 뒤에 남은 (나머지)/(분모) 형태를 부분분수 분해한다. 예를 들어,

$$G(s) = \frac{s^3 + 2s}{s(s+1)}$$

와 같은 경우 분자 차수가 3, 분모 차수가 2이므로 먼저 나눗셈을 해야 한다. 다항식 나눗셈을 수행하면

$$s^3 + 2s = (s^2 - s + 1)\cdot s(s+1) + (\text{나머지})$$

가 될 것으로 예상되지만, 실제로는 정확한 나머지를 구해봐야 한다. 나눗셈 과정은

$$s^3 \div (s^2 + s) = s,$$

이어서 $s \cdot (s^2 + s) = s^3 + s^2$. 이를 빼면 $-s^2 + 2s$. 다시

$$-s^2 \div (s^2 + s) = -1,$$

곱하면 $-1\cdot(s^2 + s) = -s^2 - s$. 이를 빼면 $3s$. 따라서

$$s^3 + 2s = (s - 1)\,s(s+1) + 3s.$$

결국,

$$G(s) = s - 1 + \frac{3s}{s(s+1)}.$$

즉,

$$G(s) = s - 1 + \frac{3s}{s(s+1)}.$$

이제 남은 항 $\frac{3s}{s(s+1)}$에 대해 부분분수 분해를 수행하면,

$$\frac{3s}{s(s+1)} = \frac{A}{s} + \frac{B}{s+1},$$

통분 후 $3s = A(s+1) + B s$에서 계수를 비교한다. $s$ 항에서는 $3 = A + B$, 상수 항에서는 $0 = A$. 따라서 $A=0$, $B=3$. 즉,

$$\frac{3s}{s(s+1)} = 3 \cdot \frac{1}{s+1}.$$

결과적으로,

$$G(s) = s - 1 + 3\cdot\frac{1}{s+1}.$$

이를 역변환하면 다항식 항 $s - 1$에 대한 역변환과 라플라스 표에 있는 $\frac{1}{s+1}$ 항의 역변환을 각각 구해 합산할 수 있다. 특히 $s$ 항은 시간영역에서의 미분 특성과 연관이 있으므로 주의 깊게 다뤄야 한다.

이차 인수에 대한 분해

허근을 포함한 이차 인수의 경우, $s^2 + 2\zeta \omega s + \omega^2$ 형태가 많이 등장한다. 이때 부분분수 분해 시에는 $\frac{As + B}{s^2 + 2\zeta \omega s + \omega^2}$와 같은 항이 반드시 존재한다. 이를 구체적으로 살펴보면, 원하는 계수를 찾은 뒤 $\zeta$와 $\omega$에 따라 사인 또는 코사인(혹은 그 결합) 형태의 역변환을 획득할 수 있다. 대표적으로,

$$\mathcal{L}^{-1}\Big\{\frac{s+a}{(s+a)^2 + b^2}\Big\} = e^{-a t}\cos(b t), \\ \mathcal{L}^{-1}\Big\{\frac{b}{(s+a)^2 + b^2}\Big\} = e^{-a t}\sin(b t).$$

위와 같은 공식들을 적절히 조합하여 실제 해석 문제를 해결하게 된다. 문제 상황에 따라 $b^2$ 대신 복소 루트 $\alpha \pm i\beta$를 직접 다루는 방법을 택하기도 한다. 그러나 라플라스 표준 변환 목록이 부호화해둔 형태를 잘 알고 있다면, 분모에 $\alpha^2 + \beta^2$가 나타나게끔 맞춰서 분해 계수를 구한 뒤, 사인∙코사인 형태 역변환 공식을 쉽게 적용할 수 있다.

복합 예시

아래와 같은 함수를 가정해 보자.

$$H(s) = \frac{s^2 + 2}{(s^2 + 1)(s+2)}.$$

분모는 $(s^2 + 1)$ 이차 인수와 $(s+2)$ 일차 인수의 곱이다. 부분분수 분해의 표준 형태는

$$\frac{s^2 + 2}{(s^2 + 1)(s+2)} = \frac{As + B}{s^2 + 1} + \frac{C}{s+2}.$$

분해 상수 $A, B, C$를 구하기 위해 통분 후 항등식을 전개하면,

$$s^2 + 2 = (As + B)(s+2) + C(s^2 + 1).$$

우선 $(As + B)(s+2) = A s^2 + 2As + B s + 2B = A s^2 + (2A + B)s + 2B$이므로,

$$s^2 + 2 = A s^2 + (2A + B)s + 2B + C s^2 + C.$$

계수별로 정리하면

$$s^2 + 2 = (A + C) s^2 + (2A + B)s + (2B + C).$$

계수 비교를 통해,

$$A + C = 1,\\2A + B = 0,\\2B + C = 2.$$

이 세 식을 풀면, 먼저 $2A + B = 0$에서 $B = -2A$. 이어 $A + C = 1$에서 $C = 1 - A$. 마지막으로 $2B + C = 2$를 대입하면 $2(-2A) + (1 - A) = 2$, 즉 $-4A + 1 - A = 2$, 그래서 $-5A = 1$, 결국 $A = -\frac{1}{5}$. 따라서 $B = -2A = \frac{2}{5}$, $C = 1 - \Big(-\frac{1}{5}\Big) = 1 + \frac{1}{5} = \frac{6}{5}$. 그러므로

$$H(s) = \frac{-\frac{1}{5}s + \frac{2}{5}}{s^2 + 1} + \frac{\frac{6}{5}}{s+2}.$$

결과적으로

$$H(s) = -\frac{1}{5}\frac{s}{s^2 + 1} + \frac{2}{5}\frac{1}{s^2 + 1} + \frac{6}{5}\frac{1}{s+2}.$$

라플라스 변환 표에서 $\frac{s}{s^2+1}$의 역변환은 $\cos t$, 그리고 $\frac{1}{s^2+1}$의 역변환은 $\sin t$ 임을 기억하면,

$$\mathcal{L}^{-1}\Big\{-\frac{1}{5}\frac{s}{s^2 + 1}\Big\} = -\frac{1}{5}\cos t, \\ \mathcal{L}^{-1}\Big\{\frac{2}{5}\frac{1}{s^2 + 1}\Big\} = \frac{2}{5}\sin t, \\ \mathcal{L}^{-1}\Big\{\frac{6}{5}\frac{1}{s+2}\Big\} = \frac{6}{5} e^{-2 t}.$$

을 각각 얻게 되므로 최종적으로

$$\mathcal{L}^{-1}\{H(s)\} = -\frac{1}{5}\cos t + \frac{2}{5}\sin t + \frac{6}{5}e^{-2t}.$$

이처럼 부분분수 분해를 통해 $s$-영역에서 복합적으로 얽혀 있던 항들이 시간 영역에서 명확한 해석이 가능한 형태(사인, 코사인, 지수 등)로 나타나게 된다.

복소적분과 잔여정리(Residue Theorem)의 활용

라플라스 역변환을 이론적으로 엄밀하게 정의하려면, 복소적분학에서 제시하는 Bromwich 적분(라플라스 역변환의 멜린-바른스(Mellin-Barnes) 형태)을 사용해야 한다. 이 적분은 일반적으로 복소평면 상에서 적절한 경로를 따라

$$f(t) = \mathcal{L}^{-1}\{F(s)\} = \frac{1}{2\pi i} \int_{\gamma - i\infty}^{\gamma + i\infty} e^{st} F(s) \, ds$$

와 같이 정의되는데, 모든 특이점(pole)이 오른쪽 반평면 내로 들어오지 않도록 수직선 $\Re(s) = \gamma$를 잡는다. 이 적분을 실제로 계산하기 위해서는 잔여정리(Residue Theorem)를 활용하게 되고, 결국 $F(s)$의 특이점들이 위치한 $s$-평면의 극(pole)을 찾아 각 극에서의 잔여를 구하는 방식으로 역변환이 이루어진다.

부분분수 분해는 이 과정을 단순화시키는 매우 강력한 기법이다. $F(s)$의 분모가 유리함수 형태인 경우, 모든 극을 일차 인수나 이차 인수 형태로 인수분해할 수 있다면, 각 극에서의 잔여를 빠르게 구할 수 있고 최종적으로 시간영역에서의 함수 $f(t)$가 얻어진다. 특히 잔여정리를 사용해 직접 계산하면 복잡해 보일 수 있으나, 부분분수 분해를 통해

$$F(s) = \sum \frac{\alpha_{ij}(s)}{(s-a_i)^{n_i}\dots}$$

와 같은 형태로 쪼개놓고 나면, 라플라스 변환 표에 등장하는 표준 역변환 식과 쉽게 대응된다.

Bromwich 적분을 통해 바라보면, 각 극 근처에서 $e^{st}$의 지수함수가 적분을 주도하는 형태가 되므로, 극의 계수(부분분수 분해에서의 분해 계수 또는 그와 동등한 잔여값)가 시간영역 해석에 직접 연결된다고 볼 수 있다. 예를 들어 서로 다른 단순근이 여러 개 있는 경우, 시간함수는 모든 극에 대한 지수 항들의 합으로 나타난다. 반복근(repeated pole)이 있다면 $t^n e^{at}$ 형태가 추가로 발생한다. 이차 인수, 혹은 그 이상의 복소근들이 포함되면 사인∙코사인 혹은 그 일반화 형태(감쇠진동 등)가 나타난다.

이렇듯 잔여정리를 통한 복소적분 계산과 부분분수 분해 기법은 본질적으로 같은 계산을 다른 언어로 표현한 것에 가깝다. 복소함수 해석학에서는 극에서의 잔여를 산출하는 과정을 통해 계수를 얻고, 실함수 영역에서는 부분분수 분해를 통해 같은 계수를 얻는다. 따라서 라플라스 변환에서 부분분수 분해가 왜 강력한지, 잔여정리 관점에서 쉽게 이해할 수 있다. 결국 유리함수로 표현된 라플라스 변환식에 대해서는 이 과정을 회피하기 어렵고, 두 방법은 궁극적으로 동일한 목적(극과 계수 식별)을 달성한다.

다변수 라플라스 변환에서의 부분분해 확장

일반적인 단일 변수 $s$에 대한 라플라스 변환 외에도, 편미분방정식을 취급할 때는 2차원 이상의 다중 라플라스 변환(two-sided Laplace transform)을 고려하기도 한다. 예를 들어 $s$와 $p$ 같은 두 개의 복소변수에 대해 동시에 라플라스 변환을 취하면, 분모가 $s, p$를 모두 포함하는 2변수 유리함수 형태가 될 수 있다. 이때 부분분수 분해는 각각의 변수에 대해 인수분해가 가능해야 하므로 훨씬 복잡해진다.

복수 변수 라플라스 변환의 역변환은 일반적으로

$$f(x, y) = \mathcal{L}^{-1}\{F(s, p)\} = \frac{1}{(2\pi i)^2}\iint_{\Gamma} e^{s x + p y} F(s, p)\, ds\, dp$$

와 같이 적분을 정의한다. 다만, 실제로는 각 변수별로 독립된 형태로 인수분해가 가능할 때만 부분분수 분해를 간단히 적용할 수 있고, 그렇지 않은 경우라면 전개 자체가 더 복잡해진다. 한편, 이론적으로는 2중 잔여정리를 이용하거나, 한 변수를 먼저 역변환하고 나머지를 다시 역변환하는 순차적인 접근 방법 등이 사용될 수 있다.

미분방정식의 응용 예

선형 미분방정식을 라플라스 변환으로 푸는 과정에서, 얻어진 $F(s)$를 시간영역으로 돌릴 때 부분분수 분해가 필수적이다. 예를 들어 2차 이상 미분방정식을 풀면 항상 $s^2, s^3$ 등 고차 항이 나타나는데, 해가 지수함수·사인·코사인·다항식 등으로 조합될 때, 바로 그 조합의 근거는 부분분수 분해 결과라고 볼 수 있다.

시스템 동역학이나 제어공학에서는 전달함수(transfer function)의 형태가

$$G(s) = \frac{Y(s)}{U(s)} = \frac{\text{(출력의 라플라스 변환)}}{\text{(입력의 라플라스 변환)}}$$

와 같은 유리함수 형태로 주어지는 경우가 많다. 이 전달함수를 시간영역에서의 임펄스 응답 함수나 계단응답 함수를 직접 구할 때 역시 부분분수 분해 과정을 거쳐야 한다. 특히 실제 계통에서 감쇠진동, 정상상태 오차 등을 해석할 때, 각 항이 지수적으로 사라지거나 사인∙코사인으로 진동함을 부분분수 항별로 구체적으로 구분해낼 수 있다.

--- 정리

라플라스 역변환은 본질적으로 복소적분을 통해 정의되지만, 실무적으로는 부분분수 분해를 이용한 표준 라플라스 변환 표 매칭을 통해 훨씬 간결하게 처리된다. 분모의 모든 근이 실수든 복소수든, 중근이든 아니든 간에, 분해 형태를 잘 갖추고 있으면 각각의 항이 표준형 역변환 식과 일대일로 대응되기 때문이다. 이러한 이유로 여러 분야에서 부분분수 분해는 사실상 라플라스 변환 계산의 필수 도구로 자리잡았다.

회로망 해석에서의 활용

라플라스 변환은 선형회로망(Linear Network) 해석에 필수적인 도구로 널리 알려져 있다. 임피던스나 어드미턴스의 라플라스 표현이 유리함수 형태를 띠는 경우가 많으므로, 해당 전달함수나 임펄스 응답을 구할 때 부분분수 분해가 직접적으로 활용된다. 예를 들어 RLC 직병렬 회로나 Op-amp를 포함한 선형 증폭기 회로의 입출력 특성을 해석할 때, 회로 방정식을 라플라스 영역으로 변환한 뒤, 부분분수 분해로 나눈 각 항을 역변환하여 시간영역에서의 전압이나 전류를 구하게 된다.

특히, 고차 필터 회로(예: 2차 또는 3차 이상의 필터)에서 전달함수가 다항식 비율로 표현될 경우, 차수가 높아질수록 분해 항이 복잡해진다. 그러나 체계적으로 단계별 분해를 수행하면, 각 항이 지수함수, 사인 또는 코사인 형태로 정리되어, 필터의 공진 주파수나 감쇠 비율을 명확히 파악할 수 있다. 이렇게 얻어진 시간영역 응답을 통해 회로 설계자가 공진 현상을 억제하거나 신호의 위상 특성을 조절할 수 있다.

부분분수 분해의 관점에서, 인덕터나 커패시터가 포함된 회로망은 대부분 다항식 분모에 $s$ 또는 $1/s$ 항이 섞인 형태로 나타난다. 기본 인덕터의 임피던스는 $sL$, 커패시터의 임피던스는 $1/(sC)$이므로, 회로방정식을 라플라스 영역에서 종합하면 분모 다항식이 $s$ 변수로만 구성되지는 않고, 분자에도 $s$ 항이 섞일 수 있다. 이때는 필요한 만큼 다항식 나눗셈을 한 뒤 부분분수 분해로 정리하여, 각 항을 시간영역에서 예측 가능한 전형적인 형태로 복원한다.

전달함수와 특성방정식

시스템 해석 관점에서, 전달함수 $G(s)$의 극(poles)은 시스템의 특성방정식에 해당하는 분모를 $0$으로 만드는 $s$ 값들이다. 부분분수 분해를 수행하면, 각 극에 대응하는 항들을 별도로 분리할 수 있고, 이는 시간영역에서 거동양식(mode)이 어떻게 구성되는지 직접 보여준다. 예컨대, 실수 극이 음수 실부를 가지면 지수감쇠 형태가 되며, 복소 극이 존재하면 사인∙코사인 등 진동모드가 반영된다.

제어 이론에서는 폐루프 전달함수(closed-loop transfer function)의 극과 영점(zeros)을 분석하여, 시스템이 안정적인지 혹은 특정 주파수 특성을 갖는지 등의 성능을 평가한다. 부분분수 분해를 통해 각 극에서 시간이 지남에 따라 어떠한 응답이 발생하는지 정확히 예측할 수 있으며, 임펄스 응답이나 계단 응답 등 다양한 입력 조건에 대한 출력 신호를 구체적으로 제시할 수 있다.

예를 들어,

$$G(s) = \frac{N(s)}{D(s)}$$

가 폐루프 전달함수라면, $D(s)=0$인 근들이 곧 특성근(characteristic roots)이 된다. 부분분수 분해 결과에서 $(s - a_i)$, $(s^2 + 2\zeta_i \omega_i s + \omega_i^2)$ 같은 인수들이 나타나면, 각각 $t$영역에서 $e^{a_i t}$ 혹은 $e^{-\zeta_i \omega_i t}\sin(\omega_d t)$, $e^{-\zeta_i \omega_i t}\cos(\omega_d t)$와 같은 응답 항이 생성된다. 따라서 라플라스 변환 표를 적극 활용하여, 부분분수 각각의 계수와 형태를 시간에 종속되는 해로 재구성할 수 있다.

미분방정식 계수의 물리적 의미

시스템 혹은 회로에 대한 미분방정식을 세울 때, 각 항의 계수는 저항, 인덕턴스, 커패시턴스, 질량, 스프링상수 등 물리적 파라미터를 반영한다. 예를 들어 2차 시스템인

$$m \frac{d^2 y}{dt^2} + c \frac{d y}{dt} + k y = u(t)$$

라는 방정식을 라플라스 영역으로 옮기면

$$m s^2 Y(s) + c s Y(s) + k Y(s) = U(s),$$

따라서

$$\frac{Y(s)}{U(s)} = \frac{1}{m s^2 + c s + k}.$$

이전까지 논의한 부분분수 분해를 통해 분모를 $s^2 + 2\zeta\omega s + \omega^2$ 등 적절한 형태로 맞추면, 그 계수에서 고유진동수 $\omega = \sqrt{\frac{k}{m}}$와 감쇠계수 $\zeta = \frac{c}{2\sqrt{km}}$를 바로 확인할 수 있다. 또, 이를 다시 시간영역으로 역변환하면 지수감쇠진동이나 임계감쇠, 또는 저감쇠 오실레이션이 나타난다. 결국 부분분수 분해는 이론적으로나 계산적으로나, 실제 물리계를 이해하는 열쇠 역할을 수행한다.

복소평면 해석과 BIBO 안정성

입력이 유한하고(유계) 출력이 유한(유계)한, 즉 BIBO(Bounded Input, Bounded Output) 안정성을 확인할 때도 분모의 실부가 음수인지, 혹은 복소근이 오른쪽 반평면에 존재하는지 여부가 중요하다. 부분분수 분해를 해보면, 각 항이 $e^{\alpha t}$ 혹은 $e^{-\alpha t}\sin(\beta t)$처럼 나타나며, 이때 지수항의 지수부 $\alpha$가 음수인지 양수인지에 따라 응답이 발산할지 감쇠할지가 결정된다. 제어공학에서 라우스-후르비츠(Routh-Hurwitz) 판별법을 사용하거나 근궤적(root locus) 해석을 진행해도, 결국에는 분모의 극을 파악한다는 점에서 부분분수 분해 원리와 맞닿아 있다.

BIBO 안정성을 직관적으로 살펴보려면, 시간을 무한히 큰 값으로 보냈을 때 해결계가 상한 없이 증가하는 항을 포함하는지 확인하면 된다. 부분분수 형태에서 $s=a$라는 실수 극이 양의 값을 갖거나, 복소수 극이 오른쪽 반평면에 위치한다면, 해당 항은 지수성분으로 인해 폭주하게 된다. 반면 모든 극이 왼쪽 반평면에 있으면, 응답은 결국 감쇠 혹은 유계 진동으로 남게 된다.

---

부분분수 분해 기법은 라플라스 역변환 과정을 단순화하여, 물리 시스템이나 회로망 해석, 제어 시스템 분석 등 다양한 분야에서 유용하게 쓰인다. 극과 영점을 분해하는 과정을 통해 실제 해석에서 필요한 모드(mode)들을 시간영역으로 빠르고 정확하게 투영할 수 있고, 시스템의 안정성과 동역학적 특성 역시 투명하게 확인된다. 이러한 이유로 라플라스 변환 과정에서 부분분수 분해는 언제나 핵심적인 해석 기법으로 취급된다.

시간영역의 특수 신호와 부분분수 분해

라플라스 역변환은 단순히 지수함수나 사인∙코사인 등 기본 해석 함수로 귀결되는 경우가 많지만, 물리적·공학적 맥락에서는 디랙 델타나 스텝함수(Heaviside step function) 등 특수 신호가 혼합된 형태의 문제도 자주 등장한다. 예컨대, 입력이 $u(t - a)$(단계 시프트)나 $\delta(t - a)$(임펄스 시프트)로 주어지면, 라플라스 영역에서

$$\mathcal{L}\{u(t - a)\} = \frac{e^{-as}}{s}, \quad \mathcal{L}\{\delta(t - a)\} = e^{-as},$$

와 같은 지연 요소 $e^{-as}$가 곱해진 형태를 얻게 된다. 이때에도 분수 형태의 일부에 $e^{-as}$가 포함된 상황을 보고, 나머지 부분을 부분분수로 분해할 수 있다. 지연 연산자 $e^{-as}$가 기저함수와 곱해진 상태에서도, 해당 기저함수가 이미 부분분수 분해로 얻어진 표준형 역변환을 갖고 있다면, 최종 해를 시간영역에서 시프트된 형태로 해석할 수 있다.

예를 들어,

$$F(s) = e^{-as} \cdot \frac{1}{s^2 + \omega^2}$$

이라면, 우선 $\frac{1}{s^2 + \omega^2}$ 자체에 대한 역변환이 $\sin(\omega t)/\omega$라는 사실을 알고 있으므로, $e^{-as}$가 곱해져 있으면 시간영역에서 $f(t - a)$ 형태가 되며, 구체적으로는

$$\mathcal{L}^{-1}\big\{e^{-as}\frac{1}{s^2 + \omega^2}\big\} = u(t-a)\,\frac{\sin[\omega(t - a)]}{\omega}.$$

이런 과정을 정당화하기 위해서도, 만약 $\frac{1}{s^2 + \omega^2}$를 부분분수 분해해야 할 상황이라면(물론 이미 표준형으로 간단하지만), 분해 결과에 곱해진 $e^{-as}$를 그대로 유지한 채, 나머지 항들을 각기 역변환하여 합산하면 된다. 결론적으로 시간지연 요소가 있는 경우에도 부분분수 분해 과정의 근본 골격은 변하지 않는다.

복합 지연 및 중첩 적분 문제

물리계나 회로망에 여러 개의 시점에서 임펄스나 스텝이 입력으로 들어오면, $F(s)$에 $e^{-a_1 s}, e^{-a_2 s}, \dots$가 서로 곱해진 형태가 나타날 수 있다. 한 번에 해결하기에는 복잡해 보이지만, 실제로는 각 지연 항마다 별도의 부분분수 계수 설정을 고려하면 된다. 물론 지연 항이 많을수록 분해 항이 기하급수적으로 늘어날 위험이 있지만, 원리는 동일하다. 최종적으로 시간영역에서도 여러 개의 시프트가 들어간 함수를 더한 형태가 된다.

복합적 외부입력을 갖는 시스템 해석에서 컨벌루션(Convolution) 연산을 사용하는 경우도 비슷하다. 라플라스 영역에서 곱셈으로 나타나면, 시간영역에서는 합성적 적분이 되는데, 이 역시 각 항에 대해 부분분수 분해를 시도하고, 그 뒤 역변환 표나 컨벌루션 정리를 통해 해를 구하면 된다. 대표적인 예는

$$Y(s) = F(s)\cdot G(s)$$

에서의 $\mathcal{L}^{-1}{Y(s)} = f(t) * g(t)$ 형태이다. 이때 $F(s)$나 $G(s)$를 각각 부분분수 분해하여 역변환하면 쉽게 $f(t)$와 $g(t)$를 얻을 수 있고, 이어서 시간영역 컨벌루션을 통해 최종 해석이 가능해진다.

고차(高次) 편미분방정식과 부분분수 분해

편미분방정식을 해석할 때도 라플라스 변환을 적용할 수 있다. 대표적으로 1차원 열방정식(Heat Equation)이나 파동방정식(Wave Equation)의 경계조건(초기조건)을 라플라스 영역에서 처리하면, 특정 부분분수 항이 다중으로 나타나는 경우가 있다. 변수가 시간 $t$ 뿐 아니라 공간변수 $x$도 함께 있을 때, 종종 $F(s, k)$ 형태로 변환(시간에 대한 라플라스, 공간에 대한 푸리에 변환 등)을 취해 이중 변환 공간에서 문제를 푸는 방식이 사용된다.

이 경우에도, 시간 변수에 대한 라플라스 변환으로 얻어진 표현이 여전히 유리함수 꼴을 가진다면, 부분분수 분해의 이점이 유효하다. 공간변수 측면에서도 마찬가지로 분해나 인수화를 시도한다. 예컨대 파동방정식에서 경계조건이 특정 형태의 임펄스이거나 지수함수로 주어지면, 그에 상응하는 라플라스 함수를 부분분수 형태로 나눠서 각각의 모드해(모드분리 방식)로 해석 가능하다. 실제 해석 과정은 복잡해질 수 있지만, 라플라스-푸리에 혼합변환을 취한 뒤의 역변환에서 부분분수를 찾는 기본 골자는 달라지지 않는다.

복소근의 쌍과 오일러 공식

허근을 갖는 이차 인수의 역변환에서 나타나는 사인∙코사인 항은, 오일러 공식에 의해 복소지수 함수로도 표현할 수 있다. 예를 들어 $e^{at}\sin(\omega t)$는 $\frac{1}{2i}\big(e^{(a + i\omega)t} - e^{(a - i\omega)t}\big)$이 되는데, 이는

$$s = -a \pm i\omega$$

라는 극들을 라플라스 영역에서 찾았을 때, 그 잔여(Residue)를 통해 얻어지는 결과와 일치한다. 실제로 복소함수 해석 관점에서 보면, $(s^2 + 2as + b^2) = 0$ 형태의 근을 $s = -a \pm i \sqrt{b^2 - a^2}$로 표현했을 때, 이들이 곧 지수진동항을 만들어낸다. 부분분수 분해에서 $s^2 + \omega^2$ 같은 이차 항을 만나면, $\frac{A s + B}{s^2 + \omega^2}$ 꼴로 분해하는 과정은 결국 오일러 공식으로 이어지는 것과 동일하다는 사실을 알 수 있다.

대수적 관점과 시스템 해석

라플라스 변환 자체는 시간 영역에서의 선형 미분연산을 대수적 곱셈으로 대체한다. 이에 따라 역변환도 “라플라스 영역에서의 인수분해”가 실제 미분방정식 해의 모드를 식별하는 것과 동등해진다. 단순근이든 중근이든 복소근이든, 각각 시간영역에서 지수감쇠나 진동, 혹은 다항식 곱 지수함수 형태를 만들어내며, 이는 결국 부분분수 항에서 얻어지는 계수에 의해 최종적으로 결정된다.

특히, 높은 차수의 시스템 해석이나 고차 필터, 진동계 문제에서는 부분분수 분해 단계가 수작업으로는 까다롭지만, 컴퓨터 대수 시스템(CAS)을 활용하면 자동으로 분해해준다. 그럼에도 원리를 제대로 이해해야, 유도된 분해 결과를 해석할 때 놓칠 수 있는 물리적 의미나 특정 극에서 발생하는 특이성 등을 명확히 파악할 수 있다.

기타 응용과 마찬가지 원리

부분분수 분해는 라플라스 변환에 국한된 기법이 아니다. 푸리에 변환, Z 변환(이산 시간), 멜린 변환 등에서도 유리함수 꼴로 표현되는 영역 변환식이 존재하면, 동일한 방식으로 부분분수를 쓸 수 있다. Z 변환에서의 역변환 역시 “극-영점(pole-zero)” 해석을 통해 부분분수 분해 후 표준 Z 변환 표와 매칭하는 과정이 핵심이다. 따라서 단일 주파수 영역 해석뿐 아니라 다양한 변환해석 방법에서도, 부분분수 분해는 변함없이 중요한 역할을 담당한다.