REMI Lab강의 노트

Signals and Systems L05

5장 푸리에 급수

실수 푸리에 급수

주기 신호를 상수와 코사인·사인의 합으로 쓰면, 각 항의 계수는 어떻게 정해질까? 앞 절의 방법대로 신호와 각 기저의 내적을 기저의 길이 제곱으로 나누면 된다.

주기가 2π2\pi인 함수 f(x)f(x)를 상수와 정수 차수의 코사인·사인의 합으로 쓴다.

f(x)=a0+∑n=1∞(ancos⁡nx+bnsin⁡nx)f(x)=a_0+\sum_{n=1}^{\infty}\left(a_n\cos nx+b_n\sin nx\right) a0=12π∫−ππf(x) dxan=1π∫−ππf(x)cos⁡nx dxbn=1π∫−ππf(x)sin⁡nx dx\begin{aligned} a_0&=\frac{1}{2\pi}\int_{-\pi}^{\pi}f(x)\,dx\\ a_n&=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\cos nx\,dx\\ b_n&=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\sin nx\,dx \end{aligned}

a0a_0는 한 주기의 평균, 곧 직류 성분이다. ana_n과 bnb_n은 ff와 cos⁡nx\cos nx, sin⁡nx\sin nx의 내적을 기저의 길이 제곱으로 나눈 값이다. 한 주기에서 ∫12 dx=2π\int 1^2\,dx=2\pi, ∫cos⁡2nx dx=∫sin⁡2nx dx=π\int\cos^2 nx\,dx=\int\sin^2 nx\,dx=\pi이므로 앞에 붙는 수가 1/(2π)1/(2\pi)와 1/π1/\pi로 다르다.

상수항을 a0/2a_0/2로 쓰고 a0=1π∫−ππf dxa_0=\frac{1}{\pi}\int_{-\pi}^{\pi}f\,dx로 정의하는 교재도 있다. 그 a0a_0는 이 장의 a0a_0의 두 배이다. 다른 자료를 읽을 때는 상수항의 정의부터 확인해야 한다.

주기가 TT인 신호 x(t)x(t)는 x=ω0tx=\omega_0 t, ω0=2π/T\omega_0=2\pi/T로 바꿔 쓰면 된다.

x(t)=a0+∑n=1∞(ancos⁡nω0t+bnsin⁡nω0t)x(t)=a_0+\sum_{n=1}^{\infty}\left(a_n\cos n\omega_0 t+b_n\sin n\omega_0 t\right) a0=1T∫Tx(t) dtan=2T∫Tx(t)cos⁡nω0t dtbn=2T∫Tx(t)sin⁡nω0t dt\begin{aligned} a_0&=\frac{1}{T}\int_T x(t)\,dt\\ a_n&=\frac{2}{T}\int_T x(t)\cos n\omega_0 t\,dt\\ b_n&=\frac{2}{T}\int_T x(t)\sin n\omega_0 t\,dt \end{aligned}

∫T\int_T는 길이 TT인 아무 한 주기에서 적분한다는 뜻이다.

같은 차수의 코사인과 사인은 하나의 정현파로 합쳐진다.

ancos⁡nx+bnsin⁡nx=Ancos⁡(nx−θn)An=an2+bn2,cos⁡θn=anAn,sin⁡θn=bnAn\begin{gathered} a_n\cos nx+b_n\sin nx=A_n\cos(nx-\theta_n)\\ A_n=\sqrt{a_n^2+b_n^2},\quad \cos\theta_n=\frac{a_n}{A_n},\quad \sin\theta_n=\frac{b_n}{A_n} \end{gathered}

그래서 (an,bn)(a_n,b_n)은 nn차 성분의 크기와 위상을 담은 2차원 좌표이다.

그림은 같은 주파수로 도는 두 벡터를 머리와 꼬리로 이어 더하고, 같은 시각의 실수부를 파형으로 그린다. A1=anA_1=a_n, A2=bnA_2=b_n, φ=−90∘\varphi=-90^\circ로 두었으므로 A1cos⁡θ+A2cos⁡(θ−90∘)=ancos⁡θ+bnsin⁡θA_1\cos\theta+A_2\cos(\theta-90^\circ)=a_n\cos\theta+b_n\sin\theta이다. 두 벡터는 직각으로 이어지고, 합은 같은 주파수의 정현파 하나이며 그 진폭은 합 벡터(초록)의 길이이다.

같은 주파수로 도는 두 벡터의 합과 그 실수부 파형을 함께 그린다. 세 파형의 주파수는 1 Hz이다.

해 보기

  1. an=bn=1a_n=b_n=1이면 합의 진폭과 위상이 얼마인지 확인하라.
    관찰

    합의 진폭은 2=1.414\sqrt2=1.414, 위상은 −45°이다. cos⁡θ+sin⁡θ=2cos⁡(θ−45∘)\cos\theta+\sin\theta=\sqrt2\cos(\theta-45^\circ)이다.

  2. an=0.6a_n=0.6, bn=0.8b_n=0.8이면 합의 진폭과 위상이 얼마인지 확인하라.
    관찰

    합의 진폭은 0.62+0.82=1\sqrt{0.6^2+0.8^2}=1이고 위상은 −53.1°이다. cos⁡θn=0.6\cos\theta_n=0.6, sin⁡θn=0.8\sin\theta_n=0.8인 θn=53.1∘\theta_n=53.1^\circ만큼 늦은 코사인이다.

  3. 두 벡터가 반대 방향이면 합이 어떻게 되는지 보고, 같은 주파수 두 정현파의 합의 진폭이 무엇으로 정해지는지 생각해 보라.
    관찰

    합 벡터의 길이가 0이고 파형도 0이다. 같은 주파수 두 정현파를 더한 진폭은 두 진폭의 합이 아니라 벡터 합의 길이이다.

예제 1: 사각파

진폭이 AA(양의 상수)이고 주기가 2π2\pi인 사각파의 계수를 구한다.

f(x)={−A,−π<x<0−A,0<x<πf(x+2π)=f(x)\begin{gathered} f(x)=\begin{cases}-A, & -\pi<x<0\\ \phantom{-}A, & 0<x<\pi\end{cases}\\ f(x+2\pi)=f(x) \end{gathered}

한 주기의 평균이 0이므로 a0=0a_0=0이다. ff는 홀함수이고 cos⁡nx\cos nx는 짝함수이므로 곱 fcos⁡nxf\cos nx는 홀함수이고, 대칭 구간 [−π,π][-\pi,\pi]에서 적분하면 0이다. 직접 계산해도 같다.

an=1π[∫−π0(−A)cos⁡nx dx+∫0πAcos⁡nx dx]=1π[−Asin⁡nxn∣−π0+Asin⁡nxn∣0π]=0\begin{aligned} a_n&=\frac{1}{\pi}\left[\int_{-\pi}^{0}(-A)\cos nx\,dx+\int_{0}^{\pi}A\cos nx\,dx\right]\\ &=\frac{1}{\pi}\left[-A\frac{\sin nx}{n}\Big|_{-\pi}^{0}+A\frac{\sin nx}{n}\Big|_{0}^{\pi}\right]=0 \end{aligned} bn=1π[∫−π0(−A)sin⁡nx dx+∫0πAsin⁡nx dx]=1π[Acos⁡nxn∣−π0−Acos⁡nxn∣0π]=Anπ[cos⁡0−cos⁡(−nπ)−cos⁡nπ+cos⁡0]=2Anπ(1−cos⁡nπ)\begin{aligned} b_n&=\frac{1}{\pi}\left[\int_{-\pi}^{0}(-A)\sin nx\,dx+\int_{0}^{\pi}A\sin nx\,dx\right]\\ &=\frac{1}{\pi}\left[A\frac{\cos nx}{n}\Big|_{-\pi}^{0}-A\frac{\cos nx}{n}\Big|_{0}^{\pi}\right]\\ &=\frac{A}{n\pi}\left[\cos 0-\cos(-n\pi)-\cos n\pi+\cos 0\right]\\ &=\frac{2A}{n\pi}\left(1-\cos n\pi\right) \end{aligned}

cos⁡nπ=(−1)n\cos n\pi=(-1)^n이므로 짝수 nn에서는 bn=0b_n=0, 홀수 nn에서는 bn=4A/(nπ)b_n=4A/(n\pi)이다. 따라서

f(x)=4Aπ(sin⁡x+13sin⁡3x+15sin⁡5x+⋯ )f(x)=\frac{4A}{\pi}\left(\sin x+\frac{1}{3}\sin 3x+\frac{1}{5}\sin 5x+\cdots\right)

항을 앞에서부터 하나씩 더한 부분합을 S1,S2,S3,…S_1,S_2,S_3,\dots라고 하자.

S1=4Aπsin⁡xS2=4Aπ(sin⁡x+13sin⁡3x)S3=4Aπ(sin⁡x+13sin⁡3x+15sin⁡5x)\begin{aligned} S_1&=\frac{4A}{\pi}\sin x\\ S_2&=\frac{4A}{\pi}\left(\sin x+\frac{1}{3}\sin 3x\right)\\ S_3&=\frac{4A}{\pi}\left(\sin x+\frac{1}{3}\sin 3x+\frac{1}{5}\sin 5x\right) \end{aligned}

S1S_1은 사각파를 흉내 내는 사인 하나이다. 여기에 4A3πsin⁡3x\frac{4A}{3\pi}\sin 3x를 더하면 피크가 눌리고 양 끝이 가팔라지며, 4A5πsin⁡5x\frac{4A}{5\pi}\sin 5x를 더하면 평평한 부분이 더 평평해진다. 더하는 성분의 진폭은 차수에 반비례해서 작아진다.

그림은 A=1A=1, 주기 T=1T=1 s인 같은 사각파(x=2πtx=2\pi t)와 그 부분합을 두 주기 그린다. NN은 합에 넣을 가장 높은 차수이므로 N=1,3,5N=1,3,5가 각각 S1,S2,S3S_1,S_2,S_3이다. 위 줄은 차수 NN에서 새로 더하는 항이고, 아래 그림의 회색은 한 차수 전(N−1N-1까지)의 부분합, 초록은 NN까지의 부분합이다. NN이 하나 오를 때마다 위 줄의 항이 회색 곡선에 더해져 초록 곡선이 되고, 평평한 부분은 펴지며 점프 쪽은 가팔라진다.

위 줄의 붉은 곡선이 차수 N의 항이고, 아래 그림의 점선이 사각파, 회색이 N − 1까지, 초록이 |k| ≤ N인 성분을 모두 더한 부분합, 회색과 초록 사이의 옅은 띠가 위 줄의 항이며, 맨 아래 줄은 지금까지 더한 계수 bₖ이다. 처음 상태 N = 1의 부분합은 (4/π)sin 2πt 하나이고, 피크 4/π ≈ 1.27이 사각파의 1보다 높다.

해 보기

  1. N=3N=3까지 더하면 부분합의 피크와 점프 쪽이 어떻게 바뀌는지 확인하라.
    관찰

    위 줄의 항 43πsin⁡6πt\frac{4}{3\pi}\sin6\pi t가 더해져(S2S_2) 피크가 눌리고 점프 쪽이 가팔라진다. 회색 곡선은 N=2N=2까지의 합인데, b2=0b_2=0이라 N=1N=1과 같은 곡선이다.

  2. N=5N=5까지 더하면 평평한 부분과 피크가 어떻게 바뀌는지 확인하라.
    관찰

    45πsin⁡10πt\frac{4}{5\pi}\sin10\pi t가 더해져(S3S_3) 평평한 부분이 더 평평해지고, 피크는 1.188로 낮아진다.

  3. 삼각파를 N=3N=3까지 더한 부분합이 삼각파에 얼마나 가까운지 확인하라.
    관찰

    두 성분만으로 부분합이 삼각파와 거의 겹친다. 꼭짓점에서도 크기가 0.90으로 1에 조금 못 미칠 뿐이다.

계수는 한 주기의 사각파에 kk차 기준 cos⁡2πkt\cos 2\pi kt, −sin⁡2πkt-\sin 2\pi kt를 곱한 두 곱의 평균으로 구한다. 그림 아래 첫 줄은 두 평균으로 쓴 복소 계수 ckc_k로, 뒤의 복소 급수 절에서 다룬다. 둘째 줄은 같은 평균에 2를 곱한 실수 급수의 계수 aka_k, bkb_k이고, 두 줄은 ck=(ak−jbk)/2c_k=(a_k-jb_k)/2로 이어진다. k=3k=3에서는 x×(−sin⁡)x\times(-\sin) 곱의 면적이 한쪽 부호로 남고, 둘째 줄의 b3b_3이 앞에서 구한 4/(3π)4/(3\pi)이다.

맨 위 그림은 한 주기의 입력(검정)과 기준 cos 2πkt(파랑 점선), −sin 2πkt(보라 점선)이고, 그 아래 두 그림은 곱과 그 부호 있는 면적, 점선은 한 주기 평균이다. 그림 밑 줄은 그 면적을 주기 T = 1 s로 나눠 평균을 구하고, 맨 아래 두 줄은 두 평균으로 쓴 cₖ와 aₖ = 2⟨x cos⟩, bₖ = −2⟨x(−sin)⟩이다(k = 0이면 a₀ = c₀가 평균이다).

해 보기

  1. k=3k=3에서 두 곱의 평균을 읽고, c3c_3, a3a_3, b3b_3을 구하라.
    관찰

    cos 곱의 평균은 0이고, −sin⁡-\sin 곱은 면적 −0.2122-0.2122를 T=1T=1 s로 나눈 평균이 −b3/2=−2/(3π)≈−0.212-b_3/2=-2/(3\pi)\approx-0.212이다. 그래서 c3=−j0.212c_3=-j0.212이고, 둘째 줄은 a3=0a_3=0, b3=0.424=4/(3π)b_3=0.424=4/(3\pi)이다.

  2. k=2k=2에서 두 곱의 평균이 어떻게 되는지 확인하라.
    관찰

    두 곱의 면적이 모두 상쇄되어 c2=0c_2=0이다. 사각파에는 짝수 차수가 없다.

  3. 삼각파의 k=3k=3 계수를 구하고, k=1k=1의 계수와 비교하라.
    관찰

    이번에는 −sin⁡-\sin 곱의 평균이 0이고 cos 곱의 평균 −0.045-0.045만 남아 a3=−0.090a_3=-0.090, b3=0b_3=0이다. k=1k=1의 −0.405-0.405와 비교하면 1/9이다.

입력을 삼각파로 바꾸면 코사인 쪽 계수만 남고, 높은 차수의 계수가 차수의 제곱에 반비례해 더 빨리 줄어든다. 삼각파의 기울기는 사각파이고, ∫sin⁡(nω0t) dt=−cos⁡(nω0t)/(nω0)\int\sin(n\omega_0t)\,dt=-\cos(n\omega_0t)/(n\omega_0)이므로 1/n1/n인 사각파 계수에 1/n1/n이 한 번 더 곱해진다. 값이 끊기지 않는 삼각파는 날카로운 모서리를 만들 높은 주파수 성분이 덜 필요하다.

nn차 성분의 계수가 an=0.6a_n=0.6, bn=−0.8b_n=-0.8이다. ancos⁡nx+bnsin⁡nx=Ancos⁡(nx−θn)a_n\cos nx+b_n\sin nx=A_n\cos(nx-\theta_n)의 AnA_n과 θn\theta_n을 고르라.

  1. (가) An=1A_n=1, θn=53.1∘\theta_n=53.1^\circ
  2. (나) An=1A_n=1, θn=−36.9∘\theta_n=-36.9^\circ
  3. (다) An=1.4A_n=1.4, θn=−53.1∘\theta_n=-53.1^\circ
  4. (라) An=1A_n=1, θn=−53.1∘\theta_n=-53.1^\circ
풀이

답 (라) An=1A_n=1, θn=−53.1∘\theta_n=-53.1^\circ

  • (가) sin⁡θn=bn/An=−0.8\sin\theta_n=b_n/A_n=-0.8이라 θn\theta_n은 음수이다.
  • (나) cos⁡θn=0.6\cos\theta_n=0.6을 만족하는 각은 ±53.1°이다. 36.9°는 cos⁡θn=0.8\cos\theta_n=0.8일 때이다.
  • (다) 진폭은 두 계수의 합이 아니라 an2+bn2\sqrt{a_n^2+b_n^2}이다.

An=0.62+0.82=1A_n=\sqrt{0.6^2+0.8^2}=1이고 cos⁡θn=0.6\cos\theta_n=0.6, sin⁡θn=−0.8\sin\theta_n=-0.8이므로 θn=−53.1∘\theta_n=-53.1^\circ이다. 따라서 0.6cos⁡nx−0.8sin⁡nx=cos⁡(nx+53.1∘)0.6\cos nx-0.8\sin nx=\cos(nx+53.1^\circ)이다.

상수항을 a0/2a_0/2로 쓰는 교재에서 어떤 주기 신호의 a0a_0가 3이다. 이 장의 규약(a0a_0가 한 주기 평균)으로 쓴 a0a_0를 구하라.

풀이

답 1.5

힌트 그 교재의 상수항 a0/2a_0/2가 곧 한 주기 평균이다.

그 교재의 a0=1π∫−ππf dxa_0=\frac1\pi\int_{-\pi}^{\pi}f\,dx는 평균의 두 배이다. 상수항 a0/2=1.5a_0/2=1.5가 평균이므로 이 장의 a0=1.5a_0=1.5이다.