[SS] Resolution of DFT

DFT의 경우 spectrum도 discrete하게 존재하기 때문에

spectrum에서의 sampling interval(실제 값을 가진 샘플의 간격)이 지나치게 넓을 경우,

Picket Fence Effect로 인한 문제점 발생. (=너무 듬성듬성하게 샘플링하면 문제가 되는 것과 같음)

• DFT의 Picket Fence Effect → 스펙트럼 샘플의 간격이 넓으면 실제 스펙트럼의 모양을 파악하기 어려움.
  → 스펙트럼의 해상도 (spectrum resolution)가 중요함!!!
• DFT에 의한 스펙트럼 해상도 → $x[n]$의 샘플 수에 의존 ( $\langle 2\pi \rangle $구간 고정됨)

picket fence에서 가로막힌 간격(샘플간의 간격)으로 인해 spectrum의 실제 모습과 다른 signal을 얻게 됨.

 

 

---

Zero Padding

참고로, 만약 time domain에서 $N_1$개의 유효 데이터가 있는데,

주파수 영역에서 $N_2$ 개의 샘플이 요구될 경우, 다음의 zero-padding으로 처리하는게 일반적임.

• $N_2 - N_1$개의 $0$을 뒤에 첨가
• 데이터 수는 증가시키되 스펙트럼 형태 변화는 생기지 않도록 처리함 : 단순히 frequency resolution만 증가시킴.
• 단, zero-padding에 의한 frequency resolution 증가는 스펙트럼 정확도(accuracy)의 향상이 아님. (accuracy는 동일)

정확도(accuracy)은 신호의 유효 샘플 수 $N_1$ 을 더 늘려야 함

유효 샘플 수 $N_1$ 의 증가 → 해상도 향상 & 주파수 중첩 감소 (=정확성 증가)

 

zero padding을 한 경우, 스펙트럼의 모양은 그대로 유지되면서 frequency resolution만 향상됨. 반면 보다 많은 샘플을 직접 얻은 경우, 보다 정확한 spectrum이 만들어짐 (위 그림에서 신호가 올라가는 기울기가 보다 정확함을 알 수 있음)

---

Example : Zero Padding

$N=4$인 signal $x[n]$의 DTFT는 다음과 같음.
$$ X(e^{j\omega n})=\displaystyle \sum^{3}_{n=0}x[n]e^{-j\omega n} $$
$x[n]$을 zero padding ($N: 4 \to 8$)수행 (신호의 전체 길이만 증가할 뿐, 샘플간의 간격은 그대로임).

DTFT관점에서 보면 다음과 같음.

$$ x_{zp}[n] = \left[ x[0], x[1], x[2], x[3], 0, 0, 0, 0\right] \\ X'(e^{j\omega n})=\displaystyle \sum^{7}_{n=0} x_{zp}[n]e^{-j\omega n} = \sum^3_{n=0}x[n]e^{-j\omega n}\\ X'(e^{j\omega n})=X(e^{j\omega n}) $$

• spectrum의 범위(sampling freq.에 의해 좌우)는 그대로이면서,
• spectrum의 샘플 수가 8개로 늘어났기 때문에 spectrum간 샘플의 간격이 1/2로 줄어들었음 (동일한 주파수대역에 2배 많은 샘플이 위치).
• 하지만, spectrum의 모양은 위 식에서 보이듯이 전혀 변화가 없음.
• 즉, 같은 모양에서 샘플이 보다 촘촘히 들어가 있게 됨.

이를 DFT관점에서 보면 다음과 같음

$$\begin{array}{clll} X'[0] &= \displaystyle \sum^{7}_{n=0} x_{zp}[n] &= \sum^{3}_{n=0} x_{zp}[n] &=\color{red}{ \sum^{3}_{n=0} x[n] =X[0]} \\ X'[1] &= \displaystyle \sum^{7}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}n} &= \sum^{3}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}n} \\ X'[2] &= \displaystyle \sum^{7}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}2n}&= \sum^{3}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}2n}&=\color{red}{ \sum^{3}_{n=0} x[n]e^{-j\frac{2\pi}{4}n}=X[1]}\\ X'[3] &= \displaystyle \sum^{7}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}3n}&= \sum^{3}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}3n}\\ X'[4] &= \displaystyle \sum^{7}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}4n}&= \sum^{3}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}4n}&=\color{red}{ \sum^{3}_{n=0} x[n]e^{-j\frac{2\pi}{4}2n}=X[2]}\\ X'[5] &= \displaystyle \sum^{7}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}5n}&= \sum^{3}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}5n}\\ X'[6] &= \displaystyle \sum^{7}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}6n}&= \sum^{3}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}6n}&=\color{red}{ \sum^{3}_{n=0} x[n]e^{-j\frac{2\pi}{4}3n}=X[3]}\\ X'[7] &= \displaystyle \sum^{7}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}7n}&= \sum^{3}_{n=0} x_{zp}[n] e^{-j\frac{2\pi}{8}7n}\\ X'[8] &= X'[0]&&=X[0]\\ \vdots \end{array}$$

• $X'[0],X'[2],X'[4],X'[6]$은 zero padding전의 $x[n]$의 $X[0],X[1],X[2],X[3]$과 같은 값을 가지며,
• 이들 값들 사이에 $X'[1],X'[3],X'[5],X'[7]$이 추가되는 형태임을 알 수 있음.

참고 : DFT

$$ X[k] = \displaystyle \sum^{N-1}_{n=0} x[n] e^{-j\frac{2\pi}{N}kn},\quad k=0,1,\dots,N-1$$

 

---

Frequency Resolution

Spectrum에서 인접한 두 frequency component를 구분해낼 수 있는 정도.

Spectrum에서 Sampling interval을 의미함.

---

Digital Frequency Resolution

$$\Delta \omega = \frac{2\pi}{N}$$ or $$\Delta F=\frac{1}{N}$$

---

Analog Frequency Resolution

$$\Delta f = \frac{f_s}{N}= \frac{\text{한 주기 구간에 해당하는 주파수,}f_s} {\text{# of samples}}=\frac{1}{t_s}=\frac{1}{NT_s}$$

 

• $t_s$ : time domain signal의 period (지속시간이라고도 불림), 주기 $T$에 해당함.
• $T_s$ : sampling interval.
• $f_s$ : spectrum에서의 한주기에 해당.

---

Example : Frequency Resolution

지속 시간이 2초($t_s$ or $T$)이고 해당 신호의 유효 대역폭($f_b$)이 10kHz 인 continuous signal에 대해

주파수 스펙트럼에서 aliasing(or 절단 왜곡)이 발생하지 않도록

DFT(Discrete Fourier Transform)을 해야 한다.

이때 DFT에서 사용할 샘플의 N과 frequency resolution Δf을 구하시오.

 

풀이

$$f_s = \frac{1}{T_s}=\frac{N}{t_s} \ge 2f_b$$

 

위 식을 통해, $N$은 다음과 같음.

$$\frac{N}{t_s} \ge 2f_b \\ N \ge t_s \times 2f_b \\ N \ge 2 (\text{sec}) \times 2 \times 10 \times 10^{3} (\text{Hz}) \\ \therefore N \ge 4 \times 10^4$$

 

frequency resolution은 다음과 같음.

$$\Delta f = \frac{f_s}{N}= \frac{1}{t_s} = \frac{1}{NT_s} = 0.5 (\text{Hz})$$