REMI Lab강의 노트

Signals and Systems L07

7장 이산시간 시스템의 푸리에 해석

이산시간 정현파와 Ω

DTFT의 Ω\Omega는 무엇을 재는 양일까? 샘플 하나마다 얼마나 도는지, 곧 회전 속도를 잰다.

샘플마다 Ω씩 도는 회전

4장에서 ejωte^{j\omega t}는 단위원 위를 각속도 ω\omega로 쉬지 않고 도는 점이었다. 수열 ejΩne^{j\Omega n}은 같은 원 위에서 샘플 하나마다 Ω\Omega rad씩 건너뛰는 점이다. n=0n=0에서 1에 있고, n=1n=1에서 각 Ω\Omega, n=2n=2에서 각 2Ω2\Omega에 놓인다.

ejΩn=cos⁡(Ωn)+jsin⁡(Ωn)e^{j\Omega n}=\cos(\Omega n)+j\sin(\Omega n)

이산시간 정현파 x[n]=Acos⁡(Ωn+θ)x[n]=A\cos(\Omega n+\theta)는 이 회전의 가로 좌표이다. 진폭 AA를 곱하고 시작 각을 θ\theta만큼 돌렸을 뿐이다. 그래서 Ω\Omega를 이산시간 정현파의 주파수라고 부르고, 단위는 rad/sample이다.

이산시간 정현파는 2장의 연속 정현파와 닮았지만 두 가지가 다르다. 항상 주기적인 것은 아니고, Ω\Omega가 커진다고 늘 빨라지지도 않는다.

주기적인 경우와 아닌 경우

x[n+N]=x[n]x[n+N]=x[n]이 되려면 NN샘플 동안 도는 각 ΩN\Omega N이 2π2\pi의 정수배여야 한다. 곧 회전이 NN샘플 만에 정확히 몇 바퀴를 돌고 제자리로 와야 한다. ΩN=2πm\Omega N=2\pi m을 풀면

Ω2π=mN\frac{\Omega}{2\pi}=\frac{m}{N}

이므로, Ω/2π\Omega/2\pi가 유리수일 때만 주기적이다. 기약분수 m/Nm/N의 분모 NN이 주기이다. cos⁡(πn/6)\cos(\pi n/6)은 Ω/2π=1/12\Omega/2\pi=1/12라 주기가 12이고, cos⁡(3πn/8)\cos(3\pi n/8)은 3/163/16이라 주기가 16이다(그동안 연속 곡선은 세 번 돈다). cos⁡(n)\cos(n)은 Ω/2π=1/2π\Omega/2\pi=1/2\pi가 무리수라 샘플이 정확히 되풀이되지 않는다. 샘플을 이은 모양은 주기 2π2\pi인 곡선처럼 보이지만 수열 자체는 주기가 없다.

Ω의 2π 주기성

Ω\Omega는 2π2\pi마다 같은 수열을 만든다. 한 바퀴를 더 돌면 같은 자리에 오기 때문이다. nn이 정수이므로 cos⁡((Ω+2π)n+θ)=cos⁡(Ωn+θ)\cos((\Omega+2\pi)n+\theta)=\cos(\Omega n+\theta)이다. 따라서 0≤Ω<2π0\le\Omega<2\pi만 보면 된다. Ω=0\Omega=0이면 점이 움직이지 않아 모든 샘플이 같은 상수이다.

가장 빠른 진동

Ω\Omega를 0에서 π\pi까지 키우면 진동이 빨라진다. Ω=π\Omega=\pi에서 점은 샘플마다 반 바퀴를 돌아 1과 −1-1을 오간다. cos⁡(πn)=(−1)n\cos(\pi n)=(-1)^n은 샘플마다 부호가 바뀌어 가장 빠르다. π\pi를 넘으면 다시 느려진다.

cos⁡((2π−Ω)n)=cos⁡(2πn−Ωn)=cos⁡(Ωn)\begin{aligned} \cos((2\pi-\Omega)n)&=\cos(2\pi n-\Omega n)\\ &=\cos(\Omega n) \end{aligned}

이기 때문이다. Ω=15π/8\Omega=15\pi/8의 샘플은 π/8\pi/8의 샘플과 똑같다. 샘플 위로 빠르게 흔들리는 연속 곡선을 그릴 수도 있지만, 같은 점을 지나는 곡선은 그것 말고도 많고 샘플만으로는 구별하지 못한다.

Ω\Omega 샘플 주기 NN
00 1,1,1,…1,1,1,\ldots 1
π/8\pi/8 cos⁡(πn/8)\cos(\pi n/8) 16
π/2=4π/8\pi/2=4\pi/8 1,0,−1,0,…1,0,-1,0,\ldots 4
π=8π/8\pi=8\pi/8 1,−1,1,−1,…1,-1,1,-1,\ldots 2
15π/815\pi/8 cos⁡(πn/8)\cos(\pi n/8)과 같음 16
2π=16π/82\pi=16\pi/8 1,1,1,…1,1,1,\ldots 1

회전 방향

복소 정현파 ejΩne^{j\Omega n}에서는 Ω\Omega와 2π−Ω2\pi-\Omega가 같지 않다. ej(2π−Ω)n=e−jΩne^{j(2\pi-\Omega)n}=e^{-j\Omega n}이므로 같은 속도로 반대 방향으로 돈다. 15π/815\pi/8은 샘플마다 −π/8-\pi/8씩 거꾸로 도는 회전이다. 코사인은 두 방향의 회전을 모두 담고 있어서 방향을 구별하지 못할 뿐이다. 아래 그림에서 Ω=π/8\Omega=\pi/8이면 샘플마다 22.5°씩 돌아 8샘플에 반 바퀴를 돌고, Ω=π\Omega=\pi이면 오른쪽과 왼쪽만 오간다. Ω=15π/8\Omega=15\pi/8은 샘플마다 337.5°를 돌지만, 샘플만 보면 시계 방향으로 22.5°씩 도는 것과 같다.

e^(jΩn)을 n = 0…7을 페이저로 그린 그림. Ω = π/8은 샘플마다 반시계로 22.5°, π/2는 90°, π는 180°씩 돌고, 15π/8은 시계 방향으로 22.5°씩 도는 것과 같다n01234567샘플당π/8+22.5°π/2+90°π180°15π/8−22.5°
한 줄이 한 Ω이고, 페이저 하나가 샘플 n = 0…7 하나의 점이다. 페이저는 ejΩn이 가리키는 방향이고, 반시계 방향이 +이다.

바퀴를 일정한 간격으로 찍은 영상에서도 같은 일이 일어난다. 1초에 7바퀴 도는 바퀴의 한 점을 1초에 10장 찍으면, 한 장 사이에 바퀴는 실제로 +252°를 돌지만 사진에서는 점이 −108°씩 옮겨 가는 것으로 보인다. 두 각은 같은 자리에서 끝난다. 한 장 사이에 도는 각은 Ω=2πf/fs\Omega=2\pi f/f_s이고, Ω\Omega가 π\pi와 2π2\pi 사이이면 사진은 반대 방향으로 2π−Ω2\pi-\Omega씩 도는 것과 같다. 그 회전의 주파수는 fa=f−fsf_a=f-f_s이고, 두 회전을 시간 축으로 펼친 두 나선은 사진을 찍는 순간마다 같은 점을 지난다.

바퀴의 한 점이 실제 f(Hz)로 돌고(+f는 반시계 방향) 사진을 1초에 fₛ장 찍으며, 사진을 한 원에 겹치고 찍은 순서대로 번호를 붙였다. 바깥 검은 화살표는 한 장 사이에 실제로 돈 각, 초록 화살표는 사진에 보이는 각이고, ②는 두 회전을 시간 축으로 펼친 것이다.

해 보기

  1. 사진 0부터 한 장씩 넘기며 사진 번호가 옮겨 가는 방향을 확인하라.
    관찰

    바퀴는 한 장마다 반시계 방향으로 Ω = 1.4π(252°)를 돌지만, 사진 번호는 시계 방향으로 0.6π(108°)씩 옮겨 간다. ②에서 +7 Hz 나선과 −3 Hz 나선이 사진마다 같은 점을 지난다(fₐ = −3 Hz).

  2. 실제 f가 3 Hz이면 사진 속 점이 어느 방향으로 도는지 확인하라.
    관찰

    Ω = 0.6π(108°)로 π보다 작아서, 사진 속 점도 실제와 같은 반시계 방향으로 돈다. 두 나선이 하나로 겹친다.

  3. 실제 f가 5 Hz(Ω = π)이면 사진만으로 도는 방향을 정할 수 있는지 확인하라.
    관찰

    한 장마다 반 바퀴씩 돌아서 사진은 두 자리만 오간다. +180°로 봐도 −180°로 봐도 같은 사진이라, 사진만으로는 도는 방향을 정할 수 없다.