2020년 10월 7일 수요일

라플라스 변환의 성질(Properties of Laplace Transform)

[경고] 아래 글을 읽지 않고 "라플라스 변환의 성질"을 보면 바보로 느껴질 수 있습니다.


[그림 1] 교류 회로 해석에 이용된 라플라스 변환(출처: wikipedia.org)

푸리에 변환(Fourier transform)복소 영역(complex domain) 일반화인 라플라스 변환(Laplace transform)의 다양한 성질은 푸리에 변환의 성질을 기반으로 쉽게 유도될 수 있다.

 
   1. 기본(basics)   

[정의]

                  (1.1)

                  (1.2)

[선형 사상(linear mapping or linearity)]

                  (1.3)

여기서 $G(s)$는 $g(t)$의 라플라스 변환이다.

[시간 및 주파수 편이(time and frequency shifting)]

                  (1.3)

[증명]
식 (1.1)에 시간 이동된 함수를 대입한다.

                  (1.4)

주파수 이동도 식 (1.4)와 비슷하게 적분할 수 있다.
______________________________

[시간 비율 조정(time scaling)]

                  (1.5)

여기서 $a > 0$이다.

[증명]
변수 $t$를 $at$로 치환하여 적분한다.

                  (1.6)
______________________________

[켤레 복소수(complex conjugate)와 대칭성(symmetry)]

                  (1.7)

                  (1.8)

여기서 $(\cdot)^*$는 켤레 복소수, 식 (1.7)과 (1.8)은 각각 $f(t)$가 복소수 및 실수인 경우이다.

[미분(differentiation)]

                  (1.9)

여기서 $f^{(n)}(t)$는 $f(t)$의 $n$번 미분, $0^{-}$는 0보다 작은 위치에서 0으로 한없이 접근한다는 표현이다.

[증명]
라플라스 변환의 정의에 대해 부분 적분(integration by parts)을 적용한다.

                  (1.10)

여기서 $s$의 실수부는 $0$보다 크다고 가정하므로 $t$가 커질 때 $e^{-st}$는 $0$이 된다. 고계 미분에 대해서는 식 (1.9)의 첫째식을 연속적으로 적용한다. 그러면 식 (1.9)의 둘째식을 얻을 수 있다.
______________________________

식 (1.9)에서 특이한 점은 $t$ = $0$의 함수값을 $f(0^-)$로 정의한 부분이다[1]. 함수 $f(t)$가 연속인 경우는 $f(0^+)$이든 $f(0^-)$이든 같아서 문제가 되지 않는다. 하지만 $f(t)$가 $t$ = $0$에서 불연속이면 식 (1.9)처럼 $f(0^-)$로 정확히 기술해야 한다. 이런 헷갈림을 피하기 위해 $t$ = $0$ 근방의 불연속 특성까지 고려해서 라플라스 변환을 다음처럼 다시 정의한다[1].

                  (1.11)

예를 들어, 식 (1.1) 대신 식 (1.11)을 사용하면, $f(t)$ = $\delta(t)$의 라플라스 변환은 논란의 여지 없이 $1$이 될 수 있다. 왜냐하면 $t$ = $0^-$ 다음에 $t$ = $0$이 나오기 때문이다. 식 (1.1)의 정의를 이용해 $\delta(t)$를 라플라스 변환할 때는 애매한 부분이 있다. 적분의 시작점이 $t$ = $0$을 포함하는지 하지 않는지가 불분명해서 라플라스 변환이 $1$이 될 수도 $0$이 될 수도 있다.
새로운 정의인 식 (1.11)의 라플라스 변환을 쓰면, 적분은 $t$ = $0^-$에서 시작하기 때문에 식 (1.9)의 첫째식은 당연하다. 혹시라도 우리의 추론에 미심쩍은 부분이 있을 수도 있다. 그래서 단위 계단 함수(unit step function) $u(t)$를 이용해 $t$ = $0$에서 불연속인 함수 $f(t)$를 $g(t)$로 다시 정의한다.   

                  (1.12)

단위 계단 함수의 성질에 의해 $u(t)$의 미분은 디랙 델타 함수(Dirac delta function) $\delta(t)$이다. 이 결과를 식 (1.11)에 대입해 정리한다.

                  (1.13)

결국 의심스러운 부분 없이 명확하게 식 (1.9)의 첫째식이 불연속의 경우에도 잘 성립함을 알 수 있다.

[적분(integration)]

                  (1.14)

[증명]
식 (1.14)에 나온 적분을 식 (4.1)의 길쌈(convolution)으로 바꾸어서 식 (4.3)과 같은 라플라스 변환을 적용한다.

                  (1.15)

여기서 $u(t)$는 단위 계단 함수(unit step function)이다.
______________________________

[특수한 변환값]

                  (1.16)

여기서 라플라스 변환은 명확한 식 (1.11)을 사용하며 $s$의 실수부는 $0$보다 크다.

[초기값 정리(initial value theorem)]

                  (1.17)

[증명]
식 (1.1)에 $s$를 곱해서 부분 적분을 한다.

                  (1.18)

식 (1.18)의 최종 결과에 $s \to \infty$를 적용하면, 적분 항은 $0$이 되므로 남는 항은 $f(0)$이다. 따라서 연속 함수인 경우는 식 (1.17)이 잘 성립한다. 만약 $t$ = $0$이라면 결과는 어떻게 될까? 함수 $f(t)$를 식 (1.12)에 있는 $g(t)$로 바꾸고 식 (1.11)을 이용해 라플라스 변환한다.

                  (1.19)

따라서 불연속 함수의 경우에도 식 (1.17)이 잘 성립한다.
______________________________

미분의 라플라스 변환인 식 (1.9)의 첫째식과 식 (1.17)을 비교하면 $t$ = $0$에서 접근하는 $f(t)$의 극한이 좌극한 $f(0^-)$와 우극한 $f(0^+)$로 서로 다르다. 이런 애매한 부분의 해결책은 명확한 라플라스 변환인 식 (1.11)이다.

[최종값 정리(final value theorem)]

                  (1.20)

여기서 $F(s)$의 모든 극점(pole)은 실수부가 $0$보다 작은 조건[혹은 좌반면(左半面, left half-plane) 조건]을 만족해야 한다.

[증명]
식 (1.19)의 마지막식에 $s \to 0$ 조건을 적용해서 정리해본다.

                  (1.21)

그러면 식 (1.21)이 쉽게 증명되지만 조심할 부분도 있다. 극한 $s \to 0$ 조건에서 식 (1.21)의 결과가 $sF(s)$라는 보장이 있는가? 즉 $s \to 0$일 때, $f(t)$와 $F(s)$는 라플라스 변환 관계를 가지고 있을까? 이를 위해 식 (1.2)에 쓴 라플라스 역변환을 고려한다.

[그림 1.1] 라플라스 역변환을 위한 닫힌 경로

라플라스 변환 $F(s)$가 함수 $f(t)$로 연결되려면 [그림 1.1]에 있는 닫힌 경로 상의 복소 적분(complex integral)이 잘 정의되어야 한다. 시간 $t \to \infty$인 경우는 복소 적분의 수렴을 위해 닫힌 경로로 $c_1 + c_2$를 택한다. 여기서 $c_1$ 상의 $\Re[s]$ = $\sigma$이며 식 (1.2)에 의해 $s \to 0$일 때 $\sigma \to 0$이 된다. 그러면 닫힌 경로 안에 $F(s)$의 모든 극점이 존재해야 $f(t)$가 라플라스 역변환이 된다. 하지만 $s \to 0$인 조건에서는 $F(s)$의 극점이 닫힌 경로 안에 없을 수도 있다. 따라서 식 (1.20)이 성립하기 위해 $F(s)$의 모든 극점은 반드시 실수부가 음수인 위치에 있어야 한다.
______________________________

최종값 정리와 다르게 식 (1.17)에 있는 초기값 정리는 $s \to \infty$로 간다. 그래서 $F(s)$의 극점은 어디에 있든지 관계가 없다.

[주기 함수(periodic function)]

                  (1.22)

여기서 $\Re[s] > 0$, $f(t)$는 주기 $T$를 가진 주기 함수이다.

[증명]
주기 함수의 특성인 $f(t)$ = $f(t+T)$를 이용해 적분을 $T$로 나누어 다시 적분한다.

                  (1.23)

여기서 공비 $e^{-sT}$의 크기가 $1$보다 작기 때문에 무한 등비 급수는 수렴한다.
______________________________

[거듭제곱 함수(power function)]

                  (1.24)

여기서 $n$은 0 혹은 자연수이다.

[증명]
라플라스 함수의 정의인 식 (1.1)을 $s$에 대해 미분해서 유도한다.
______________________________

단순한 거듭제곱 함수 $t^n$의 라플라스 변환은 식 (2.1)에 따라 $1/s$를 $n$번 미분하고 부호를 $(-1)^n$으로 결정하면 된다. 이 결과는 $\nu$ = $n$인 식 (2.4)와 동일하다.


   2. 초등 함수의 변환(transform of elementary functions)   

[단위 계단 함수(unit step function)]

                  (2.1)

여기서 $\Re[s] > 0$이다. 라플라스 변환 $F(s)$의 극점이 $s$ = $0$에 있지만, 예외적으로 식 (1.20)의 최종값 정리가 성립한다. 즉 $\lim_{s \to 0} s F(s)$ = $1$은 $u(t)$의 최종값이다.

[지수 함수(exponential function)]

                  (2.2)

여기서 $\Re[s] > -a$이다. 만약 $a > 0$이면, $F(s)$의 극점이 좌반면[$\Re[s] < 0$]에 있어서 $e^{-at}$의 최종값은 $0$이 된다. 하지만 $a < 0$이라면, 우반면[$\Re[s] > 0$]에 $F(s)$의 극점이 있어서 최종값 정리가 성립하지 않는다.

[삼각 함수(trigonometric function)]

                  (2.3)

여기서 $\Re[s] > 0$이다. 라플라스 변환 $F(s)$의 극점이 허수축에 있어서 최종값 정리가 성립하지 않는다.

[거듭제곱 함수(power function)]

                  (2.4)

여기서 $\Re[s] > 0$이다.

[증명]
감마 함수(gamma function)의 정의를 이용하여 식 (1.1)의 적분을 정리한다.
______________________________


   3. 특수 함수의 변환(transform of special functions)   

[베셀 함수(Bessel function)] [2]

                  (3.1)

여기서 $\Re[\nu] > -1$이다.

[증명]
제1종 베셀 함수를 일반적으로 정의하는 쉴레플리–좀머펠트 적분(Schläfli–Sommerfeld integral)부터 출발한다.

                      (3.2)

여기서 $\mathcal{R}$은 복소 평면에 정의된 사각형 경로이다. 식 (3.2)의 우변에 라플라스 변환을 적용해서 정리한다.

                      (3.3)

여기서 사각형 경로상에서 $\sinh u$의 실수부는 항상 음수, 무한대에서 발산하지 않도록 $e^{-(\nu+1)u}$의 지수는 $\Re[\nu] > -1$인 조건을 가져야 한다. 식 (3.3)의 피적분 함수의 유수(residue)는 $u$ = $\sinh^{-1} s$이므로, 식 (3.3)은 유수 정리(residue theorem)를 이용해 계산한다.

                      (3.4)

식 (3.4)에 역쌍곡 사인 함수(inverse hyperbolic sine function)의 로그 함수 표현식 $\sinh^{-1} s$ = $\log (\sqrt{s^2 + 1} + s)$를 대입해서 예쁘게 만든다.

                      (3.5)

식 (3.5)에 다시 식 (1.5)를 사용해서 증명을 완성한다.
______________________________

식 (3.5)에서 $s$ = $0$으로 놓으면 베셀 함수의 적분은 항상 1로 수렴한다.

                  (3.6)

통상적으로 베셀 함수의 정적분은 구하기 어렵지만, 라플라스 변환 관점의 적분은 조금 수월하다.


   4. 길쌈(convolution)   

[그림 4.1] 감쇠 지수 함수의 길쌈 예(출처: wikipedia.org)

[정의]

                  (4.1)

길쌈(convolution)은 [그림 4.1]처럼 한 신호를 뒤집어서 다른 신호와 곱하기 때문에, 신호의 응답(response)을 모형화할 때 매우 유용하게 사용된다.

[대수적 성질(algebraic properties)]

                  (4.2)

[길쌈 정리(convolution theorem)]

                  (4.3)

여기서 둘째식에 나오는 적분은 $\Re[s] > 2 \sigma$, $\sigma$는 $F(s)$와 $G(s)$의 모든 극점을 포함하도록 정하며, $u$의 적분 경로 상에서 $F(s)$와 $G(s)$는 잘 정의되며, 복소 영역으로 확장한 $u$의 적분 경로는 반시계 방향으로 설정한다.

[증명]
길쌈인 식 (4.1)의 둘째 줄을 식 (1.1)에 대입해서 차례로 정리한다.

                  (4.4)

그러면 식 (4.3)의 첫째식이 증명된다. 식 (4.3)의 둘째식을 위해서는 식 (1.2)를 이용해 $f(t), g(t)$를 바꾸어 쓰면 된다.

                  (4.5)

여기서 지수 함수의 수렴을 위해 $\Re[s] > 2 \sigma$로 선택하며, $t \ge 0$이라서 $u, v$의 적분 경로는 반시계 방향으로 돈다. 라플라스 변환의 수렴 조건에 의해 다음 관계도 성립한다.

                  (4.6)

식 (4.6)을 이용해 식 (4.5)를 정리하면, $u$에 대한 피적분 함수는 $F(u) G(s-u)$가 된다.
______________________________

[파르세발의 정리(Parseval's theorem)]

                  (4.7)

[증명]
식 (1.7)을 식 (4.3)의 둘째식에 넣고 $s$ = $0$으로 고정해서 정리하면 식 (4.7)을 유도할 수 있다.
______________________________

식 (4.7)을 더 간단히 하기 위해 $g(t)$ = $f(t)$로 가정하면, 새로운 파르세발의 정리가 나온다.

                  (4.8)

함수 $f(t)$가 실수인 경우는 다음처럼 더 간략화할 수도 있다.

                  (4.9)

예를 들어, $f(t)$ = $u(t)$, $g(t)$ = $e^{-at}$로 가정해서 파르세발의 정리를 적용한다. 먼저 식 (4.7)의 좌변은 간단히 계산된다.

                  (4.10)

여기서 $a > 0$이다.

[그림 4.2] 파르세발의 정리를 적용하기 위한 닫힌 경로

라플라스 변환인 $F(s)$ = $1/s$, $G(s)$ = $1/(s+a)$를 식 (4.7)의 우변에 대입해서 [그림 4.2]와 같은 복소 영역에서 적분을 한다.

                  (4.11)

여기서 극점(pole)은 $s$ = $0$과 $a$, 적분 경로 $c_1 + c_2$에 포함되는 극점은 $s$ = $0$이다.


[참고문헌]
[1] K. H. Lundberg, H. R. Miller, and D. L. Trumper, "Initial conditions, generalized functions, and the Laplace transform troubles at the origin," IEEE Control Systems Magazine, vol. 27, no. 1, pp. 22–35, Feb. 2007. (방문일 2020-10-09)
[2] "Laplace transform of Bessel function of the first kind," ProofWiki. (방문일 2023-11-12)

[다음 읽을거리]

2020년 10월 5일 월요일

라플라스 변환(Laplace Transform)

[경고] 아래 글을 읽지 않고 "라플라스 변환"을 보면 바보로 느껴질 수 있습니다.


[라플라스 변환 소개]

푸리에 변환(Fourier transform)은 단순 정적분(definite integral)에 해당하는 유한 구간 적분에 강하다. 만약 $t$ 혹은 $\omega$에 대한 적분 구간이 무한(infinite)[$-\infty < t < \infty$]이나 반무한(semi-infinite)[$t \le c$ 혹은 $t \ge c$]으로 가면, 푸리에 변환 혹은 역변환은 적분이 잘 구해지지 않는다.

                  (1)

무한한 적분 구간에서 식 (1)의 결과를 안정되게 구할 때는 감쇠를 이용한 극한(limit)이나 복소 함수론(complex analysis)에 기반한 경로 적분을 적극적으로 써야 한다. 하지만 이러한 방법론은 너무 수학적이라서 쉽게 쓰려고 도입한 적분 변환(integral transform)의 취지에 적합하지 않다. 푸리에 변환이 유용해서 버리기는 아깝지만 무한 영역의 적분에서 문제가 분명히 생긴다. 이러한 애매함을 감쇠라는 현실적인 양으로 해결한 엄청난 도구가 바로 라플라스 변환(Laplace transform)이다.

                  (2)

여기서 $s$ = $\sigma + i \omega$, $\sigma$는 감쇠(減衰, attenuation)이다. 라플라스 변환은 천제 현상을 기술하는 미분 방정식(differential equation)을 풀기 위한 라플라스Pierre-Simon Laplace(1749–1827) 노력의 결정체이다[2]. 후에 헤비사이드Oliver Heaviside(1850–1925)라플라스 변환을 재발견하고 현재와 같은 쓰임새를 만들어냈다.
식 (2)에 도입한 $\sigma$가 너무 작위적이라 생각할지도 모른다. 그러나 현존하는 모든 시스템은 크던 작던 감쇠가 있기 때문에, 감쇠없이 항상 진동한다는 푸리에 변환이 오히려 현실에 잘 맞지 않고, 시간에 따라 감쇠한다는 라플라스 변환이 더 현실적이다.

[그림 1] 시스템 과도 응답의 예(출처: wikipedia.org)

라플라스 변환은 적분이 $0$부터 시작하므로 [그림 1]과 같은 과도 응답(transient response)의 분석에 적합하다. 푸리에 변환은 적분의 시작과 끝이 무한대라서 사실 응답의 시작과 끝이 없다. 그래서 푸리에 변환은 정현파처럼 끊임없이 반복적인 정상 상태 응답(steady-state response)의 분석에 좋다. 라플라스 변환과 푸리에 변환의 관계는 다음과 같다.

                  (3)

여기서 $u(t)$는 단위 계단 함수(unit step function)이다. 푸리에 변환이 실패한 무한 적분 유도가 라플라스 변환에서는 연습 문제 수준인 이유를 식 (3)에서 명확히 볼 수 있다. 라플라스 변환에는 적분을 강제로 수렴시키는 감쇠 $\sigma$와 함께 단방향 적분을 만드는 단위 계단 함수 $u(t)$가 필수 요소이다. 그래서 단위 계단 함수의 또 다른 이름인 헤비사이드 계단 함수(Heaviside step function)에 라플라스 변환의 제안자인 헤비사이드Oliver Heaviside(1850–1925) 이름이 들어간다.
라플라스 변환의 위력을 느끼기 위해 $u(t)$의 라플라스 변환을 구해보자.

                  (4)

반면에 $u(t)$에 대한 푸리에 변환을 구하려면, 로렌츠–코쉬 함수(Lorentz–Cauchy function)의 극한을 쓰거나 복소 함수론에 나오는 경로 적분(contour integral)을 사용해야 한다.

                  (5)

복소수 $s$ = $i \omega$로 두면 두 결과는 매우 비슷하다. 하지만 무한 적분을 계산할 때 사용한 수학 수준은 하늘과 땅만큼 차이난다.
식 (2)에 소개한 라플라스 변환은 푸리에 변환만큼 흔하게 볼 수 있는 적분 변환이다. 하지만 그 역인 라플라스 역변환(inverse Laplace transform)은 어떻게 표현할까? 아마 라플라스 역변환은 거의 본 적이 없을 것이다. 그러나 비슷하지만 다른 푸리에 역변환은 푸리에 변환만큼 자주 사용된다. 이런 차이는 어디서 오는 것일까? 라플라스 역변환과 관련된 멱급수(power series)의 수렴 조건과 원래 함수로의 일대일 대응 문제는 헤비사이드가 제안한 라플라스 변환이 수학자들에게 공격받은 주요 약점이었다. 이런 상황에서도 공학자 출신인 헤비사이드는 자신만의 변환 개념을 그대로 밀어붙여 미분 방정식의 놀라운 해법인 연산 미적분학(operational calculus)을 공학적으로 완성시켜 버렸다. 연산 미적분학과 복소 함수론(complex analysis)이 결합하면 자연스럽게 라플라스 변환이 만들어진다. 즉 복소 함수론으로 연산 미적분학을 완벽하게 개선한 결과물이 우리가 쓰는 라플라스 변환이다. 연산 미적분학의 수학적 엄밀성과 라플라스 역변환의 존재성을 헤비사이드 자신도 증명하지 못했지만 괜찮았다. 수학적 기반이 부족해도 헤비사이드의 연산 미적분학은 실무에서 완벽하게 잘 동작했다. 그래서 자신이 만든 연산 미적분학이 미분 방정식을 쉽게 푸는 새로운 세상을 열었다고 믿었다. 실제로도 공학자의 공부 분량을 현저하게 줄인 주요 개념이 연산 미적분학과 라플라스 변환이다. 헤비사이드가 그토록 바라던 연산 미적분학의 수렴성과 라플라스 역변환은 브롬위치Thomas John I'Anson Bromwich(1875–1929)가 복소 함수론을 도입하여 완전하게 해결했다. 이런 이유로 라플라스 역변환을 브롬위치 적분(Bromwich integral)으로도 부른다. 결국 복소 함수론에 의해 완전체가 된 연산 미적분학과 라플라스 변환은 수학자들의 뜨거운 찬사를 받으며 19세기 수학의 전설이 되었다.

[라플라스 역변환 혹은 브롬위치 적분]

                  (6)

여기서 $\sigma$는 복소 적분이 수렴하는 위치로 정한다.

[증명]
식 (6)의 우변에 식 (2)를 넣어서 적분을 정리한다.

                  (7)

[그림 2] 브롬위치 적분을 위한 닫힌 경로

만약 식 (7)에서 $t \ne t'$이라면, 식 (7)의 복소 적분은 [그림 2]와 같은 닫힌 경로의 내부에 극점(pole)이 없다. 그래서 다음 적분은 코쉬의 적분 정리(Cauchy's integral theorem)에 의해 항상 $0$이다.

                  (8)

여기서 반원 경로의 반지름 $R$을 무한대로 보낸다. 시간이 $t > t'$와 $t < t'$일 때는 경로를 각각 $c_2$와 $c_3$로 택해서 닫힌 경로의 내부 영역이 잘 수렴하도록 한다. 이러한 조건에서 $R$을 계속 키우면, 푸리에 변환과 디랙 델타 함수의 증명과 동일하게 $c_2$와 $c_3$ 경로 상의 복소 적분은 각각 $0$으로 수렴한다. 따라서 $t \ne t'$이라면, 식 (8)에 의해 경로 $c_1$에 대한 적분도 당연히 $0$이 된다. 이 결과를 사용해 식 (7)에 있는 $t'$의 적분 영역을 $t - \Delta \le t \le t + \Delta$로 제한해도 된다. 시간 $t'$의 적분 영역을 매우 축소한 경우, 실수 $t'$와 복소수 $s$에 대한 적분을 다시 해본다.

                  (9)

닫힌 경로 $c_1 + c_2$ 내부에는 극점이 $z = 0$에 있어서 유수 정리(residue theorem)를 사용할 수 있다. 반대로 $c_1 + c_3$ 내부는 해석적이라서 극점이 없다. 그래서 이 경로를 사용한 복소 적분은 항상 $0$이다. 식 (8)과 (9)를 종합해서 다음과 같은 디랙 델타 함수(Dirac delta function)의 관계식을 얻는다.

                  (10)

식 (10)을 식 (7)에 대입하면, 쉽게 식 (2)의 결과를 증명할 수 있다. 
______________________________

브롬위치 적분의 증명을 고찰하면, 브롬위치 이전 수학자들이 라플라스 역변환을 얻지 못했던 이유를 분명히 볼 수 있다. 푸리에 변환과 역변환의 적분은 실수 영역에서 둘 다 정의되어서 굉장히 상식적이다. 하지만 라플라스 변환은 실수 영역에서 1차원적 정적분(definite integral)을 하고 라플라스 역변환은 복소 평면에서 2차원적 선 적분(line integral)을 한다. 변환과 역변환이 서로 다른 영역과 적분 방식을 가져도 되는가? 라플라스 변환과 역변환이 가진 이러한 비대칭적 특성은 우리 상식에 반하지만, 식 (6)처럼 라플라스 역변환도 잘 정의된다. 그래서 헤비사이드가 처음에 가졌던 확신이 타당함을 브롬위치가 복소 함수론으로 확실하게 증명했다.
간단한 처리가 가능한 라플라스 변환과 달리 수치 계산이 어려운 라플라스 역변환은 의외로 다양한 응용을 가지고 있다[3]. 그래서 요즘은 적분 핵심(integral kernel) $e^{st}$를 쌍곡선 함수(hyperbolic function)로 근사화해 계산을 가속하는 고속 라플라스 역변환(fast inverse Laplace transform, FILT)을 많이 연구하고 있다[4]. FILT의 연구 관점은 단순화[계산 공식의 간략화]와 병렬화[동시 계산으로 시간 절약]가 중요하다.

브롬위치 적분이 라플라스 역변환임을 더 극적으로 보기 위해, 평행 이동된 단위 계단 함수 $u(t-t_0)$의 라플라스 변환을 생각한다.

                  (11)

여기서 수렴을 위한 조건은 $\Re[s] > 0$이다. 식 (11)의 결과를 식 (6)에 대입해서 복소 적분을 계산한다. 만약 $t > t_0$이라면, [그림 2]에서 닫힌 경로를 $c_1+c_2$로 택해서 유수(residue)를 구한다. 

                  (12)

반대로 $t < t_0$인 경우, 피적분 함수를 수렴시키기 위해 닫힌 경로는 $c_1+c_3$가 되어야 한다. 그러면 닫힌 경로 내부에 극점이 없어서 복소 적분은 $0$이 된다. 이 두 결과를 종합하면, $e^{-s t_0}/s$에 대한 식 (2)의 적분은 $u(t-t_0)$이다. 따라서 브롬위치 적분은 정확히 라플라스 역변환이다. 이상의 결과를 명확하게 명제로 표현하면 다음과 같다.

[라플라스 변환의 일대일 대응(one-to-one mapping)]
두 함수의 라플라스 변환이 같으면 두 함수는 유한한 점에서만 다르고 나머지 영역에서는 같다.

                  (13)

[증명]
라플라스 변환의 역변환이 존재하기 때문에 일대일 대응은 당연히 성립한다. 다만 라플라스 변환은 리만 적분(Riemann integral)으로 정의하므로, 유한한 점의 적분은 기여가 없다. 그래서 두 함수는 유한한 점에서 다를 수 있다.
______________________________

두 함수가 무한한 점에서 다를 경우는 디랙 델타 함수(Dirac delta function)가 있다는 뜻이다. 그래서 두 함수의 라플라스 변환은 같을 수 없다. 다시 말해 무한한 점에서 다르면 라플라스 변환의 일대일 대응이 잘 성립한다. 다만 이런 조건이 성립하려면 라플라스 변환이 잘 존재해야 한다. 식 (11)의 예에서 보듯이, 라플라스 변환을 구성하는 $s$에 따라 적분이 존재할 수도 혹은 없을 수도 있다. 식 (6)에 정의한 브롬위치 적분을 잘 조합하면 라플라스 변환에 대한 정확한 수렴 조건(convergence condition)을 얻을 수 있다.

[라플라스 변환의 수렴 조건]
라플라스 변환 $F(s)$의 극점보다 우반면(right half-plane)에 있는 $s$를 택해야 라플라스 변환이 수렴한다.

                  (14)

여기서 $\Re[s] > \sigma$이다.

[증명]
식 (6)을 식 (2)에 대입해서 정리한다.

                  (15)

여기서 지수 함수를 수렴시키기 위해 $\Re[s] > \sigma$인 조건을 택한다. 식 (15)의 마지막식은 아래에 있는 코쉬의 적분 공식(Cauchy's integral formula)과 매우 비슷하다.

                  (16)

식 (16)이 성립하기 위해서는 닫힌 경로 내부에서 $f(z)$가 해석적이어야 한다. 그래서 식 (15)에 있는 경로를 [그림 2]에 있는 $c_1 + c_3$으로 잡는다. 또한 $c_1 + c_3$ 내부에 $F(u)$의 극점이 있으면 안되므로, $F(u)$의 극점은 $s$의 좌반면(left half-plane)에 있도록 라플라스 변환의 $s$를 정한다. 그러면 코쉬의 적분 공식에 의해 식 (15)는 $F(s)$가 된다.
______________________________

푸리에 변환과 역변환은 동일한 실수 영역에서 적분이 이루어져서 증명 과정이 비슷하다. 하지만 라플라스 변환은 역변환과 적분 영역이 아주 달라서 상호간의 증명 과정도 서로 차이가 많이 난다. 라플라스 역변환에 변환을 대입한 식 (7)은 디랙 델타 함수(Dirac delta function)를 이용해 증명한다. 반면에 라플라스 변환에 역변환을 넣은 식 (15)의 경우에는 코쉬의 적분 공식(Cauchy's integral formula)이 필수적이다.
식 (2)에 있는 통상적인 라플라스 변환은 적분이 $0$에서 시작하지만, 푸리에 변환처럼 $t$에 대한 적분 구간이 모든 실수 영역이 되게 라플라스 변환을 양방향으로 정의할 수도 있다.

                  (17)

식 (17)과 같은 적분 변환은 양방향 라플라스 변환(bilateral Laplace transform or two-sided Laplace transform)이라 한다. 식 (17)에 대비되는 식 (2)는 단방향 라플라스 변환(unilateral Laplace transform or one-sided Laplace transform)으로 부르기도 한다. 우리가 $t$에 대한 적분 구간을 모든 실수로 늘리기는 했지만, 양방향 라플라스 역변환(inverse bilateral Laplace transform)은 식 (6)과 동일하다. 왜냐하면 식 (10)처럼 $t$가 아닌 $s$에 대한 적분이 디랙 델타 함수를 생성하기 때문이다.

                  (18)

양방향 라플라스 변환은 다른 적분 변환과 다음 관계를 가지고 있다.

                  (19)

                  (20)

여기서 $\mathfrak{F}[\cdot]$는 푸리에 변환을 의미한다. 따라서 양방향 라플라스 변환은 푸리에 변환을 포함한 다양한 적분 변환을 만들 수 있는 범용성이 있다.

[참고문헌]
[1] H. Jeffreys, "Bromwich's work on operational methods," J. London Math. Soc., vol. 3,  220–223, Jul. 1930.
[2] 김계환, 김성숙, "라플라스의 생애와 현대과학에 미친 영향", 한국수학사학회지, 제32권, 제6호, pp. 271–279, 2019년 12월.
[3] T. Hosono, "Numerical inversion of Laplace transform and some applications to wave optics", Radio Sci., vol. 16, no. 6, pp. 1015–1019, Nov./Dec. 1981.
[4] D. Wu, S. Kishimoto, and S. Ohnuki, "Optimal parallel algorithm of fast inverse Laplace transform for electromagnetic analysis," IEEE Antennas Wirel. Propag. Lett., vol. 19, no. 12, pp. 2018–2022, Dec. 2020.

[다음 읽을거리]