跳到主要内容

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

Author

diohabara, 祭音Myyura

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=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πit/T,ci=1TT0T0+Tx(t)ej2πit/Tdt.x(t)=\sum_{i=-\infty}^{\infty}c_i e^{j2\pi it/T},\qquad c_i=\frac1T\int_{T_0}^{T_0+T}x(t)e^{-j2\pi it/T}\,dt.

信号 x(t)x(t) のフーリエ変換 X(ω)X(\omega)

X(ω)=x(t)ejωtdt.X(\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)dt.x_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)X_2(\omega)

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

题目描述

给定信号 f(t)f(t),以时间间隔 tst_s 对其采样。把每隔 tst_s 排列的狄拉克 δ\delta 函数组成的信号称为单位冲激列 δs(t)\delta_s(t)

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

采样后的信号为

fs(t)=f(t)δs(t).f_s(t)=f(t)\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=2\pi/t_s

(3) 设 f(t)f(t)fs(t)f_s(t) 的傅里叶变换分别为 F(ω)F(\omega)Fs(ω)F_s(\omega),用 F(ω)F(\omega) 表示 Fs(ω)F_s(\omega)

(4) 给出混叠(折叠失真)的定义,说明它在 (3) 的结果中表现为何种现象,并用 ωs\omega_s 写出为避免混叠而要求 F(ω)F(\omega) 满足的条件。

必要时可使用以下公式。周期为 TT 的信号 x(t)x(t) 的傅里叶级数为

x(t)=i=ciej2πit/T,ci=1TT0T0+Tx(t)ej2πit/Tdt.\begin{aligned} x(t)&=\sum_{i=-\infty}^{\infty}c_i e^{j2\pi it/T},\\ c_i&=\frac1T\int_{T_0}^{T_0+T}x(t)e^{-j2\pi it/T}\,\mathrm dt. \end{aligned}

信号 x(t)x(t) 的傅里叶变换定义为

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

且常数信号 x(t)=1x(t)=1 的傅里叶变换为 2πδ(ω)2\pi\delta(\omega)。卷积定义为

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

x1(t),x2(t)x_1(t),x_2(t) 的傅里叶变换分别为 X1(ω),X2(ω)X_1(\omega),X_2(\omega),则 x1x2x_1*x_2 的傅里叶变换为 X1X2X_1X_2,而 x1x2x_1x_2 的傅里叶变换为 12πX1X2\frac1{2\pi}X_1*X_2

Kai

(1)

一周期 [ts/2,ts/2)[-t_s/2,t_s/2) に含まれるインパルスは t=0t=0 の一つなので、

ci=1tsts/2ts/2δ(t)ejiωstdt=1ts.c_i=\frac1{t_s}\int_{-t_s/2}^{t_s/2}\delta(t)e^{-ji\omega_s t}\,dt =\frac1{t_s}.

したがって、

δs(t)=1tsi=ejiωst.\boxed{\delta_s(t)=\frac1{t_s}\sum_{i=-\infty}^{\infty}e^{ji\omega_s t}}.

(2)

ejiωste^{ji\omega_s t} のフーリエ変換は 2πδ(ωiωs)2\pi\delta(\omega-i\omega_s) だから、

Δs(ω)=2πtsi=δ(ωiωs)=ωsi=δ(ωiωs).\boxed{\Delta_s(\omega)=\frac{2\pi}{t_s}\sum_{i=-\infty}^{\infty}\delta(\omega-i\omega_s) =\omega_s\sum_{i=-\infty}^{\infty}\delta(\omega-i\omega_s)}.

(3)

時間領域の積は周波数領域の畳み込みになるので、

Fs(ω)=12π(FΔs)(ω)=1tsi=F(ωiωs).\begin{aligned} F_s(\omega)&=\frac1{2\pi}(F*\Delta_s)(\omega)\\ &=\boxed{\frac1{t_s}\sum_{i=-\infty}^{\infty}F(\omega-i\omega_s)}. \end{aligned}

(4)

エイリアシングとは、標本化によって異なる周波数の成分が区別できなくなり、別の周波数成分として現れる現象である。実際、ω\omegaω+kωs\omega+k\omega_s は標本時刻で同じ複素指数列を与える。

(3) では、ωs\omega_s 間隔で複製されたスペクトルが重なり、加算されることに相当する。通常の低域信号でこれを防ぐ条件は

F(ω)=0(ωωs2).\boxed{F(\omega)=0\qquad\left(|\omega|\ge\frac{\omega_s}{2}\right)}.

すなわち、信号の最高角周波数を ωm\omega_m とすれば、ωs>2ωm\omega_s>2\omega_m とすればよい。