# 물리학적 예제 (진동, 전기회로 등)

#### 감쇠 및 외력에 의한 진동의 해석

질량이 m, 감쇠계수가 c, 스프링 상수가 k인 1차원 계에서 변위 x(t)가 외부 힘 $F(t)$에 의해 작용받을 때, 운동방정식은

$$
m \frac{d^2 x}{dt^2} + c \frac{dx}{dt} + k x = F(t)
$$

형태로 정리된다. 라플라스 변환을 적용하기 위해 초기조건을 $x(0) = x\_0, \frac{dx}{dt}(0) = v\_0$라 할 때, 해당 식의 라플라스 변환은

$$
m \left\[ s^2 X(s) - s x\_0 - v\_0 \right] + c \left\[ s X(s) - x\_0 \right] + k X(s) = F(s)
$$

가 된다. 여기서 $X(s)$는 $x(t)$의 라플라스 변환, $F(s)$는 $F(t)$의 라플라스 변환이다. 이를 정리하면

$$
\left( m s^2 + c s + k \right) X(s) - m \left( s x\_0 + v\_0 \right) - c x\_0 = F(s)
$$

이 되고, 따라서

$$
X(s) = \frac{m \left( s x\_0 + v\_0 \right) + c x\_0 + F(s)}{m s^2 + c s + k}
$$

형태를 얻는다. 이후 역라플라스 변환을 통해 $x(t)$를 구한다.

감쇠계수를 고려하지 않는 단순 조화진동자의 경우 $c=0$이므로

$$
m \frac{d^2 x}{dt^2} + k x = F(t)
$$

가 되며, 라플라스 변환 후에는

$$
m \left\[ s^2 X(s) - s x\_0 - v\_0 \right] + k X(s) = F(s)
$$

가 되고

$$
X(s) = \frac{m \left( s x\_0 + v\_0 \right) + F(s)}{m s^2 + k}
$$

로 단순화된다. 적절한 $F(t)$에 대한 $F(s)$를 대입하고 역변환하면 해가 구해진다.

감쇠 항이 존재하는 경우 $c \ne 0$이 되면 계의 해는 지수 감쇠 인자를 갖는 부분과 스프링에 의한 진동 항을 결합한 형태로 나타나게 된다. 예를 들어 외력이 $F(t) = F\_0 \cos(\omega t)$일 때, 라플라스 변환은

$$
F(s) = \frac{F\_0 , s}{s^2 + \omega^2}
$$

가 되므로

$$
X(s) = \frac{m \left( s x\_0 + v\_0 \right) + c x\_0 + \frac{F\_0 , s}{s^2 + \omega^2}}{m s^2 + c s + k}
$$

를 얻게 된다. 이를 부분분수 전개 후 역라플라스 변환하면, 동특성과 강제진동특성을 포함한 물리적인 해석이 가능하다. 계수가 특정 관계를 만족하면 공진 현상이나 임계감쇠 등의 조건도 이 해에서 유도할 수 있다.

#### RLC 직렬 회로의 해석

회로이론에서 저항 R, 인덕터 L, 커패시터 C가 직렬로 연결되고 이에 걸리는 회로전압을 $v(t)$라고 할 때, 회로에 흐르는 전류 $i(t)$를 구하기 위한 기본방정식은

$$
L \frac{di}{dt} + R i + \frac{1}{C} \int i(t) , dt = v(t)
$$

가 된다. 라플라스 변환을 적용하기 전, 초기 커패시터 전압을 $v\_C(0) = V\_0$라 두면

$$
\frac{1}{C} \int i(t) , dt = \frac{1}{C} q(t)
$$

에서 $q(t)$가 커패시터에 축적된 전하이고, $q'(t) = i(t)$ 관계가 성립한다. 초기 전하가 $q(0) = C V\_0$임을 사용하면, 적분항은 라플라스 영역에서

$$
\frac{1}{C} \left\[ \frac{I(s)}{s} + q(0) \frac{1}{s} \bigg|\_{q(0)=C V\_0} \right]
$$

로 나타난다. 식을 더 체계적으로 쓰기 위해 미분항과 적분항을 모두 라플라스 변환하면

$$
L \left\[ s I(s) - i(0) \right] + R I(s) + \frac{1}{C} \left\[ \frac{I(s)}{s} + \frac{C V\_0}{s} \right] = V(s)
$$

가 되며, 이를 정리하면

$$
\left( L s + R + \frac{1}{sC} \right) I(s) - L i(0) + \frac{V\_0}{s} = V(s)
$$

가 된다. 보통은 초기전류 $i(0) = 0$이나 초기전압 $V\_0 = 0$ 같은 단순 초기조건으로 가정하는 경우가 많지만, 라플라스 변환은 이 두 초기조건이 동시에 비영(非零)일 때에도 문제를 정교하게 다룰 수 있게 해준다.

시스템이 단순화된 경우, 예를 들어 $V\_0 = 0, i(0) = 0$이라고 가정하고 정현파 전압원 $v(t) = V\_0 \cos(\omega t)$를 인가한다면

$$
V(s) = \frac{V\_0 , s}{s^2 + \omega^2}
$$

이며, 회로의 임피던스에 해당하는 라플라스 표현은

$$
Z(s) = L s + R + \frac{1}{s C}
$$

이므로

$$
I(s) = \frac{\frac{V\_0 , s}{s^2 + \omega^2}}{L s + R + \frac{1}{s C}}
$$

형태가 되고, 다시 부분분수 전개 후 역라플라스 변환을 취하면 시간영역에서 $i(t)$를 구할 수 있다. 이 해에는 지수함수 감쇠, 공진 등에 의한 진동 특성이 모두 반영되며, 특정 주파수 구간에서의 위상 천이와 진폭 변화를 체계적으로 이해할 수 있다. 실제 회로의 과도응답이나 정상상태 해석은 이처럼 라플라스 변환을 사용하면 미분방정식을 풀 때 생기는 불편한 초기조건 처리를 훨씬 쉽게 진행할 수 있게 된다.

#### 회로 해석의 다이어그램

단순 RLC 직렬회로를 도식으로 나타낼 때는 다음과 같은 mermaid 예시로도 표현할 수 있다

{% @mermaid/diagram content="flowchart LR
A\["전압원 v(t)"] --> B\[R] --> C\[L] --> D\[C] --> E\[접지]
D -->|병렬| E" %}

여기서 A에서 공급되는 전압 $v(t)$가 R, L, C를 차례로 거쳐 흐른 뒤 접지에 도달한다. 이 회로에 대한 미분방정식을 라플라스 변환으로 변환해 해를 구하는 기본 절차는 앞에서 제시한 것과 같다.

#### RLC 병렬 회로의 해석

전압원 $v(t)$가 병렬로 구성된 R, L, C를 동시에 연결하는 경우, 회로방정식은 전류의 합이 성분별 전류를 모두 합한 것과 같다는 형태가 된다. 병렬 회로는 각 소자에 걸리는 전압이 동일하고 소자 전류가 서로 다르게 흐른다는 점이 직렬 회로와 다르다. 이때 모든 소자 양단 전압을 $v(t)$, 전체 전류를 $i(t)$, 저항을 R, 인덕턴스를 L, 커패시턴스를 C라고 하면

$$
i(t) = i\_R(t) + i\_L(t) + i\_C(t)
$$

가 된다. 여기서

$$
i\_R(t) = \frac{v(t)}{R}
$$

$$
i\_L(t) = \frac{1}{L} \int\_0^t v(\tau),d\tau \quad \text{(초기 전류가 없다고 가정할 때)}
$$

$$
i\_C(t) = C \frac{dv}{dt} \quad \text{(초기 전하가 없다고 가정할 때)}
$$

등의 식이 성립한다. 직접 적분과 미분 형태로 나타나므로, 이를 시간영역에서 바로 푸는 것은 까다로울 수 있다. 라플라스 변환을 사용하면 이들의 합으로 구성된 전체 전류 $I(s)$에 대한 표현을 얻을 수 있다.

라플라스 변환에서

$$
\mathcal{L} { i\_R(t) } = \frac{V(s)}{R}
$$

$$
\mathcal{L} { i\_L(t) } = \frac{V(s)}{s L}
$$

$$
\mathcal{L} { i\_C(t) } = C \left\[ s V(s) - v(0) \right]
$$

가 되므로, 초기 전압 $v(0)$이 0이라면

$$
I(s) = \frac{V(s)}{R} + \frac{V(s)}{sL} + C , s V(s)
$$

이 되고, 이를 정리하면

$$
I(s) = V(s) \left( \frac{1}{R} + \frac{1}{sL} + sC \right).
$$

주어진 회로가 전압원 $v(t)$에 의해 구동되고 있다면, $i(t)$와 $v(t)$의 관계는 전류-전압 비의 라플라스 표현을 통해 손쉽게 구한다. 예컨대 $v(t) = V\_0 \cos(\omega t)$인 경우

$$
V(s) = \frac{V\_0 , s}{s^2 + \omega^2}
$$

이므로

$$
I(s) = \frac{V\_0 , s}{s^2 + \omega^2} \left( \frac{1}{R} + \frac{1}{sL} + sC \right).
$$

역라플라스 변환을 수행하면 시간영역에서의 $i(t)$가 유도되며, 그 해는 지수 응답(과도성분)과 정상상태에서의 정현 응답이 결합된 형태로 표현된다. 이 계의 공진 주파수는 분모의 극점으로부터 산출할 수 있으며, 병렬 공진 시점에서는 유도성 임피던스와 커패시턴스 임피던스가 상쇄되어 회로가 저항의 역할만 하거나, 특정 조건에서는 반대로 급격한 임피던스 증가 현상이 발생한다.

#### 라플라스 변환을 통한 고차 미분방정식 해석

계의 차수가 높아질수록, 예컨대 2차 미분이 아닌 3차, 4차 미분으로 이루어진 운동방정식이나 회로방정식을 다룰 때 라플라스 변환은 더욱 유용한 도구가 된다. 예를 들어, 2자유도 진동계에서 질량이 각각 $m\_1, m\_2$이고, 스프링 상수가 각각 $k\_1, k\_2, k\_3$이며, 감쇠 요소가 $c\_1, c\_2, c\_3$인 계를 생각하자. 질량 $m\_1$의 변위를 $x\_1(t)$, 질량 $m\_2$의 변위를 $x\_2(t)$라 하면 운동방정식은

$$
m\_1 \frac{d^2 x\_1}{dt^2} + (c\_1 + c\_2) \frac{d x\_1}{dt} - c\_2 \frac{d x\_2}{dt} + (k\_1 + k\_2) x\_1 - k\_2 x\_2 = F\_1(t)
$$

$$
m\_2 \frac{d^2 x\_2}{dt^2} - c\_2 \frac{d x\_1}{dt} + (c\_2 + c\_3) \frac{d x\_2}{dt} - k\_2 x\_1 + (k\_2 + k\_3) x\_2 = F\_2(t)
$$

형태가 된다. 이를 라플라스 변환하면

$$
m\_1 \left\[s^2 X\_1(s) - s x\_1(0) - \dot{x}\_1(0)\right]  + (c\_1 + c\_2)\left\[s X\_1(s) - x\_1(0)\right]  - c\_2\left\[s X\_2(s) - x\_2(0)\right]  + (k\_1 + k\_2) X\_1(s)  - k\_2 X\_2(s) = F\_1(s)
$$

$$
m\_2 \left\[s^2 X\_2(s) - s x\_2(0) - \dot{x}\_2(0)\right] - c\_2 \left\[s X\_1(s) - x\_1(0)\right] + (c\_2 + c\_3)\left\[s X\_2(s) - x\_2(0)\right] - k\_2 X\_1(s)  + (k\_2 + k\_3) X\_2(s) = F\_2(s)
$$

와 같이 표현된다. 이를 정리하면

$$
\begin{pmatrix} m\_1 s^2 + (c\_1 + c\_2) s + (k\_1 + k\_2) & -c\_2 s - k\_2 \ - c\_2 s - k\_2 & m\_2 s^2 + (c\_2 + c\_3) s + (k\_2 + k\_3) \end{pmatrix} \begin{pmatrix} X\_1(s) \ X\_2(s) \end{pmatrix} = \begin{pmatrix} \tilde{G}\_1 \ \tilde{G}\_2 \end{pmatrix}
$$

꼴이 되며, 우변의 $\tilde{G}\_1, \tilde{G}\_2$는 $F\_1(s)$, $F\_2(s)$ 및 초기조건 항들이 모두 포함된 형태다. 이 연립방정식을 풀면 $X\_1(s)$와 $X\_2(s)$가 주어진다. 다시 역라플라스 변환을 취하면 $x\_1(t), x\_2(t)$가 완전한 형태로 유도된다.

이처럼 라플라스 변환을 사용하면 고차 미분항을 갖는 다자유도계의 해석도 표준화된 연립 선형대수 형태로 바꿀 수 있고, 각 변수의 시간응답을 전개하기에 간결하다. 더 나아가 해가 갖는 극점(폴) 구조를 통해 계의 고유진동수 및 감쇠비를 일관성 있게 해석할 수도 있다.

#### 실제 회로의 예: OP 앰프를 포함한 라플라스 해석

연산증폭기(OP 앰프)를 포함하는 회로를 라플라스 영역에서 해석하면, 가상의 무한대 게인(ideal OP amp) 가정에 따라 입력단자 전압이 같아야 한다는 조건이나 입력단 전류가 0이라는 조건 등을 쉽게 반영할 수 있다. 예를 들어, 간단한 미분회로(ideal differentiator)는 OP 앰프의 반전 입력 단자에 커패시터 C를 직렬 연결하고, 피드백 저항 R을 출력에서 반전 입력으로 연결한 구조다. 입력 전압을 $v\_\mathrm{in}(t)$, 출력 전압을 $v\_\mathrm{out}(t)$라고 하면, 이상 OP 앰프 가정 하에 라플라스 영역에서

$$
V\_\mathrm{out}(s) = - \frac{1}{R C} \left\[ s V\_\mathrm{in}(s) - v\_C(0) \right]
$$

와 같이 나타날 수 있다. 여기서 $v\_C(0)$는 커패시터의 초기 전압이다. 초기조건이 0이고, $v\_\mathrm{in}(t)$가 0 이전에 존재하지 않으면

$$
V\_\mathrm{out}(s) = - \frac{s}{R C} V\_\mathrm{in}(s).
$$

시간영역에서의 이상적인 표현은

$$
v\_\mathrm{out}(t) = - \frac{1}{R C} \frac{dv\_\mathrm{in}}{dt}
$$

이지만, 실제 회로에서는 이상적인 OP 앰프 특성이나 무시할 수 없는 기생요소, 잡음 등이 존재하여 고주파 대역에서 정확히 미분동작을 하기가 어렵다. 그럼에도 라플라스 변환을 통해 이론적으로 설계 개념을 잡는 것이 효과적이다. 필요하면 보상회로를 추가하여 고주파 응답을 제한하거나, 적정한 주파수 스펙에서 동작하도록 할 수 있다. 이런 전체 동작 구조를 라플라스 영역에서 분석하면, 단순 미분 또는 적분 기능뿐 아니라 저역통과, 고역통과, 대역통과, 대역저지 등 다양한 2차, 3차 이상의 복합 필터 설계에도 적용 가능하다.

#### 블록선도로 표현되는 시스템 해석

연산증폭기나 여러 회로 요소가 상호작용하는 복잡한 시스템의 경우, 블록선도를 그리고 그 블록에 해당하는 라플라스 전달함수를 부여하여 해석할 수 있다. 예를 들어

{% @mermaid/diagram content="flowchart LR
A\["v\_in(t)"] -->|입력| B\["블록 G1(s)"]
B --> C\["블록 G2(s)"]
C --> D\["v\_out(t)"]
C -->|피드백| B" %}

이런 구조는

$$
v\_\mathrm{out}(s) = \frac{G\_1(s) G\_2(s)}{1 + G\_1(s) G\_2(s)} , v\_\mathrm{in}(s)
$$

와 같은 폐루프 전달함수를 가진다. 라플라스 변환은 이렇게 시스템 전반의 동특성을 간결하게 나타내고, 고루프 이득(limit of large open-loop gain)에서의 동작, 공진 주파수, 위상 여유 등을 음미하기에 효과적이다.

#### 2포트 네트워크의 라플라스 해석

여러 전기·전자 시스템을 서브 블록 단위로 묶어 해석할 때, 2포트 회로망(two-port network)은 특히 중요한 역할을 한다. 입력단과 출력단이 각각 두 개씩의 단자를 공유하는 구조에서, 이를 라플라스 영역에서 간단히 표현할 수 있다. 가령 한 쌍의 단자를 입력, 다른 쌍의 단자를 출력으로 삼는 상황을 생각하면, 2포트 파라미터로 임피던스 형태(Z-파라미터), 어드미턴스 형태(Y-파라미터), 전달파라미터(ABCD-파라미터), 혼합파라미터(h-파라미터) 등을 사용할 수 있다.

예를 들어 Z-파라미터로 표현하면

$$
\begin{pmatrix}  V\_1(s) \  V\_2(s) \end{pmatrix} = \begin{pmatrix} Z\_{11}(s) & Z\_{12}(s) \ Z\_{21}(s) & Z\_{22}(s) \end{pmatrix} \begin{pmatrix} I\_1(s) \ I\_2(s) \end{pmatrix}
$$

꼴이 되며, 각 $Z\_{ij}(s)$는 주파수(또는 s)에 따라 달라지는 일반적인 임피던스를 나타낸다. 라플라스 영역에서 2포트 네트워크가 직렬, 병렬, 캐스케이드, 브리지 등으로 결합될 때, 대응하는 네트워크 파라미터를 적절히 합성하여 전체 시스템 파라미터를 구할 수 있다.

2포트 해석에서 가장 흔히 쓰이는 것은 ABCD(전달 파라미터) 형태인데,

$$
\begin{pmatrix} V\_1(s) \  I\_1(s) \end{pmatrix} = \begin{pmatrix} A(s) & B(s) \ C(s) & D(s) \end{pmatrix} \begin{pmatrix} V\_2(s) \ - I\_2(s) \end{pmatrix}
$$

라고 정의한다. 여기서 $A(s), B(s), C(s), D(s)$가 각각 2포트의 특성을 라플라스 영역에서 나타내며, 직렬 연결(캐스케이드) 시에는 행렬 곱셈으로 간단히 전체 전달 파라미터를 구할 수 있다. 필터, 증폭기, 전송회로 등을 이런 방식으로 블록 모듈화하면, 복잡한 시스템을 단계적으로 설계하고 해석하기가 한결 용이해진다.

#### 텔레그래퍼 방정식을 통한 전송선로 해석

고주파 영역에서 전송선로 효과가 무시될 수 없는 경우, 전송선로에 대한 해석이 중요해진다. 전송선로는 일반적으로 분포정수(직렬 인덕턴스와 저항, 병렬 커패시턴스와 컨덕턴스)가 선로를 따라 연속적으로 존재하는 구조로서, 이를 미분방정식 형태로 나타내면 텔레그래퍼 방정식을 얻게 된다. 길이 좌표를 x, 시간은 t, 선로 위의 전압과 전류를 각각 $v(x,t)$, $i(x,t)$라고 하면, 단위 길이당 매개변수를 L(인덕턴스), R(저항), C(커패시턴스), G(컨덕턴스)로 두었을 때

$$
\frac{\partial v(x,t)}{\partial x} = -L \frac{\partial i(x,t)}{\partial t} - R , i(x,t)
$$

$$
\frac{\partial i(x,t)}{\partial x} = -C \frac{\partial v(x,t)}{\partial t} - G , v(x,t)
$$

가 텔레그래퍼 방정식이다. 이를 라플라스 변환(시간 t에 대하여)하면, s-영역에서

$$
\frac{\partial V(x,s)}{\partial x} = - \bigl\[ s L + R \bigr] I(x,s)∂I(x,s)
$$

$$
\frac{\partial I(x,s)}{\partial x} = - \bigl\[ s C + G \bigr] V(x,s)
$$

형태가 되고, x에 대한 1차 미분방정식을 결합하면

$$
\frac{\partial^2 V(x,s)}{\partial x^2} = \gamma^2(s) , V(x,s),
\quad
\gamma(s) = \sqrt{\bigl\[s L + R\bigr]\bigl\[s C + G\bigr]}
$$

와 같은 형태로 묶여진다. 초기 및 경계조건(예: 선로 길이 끝단에 연결된 임피던스 등)을 적절히 설정하여 이 선로 문제를 풀면, 반사와 전파, 감쇠 등에 대한 상세한 해석을 할 수 있다. 이때 라플라스 해석은 과도응답을 포함한 시간영역 거동을 구하는 데 편리하다. 포인트-매칭이나 모드 전개 방식으로 역라플라스 변환을 시도하기도 하며, 단순화된 경우(무손실 전송선로, R=G=0)는 보다 간단한 폐형해를 얻을 수 있다.

병렬 또는 직렬로 가지를 이루는 복합 전송선로 구조나, 여러개의 분포정수를 결합한 파동가이드형 구조가 나타나면, 적절한 경계조건을 매 구간마다 부여하고 구역해를 붙여나가는 식으로 라플라스 영역에서 전체 해를 구성한다. 이렇게 구한 해는 고주파 안테나 급전선 설계나 고속 디지털 회로 기판(PCB) 배선 해석에도 폭넓게 응용된다.

#### 기계-전기 아날로지를 통한 라플라스 변환 적용

물리학적 예제를 확장해보면, 기계계(질량-스프링-댐퍼)와 전기계(RLC)의 상호 대응관계를 정립하여, 한쪽의 문제를 다른 쪽의 형태로 치환해 푸는 기법이 자주 사용된다. 예를 들어 변위-속도-가속도와 전하-전류-전압의 대응, 스프링 상수와 커패시턴스, 질량과 인덕턴스, 감쇠와 저항의 대응 등이다. 이를 기계-전기 아날로지라고 하며, 라플라스 변환을 적용할 경우 이런 대응관계를 더욱 일관성 있게 파악할 수 있다.

예를 들어 질량 m, 감쇠계수 c, 스프링상수 k가 있는 단자에서, 힘 F(t)와 변위 x(t)에 관한 방정식

$$
m \frac{d^2 x}{dt^2} + c \frac{dx}{dt} + k x = F(t)
$$

이 있다면, 이를 전기회로에 대응시키면 인덕턴스 L, 저항 R, 커패시턴스 C에 의해

$$
L \frac{d i}{dt} + R i + \frac{1}{C} \int i(t),dt = v(t)
$$

같은 꼴의 방정식에 맵핑된다. 여기서 변위 x(t)는 커패시터 전압 v(t), 힘 F(t)는 전류원 혹은 전압원 형태로 해석되며, 미분연산과 적분연산의 대응이 라플라스 변환에서 간단히 정립된다. 이러한 유비적 해석을 통해 전자·기계·유체 등 다양한 물리계를 통합적인 시각으로 분석할 수 있다는 장점이 있다.

#### 비선형 진동계의 라플라스 해석 접근

실제 물리계나 공학적 시스템은 종종 스프링이나 회로 소자가 선형 범위에서 벗어나, $k x$ 관계 대신 $k x + \alpha x^3$ 같은 비선형 항을 포함하기도 한다. 예를 들어 더핑(Duffing) 방정식은

$$
m \frac{d^2 x}{dt^2} + c \frac{dx}{dt} + k x + \alpha x^3 = F(t)
$$

형태를 가지며, 스프링의 강성이 선형 항 $k x$와 비선형 항 $\alpha x^3$가 결합된 구조다. 순수 선형계의 경우와 달리, 이러한 비선형성이 존재하면 직접 라플라스 변환을 적용하여 폐형해(Closed-form solution)를 구하기가 쉽지 않다.

그럼에도 불구하고, 비선형 항이 충분히 작아 선형계 해를 기준으로 근사할 수 있는 상황이라면, 라플라스 변환을 응용한 반복법(perturbation)을 사용할 수 있다. 예를 들어

$$
m \frac{d^2 x}{dt^2} + c \frac{dx}{dt} + k x = F(t) - \alpha x^3
$$

로 재정리하여

$$
X(s) = \frac{m \left( s x(0) + \dot{x}(0) \right) + c, x(0) + F(s) - \mathcal{L}{\alpha x^3}}{m s^2 + c s + k}
$$

형태를 얻은 뒤, $x^3(t)$를 이전 단계에서 구한 해의 거동으로 근사화해가면서 반복하는 식이다. 이 근사법은 비선형성이 크지 않은 강제진동계에서 준선형(jump phenomenon, subharmonic 응답 등) 해를 찾는 데 유용할 수 있다.

심하게 비선형인 경우, 고전적인 라플라스 변환만으로 해를 구하기보다는 수치적 시간영역 해석과 함께, 부분적으로 라플라스 변환을 활용하여 과도응답이나 초기조건 영향을 분석하는 방법이 선택된다. 예를 들어 모수여기나 카오스 영역에서의 진동해를 다룰 때는 시스템의 위상공간(dynamical system 관점) 해석이 더 적절할 수 있다. 하지만 낮은 차수의 연산 증폭기 회로나 스프링 계에서 발생하는 약한 비선형 효과는, 위에서 제시한 근사 라플라스 기법을 통해 일정 수준까지는 간단히 다룰 수 있다.

#### 비선형 전기회로 해석 (Chua 회로 등)

비선형 소자를 포함하는 전기회로도 마찬가지로 라플라스 변환을 직접 적용하기 어렵다. 대표적인 예가 Chua 회로로 알려진 혼돈(Chaos) 현상 발현 회로다. Chua 다이오드라고 불리는 특수한 비선형 소자를 사용하면, 다음과 같은 형태의 3차 미분방정식이 나온다.

단순화된 Chua 회로를 표현하면, 커패시터 두 개와 인덕터, 그리고 비선형 소자(Chua 다이오드) 및 저항 등이 포함된다. 각 축 전압 및 전류를 $v\_1(t), v\_2(t), i\_L(t)$로 정의할 때, 미분방정식 계는

$$
C\_1 \frac{dv\_1}{dt} = G \left( v\_2 - v\_1 \right) - f\bigl( v\_1 \bigr)
\\
C\_2 \frac{dv\_2}{dt} = G \left( v\_1 - v\_2 \right)
\\
L \frac{d i\_L}{dt} = -v\_2
$$

이 된다. 여기서 $f(v\_1)$가 비선형 특성을 나타내는 항이다. 가령

$$
f(v\_1) = m\_0 v\_1 + \tfrac{1}{2} (m\_1 - m\_0) \left\[ \lvert v\_1 + B\_p \rvert - \lvert v\_1 - B\_p \rvert \right]
$$

같은 구간별 선형 함수로 정의되기도 한다. 이를 라플라스 변환해보면, $f(v\_1)$ 항 자체가 선형이 아니므로 직접적인 선형 해법이 막히게 된다. $\mathcal{L} { f(v\_1(t)) } = F\bigl(V\_1(s)\bigr)$ 식으로 단순화되지 않기 때문이다.

따라서 이런 비선형 전기회로를 해석할 때는, 라플라스 변환을 통한 해석을 부분적으로 이용하되 전체 해는 수치적(예: Runge-Kutta 방법)으로 풀어야 한다. 순수 정상상태 해석이나 소신호선형화(small-signal linearization) 관점에서는 $f(v\_1)$을 특정 평형점 근처에서 선형화하여 라플라스 표현을 얻는 방법도 쓸 수 있다. 즉 평형점 $v\_{1e}, v\_{2e}, i\_{Le}$ 근방에서 테일러 전개하면

$$
f(v\_1) \approx f(v\_{1e}) + f'(v\_{1e}) \left( v\_1 - v\_{1e} \right)
$$

가 되어 선형화된 전달함수를 구할 수 있다. 이 과정에서 라플라스 변환이 다시 빛을 발하며, 계의 안정성(특특값, 극점 위치), 공진특성, 임계파라미터 등을 얻을 수 있다.

#### 선형 시스템과의 접합부 해석

복합시스템에서 한 부분은 선형(라플라스 변환 적용 용이), 다른 부분은 비선형(수치적 해석 필요)일 때, 이를 하이브리드 기법으로 연결하여 해결하기도 한다. 예를 들어 입력에서 일부 필터링을 담당하는 RLC 블록은 라플라스 해석을 통해 정확한 시간응답 혹은 전달함수를 얻고, 출력단에서 나타나는 비선형 소자(예: 다이오드, 트랜지스터, 포화증폭기 등)에 대해서는 시간영역 미분방정식을 수치적으로 풀어 병합한다. 간단한 예로 OP 앰프가 선형 구간에선 일정 이득을 보이다가 입력이 클 때 포화(saturation)에 진입하는 비선형을 보이는 경우, 포화 전후의 동작을 분리하여 라플라스 해석 부분과 수치분해 부분을 연결해 최종 해를 구할 수 있다.

시뮬레이션 툴(예: SPICE)들은 본질적으로 KCL, KVL 방정식을 시간영역에서 수치적 방법으로 풀면서도, 선형 서브회로는 행렬화하여 빠르게 풀어낸다. 라플라스 변환의 이론적 토대는 이러한 시뮬레이터 구조에도 중요한 개념적 근간이 된다.

#### FFT 및 z-변환과의 관계

시스템 응답을 주파수 영역에서 해석할 때, 라플라스 변환(연속-시간)과 푸리에 변환(연속-시간) 사이의 관계는

$$
\mathcal{F}{ x(t)} = X(j \omega) = X(s) \big|\_{s = j \omega}
$$

이런 형태로 이어진다. 여기서 $X(s)$는 라플라스 변환, $X(j\omega)$는 푸리에 변환 또는 주파수응답이라고도 한다. 실제 디지털 신호처리에서는 z-변환이 사용되는데, 이는 시간을 이산화(discretization)했을 때 라플라스 변환과 유사한 역할을 수행한다. z-변환의 양변 환산

$$
z = e^{sT}
$$

($T$는 샘플링 주기)을 통해, 디지털 필터 설계를 라플라스 영역에서의 연속 필터 설계와 대비시키는 기법이 잘 알려져 있다.

이처럼 라플라스 변환은 연속-시간 영역에서 물리계 전반을 다루는 핵심 도구일 뿐 아니라, 이산화와 주파수 영역 해석으로 자연스럽게 확장되어, 실제 공학 시스템의 아날로그-디지털 접합부를 분석하는 데도 유용하다.

#### 편미분방정식에서의 라플라스 변환

진동이나 전기회로 문제는 대부분 상미분방정식(ODE) 형태로 나타나지만, 더 높은 차원의 물리 현상은 편미분방정식(PDE)으로 표현되기도 한다. 예를 들어 2차원 혹은 3차원 계에서 파동, 열, 전자기적 현상을 기술할 때, 시간 t와 공간 좌표 x, y, z에 대한 종속변수로 구성된 PDE가 등장한다. 라플라스 변환은 시간 변수에 대해서만 변환을 취하므로, 공간 좌표는 그대로 보존된다. 이를 통해 복잡한 PDE를 s-영역에서 공간 좌표에 대한 편미분방정식(또는 상미분방정식) 형태로 변환할 수 있다.

대표적으로 열전달(Heat Equation) 문제를 고려할 수 있다. 1차원 열 방정식

$$
\frac{\partial u(x,t)}{\partial t} = \alpha \frac{\partial^2 u(x,t)}{\partial x^2}
$$

에서 초기조건 $u(x,0) = \phi(x)$와 경계조건(예: $u(0,t) = 0, u(L,t) = 0$)이 주어졌다고 하자. 시간 t에 대해 라플라스 변환을 취하면

$$
s U(x,s) - \phi(x) = \alpha \frac{\partial^2 U(x,s)}{\partial x^2}
$$

가 된다. 여기서 $U(x,s)$는 $u(x,t)$의 라플라스 변환이다. 그 결과, x에 대한 2차 상미분방정식

$$
\frac{\partial^2 U(x,s)}{\partial x^2} - \frac{s}{\alpha} U(x,s) = - \frac{\phi(x)}{\alpha}
$$

가 나타난다. 경계조건은 라플라스 영역에서도 동일한 x=0, x=L에서 $U(0,s)=0, U(L,s)=0$ 형태로 적용된다. 이 선형의 2차 미분방정식을 풀면 $U(x,s)$를 구할 수 있고, 최종적으로 시간영역 해 $u(x,t)$는 역라플라스 변환으로 얻게 된다.

이는 라플라스 변환이 ‘시간에 대한 적분연산’을 대수적으로 단순화해주기 때문에 가능한 접근이다. 열 방정식뿐 아니라, 물리학이나 공학에서 등장하는 전자기파 전파(맥스웰 방정식 계열을 단순화한 형태), 음향학(음파 방정식), 탄성파 해석(탄성체 파동 방정식) 등에서도 유사하게 시간 변수만 라플라스 변환을 취하고, 공간방정식을 풀어내는 절차가 흔히 사용된다.

#### 예시: 1차원 파동방정식의 과도응답

1차원 파동방정식

$$
\frac{\partial^2 u(x,t)}{\partial t^2} = c^2 \frac{\partial^2 u(x,t)}{\partial x^2}
$$

가 있다면, 이 계에도 라플라스 변환을 적용할 수 있다. 초기조건 $u(x,0) = \phi(x), \frac{\partial u}{\partial t}(x,0) = \psi(x)$가 주어졌을 때, 시간 t에 대하여 변환하면

$$
s^2 U(x,s) - s ,\phi(x) - \psi(x) = c^2 \frac{\partial^2 U(x,s)}{\partial x^2}.
$$

이는 공간좌표 x에 대한 2차 상미분방정식

$$
\frac{\partial^2 U(x,s)}{\partial x^2} - \frac{s^2}{c^2} U(x,s) = -\frac{s}{c^2} , \phi(x) - \frac{1}{c^2},\psi(x)
$$

로 정리된다. 경계조건(예: 유한한 막의 양 끝에서의 고정단, 자유단 등)을 라플라스 영역에서 반영하여 $U(x,s)$를 구한 뒤 역라플라스 변환을 통해 $u(x,t)$를 얻는다.

실제로는 파동방정식의 그린 함수나 d’Alembert 해법이 더 직접적일 수도 있으나, 라플라스 변환을 사용하면 과도응답(Initial transient response)을 보다 체계적으로 추적하기 쉽다는 장점이 있다. 경계가 복잡하거나 외력이 주어졌을 때(예: 분포하중, 시간에 따라 변하는 외부 경계입력 등), 라플라스 도메인에서 조건을 반영하고, 부분구간 해를 이어 붙여 전 영역 해석을 수행한다. 이 방법은 구간별 정의된 PDE를 풀어야 하는 복잡 구조물(다층 매질, 분절된 막)에서도 응용할 수 있다.

#### 다자유도 진동계의 PDE로의 확장

단순히 물리량이 공간적으로 연속이 아닌 개별 질점과 연결 스프링으로 구성된 고차원 진동계도 충분히 복잡하지만, 실제 구조물(예: 빔, 판, 쉘)처럼 공간에 연속적으로 분포된 탄성계는 PDE로 표현된다. 예를 들어 오일러-베르누이 빔 방정식

$$
\rho A \frac{\partial^2 w(x,t)}{\partial t^2} + E I \frac{\partial^4 w(x,t)}{\partial x^4} = 0
$$

에서 $w(x,t)$는 빔의 처짐(displacement), $\rho$는 밀도, A는 단면적, E는 탄성계수, I는 단면 2차 모멘트 등을 나타낸다. 진동해석 시 초기조건과 경계조건(캔틸레버, 양단 고정, 양단 힌지 등)을 부여하고 라플라스 변환을 t에 대해 취하면

$$
s^2 \rho A , W(x,s) + E I \frac{\partial^4 W(x,s)}{\partial x^4} = - \rho A \left\[ s w(x,0) + \frac{\partial w}{\partial t}(x,0) \right]
$$

같은 상미분방정식을 x에 대해 푸는 형태가 된다. 여기서 $W(x,s)$는 $w(x,t)$의 라플라스 변환이다. 경계조건은 빔 단부에서의 회전, 변위, 모멘트, 전단력을 어떻게 처리할지에 따라 라플라스 영역에서도 대응되는 조건을 세울 수 있다. 이를 해석해 얻은 $W(x,s)$에 역라플라스 변환을 취해 시간영역으로 복원하면, 빔의 모든 위치 x에서 시간 t에 따른 처짐 $w(x,t)$를 구할 수 있다. 이 방법은 필연적으로 복잡한 부분분수 전개나 특이함수(베셀함수, 힌지함수 등)를 사용할 수도 있으나, 특정 단순 경계나 주파수 대역에 집중할 경우 유용하다.

#### 변분법, 모드해석과의 결합

편미분방정식을 라플라스로 접근할 때, 직접 미분방정식을 푸는 대신 고유치 문제(eigenvalue problem)로 변환하고 모드 전개를 취하는 방식도 널리 사용된다. 예를 들어 빔이나 막 방정식에서, 공간모드 $\phi\_n(x)$들을 미리 구해두고, 변위를

$$
w(x,t) = \sum\_n q\_n(t) , \phi\_n(x)
$$

로 전개하면, 각 모드별로 상미분방정식(시간 t에 대한)이 등장한다. 이 상미분방정식에 라플라스 변환을 적용해 $Q\_n(s)$를 구한 뒤, 다시

$$
W(x,s) = \sum\_n Q\_n(s) , \phi\_n(x)
$$

로 표현한다. 이 과정을 거치면 PDE가 모드별 ODE들의 집합으로 환원되므로, 해석이 훨씬 단순화된다. 다만 모드 수가 많거나 공진특성이 복합적인 대형 구조에서는, 이 모드합도 수치화가 필수적이 된다.

모드해석은 전자기적 문제(가우스 방정식, 맥스웰 방정식)에서도 유사하게 나타난다. 예를 들어 도파관 파동해석에서 모드함수들을 찾고, 시간종속형식에 라플라스 변환을 취하거나, 고정 주파수에서의 정상파 해석을 시간영역 해석으로 확장할 때 적용이 가능하다. 실제로는 s가 복소 주파수로 확장된 해석이므로, 소실계수(attenuation constant)나 위상천이(phase shift) 등이 자연스럽게 반영된다.

#### 복합 물리계에서의 적분방정식 접근

라플라스 변환은 PDE뿐 아니라 적분방정식이나 적분-미분 방정식 형태에도 확장 적용 가능하다. 예컨대 역학계의 지연요소나 열 복사(radiative heat transfer)에 의한 복합 연산이 추가되면, 시간 영역에서의 지연 커널(convolution integral)이 등장하게 된다. 이를

$$
u(t) = \int\_0^t K(\tau) x(t-\tau), d\tau
$$

같은 형태로 쓰면, 라플라스 영역에서 곱셈 형태 $U(s) = K(s) X(s)$로 간단히 표현된다. PDE에서 시간 경과에 따라 누적되는 어떤 비선형 손실이나, 기억작용(memory effect)이 있는 재료를 다룰 때도 마찬가지다. 재료가 선형 점탄성(viscoelasticity)일 경우 응력과 변형률의 관계가 적분 연산으로 주어지지만, 라플라스 변환을 취하면 모듈러스가 $E(s)$ 꼴로 표현되어 상미분방정식 수준으로 바뀐다.

이처럼 복합 물리계에서 라플라스 변환은 선형 영역의 지연효과, 메모리 효과, 경계조건 처리를 모두 깔끔하게 정리해주는 도구로 널리 응용된다. 실제 공학 설계에서는 이를 전용 소프트웨어(유한요소 해석, 유한차분 해석 등)로 구현할 때 내부적으로 라플라스 변환이 아닌 다른 알고리즘을 쓸 수도 있다. 하지만 이론적 기초를 학습하고, 간단한 모델에서 문제의 핵심 역학을 파악하는 데에는 여전히 라플라스 변환이 유용하다.

#### 특수함수 입력(단위계단, 델타함수)에 대한 물리계 응답

라플라스 변환은 단위계단 함수나 델타 함수 같은 특수 입력이 주어졌을 때, 계의 과도응답을 분석하기에 매우 유리하다. 실제 물리학 혹은 공학에서는 “갑자기 스위치를 켠 전압” (단위계단 입력)이나 “매우 짧은 충격력을 가한 상황” (델타 입력)을 쉽게 모델링할 수 있다는 점에서 이러한 특수함수의 활용도가 높다.

단위계단 함수 $u(t)$(Heaviside step function)의 라플라스 변환은

$$
\mathcal{L}{ u(t) } = \frac{1}{s}
$$

이며, 델타 함수 $\delta(t)$(Dirac delta)의 라플라스 변환은

$$
\mathcal{L}{\delta(t)} = 1
$$

이다. 예를 들어 전기회로에서 스위치를 켰을 때, 전압원이 $v(t) = V\_0 u(t)$로 주어진다고 하면 $V(s) = \frac{V\_0}{s}$가 된다. 이 전압을 직렬 RLC 회로에 가하면,

$$
I(s) = \frac{V(s)}{L s + R + \frac{1}{s C}} = \frac{V\_0/s}{L s + R + \frac{1}{s C}}
$$

가 되므로, 이후 부분분수 전개와 역라플라스 변환을 통해 시정수 및 고유진동수, 감쇠 등 과도해를 명확히 파악할 수 있다.

델타 함수 입력은 순수한 충격 자극으로서, 물리계가 초기조건을 갖도록 만드는 역할을 한다. 예를 들어 기계계에서 $\delta(t)$만큼의 외력을 인가하면 초기 순간에만 힘이 작용하므로, 이후에는 자유진동(자유응답)만 남게 된다. 따라서 $F(t) = F\_0 \delta(t)$의 라플라스 변환은 $F\_0$가 되고, 2차 진동계

$$
m \frac{d^2 x}{dt^2} + c \frac{dx}{dt} + k x = F\_0 \delta(t)
$$

에 대해 라플라스 변환을 하면

$$
m (s^2 X(s) - s x\_0 - \dot{x}(0)) + c (s X(s) - x\_0) + k X(s) = F\_0.
$$

초기 변위, 속도가 0이었다고 가정하면

$$
(m s^2 + c s + k) X(s) = F\_0,
\\
X(s) = \frac{F\_0}{m s^2 + c s + k}.
$$

이 식의 역변환을 통해 $x(t)$가 실제로 어떤 자유진동 모드를 갖는지, 충격 후에 어떻게 감쇠해가는지를 바로 알 수 있다.

#### 초기조건으로서의 충격 및 스위치 동작

실제 물리계나 회로는 초기 변위, 초기 전압, 초기 전류 등 “초기조건”을 가질 수 있다. 이를 $\delta(t)$나 $u(t)$ 형태의 특수함수 입력으로 치환하여 해석하기도 한다. 예컨대 RC 회로에서 커패시터에 초기전압 $V\_0$가 걸려 있다면, t=0 직전에 $\delta(t)$ 유형의 순간 충전이 발생했다고 볼 수도 있다. 이런 식으로 초기조건을 특수신호 입력 형태로 등가 치환하면, 회로방정식을 “무(無) 초기조건” + “특수신호 입력” 형태로 재해석하는 것이 가능하다.

라플라스 변환은 초기조건 처리에 편리한 이유가, 모든 미분항에 대해 해당 초기값이 명시적으로 등장하기 때문이다. 예를 들어 2차 미분에 대한 라플라스 변환 $\mathcal{L}{\ddot{x}(t)} = s^2 X(s) - s x(0) - \dot{x}(0)$ 꼴로, $x(0), \dot{x}(0)$이 자동으로 분리된다. 그 결과, 강제신호(외부입력)만 변환해도, 자연스럽게 초기조건이 라플라스 영역 식에 포함되어 해석이 일관되게 진행된다.

#### 과도응답과 정상상태응답 분리

선형시스템에서 입력이 정현파(또는 복소지수)인 경우, 장시간이 흐른 뒤에는 결국 정상상태응답(steady-state response)이 남고, 초기에 나타나는 과도응답(transient response)은 사라진다. 라플라스 변환으로 해를 전개하면, 보통

$$
X(s) = \underbrace{\frac{\text{초기조건 항}}{m s^2 + c s + k}}*{\text{과도응답 성분}}  + \underbrace{\frac{\text{외력항}}{m s^2 + c s + k}}*{\text{강제응답 성분}}
$$

같이 분할된 구조를 바로 확인할 수 있다. 여기서 과도응답에 해당하는 항들은 보통 $s$ 평면에서 음의 실부분을 갖는 극점을 지니므로, 시간영역에서는 $e^{-\alpha t}$ 형태의 소멸항(감쇠항)으로 나타난다. 한편 강제응답은 안정계라면 영구진동(정현응답) 또는 일정값(DC 응답) 형태로 남는다.

특히 라플라스 영역에서 $s=j\omega$로 대입했을 때, 전달함수의 주파수응답을 구해볼 수 있다. 이는 비단 진동계 분석뿐 아니라 회로, 제어 시스템에서도 중요한 해석 기법이다. 시스템이 특정 주파수에서 공진(레조넌스)을 일으키는지, 얼마나 큰 진폭 증폭을 보이는지를 손쉽게 파악할 수 있다. 과도영역에서는 폴 위치가 (Re(s)<0)인지 (Re(s)=0)인지 등에 따라 감쇠, 발산, 한계주기운동 등이 결정된다.

#### 라플라스 영역에서의 폴과 제로 해석

전달함수

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

가 주어졌을 때, $D(s)=0$을 만족하는 해들을 폴(pole)이라 하고 $N(s)=0$인 해들을 제로(zero)라 한다. 물리학적 시스템에서는 폴이 바로 자연진동 고유값, 혹은 자유응답 모드에 대응한다. 예를 들어 2차 시스템

$$
G(s) = \frac{\omega\_0^2}{s^2 + 2\zeta \omega\_0 , s + \omega\_0^2}
$$

에서 분모 $s^2 + 2\zeta \omega\_0 , s + \omega\_0^2 =0$을 풀면

$$
s = - \zeta \omega\_0 \pm \omega\_0 \sqrt{\zeta^2 - 1}
$$

가 폴이 된다. 감쇠비 $\zeta<1$이면 (복소수) 음의 실부분을 갖는 복소폴이며, $\zeta=1$이면 중근(임계감쇠), $\zeta>1$이면 서로 다른 음의 실근(과감쇠)이다. 제로가 포함된 시스템이라면, 그 위치에 따라 동특성이 달라진다. 예를 들어 어떤 회로에서는 제로가 오른쪽 반평면(Re(s)>0)에 존재하면 초과위상(lead effect)이나 불안정 모드가 관여하기도 한다.

물리계에서 폴과 제로의 배치는 가장 기본적인 동특성 및 주파수응답 특성을 결정짓는다. 라플라스 변환은 이를 간결한 다항식 형태로 추출하므로, 복잡한 미분연산을 ‘유리함수로의 단순화’로 연결해준다.

#### 시스템 안정성 판별

미분방정식으로 표현되는 계가 주어진 초기조건이나 입력에 대해 해가 발산하지 않고 유계(bounded)로 유지되는지 판별하는 것은 공학에서 매우 중요하다. 라플라스 영역에서의 시스템 안정성은 분모 다항식(특성방정식)의 폴 위치로 판별할 수 있다. 모든 폴이 왼쪽 반평면(Re(s)<0)에 존재하면 계가 안정(Asymptotically stable), 오른쪽 반평면에 하나라도 폴이 있으면 불안정, 허수축(j\omega 축)에 있으면 경계 안정(Marginally stable) 또는 준안정으로 구분한다.

기계 진동계에서 $\zeta>0$이 보장되면 폴이 왼편으로 오므로 감쇠성 진동이 되고, $c=0$이라면 $\zeta=0$이 되어 폴이 허수축에 놓인 순수진동이 이루어진다. 전자회로에서도 RLC 회로의 저항이 충분하면 왼편 폴로 안정동작하나, 증폭기 회로처럼 능동소자가 있으면 양의 실부분 폴이 생길 수 있어 발진을 일으키기도 한다.

라플라스 관점에서 안정성 판별은 Routh-Hurwitz 판별법이나 Nyquist 판별법 등을 통해 알기 쉽게 구현된다. 예를 들어 2차 방정식 $s^2 + 2\zeta \omega\_0 s + \omega\_0^2=0$에서 $\zeta>0$이면 안정, $\zeta=0$이면 경계 안정, $\zeta<0$이면 불안정인 것을 바로 알 수 있다.

#### 제어 시스템 설계에서의 활용

제어 이론에서 라플라스 변환은 시스템 모델링 및 폐루프 전달함수 구성에 기반이 된다. 오차율, 감쇠비, 오버슈트, 상승시간 등의 성능 지표를 분석할 때도 라플라스 영역에서의 극점, 영점, 폐루프 전달함수 형태가 핵심 지표로 작용한다. 단순 PD, PI, PID 제어부터, 고차 보상기(lead-lag, 상태피드백 등)에 이르기까지, 선형 제어 설계는 대부분 라플라스 기반으로 수행된다.

예를 들어 PID 제어기는

$$
G\_c(s) = K\_p + \frac{K\_i}{s} + K\_d s
$$

같은 꼴로 표현된다. 프로세스(plant) $G\_p(s)$와 직렬 연결(병렬이 아니라 신호적 측면에서 합성된 전송함수로서)하면, 폐루프 전달함수는

$$
T(s) = \frac{G\_p(s) G\_c(s)}{1 + G\_p(s) G\_c(s)}.
$$

여기서 $T(s)$의 분모를 특성방정식이라 부르고, 이 극점들의 분포를 조절하여 원하는 동특성(빠른 응답, 작은 오버슈트, 제한된 스테디스테이트 오차 등)을 구현한다. 모두 라플라스 변환의 ‘유리함수 표현’을 통해 얻을 수 있는 구조다.

물리학적 관점에서 본다면, 기계 진동계를 능동 제어하는 경우(예: 액추에이터로 힘을 가해 진동을 저감)나, 전자회로에서 주파수 응답 특성을 보상(Compensation)하여 발진 혹은 과도 진동을 줄이고자 할 때, 이런 라플라스 해석과 제어 설계가 자연스럽게 결합된다.

#### 추가적 응용: 네트워크 합성, 필터 설계

회로이론에서 임의의 양단 임피던스나 전달함수를 만족하도록 R, L, C 소자를 조합해내는 문제를 네트워크 합성(network synthesis)이라고 한다. 이때도 라플라스 변환이나 푸리에 변환이 필수적이다. 전달함수를

$$
H(s) = \frac{V\_\mathrm{out}(s)}{V\_\mathrm{in}(s)}
$$

로 정의했을 때, 여기에 필요한 영점/극점의 배치를 RLC 소자들이 어떻게 실현할 수 있는지(부분분수 전개, Foster-형, Cauer-형 등) 체계적으로 해석한다. 무리없이 수동소자로만 구성 가능한가, 능동소자가 필요한가, 특정 주파수대역에서 얼마나 높은 차수의 roll-off(감쇠)를 구현해야 하는가 등의 문제를 다룰 수 있다. 이는 고성능 필터, 주파수 분할 장치, 임피던스 정합 회로 등에 자주 적용된다.

마찬가지로, 전기음향 분야나 디지털 오디오 시스템에서 요구되는 등화기(equalizer)와 크로스오버 네트워크 또한 라플라스(또는 z-변환) 해석 하에 설계된다. 임의의 크로스오버 주파수를 설정하고, 각 대역에서 위상왜곡을 최소화하도록 2차, 3차 이상의 필터를 직렬-병렬 조합하여 구현하기 위해선, 전부 전달함수 설계가 선행되어야 한다.

#### 디지털 구현을 위한 변환

아날로그 방식으로 설명된 물리계나 회로모델을 디지털 제어기로 구현하려면, 샘플링과 양자화 과정을 거쳐야 한다. 이때 라플라스 변환 모델을 적절히 이산화(discretization)하는 과정에서

$$
z = e^{s T}
$$

같은 대응관계가 활용된다. T는 샘플링 주기이며,

$$
G(z) = G\bigl(e^{sT}\bigr)
$$

로 치환해 디지털 필터나 디지털 제어 알고리즘을 설계하는 것이다. 이 과정을 Tustin의 방법(양자화 전쌍), Zero-Order Hold 등 다양한 방식으로 구체화할 수 있으며, 실제 설계 시엔 라플라스 변환을 기반으로 얻은 아날로그 전달함수를 적절히 변환해오거나, 반대로 z-변환에서 출발해 라플라스 영역 해석을 대조하기도 한다.

최종적으로는 아날로그-디지털 혼합환경(예: 센서로부터 측정한 연속신호를 A/D로 변환해 DSP나 마이크로컨트롤러에서 연산하고, D/A로 변환해 액추에이터를 구동)에서 라플라스 변환과 z-변환을 자유롭게 넘나들어야 한다. 이 또한 물리계(연속시간, 미분방정식)와 디지털 처리(이산시간, 차분방정식)를 연결하는 핵심 이론적 고리가 되어 준다.

#### 라플라스 변환과 분산(distributed) 시스템

공학 현장에서 하나의 기기나 회로만 독립적으로 있는 경우도 있지만, 대규모 네트워크 또는 여러 장치가 상호작용하는 분산(distributed) 환경이 등장하기도 한다. 가령 전력망, 통신망, 대규모 메카트로닉스 시스템 등에서 수많은 서브시스템을 연결해 전체 시스템의 동작을 분석해야 하는 경우가 있다. 이때 각각의 서브시스템이 일정 부분 선형성을 가진다면, 라플라스 변환을 통해 서브시스템별 전달함수를 구한 뒤 네트워크 형태로 합성할 수 있다. 전력계통에서 선로 임피던스, 발전기 모델, 부하 모델 등을 각각 s-영역 표현으로 얻고, 이를 연결하여 전체 계통의 과도안정도나 단락사고 시 응답 등을 해석하는 식이다.

통신망 분석에서는 선형 필터나 증폭기, 변조기 등을 라플라스(또는 푸리에) 도메인에서 연결해 링크의 응답, 신호대잡음비(SNR), 왜곡 등을 계산한다. 고주파 대역(마이크로파)에서 전송선로 모형이 중요해질 때도, 앞서 말한 텔레그래퍼 방정식을 라플라스 도메인으로 변환해 각 구간별 특성을 블록으로 취급하면 망 전체 해석이 가능하다.

#### 다중물리(Multiphysics) 시스템에서의 라플라스 변환

현대 공학은 기계, 전기, 열, 유체, 화학 등의 여러 물리 현상을 동시에 고려해야 하는 복합시스템을 다룬다. 예를 들어 전기모터를 열 방출까지 고려하며 해석하거나, 유체의 흐름에 의해 기계구조물 진동이 발생하고 그로 인해 전기신호에 피드백이 생기는 등 복합적 상호작용이 무시할 수 없을 정도로 커지는 사례가 늘어난다. 이럴 때도 라플라스 변환 접근은 개별적인 물리 도메인에서의 미분방정식을 블록화하여 연결하기 용이하다.

예를 들어 열 방정식을 통해 온도 분포가 시간에 따라 변화하는 것을 파악한 다음, 그 온도가 전기회로나 기계 부품의 파라미터(R, L, C, k, c 등)에 변화를 일으키는 식으로 연결할 수 있다. 필요하면 이 변화가 선형 근사로 기술될 수 있는 구간에서만 라플라스 변환을 적용하고, 큰 온도 변화 시에는 그 구간을 재설정하여 해석을 반복하기도 한다. 비선형, 시간변화 파라미터가 존재해도 구간별 선형화(piecewise linearization)를 하면 라플라스 해석을 단계적으로 이어붙일 수 있다.

#### 그린 함수와 라플라스 변환

편미분방정식이나 복잡한 경계조건 문제를 해석할 때, 그린 함수(Green’s function)를 구해 점근해(積分核)로 활용하는 방법이 있다. 어떤 선형 연산자 L에 대해

$$
L {G(x,\tau)} = \delta(x-\tau)
$$

형태의 근본해를 구하면, 임의의 외력(혹은 외부자극) 함수가 주어졌을 때 그 해를 적분 연산으로 나타낼 수 있다. 시간과 공간을 모두 포함한 문제에서도, 시간에 대해서만 라플라스 변환을 취하고 공간 변수에 대해 그린 함수를 구하면, 한 번 계산한 그린 함수를 가지고 여러 종류의 외력 입력에 대한 응답을 간단히 적분 형태로 표현할 수 있다.

대표적으로, 1차원 열방정식에서의 그린 함수는 보통 가우시안 꼴(열핵 함수)로 나타나며, 이에 대한 라플라스 변환 표현은 $e^{- \sqrt{s/\alpha},(x-\tau)}$ 꼴로 이어질 수 있다. 파동방정식, 탄성체 방정식 등도 마찬가지로 그린 함수를 구한 뒤, 라플라스 변환을 적용하면 특정 경계조건에서의 과도파 전달이나 구조물 내부 응력 분포를 체계적으로 계산할 수 있다.

#### 라플라스 변환의 수치적 기법

실제 복잡한 시스템에서 역라플라스 변환의 폐형해를 얻기는 어렵거나 불가능한 경우가 많다. 이때는 직접 수치적 알고리즘을 이용해 $X(s)$로부터 $x(t)$를 복원한다. 대표적인 기법으로는 Bromwich 적분을 근사하는 방법, Talbot’s 방법, Durbin’s 방법 등이 있다. Bromwich 적분은

$$
x(t) = \frac{1}{2\pi j} \int\_{\gamma - j\infty}^{\gamma + j\infty} e^{s t} X(s), ds
$$

형태인데, 이를 직접 적분하기엔 어려움이 있으므로 적절한 수치积分 기법과 경로 변형(Deformation of contour)을 결합한다. Talbot’s 방법은 복소평면에서의 적분경로를 곡선 형태로 잡고 수치 근사를 수행한다.

수치 라플라스 역변환은 ODE 솔버와는 달리, 전 구간 시간응답을 한 번의 적분(혹은 근사합)으로 얻을 수 있다는 이점이 있다. 그러나 $X(s)$가 복잡하면 복소평면 상의 극점 구조, 특이점 등이 섬세하게 반영되어야 하므로 계산이 만만치 않을 수 있다. 특히 Re(s) < 0 방향으로 깊이 들어가 있는 폴이나 브랜치 컷(Branch cut)이 존재하면, 경로 배치와 부분분수 전개가 까다로워진다.

또 다른 방법으로는, 시간영역에서 적분-미분 방정식을 수치해석하는 (예: Runge-Kutta) 방식과, 라플라스 역변환을 수치로 시도하는 방식을 혼합하여, 주파수 응답 혹은 특정 t에 대한 응답값만 계산하는 식이다. 이 과정에서 FFT 계열 알고리즘을 사용하여, 라플라스 영역에서 jω축 부근 정보를 집중적으로 다루고, 필요한 경우 s 평면의 실제부(Re(s))에 대한 이동(감쇠 인자)을 포함하도록 설계하기도 한다.

#### 고차, 고차원 시스템에서 라플라스 해석의 한계와 대안

계의 차수와 차원이 매우 커지면(수백, 수천 차수 이상의 대규모 시스템), 직접 전형적인 라플라스 변환이나 부분분수 전개로 해를 기술하는 것은 사실상 불가능해진다. 이런 대규모 문제는 보통 행렬 미분방정식 형태로 정리하고, 수치선형대수를 이용해 해석한다. 시스템 매트릭스 \mathbf{A}, 입력 매트릭스 \mathbf{B}, 출력 매트릭스 \mathbf{C} 등을 정의하면,

$$
\dot{\mathbf{x}}(t) = \mathbf{A},\mathbf{x}(t) + \mathbf{B},\mathbf{u}(t)
$$

$$
\mathbf{y}(t) = \mathbf{C},\mathbf{x}(t) + \mathbf{D},\mathbf{u}(t)
$$

같은 상태방정식이 되고, 라플라스 변환에서

$$
s \mathbf{X}(s) - \mathbf{x}(0) = \mathbf{A},\mathbf{X}(s) + \mathbf{B},\mathbf{U}(s)
$$

로 표현되므로

$$
\mathbf{X}(s) = (s \mathbf{I} - \mathbf{A})^{-1} \bigl\[ \mathbf{x}(0) + \mathbf{B},\mathbf{U}(s)\bigr].
$$

이때 $(s \mathbf{I} - \mathbf{A})^{-1}$를 전개하면, 고차 다항식의 특성방정식에 해당하는 행렬식(det(s \mathbf{I}-\mathbf{A}))이 나와 시스템 극점들이 결정된다. 극점 수(계의 차수)가 매우 많아지면, 부분분수 전개가 복잡해지고 해가 상징적으로 간단히 표현되지 않는다. 그래서 수치선형대수를 통해 고유치-고유벡터(eigenvalue-eigenvector)로 분해하거나, 모드 축소(Modal truncation), 균형 축소(Balanced truncation) 등으로 차수를 줄이는 방법을 쓴다.

결국, 라플라스 변환 자체는 이론적으로는 모든 선형계 해석에 통용되는 통일된 언어지만, 대규모 현실 문제에선 수치선형대수 기법과 결합하여 부분적으로만 사용하거나, 주파수응답 해석에 집중하는 방식으로 사용하게 된다. 그럼에도 불구하고 계통의 극점, 영점을 파악하고 안정성, 공진 특성 등을 분석하는 근본 틀은 여전히 라플라스 변환 언어에 토대를 두고 있다고 할 수 있다.