REMI Lab강의 노트

Signals and Systems L03

3장 연속시간 시스템과 컨볼루션

임펄스 응답으로 본 LTI 시스템

임펄스 응답 하나만 보고 LTI 시스템의 출력과 성질을 모두 알 수 있을까?

LTI 시스템은 임펄스 응답 하나로 정해진다. 입력에서 출력까지 가는 흐름은 네 단계이다.

  1. 입력을 가중 임펄스로 쪼갠다: x(t)=∫x(τ) δ(t−τ) dτx(t)=\int x(\tau)\,\delta(t-\tau)\,d\tau.
  2. 시스템에 넣는다: y(t)=H{∫x(τ) δ(t−τ) dτ}y(t)=\mathcal H\left\{\int x(\tau)\,\delta(t-\tau)\,d\tau\right\}.
  3. LTI이면 시스템과 적분의 순서를 바꿀 수 있다: y(t)=∫x(τ) H{δ(t−τ)} dτy(t)=\int x(\tau)\,\mathcal H\{\delta(t-\tau)\}\,d\tau.
  4. 따라서 출력은 입력과 임펄스 응답의 컨볼루션이다: y=x∗hy=x*h.

LTI 시스템의 성질도 hh에서 바로 알 수 있다. t<0t<0에서 h(t)=0h(t)=0이면 컨볼루션에 현재와 과거의 입력만 들어가므로 인과이다(예제 2의 hh는 인과가 아니다). h(t)=c δ(t)h(t)=c\,\delta(t) 꼴이면 y(t)=c x(t)y(t)=c\,x(t)로 현재 입력만 쓰므로 무기억이다. 입력이 ∣x(t)∣≤M<∞|x(t)|\le M<\infty이면

∣y(t)∣≤∫−∞∞∣h(τ)∣ ∣x(t−τ)∣ dτ≤M∫−∞∞∣h(τ)∣ dτ|y(t)|\le\int_{-\infty}^{\infty}|h(\tau)|\,|x(t-\tau)|\,d\tau \le M\int_{-\infty}^{\infty}|h(\tau)|\,d\tau

이므로 ∫∣h(τ)∣ dτ\int|h(\tau)|\,d\tau가 유한하면 출력에도 시간과 무관한 유한한 상한이 있어 BIBO 안정이다.

고에너지 광자 검출기가 LTI 시스템의 예이다. 광자 하나가 들어오면 검출기는 빠르게 켜졌다가 지수 함수처럼 줄어드는 빛을 낸다. 이것이 검출기의 임펄스 응답이다. 광자가 시각 t1,t2,…t_1, t_2, \dots에 들어오면 입력은 임펄스의 합 ∑kδ(t−tk)\sum_k\delta(t-t_k)이고, 출력은 같은 응답을 도착 시각만큼 옮겨 더한 ∑kh(t−tk)\sum_k h(t-t_k)이다. 광자가 드물면 응답이 하나씩 따로 보이고, 잦으면 앞 응답이 다 줄어들기 전에 다음 응답이 얹혀 쌓인다. 검출기가 포화하지 않을 만큼 광자가 적을 때 이 LTI 설명이 맞는다.

광자가 3개 들어올 때와 9개 들어올 때를 비교하면, 광자(화살표)마다 같은 모양의 응답이 하나씩 생긴다. 9개일 때는 점선으로 그린 응답들이 겹친 곳에서 실선 출력이 그 합으로 쌓여 올라간다.

광자가 도착할 때마다 검출기 출력에 같은 모양의 지수 감쇠 응답이 하나씩 생기고, 출력은 그 응답들의 합이다. 광자가 드물면 응답이 따로 보이고, 잦으면 겹쳐 쌓인다.광자 3개t광자 9개t
회색 화살표: 광자 도착(입력 임펄스). 점선: 광자 하나의 응답. 실선: 검출기 출력(응답의 합). 감쇠 시정수와 도착 시각은 설명용 값이다.

RC 회로도 같은 생각으로 설명한다. 임펄스 응답이 h(t)=1τe−t/τu(t)h(t)=\frac1\tau e^{-t/\tau}u(t)인 시스템이 있다(여기서 τ=RC>0\tau=RC>0은 적분변수가 아니라 시정수이다). 단위 계단을 넣은 출력은 t≥0t\ge0에서 s(t)=∫0t1τe−v/τ dv=1−e−t/τs(t)=\int_0^t\frac1\tau e^{-v/\tau}\,dv=1-e^{-t/\tau}이고, t<0t<0에서는 0이다(적분변수는 vv). 이 시스템에 펄스를 넣으면, 펄스를 두 계단의 차 u(t)−u(t−T)u(t)-u(t-T)로 보고 각 계단의 응답을 빼서 출력을 얻는다. 뒤의 계단이 늦게 들어온 만큼 그 응답도 늦게 나오는 것이 시불변성이고, 부호를 바꿔 더하는 것이 선형성이다.

이 계산을 세 줄로 나누어 그리면, 위 줄은 입력 펄스와 그것을 이루는 두 계단(점선), 가운데 줄은 계단 응답 s(t)s(t)와 임펄스 응답 h(t)h(t), 아래 줄은 두 계단의 응답(점선)과 그 합인 출력이다. s(t)s(t)의 값은 h(t)h(t)를 0부터 tt까지 적분한 면적(보라색)과 같고, 출력 y=s(t)−s(t−T)y=s(t)-s(t-T)는 두 계단 응답의 높이(점선 위 두 점)를 뺀 값이다.

처음 상태는 1 V, 폭 T = 20 ms인 펄스를 τ = RC = 20 ms인 시스템에 넣은 것이고, 시스템의 초기 상태는 0이다. 계단 응답 그림의 주황 점은 t − T에서의 값 s(t − T)이고, 출력의 합에서 빼는 값이 이것이다(처음 상태 0.393).

해 보기

  1. 입력이 계단 하나이면 출력이 어떻게 되는지 확인하라.
    관찰

    늦은 계단이 없어 출력이 계단 응답 s(t)=1−e−t/τs(t)=1-e^{-t/\tau} 하나이고, 1 V로 다가간다.

  2. 펄스 폭이 80 ms이면 펄스가 끝나는 순간의 출력이 얼마인지 확인하라.
    관찰

    처음 계단의 응답이 거의 1 V까지 오른 뒤에 늦은 계단의 응답이 빠지기 시작한다. 펄스가 끝나는 순간 y=1−e−4≈0.98y=1-e^{-4}\approx0.98 V이다.

  3. 펄스 폭이 5 ms이면 출력이 어디까지 오르고 그 뒤 어떻게 줄어드는지 확인하라.
    관찰

    출력이 1−e−0.25≈0.221-e^{-0.25}\approx0.22 V까지만 오르고, 그 뒤로는 h(t)h(t)와 같은 시정수 20 ms로 지수처럼 줄어든다.

인과이고 BIBO 안정인 LTI 시스템의 임펄스 응답을 모두 고르라.

  1. (가) h(t)=e−2∣t∣h(t)=e^{-2|t|}
  2. (나) h(t)=u(t)h(t)=u(t)
  3. (다) h(t)=δ(t−1)h(t)=\delta(t-1)
  4. (라) h(t)=e−2tu(t)h(t)=e^{-2t}u(t)
풀이

답 (다) h(t)=δ(t−1)h(t)=\delta(t-1), (라) h(t)=e−2tu(t)h(t)=e^{-2t}u(t)

  • (가) ∫∣h∣ dt=1\int|h|\,dt=1이라 안정이지만 t<0t<0에서도 값이 있어 인과가 아니다.
  • (나) 인과이지만 ∫0∞1 dt\int_0^\infty1\,dt가 무한해서 안정이 아니다. 계단을 넣으면 출력이 램프로 끝없이 커진다.

인과는 t<0t<0에서 h(t)=0h(t)=0인지로, BIBO 안정은 ∫∣h(t)∣ dt\int|h(t)|\,dt가 유한한지로 판정한다. e−2tu(t)e^{-2t}u(t)는 면적이 12\tfrac12이고 δ(t−1)\delta(t-1)은 면적이 1이며, 둘 다 t<0t<0에서 0이다. δ(t−1)\delta(t-1)은 1초 지연기이다.

1 V, 폭 T=20T=20 ms인 펄스를 τ=20\tau=20 ms인 시스템에 넣었다. 펄스가 끝나는 순간의 출력 y(T)y(T)를 구하라.

풀이

답 0.632 V

y(t)=s(t)−s(t−T)y(t)=s(t)-s(t-T)에서 t=Tt=T이면 s(0)=0s(0)=0이라 y(T)=s(T)=1−e−T/τ=1−e−1≈0.632y(T)=s(T)=1-e^{-T/\tau}=1-e^{-1}\approx0.632 V이다.