REMI Lab강의 노트

Signals and Systems L02

2장 신호의 표현

기본 신호

복잡한 신호를 식으로 쓸 때 부품으로 쓰는 신호는 무엇일까?

단위 계단

단위 계단(unit step function) u(t)u(t)는 t=0t=0에서 0에서 1로 올라서는 신호이다. 단위 계단 함수 또는 Heaviside 함수라고도 한다.

u(t)={1,t≥00,otherwiseu(t)=\begin{cases}1, & t\ge0\\ 0, & \text{otherwise}\end{cases}

t=0t=0에서의 값은 교재마다 다르다. 1로 두기도 하고, ½로 두기도 하고, 정의하지 않기도 한다. 이 책에서는 u(0)=1u(0)=1로 둔다. 계단은 스위치를 t=0t=0에 켜는 순간을 나타내고, 다른 신호에 곱하면 t<0t<0 부분을 지운다. 예를 들어 e−tu(t)e^{-t}u(t)는 t≥0t\ge0에서만 e−te^{-t}이고 그 전에는 0이다.

사각 펄스

사각 펄스(rectangular function) rect(t)\mathrm{rect}(t)는 폭 1, 높이 1인 펄스이다. 사각 함수 또는 펄스 신호라고도 부른다.

rect(t)={1,−12≤t<120,otherwise\mathrm{rect}(t)=\begin{cases}1, & -\tfrac12\le t<\tfrac12\\ 0, & \text{otherwise}\end{cases}
단위 계단 u(t)tu(t)⋯−111사각 펄스 rect(t)trect(t)−½½1
채운 점이 그 시각의 값이다. u(0) = 1이고, rect(t)는 t = −½를 포함하고 t = ½는 포함하지 않는다.

폭을 TT로 바꾸려면 시간 스케일링을 쓴다. rect(t/T)\mathrm{rect}(t/T)는 −T/2≤t<T/2-T/2\le t<T/2에서 1이다. 면적은 rect(t)\mathrm{rect}(t)가 1, rect(t/T)\mathrm{rect}(t/T)가 TT이다.

디랙 델타

디랙 델타 함수(Dirac delta function) δ(t)\delta(t)는 단위 임펄스(unit impulse)라고도 부르며, 줄여서 임펄스라고 한다. 다음 두 성질로 정의한다.

δ(t)=0(t≠0),∫−∞∞δ(t) dt=1\delta(t)=0\quad(t\neq0),\qquad\int_{-\infty}^{\infty}\delta(t)\,dt=1

t=0t=0 한 점을 빼면 0인데 면적은 1이어야 하므로 δ(0)\delta(0)은 유한한 값일 수 없다(δ(0)→∞\delta(0)\to\infty). 그래서 δ(t)\delta(t)는 보통의 함수가 아니다. 한 점에서만 0이 아닌 보통의 함수는 적분하면 0이기 때문이다. δ(t)\delta(t)는 적분 안에서 무엇을 하는지로 정의되는 대상이고, 수학에서는 일반화 함수(generalized function) 또는 분포(distribution)라고 부른다. 그림에서는 화살표로 그리고, 화살표 옆에 적은 수는 높이가 아니라 면적이다.

다음 문장이 참인지 거짓인지 고르라. "t=0t=0에 그린 화살표 옆에 2라고 적힌 2δ(t)2\delta(t)는 t=0t=0에서 값이 2인 신호이다."

  1. 참
  2. 거짓
풀이

답 거짓

  • 참 δ(t)\delta(t)는 t=0t=0에서 유한한 값을 갖지 않는다. 화살표 옆의 수가 무엇을 나타내는지 다시 본다.

2δ(t)2\delta(t)는 t≠0t\neq0에서 0이고, 적분하면 2인 임펄스이다. 화살표 옆의 2는 높이가 아니라 면적이다. t=0t=0에서의 값은 유한한 수로 정해지지 않는다.

펄스의 극한으로 본 δ(t)

δ(t)\delta(t)는 폭이 좁아지는 펄스로 생각하면 머릿속에 그리기 쉽다.

gε(t)={1∣ε∣,∣t∣<∣ε∣20,otherwiseg_\varepsilon(t)=\begin{cases}\dfrac{1}{|\varepsilon|}, & |t|<\dfrac{|\varepsilon|}{2}\\[6pt] 0, & \text{otherwise}\end{cases}

gεg_\varepsilon는 폭 ∣ε∣|\varepsilon|, 높이 1/∣ε∣1/|\varepsilon|인 사각 펄스라서 ε\varepsilon가 얼마이든 면적이 1이다. ε\varepsilon를 줄이면 폭은 0으로, 높이는 무한대로 가고, 면적은 1로 남는다. 이 극한이 δ(t)\delta(t)이다.

δ(t)=lim⁡ε→0gε(t)\delta(t)=\lim_{\varepsilon\to0}g_\varepsilon(t)
폭 ε, 높이 1/ε인 펄스와 ε → 0ε = 2, 1, 1/2, 1/4의 펄스를 겹쳐 그림. 폭이 반으로 줄 때마다 높이는 두 배가 되어 면적은 1로 같다.tgε(t)ε = 2, 높이 ½ε = 1, 높이 1ε = ½, 높이 2ε = ¼, 높이 4−1−0.50.51
ε = 2, 1, ½, ¼인 gε(t)를 겹쳐 그렸다. 폭이 반으로 줄 때마다 높이가 두 배가 되어 면적은 모두 1이고, ε를 ⅛, 1/16로 줄이면 높이는 8, 16이 된다.

이 극한은 적분 안에서 뜻을 갖는다. t=0t=0에서 연속인 신호 x(t)x(t)에 gεg_\varepsilon를 곱해 적분하면

∫−∞∞x(t) gε(t) dt=1∣ε∣∫−∣ε∣/2∣ε∣/2x(t) dt\int_{-\infty}^{\infty}x(t)\,g_\varepsilon(t)\,dt=\frac{1}{|\varepsilon|}\int_{-|\varepsilon|/2}^{|\varepsilon|/2}x(t)\,dt

이고, 이것은 폭 ∣ε∣|\varepsilon|인 구간에서 xx의 평균이다. 구간이 0으로 줄면 평균은 x(0)x(0)으로 간다. 즉 δ(t)\delta(t)를 곱해 적분하면 x(0)x(0)만 남는다. 다음 절의 sifting 성질이 이것이다.

ε=0.1\varepsilon=0.1일 때 gε(t)g_\varepsilon(t)의 높이와 면적을 고르라.

  1. (가) 높이 10, 면적 10
  2. (나) 높이 10, 면적 1
  3. (다) 높이 0.1, 면적 0.01
  4. (라) 높이 1, 면적 0.1
풀이

답 (나) 높이 10, 면적 1

  • (가) 면적은 폭 × 높이 =0.1×10=0.1\times10이다.
  • (다) 폭과 높이를 바꿔 읽었다. 폭이 ∣ε∣=0.1|\varepsilon|=0.1이고 높이는 1/∣ε∣1/|\varepsilon|이다.
  • (라) 높이를 1로 두면 면적이 폭만큼 줄어든다. gεg_\varepsilon는 높이를 1/∣ε∣1/|\varepsilon|로 키워 면적을 1로 지킨다.

폭은 ∣ε∣=0.1|\varepsilon|=0.1, 높이는 1/∣ε∣=101/|\varepsilon|=10이고 면적은 0.1×10=10.1\times10=1이다. ε\varepsilon를 더 줄여도 면적은 1로 남는다.

δ와의 곱과 sifting 성질

t0t_0에서 연속인 x(t)x(t)에 δ(t−t0)\delta(t-t_0)를 곱하면, δ\delta가 0이 아닌 곳은 t=t0t=t_0 한 점뿐이다. 따라서 x(t)x(t)를 그 점의 값 x(t0)x(t_0)로 바꿔 써도 된다.

x(t) δ(t−t0)=x(t0) δ(t−t0)x(t)\,\delta(t-t_0)=x(t_0)\,\delta(t-t_0)

δ 자리에 앞 절의 펄스 gεg_\varepsilon를 t0t_0로 옮겨 놓으면 이 등식이 극한으로 보인다. gε(t−t0)g_\varepsilon(t-t_0)는 t0t_0를 가운데 둔 폭 ε\varepsilon, 높이 1/ε1/\varepsilon의 펄스라서, 곱 x(t) gε(t−t0)x(t)\,g_\varepsilon(t-t_0)의 면적은 t0t_0 둘레 폭 ε\varepsilon 구간에서 xx의 평균이다. 그림의 처음 상태 t0=0t_0=0, ε=2\varepsilon=2에서는 ①의 색칠한 윈도우 안의 평균(가로 선분)이 x(t0)x(t_0) 점보다 낮다. ε\varepsilon를 반씩 줄이면 ②의 곱은 좁고 높아지고, ③의 면적은 x(t0)x(t_0)로 다가간다.

x(t) = 1 + 0.6 cos 2t + 0.3 sin t이고 x(0) = 1.6이다. 면적은 ε = 2에서 1.273, ε = 1에서 1.505, ε = ½에서 1.575이다.

해 보기

  1. ε=14\varepsilon=\tfrac14이면 곱의 면적이 x(0)x(0)에 얼마나 가까운지 확인하라.
    관찰

    곱은 폭 0.25, 높이 약 6.5로 솟고 면적은 1.594이다. x(0)=1.6x(0)=1.6과의 차이는 0.006으로, ε=2\varepsilon=2일 때의 0.327보다 50배 넘게 작다.

  2. 윈도우가 골짜기 t0=1.5t_0=1.5에 있을 때 윈도우 안의 평균을 x(1.5)x(1.5)와 비교하라.
    관찰

    윈도우 안의 평균 0.982가 x(1.5)=0.705x(1.5)=0.705보다 높다. 골짜기에서는 윈도우 양쪽의 값이 가운데보다 커서 평균이 위로 벗어난다. 피크 t0=0t_0=0에서 평균이 x(t0)x(t_0)보다 낮았던 것과 반대이다.

  3. t0=1.5t_0=1.5에서 ε\varepsilon를 ⅛까지 줄여 가며 면적이 어느 쪽에서 x(t0)x(t_0)에 다가가는지 확인하라.
    관찰

    면적이 0.982, 0.787, 0.727, 0.711을 거쳐 0.707로 내려온다. 위에서 다가가지만 극한은 역시 x(t0)=0.705x(t_0)=0.705이다.

ε\varepsilon를 반으로 줄일 때마다 면적과 x(0)x(0)의 차이는 약 ¼로 준다. 윈도우가 좁아지면 윈도우 안의 xx가 직선에 가까워지고, 직선의 평균은 가운데 값과 같기 때문이다. 곱은 ε\varepsilon가 작아질수록 높이 x(t0)/εx(t_0)/\varepsilon인 사각형 x(t0) gε(t−t0)x(t_0)\,g_\varepsilon(t-t_0)에 가까워지고, 극한에서 위 식의 x(t0) δ(t−t0)x(t_0)\,\delta(t-t_0)가 된다.

곱은 면적이 x(t0)x(t_0)인 임펄스이다. 이것을 적분하면 δ\delta의 면적이 1이므로 다음을 얻는다.

∫−∞∞x(t) δ(t−t0) dt=x(t0)\int_{-\infty}^{\infty}x(t)\,\delta(t-t_0)\,dt=x(t_0)

신호 전체에서 t0t_0의 값 하나만 걸러 낸다는 뜻에서 sifting(체로 거르기) 성질이라고 부른다. 이름이 비슷한 shifting(이동)과는 다르다. 이 식은 x(t)x(t)가 t0t_0에서 연속일 때 성립한다.

가중 임펄스의 합

sifting 식에서 t0t_0 자리에 tt를, 적분변수 자리에 τ\tau를 넣고 δ\delta가 짝함수임(δ(τ−t)=δ(t−τ)\delta(\tau-t)=\delta(t-\tau))을 쓰면 다음 식이 된다.

x(t)=∫−∞∞x(τ) δ(t−τ) dτx(t)=\int_{-\infty}^{\infty}x(\tau)\,\delta(t-\tau)\,d\tau

sifting 식과 같은 식이지만 뜻은 반대 방향이다. sifting 식은 "임펄스로 한 점의 값을 뽑는다"이고, 이 식은 "신호 x(t)x(t)는 시각 τ\tau마다 놓인 임펄스 δ(t−τ)\delta(t-\tau)에 그 시각의 값 x(τ)x(\tau)를 가중치로 곱해 모두 더한 것"이다. 3장에서 시스템의 출력을 컨볼루션으로 구할 때 바로 이 분해에서 출발한다.

적분 안의 δ\delta에 배율이 붙어 있으면 먼저 풀어야 한다. u=atu=at로 치환할 때 a>0a>0이면 dt=du/adt=du/a이고 적분 한계가 그대로라 ∫−∞∞x(t) δ(at) dt=x(0)/a\int_{-\infty}^{\infty}x(t)\,\delta(at)\,dt=x(0)/a이다. a<0a<0이면 uu의 적분 한계가 ∞\infty에서 −∞-\infty로 바뀌므로, 이를 되돌리는 음의 부호가 더해져 적분값은 −x(0)/a=x(0)/∣a∣-x(0)/a=x(0)/|a|이다. 따라서 a≠0a\neq0인 두 경우 모두

δ(at)=1∣a∣ δ(t)\delta(at)=\frac{1}{|a|}\,\delta(t)

이다. 면적이 1/∣a∣1/|a|인 임펄스라는 뜻이다.

∫−∞∞cos⁡t δ(−2t) dt\displaystyle\int_{-\infty}^{\infty}\cos t\,\delta(-2t)\,dt를 구하라.

풀이

답 0.5

힌트 δ(at)=1∣a∣δ(t)\delta(at)=\frac{1}{|a|}\delta(t)에서 a=−2a=-2이다. 절댓값에 주의한다.

δ(−2t)=1∣−2∣δ(t)=12δ(t)\delta(-2t)=\frac{1}{|-2|}\delta(t)=\tfrac12\delta(t)이므로 적분은 12cos⁡0=0.5\tfrac12\cos0=0.5이다. 절댓값 없이 1a\frac{1}{a}를 쓰면 −0.5-0.5가 나오는데, 면적이 음수인 임펄스가 되어 틀린다.

직접 계산해 보기. 다음 적분값을 구하라.

(1) ∫−21(t+t2) δ(t−3) dt\displaystyle\int_{-2}^{1}(t+t^2)\,\delta(t-3)\,dt

풀이

0이다. 임펄스는 t=3t=3에 있는데 적분 구간 [−2,1][-2,1] 밖이다. 구간 안에서 δ(t−3)\delta(t-3)은 어디서나 0이다.

적분 구간과 임펄스 위치tδ(t − 3)[−2, 1][−2, 4]−20134
임펄스가 적분 구간 안에 있을 때만 적분값이 생긴다. t = 3은 [−2, 1] 밖, [−2, 4] 안에 있다.

(2) ∫−24(t+t2) δ(t−3) dt\displaystyle\int_{-2}^{4}(t+t^2)\,\delta(t-3)\,dt

풀이

12이다. t=3t=3이 구간 [−2,4][-2,4] 안에 있으므로 sifting 성질에 따라 3+32=123+3^2=12이다.

(3) ∫03et−2 δ(2t−4) dt\displaystyle\int_{0}^{3}e^{t-2}\,\delta(2t-4)\,dt

풀이

½이다. δ(2t−4)=δ(2(t−2))=12 δ(t−2)\delta(2t-4)=\delta\bigl(2(t-2)\bigr)=\tfrac12\,\delta(t-2)이고 t=2t=2는 구간 [0,3][0,3] 안에 있으므로 12e0=12\tfrac12e^{0}=\tfrac12이다. 배율을 풀지 않고 t=2t=2만 대입하면 틀린 답 1이 나온다.

(4) ∫−∞tδ(τ) dτ\displaystyle\int_{-\infty}^{t}\delta(\tau)\,d\tau

풀이

u(t)u(t). t<0t<0이면 구간 (−∞,t](-\infty,t]에 임펄스 위치 0이 들어가지 않아 0이고, t>0t>0이면 들어가므로 1이다. 임펄스를 처음부터 누적 적분하면 단위 계단이 된다. 거꾸로 u(t)u(t)를 미분하면 δ(t)\delta(t)이다. t=0t=0에서의 값은 약속으로 정하며, 이 책에서는 u(0)=1u(0)=1을 쓴다.