━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━計測コラム

emm136 号用━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━

基礎からの周波数分析（7）－「時間信号のサンプリング」
今回は、今までお話しした内容（デルタ関数や畳み込みなど）を基にして、時系列信号のサンプリングに関してお話しします。
センサからなどの時間信号をディジタル信号処理するためには、そのアナログ信号を一定間隔毎にディジタル値（数値データ）に変換する必要があります。この処理をサンプリング（標本化）、その一定間隔をサンプリング周期（その逆数をサンプリング周波数）と呼びます。

YY

5

図 1 時間信号とそのサンプリング
図 1 は、左側のアナログ時間信号をサンプリング周期 τ (s)でサンプルした例で、右側のドット（・）の点における数値列だけが残ることになります。この場合ではサンプリング周波数は、1/τ (Hz)です。
さて、元のアナログ信号に含まれる情報を欠落することなく、サンプルするためにはサンプリング周波数をどれくらいとすればよいのでしょうか？

直感的に考えれば、大きなサンプリング周波数（すなわち細かい間隔）でサンプルすればよいと思われます
が、処理するデータ数がむやみに多くなり、好ましくないでしょう。今回のテーマは、サンプリング周波数はどのように決めるのがよいかという点に関する話です。
時間信号のサンプリングを説明する準備として、デルタ関数（インパルス関数）の性質について１つ追加します。
「等間隔のインパルス列のフーリエ変換は同じように等間隔のインパルス列となる」
ことです。
-1-

単位インパルス関数 δ (t)を周期 τ 毎に並べたインパルス列を；  t  n  h  t  .................................（1）n  とすると、そのフーリエ変換 H  f  は；H  f  

1

n

 f   n    .................................（2）.= \0]+= \-]PO=P

NT= NS= NR= NQ= N=

=

Q=

R= S=

T=02 ;4<:J?H6

NRO=

NQO=

NPO=

PO=

QO=

RO=-図 2 等間隔 τ のインパルス列とそのフーリエ変換
この性質の証明をします。
式（1）の h  t  は  t  が周期 τ で繰り返す周期関数となり、その複素フーリエ係数 c nは（インパルス関数の性質により）;

cn 1

 2

   t  e

 j n2   t

 2dt 1.................................（3）となり；h  t  1

ej n 2  t.................................（4） n  と書き換えることができます。

-2-

これから、式（4）のフーリエ変換は；H  f   11
e j n 2  t  e  j 2  f t dt e  j 2   f  n   t dt n   n  1n f     n  と式（2）を導きだすことができます。
図 2 にあるように、周期が τ 間隔で続くインパルス列時系列のフーリエ変換は、周波数間隔が 1/τ で続く周波数軸上のインパルス列となります。
図 3 は、サンプリング周期 τ で任意の時間信号 x ( t ) をサンプルした生成したサンプル化時系
列信号 x  t  のフーリエ変換を前回話した畳み込みにより導出した説明図です。図 3 の左側の図は時間領域で、右側の図は周波数領域となっています。

Z \Y]

U \V]

DMEDNE* 3,42B7@.
SVXY-/.

T5 \V]DPE
HG5H

FL5 FK5 FJ5 FI5 F5V-0.

W5 \Y]DOE

VX5I5J5 K5L5Y* 3,42B7@.

FJG5FIG5HG5FHG5?

IG5JG5=191

Z5\Y]

U 5\V]DQE

DRE

FL5 FK5 FJ5 FI5 F55I5J5 K5L5Y* 3,42B7@.

VX

HG5+ VXV

HG5 VX* :8;C6.

* <A>C6.

図 3 周期 τ でサンプルされた時系列信号のフーリエ変換を畳み込みにより求める説明図-3-以下、この図について、説明します。
ここで、時間信号 x ( t ) （ 図 3－(a)）の周波数成分は、最大 f m (Hz) まで（すなわち帯域制限されている）でそれ以上含まないと仮定します。（図 3－(b)）この仮定は、理論的には必要ないですが、現実的に扱う時間信号は有限とみなしてよいので、このように仮定します。
時間信号 x ( t ) をサンプリング周期 τ でサンプルすると、その時系列は；
x (τ) 、 x (2τ) 、 x (3τ) 、・ ・ ・ ・ x ( n τ) 、・ ・ ・となります。この時系列信号 x  t  は；x  t   x n   t  n  .................................（5）n  と表すことができます。すなわち、x ( t ) のサンプル時系列 x  t  は、等間隔 τ で無限に続くインパルス列でその振幅は同時刻の x ( t ) の値で与えられます。図 3 では、(a) と(c) との積である(e) となります。
次に、周波数領域に移ります。
X ( f ) は時間信号 x ( t ) のフーリエ変換（図 3－(b)）で、インパルス列（図 3－(c)）のフーリエ変換も上記で説明したように周波数軸上でのインパルス列（図.3－(d)）となります。
前回説明したように、時間軸上での積は周波数軸上では畳み込み演算（周波数畳み込み定理）に対応します。それ故、サンプル時系列信号 のフーリエ変換は図 3 の(f) のように X ( f )と周波数軸上のインパルス列との畳み込みとなり、サンプリング周波数 1/τ(Hz) で繰り返す周期的関数となります。
この図 3 の例では、時間信号の周波数帯域 f m に比べてサンプリング周波数（1/τ ）が充分大きい（サンプリング周期 τ が充分小さい）場合でしたが、そうでない場合を考えて見ます。
サンプリング周期が大きくなってくると周波数インパルスの間隔が狭くなってくるので、図4 のように、畳み込み同士が重なり合う現象がおき、フーリエ変換の結果に歪みが生じます。
このような現象をエイリアシング（折り返し）、そしてこの歪みをエイリアシング歪みと呼びます。

-4-図 3 における(f) から、重ならないようにするためには、 fm  1   fm すなわち；

1 2 fm.................................（6）となるようにサンプリング周波数 1/τ を選ぶ必要があります
一般に、任意の時間関数 x ( t ) が f m で帯域制限されていてかつ式（6）を満たすサンプリング周波数 1/τ でサンプルするならば、そのサンプル時系列信号 x  t  は式（5）だけで情報を欠落することなく一意的に決定することができます。そしてそのフーリエ変換の結果X ( f ) も正確に求めることができます。
この時、下記の式でサンプルされた数値列 x ( n τ) から、時間連続関数 x ( t ) が求められます。

x t  n  x n   t  n   sin .................................（7） t  n  

この関係を、サンプリング定理と呼びます。
通常 FFT アナライザで周波数分析したい時間信号は、未知の信号でその上限の周波数 f m は不明です。そこで、サンプリング周波数 f s を適当に決めるしかありません。この場合では、分析できる最大の周波数は f s /2 となりこの周波数をナイキスト周波数と呼びます。
もし、信号の上限周波数f m よりナイキスト周波数f s /2 が小さいならば折り返し歪みが生じます。
この歪みが生じないように FFT アナライザでは、ナイキスト周波数 f s /2 で帯域制限するローパスフィルタを装備しています。これを、アンチエイリアシングフィルタ（折り返し防止フィルタ）と呼びます。

, =\-]

NRO=

NQO=

NPO=

PO=

QO=

RO=

-

図 4 エイリアシング現象が生じた例-5-

最後に、まとめです。
(1) アナログ信号をディジタル信号処理するためには、一定時間間隔でディジタル値に変換する必要がありその処理をサンプリング（標本化）、その一定間隔をサンプリング周期（逆数をサンプリング周波数）と呼びます。

(2) 等間隔のインパルス時間列のフーリエ変換は、周波数軸上で同じように等間隔のインパルス列となります。

(3) 周波数帯域 f m で帯域制限されている時間信号を正しくサンプリングするためには、そのサンプリング周波数は、1/2 f m 以上でなければならない。これをサンプリング定理と呼びます。

(4) サンプリング定理を守らないでサンプルすると、エイリアシング歪みが生じて正しくそのスペクトルを求めることができなくなります。

(5) こ の エ イ リ ア シ ン グ 現 象 を 防 止 す る た め に 、FFT ア ナ ラ イ ザ で は サ ン プ リ
ング周波数に連動した折り返し防止フィルタを装備しています。

【キーワード】サンプリング、標本化、サンプリング周期、サンプリング周波数、デルタ関数、インパルス関数、畳み込み、周波数畳み込み定理、サンプリング定理、ナイキスト周波数、アンチエイリアシングフィルタ、折り返し防止フィルタ-----------------------------------------【参考資料】
1.

「高速フーリエ変換」E. ORAN BRIGHAM 著

2.

「ディジタルフーリエ解析（I）－基礎編－」城戸健一

3.

「スペクトル解析」日野幹雄著

科学技術出版社
コロナ社

朝倉書店
以上
（Hima）-6-