跳到主要内容

東京大学 情報理工学系研究科 電子情報学専攻 2018年8月実施 専門 第5問

Author

diohabara

Description

信号 f(t)f(t) が与えられたとき, f(t)f(t) を時間間隔 tst_s で標本化することを考える。時間間隔 tst_s にデルタ関数 δ(t)\delta(t) が並ぶ信号を単位インパルス列 δs(t)\delta_s(t) と呼ぶ。すなわち,

δs(t)=i=δ(tits)\delta_s(t) = \sum_{i = -\infty}^{\infty}\delta(t - it_s)

この時, f(t)f(t) の標本化された信号 fs(t)f_s(t) は, fs(t)=f(t)δs(t)f_s(t) = f(t)\cdot\delta_s(t) と表される。

  以下の問いに答えよ。

(1) δs(t)\delta_s(t) を周期 tst_s の周期信号と考え, フーリエ級数展開せよ。

(2) δs(t)\delta_s(t) のフーリエ変換 Δs(ω)\Delta_s(\omega) を求めよ。ただし ωs=2πts\omega_s = \frac{2\pi}{t_s} とせよ。

(3) f(t)f(t) 並びに fs(t)f_s(t) のフーリエ変換をそれぞれ F(ω)F(\omega)Fs(ω)F_s(\omega) とする。 Fs(ω)F_s(\omega)F(ω)F(\omega) を用いて表せ。

(4) 折り返し歪 (エイリアシング) の定義を述べよ。また, (3) の結果においてどのような現象となるのか説明せよ。さらに, 折り返し歪が起きないために F(ω)F(\omega) が満足すべき条件を ωs\omega_s を用いて述べよ。

必要に応じて以下を使ってよい。

周期 TT の信号 x(t)x(t) のフーリエ級数展開:

x(t)=i=ciej2πitTci=1TT0T0+Tx(t)ej2πitTdt\begin{aligned} x(t) &= \sum_{i = -\infty}^{\infty}c_i e^{j\frac{2\pi it}{T}} \\ c_i &= \frac{1}{T}\int_{T_0}^{T_0 + T} x(t)e^{-j\frac{2\pi it}{T}}dt \end{aligned}

信号 xx のフーリエ変換 X(ω)X(\omega) :

X(ω)=x(t)ejωtdtX(\omega) = \int_{-\infty}^{\infty}x(t)e^{-j\omega t}dt

信号 x(t)=1x(t) = 1 のフーリエ変換は 2πδ(ω)2\pi\delta(\omega) . 

x1(t)x_1(t)x2(t)x_2(t) の畳み込み:

x1(t)x2(t)=x1(t)x2(tt)dtx_1(t)*x_2(t) = \int_{-\infty}^{\infty}x_1(t')x_2(t - t')dt'

信号 x1(t)x_1(t)x2(t)x_2(t) ののフーリエ変換をそれぞれ X1(ω)X_1(\omega)X2(ω)X_2(\omega) とすると, x1(t)x2(t)x_1(t)*x_2(t) のフーリエ変換は , X1(ω)X2(ω)X_1(\omega) \cdot X_2(\omega).

同様に x1(t)x2(t)x_1(t) \cdot x_2(t) のフーリエ変換は 12πX1(ω)X2(ω)\frac{1}{2\pi}X_1(\omega) * X_2(\omega).

Kai

(1)

δs(t)=i=δ(tits)\delta_{s}(t) = \sum_{i = -\infty}^{\infty}\delta(t - it_s)11 周期分 (ts2t<ts2)(-\frac{t_s}{2} \le t < \frac{t_s}{2}) 切り出して、係数を求める。問題文にもあるように係数を求める式は次の通り。

ci=1TT/2T/2f(t)ej2πitTdtc_i = \frac{1}{T}\int_{-T/2}^{T/2}f(t)e^{-j\frac{2\pi it}{T}dt}

この場合、T=ts,f(t)=δ(t)T = t_s ,f(t) = \delta(t) なので、これを代入して

ci=1tsts/2ts/2f(t)ej2πittsdt=1ts1=1ts\begin{aligned} c_i &= \frac{1}{t_s}\int_{-t_s/2}^{t_s/2}f(t)e^{-j\frac{2\pi it}{t_s}dt} \\ &= \frac{1}{t_s} \cdot 1 = \frac{1}{t_s} \end{aligned}

よって、係数は等しく 1ts\frac{1}{t_s} であり、フーリエ級数展開の結果は

δs(t)=i=1tsej2πitts\delta_s(t) = \sum_{i = -\infty}^{\infty}\frac{1}{t_s}e^{j\frac{2\pi it}{t_s}}

(2)

(1) の結果をフーリエ変換の公式に代入する。ただし、積分と和の交換、フーリエ変換の公式 (12πδ(ω),f(t)ejω0tF(ωω0))(1 \rightarrow 2\pi \delta(\omega),f(t) \rightarrow e^{j\omega_0 t} \rightarrow F(\omega - \omega_0))

Δs(ω)=δs(t)ejωtdt={i=1tsej2πitts}=1ts{ej(ωiωs)t}=1tsi=2πδ(ωiωs)=ωsi=δ(ωiωs)\begin{aligned} \Delta_s(\omega) &= \int_{-\infty}^{\infty}\delta_s(t)e^{-j\omega t}dt \\ &= \int_{-\infty}^{\infty}\bigg\{\sum_{i = -\infty}^{\infty}\frac{1}{t_s}e^{j\frac{2\pi it}{t_s}}\bigg\} \\ &= \frac{1}{t_s}\sum_{\infty}^{\infty}\bigg\{\int_{-\infty}^{\infty}e^{-j(\omega - i\omega_s)t}\bigg\} \\ &= \frac{1}{t_s}\sum_{i = -\infty}^{\infty}2\pi\delta(\omega - i\omega_s) \\ &= \omega_s \sum_{i = -\infty}^{\infty}\delta(\omega - i\omega_s) \end{aligned}

(3)

fs(t)=f(t)δs(t)f_s(t) = f(t) \cdot \delta_s(t) の両辺をフーリエ変換する。

下の計算では和と積分の入れ替えや、(2) で得た結論や、畳込みの定義式やフーリエ変換 を使った。

Fs(ω)=12πF(ω)Δs(ω)=1tsF(ω)(i=δ(ωiωs))=1tsi=(F(ω)δ(ωiωsω)dω)=1tsi=F(ωiωs)\begin{aligned} F_s(\omega) &= \frac{1}{2\pi}F(\omega) * \Delta_s(\omega) \\ &= \frac{1}{t_s}F(\omega) * \bigg(\sum_{i = -\infty}^{\infty}\delta(\omega - i\omega_s)\bigg) \\ &= \frac{1}{t_s}\sum_{i = -\infty}^{\infty}\bigg(\int_{-\infty}^{\infty}F(\omega')\delta(\omega - i\omega_s - \omega')d\omega'\bigg) \\ &= \frac{1}{t_s}\sum_{i = -\infty}^{\infty}F(\omega - i\omega_s) \end{aligned}

(4)

  • エイリアシングの定義

エイリアシングとは、サンプリングに従って信号の一部が本来の周波数とは異なる周波 数の成分として混入してしまい、波形に歪みが生じることを言う。

  • (3) においてどのような減少となるの

サンプリング周波数 ωs\omega_s に対して周波数 ω,ω±ωs,ω±2ω,\omega,\omega \pm \omega_s ,\omega \pm 2\omega ,\dotsb の成分がすべて Fs(ω)F_s(\omega) 上の同じ点に重なってしまうため、サンプリング後の信号 Fs(ω)F_s(\omega) を見たときに F(ω)F(\omega) のど この周波数由来なのか判別不可能になる。

  • F(ω)F(\omega) が満足すべき条件 ω>ωs2\omega > \frac{\omega_s}{2} において F(ω)=0F(\omega) = 0 を満たすこと。