離散時間フーリエ変換のソースを表示

{{脚注の不足|date=2022年12月}}
{{読み仮名|'''離散時間フーリエ変換'''|りさんじかんフーリエへんかん|{{lang-en-short|discrete-time Fourier transform}}, '''DTFT'''}}は無限長の[[離散時間信号]]を[[フーリエ変換]]様に[[周波数領域]]へ変換する操作である<ref name=":0">
{{Ilq|離散時間信号 ... の離散時間フーリエ変換 ... <math>X(e^{j\omega})</math> は ... 周波数スペクトルともよばれる．}}{{Harv|越田|2008|p=5}}</ref>。

DTFTの周波数領域の表現は常に周期的関数である。したがって1つの周期に必要な情報が全て含まれるため、DTFTを「有限な」周波数領域への変換であるということもある。

== 定義 ==
<math>\omega \in \R^1</math> を変数とし、無限長である[[離散時間信号]] <math>x[n]</math> の離散時間フーリエ変換 <math>X(\omega)</math> は次式で定義される<ref>{{Ilq|離散時間信号 ... の離散時間フーリエ変換 ... は，次のように定義される．}}{{Harv|越田|2008|p=5}}</ref><ref name=":1">{{Ilq|周波数ωが連続的な変数となっており}}{{Harv|越田|2008|p=5}}</ref>：

:<math>X(\omega) = \sum_{n=-\infty}^{\infty} x[n] \,e^{-i \omega n}</math>
<math>\omega</math> は[[正規化角周波数]]と呼ばれ（⇒ [[#周期性]]）、<math>X(\omega)</math> は[[周波数スペクトル]]とも呼ばれる<ref name=":0" />。

== 表記 ==
DTFT <math>X(\omega)\,</math> は別の表記として <math>X(e^{i \omega})\,</math> とも書かれる。

この表記には次の特徴がある。

* 周期性の強調（⇒ [[#周期性]]）
* DTFT とその元になっている <math>x(t)\,</math> のフーリエ変換 <math>X(f)\,</math>（または <math>X(\omega)\,</math>）との違いを明確化
* DTFT と[[Z変換]]との関係を強調（⇒ [[#Z変換との関係]]）

== 性質 ==

=== 周期性 ===
DTFT <math>X(\omega)</math> は[[周期関数]]である<ref>{{Ilq|離散化された時間関数のフーリエ変換は、周期的な連続周波数関数となる（離散時間フーリエ変換）}}p.6 より引用。{{Cite journal|author=Hima|year=2013|date=2013-03-22|title=基礎からの周波数分析（8）－「離散フーリエ変換（DFT）」|url=https://www.onosokki.co.jp/HP-WK/eMM_back/emm138.pdf|journal=ONO SOKKI -- info channel|issue=138|pages=1-6|publisher=小野測器}}</ref>。

変数 <math>\omega \in \R^1</math> をとる DTFT <math>X(\omega)</math> は周期 <math>2\pi</math> の周期関数である。これは[[複素指数関数]]の周期性と <math>n</math> が整数であることを用いて以下で示される：

:<math>X(\omega + 2\pi)
= \sum_{n=-\infty}^{\infty} x[n] \ e^{-i (\omega + 2\pi) n}
= \sum_{n=-\infty}^{\infty} x[n] \ e^{-i \omega n} e^{-i 2\pi n}
= \sum_{n=-\infty}^{\infty} x[n] \ e^{-i \omega n} \cdot 1
= X(\omega)</math>

よって DTFT の性質を示すには <math>-\pi \leq \omega \leq \pi</math> で十分であり、<math>\omega</math> は[[正規化角周波数]]と解釈される。

この性質は信号 <math>x[n]</math> が[[離散時間信号]]であることと[[複素指数関数]]の周期性に由来する。

=== 連続性 ===
DTFT <math>X(\omega)</math> は連続的である。

DTFTの変数 <math>\omega \in \R^1</math> は連続的な変数である<ref name=":1" />。

この性質は信号 <math>x[n]</math> が時間領域で非周期的であることと[[複素指数関数]]の周期性に由来する。

=== 対称性 ===
フーリエ変換は、実数成分と虚数成分に分離できる。

:<math>X(e^{i \omega}) = X_R(e^{i \omega}) + iX_I(e^{i \omega}) \!</math>

また、偶数成分と奇数成分に分離できる。

:<math>X(e^{i \omega}) = X_E(e^{i \omega}) + X_O(e^{i \omega}) \!</math>

{| class="wikitable"
|-
!時間領域 <br><math>x[n] \!</math>
!周波数領域 <br><math>X(e^{i \omega}) \!</math>
|-
| <math>x^*[n]\!</math>
| <math>X^*(e^{-i \omega}) \!</math>
|-
| <math>x^*[-n]\!</math>
| <math>X^*(e^{i \omega}) \!</math>
|}

=== その他の性質 ===
以下の表は、一般的な離散時間フーリエ変換を示したものである。以下のような記法を用いている。
* <math>*\!</math> は、2つの信号の[[畳み込み]]を意味する。
* <math>x[n]^*\!</math> は、関数 ''x[n]'' の[[共役複素数|複素共役]]である。
* <math>\rho_{xy} [n]\!</math> は、 ''x[n]'' と ''y[n]'' の[[相関]]を表す。

最初の列は属性の説明、第二列は時間領域での関数表現、第三列は周波数領域でのスペクトル表現である。

{| class="wikitable"
|-
! 特性
! 時間領域 <math>x[n] \!</math>
! 周波数領域  <math>X(\omega) \!</math>
! 備考
|-
| 線形性
| <math>a x[n] + b y[n] \!</math>
| <math> a X(e^{i \omega}) + b Y(e^{i \omega}) \!</math>
|
|-
| 時間におけるシフト
| <math>x[n - k] \!</math>
| <math>X(e^{i \omega}) e^{-i \omega k} \!</math>
| ''k'' は整数
|-
| 周波数におけるシフト（変調）
| <math>x[n]e^{ian} \!</math>
| <math>X(e^{i (\omega-a)}) \!</math>
| ''a'' は実数
|-
| 時間逆転
| <math>x[- n] \!</math>
| <math>X(e^{-i \omega}) \!</math> 
|
|-
| 時間共役
| <math>x[n]^* \!</math>
| <math>X(e^{-i \omega})^* \!</math>
|
|-
| 時間逆転と共役
| <math>x[-n]^* \!</math>
| <math>X(e^{i \omega})^* \!</math>
|
|-
| 周波数における微分
| <math>\frac{n}{i} x[n] \!</math>
| <math>\frac{d X(e^{i \omega})}{d \omega} \!</math>  
|
|-
| 周波数における積分
| <math>\frac{i}{n} x[n] \!</math> 
| <math>\int_{-\pi}^{\omega} X(e^{i \vartheta}) d \vartheta \!</math>
|
|-
| 時間における畳み込み
| <math>x[n] * y[n] \!</math>
| <math>X(e^{i \omega}) \cdot Y(e^{i \omega}) \!</math>
|
|-
| 時間における乗算
| <math>x[n] \cdot y[n] \!</math>
| <math>\frac{1}{2 \pi} X(e^{i \omega}) * Y(e^{i \omega}) \!</math>
|
|-
| 相関
| <math>\rho_{xy} [n] = x[-n]^* * y[n] \!</math>
| <math>R_{xy} (\omega) = X(e^{i \omega})^* \cdot Y(e^{i \omega}) \!</math> 
|
|}

=== 特定条件下での性質 ===

==== 周期信号のDTFT ====
[[離散信号#周期性|周期性]]をもち無限長である[[離散時間信号]]のDTFTは離散的である。

無限長の[[複素関数|複素]]離散時間信号 <math>x[n]</math> を[[離散信号#周期|周期]] <math>N \in \Z^+</math> の周期信号、つまり <math>x[n] = x[n - N]</math> とする。 DTFTの定義式は <math>x[n]</math> と <math>e^{-i \omega n}</math> の内積であるため、両者を同じだけシフトしても値は変わらない。そこで両者を <math>N</math> だけシフトすると周期性により以下が成り立つ：
:<math>X(\omega)
= \sum_{n=-\infty}^{\infty} x[n-N] \,e^{-i \omega (n-N)}
= \sum_{n=-\infty}^{\infty} x[n] \,e^{-i \omega n} e^{i \omega N}
= e^{i \omega N} X(\omega)

</math>
この恒等式により各 <math>\omega</math> で <math>e^{i \omega N} = 1


</math> あるいは <math>X(\omega)=0
</math> となる必要がある。[[複素指数]]が <math>1


</math> になるにはその[[偏角 (数学)|偏角]] <math>\omega N


</math> が <math>2\pi
</math> の整数倍であることが求められるため、上記の恒等式条件により周期信号のDTFTは以下の形になる：
:<math>X(\omega) =
\begin{cases}
\sum_{n=-\infty}^{\infty} x[n] \,e^{-i \omega n}, & \text{if } \omega = \frac{2\pi k}{N}, \quad k \in \Z \\
0,                                                & \text{otherwise}
\end{cases}

</math>
つまり周期性をもち無限長である離散時間信号のDTFTは離散的（かつ周期的）である。

また上記式の複素指数部を <math>g_k(n) = \exp{(-i 2\pi k n / N)}, \, k \in \Z

</math> とし、<math>g_k(n)
</math> を <math>N</math> だけシフトすると、
:<math>g_k(n-N)
= e^{-i \frac{2\pi k}{N} (n-N)}
= e^{-i \frac{2\pi k}{N} n} e^{i 2\pi k}
= e^{-i \frac{2\pi k}{N} n} \cdot 1
= g_k(n)</math>
となり、<math>g_k(n)
</math> も <math>x[n]</math> と同じ周期 <math>N</math> の周期関数と見做せる。そのため「無限長の[[総和]]」が「同じ値をもつ『長さ <math>N</math> の部分和』の無限倍」と見做せる。すなわち周期信号のDTFTは以下の形になる：
:<math>X(\omega) =
\begin{cases}
\infty \cdot \sum_{n=0}^{N-1} x[n] \,e^{-i \omega n}, & \text{if } \omega = \frac{2\pi k}{N}, \quad k \in \Z \\
0,                                                & \text{otherwise}
\end{cases}

</math>
これは[[離散フーリエ変換]]（DFT）としばしば対比される。

== 逆変換 ==
{{読み仮名|'''離散時間フーリエ逆変換'''|りさんじかんフーリエぎゃくへんかん|{{lang-en-short|inverse discrete-time Fourier transform}}, '''IDTFT'''}}は離散時間フーリエ変換表現から[[離散時間信号]]を求める演算である<ref>{{Ilq|離散時間フーリエ変換 ... から信号 ... を求める演算は離散時間フーリエ逆変換とよばれ}}{{Harv|越田|2008|p=5}}</ref>。

離散時間信号 <math>x[n]</math>、[[正規化角周波数]] <math>\omega</math>、DTFT表現 <math>X(\omega)</math> を用いて、IDTFT は次式で定義される<ref>{{Ilq|離散時間フーリエ逆変換 ... 次式で与えられる．}}{{Harv|越田|2008|p=5}}</ref>：

:<math>\begin{align}
x[n] &= \frac{1}{2 \pi}\int_{-\pi}^{\pi} X(\omega)\cdot e^{i \omega n} \, d \omega \\
     &= T \int_{-\frac{1}{2T}}^{\frac{1}{2T}} X_T(f)\cdot e^{i 2 \pi f nT}\, df \\
\end{align}
</math>

積分区間はDTFTの一周期全体であり、これは {x[n]} の標本群がDTFTの[[フーリエ級数]]展開の係数でもあることを示している。無限区間の積分では、この変換が通常の[[フーリエ変換]]の逆変換となり、ディラックのインパルスも復元する。すなわち次のようになる。

:<math>\int_{-\infty}^\infty X_T(f)\cdot e^{i 2 \pi f t}\, df \ = \ x_T(t) \ =\ \sum_{n=-\infty}^{\infty} x[n]\cdot \delta(t - n T)</math>

== 他の変換との関係 ==

=== 他のフーリエ変換との関係 ===
基本的にDTFTは[[フーリエ級数]]の逆であり、後者は継続的だが、周期的入力と離散スペクトルを持っている。これら2つの変換の応用は全く異なる。

DFT と DTFT は、標準の連続フーリエ変換を離散的データに適用しようとして自然に生まれたと見ることもできる。そういった観点では、単に入力形式が異なるだけで、変換そのものは同じである。
* 入力が離散的なら、フーリエ変換は DTFT となる。
* 入力が周期的なら、フーリエ変換はフーリエ級数となる。
* 入力が離散的かつ周期的なら、フーリエ変換は DFT となる。

==== DFTとの関係 ====
{{See also|離散フーリエ変換}}
DTFTの数値的評価では、有限長のシーケンスが明らかに必要とされる。実際、長いシーケンスは矩形[[窓関数]]で修正され、次のようになる。

:<math>X(\omega) = \sum_{n=0}^{L-1} x[n] \,e^{-i \omega n}\,</math>, &nbsp; ここで <math>L\,</math> は修正されたシーケンス長である。

これは、修正前のシーケンスのスペクトルの便利な近似として使われる。これによって解像度が悪くなるが、L を増やすことで改善される。

<math>X(\omega)</math> を (2π) の一周期上に一様に分布する任意の <math>(N)</math> 個の周波数で評価するのが一般的である。

:<math>\omega_k = \frac{2 \pi }{N} k\,</math>, &nbsp; &nbsp; ここで <math>k =  0, 1, \dots, N-1 \,</math>

これにより、次が得られる。

:<math>X[k] = X(\omega_k) = \sum_{n=0}^{L-1} x[n] \,e^{-i 2 \pi \frac{k}{N} n}</math>

<math>N \ge L\,</math> であるとき、次のようにも表せる。

<math>X[k] = \sum_{n=0}^{N-1} x[n] \,e^{-i 2 \pi \frac{k}{N} n}</math>, &nbsp; 何故なら <math>n \ge L\,</math> について <math>x[n] = 0\,</math> と定義するため。

このように変形すると、<math>X[k]\,</math> のシーケンスは[[離散フーリエ変換]]（DFT）となる。<math>N</math> はDTFTを標本化する際の解像度と定義され、<math>L</math> はDTFT自体の固有解像度である。したがって、通常これらはほぼ同じ値である。<math>N > L</math> を選択するのが一般的だが、値がゼロの項を総和に含める理由は、DFTを計算する[[高速フーリエ変換]]アルゴリズムを利用できるためである。そのことを強調する場合、「ゼロパディングDFT」あるいは「内挿DFT」と呼ぶ。しかし、値がゼロの項を使わずに単純に計算しても全く同じDFTが得られる。<math>N < L</math> の場合のDTFTも計算でき、その場合はDFTとは等価ではない。

<math>N > L</math> が一般的であることを示すため、次のシーケンスを考える。

:<math>x[n] = e^{i 2\pi \frac{1}{8} n}</math>, ここで <math>L=64</math>

下に示した2つの図は、ラベルで示される通り、異なるサイズのDFTを図示したものである。どちらの場合も支配的な周波数成分は <math>f = \begin{matrix} \frac{1}{8}\end{matrix} = 0.125\,</math> である。右の図に表れているパターンは、<math>L=64</math> の矩形窓関数の[[スペクトル漏れ]]である。左側の図がこのようになっているのは、右の図のゼロと交差している点と標本化した点が重なっている結果である。これは、有限長シーケンスのDTFTというよりも、無限に続く正弦波のような印象を与える。このような図になる原因は、矩形窓関数の使用と、64個の標本あたり8個という整数個の周期になるような周波数を選択しているためである（<math>\begin{matrix} \frac{1}{8}\end{matrix} = \begin{matrix} \frac{8}{64}\end{matrix}</math>）。

{|
|-
| [[ファイル:no-zeropad.png|none|thumb|350px|L = 64 および N = 64 についての DFT]]
| [[ファイル:zeropad.png|none|thumb|350px|L = 64 および N = 256 についての DFT]]
|}

=== Z変換との関係 ===
DTFT は[[Z変換]]の特殊ケースである。両側Z変換は次のように定義される。

:<math>X(z) = \sum_{n=-\infty}^{\infty} x[n] \,z^{-n}</math>

DTFT は <math>z = e^{i \omega}\,</math> の場合である。このとき <math>|e^{i \omega}| = 1\,</math> なので、これは[[複素数|複素平面]]での[[単位円]]付近での[[Z変換]]の評価である。

=== 標本化変換との関係 ===
名称が暗に示している通り、{x[n]} は連続時間関数 <math>x(t)\,</math> の値（標本）を表している。このときの標本化間隔を <math>T\,</math> としたとき、各標本の採取時刻は <math>t = nT\quad</math> であり、<math>1/T = f_s\,</math> が[[サンプリング周波数]]となる。DTFTは次の連続時間[[フーリエ変換]]の近似である。

:<math> X(f) = \int_{-\infty}^\infty x(t)\cdot e^{- i 2\pi f t}\,dt </math>

[[標本化定理]]で示されるように、次の[[くし型関数]]の変調に <math>x(nT)\,</math> の値を使用すると見ることもできる。

:<math>\Delta_T(t) = T\sum_{n=-\infty}^{\infty} \delta(t - nT) \ </math>

その場合得られる関数のフーリエ変換は、<math>f_s \,</math> の間隔で重ね合わせられた <math>X(f)\,</math> のコピーの総和である。

:<math>X_\mathrm{T}(f) = \sum_{k = -\infty}^{\infty} X(f - {k f_s}) </math>

以下で示すように、これは周期関数のDTFTである。そして、ある明白な条件下で、k=0 の項はほとんど全く他の項からの歪み（[[折り返し雑音]]）が観測されない。変調されたくし型関数は次の通りである。

:<math>x_\mathrm{T}(t) = T \sum_{n=-\infty}^{\infty} x(nT)\, \delta(t - n T)</math>

したがって、

:<math> \begin{align}
X_\mathrm{T}(f)
& = \int_{-\infty}^\infty \left[T \sum_{n=-\infty}^{\infty} x(nT)\, \delta(t - n T)\right] e^{- i 2\pi f t}\,dt \\
& = \sum_{n=-\infty}^{\infty} T\cdot x(nT) \int_{-\infty}^\infty \left[\delta(t - n T)\cdot e^{- i 2\pi f t}\right]\,dt \\
& = \sum_{n=-\infty}^{\infty} T\cdot x(nT) \cdot e^{- i 2\pi f n T}
\end{align}
</math>

このとき次が成り立つ。

:<math>x[n] = T \cdot x(nT)\,</math>
:<math>\omega = 2\pi f T = 2 \pi \frac{f}{f_s}</math>

つまり <math>X_\mathrm{T}(f)\,</math> は <math>X(\omega)\,</math> と同じである。

ここで、<math>f\,</math> は通常の周波数（単位時間当たりの周期数）であり、<math>f_s\,</math> はサンプリング周波数（単位時間当たりの標本数）であるから、<math>f / f_s\,</math> は「標本当たりの周期数」を意味する。これを正規化周波数（normalized frequency）と呼ぶ。上で定義されている <math>\omega\,</math> も正規化周波数だが、こちらの単位は「標本当たりのラジアン」である。正規化周波数は、期間 <math>2\pi</math> の周期を持つ関数 <math>X(\omega)</math> で表されるという特徴がある。そのため、逆変換では <math>2\pi</math> の期間のみを評価すればよい。

== 例 ==
下表は典型的な変換を示したものである。
* <math>n \!</math> は離散時間領域（標本）を表現する整数である。
* <math>\omega \!</math> は <math>(-\pi,\ \pi)</math> の範囲内の実数であり、連続角周波数（標本当たりのラジアン）を表す。
** それ以外 <math>(|\omega| > \pi \,)</math> の変換は、<math>X(\omega + 2\pi k) = X(\omega)\,</math> で定義される。
* <math>u[n] \!</math> は離散時間[[ヘヴィサイドの階段関数|単位ステップ関数]]である。
* <math>\operatorname{sinc}(t) \!</math> は正規化[[Sinc関数]]である。
* <math>\delta (\omega) \!</math> は[[ディラックのデルタ関数]]である。
* <math>\delta [n] \!</math> は[[クロネッカーのデルタ]] <math>\delta_{n,0} \!</math> である。
* <math> \operatorname{rect}(t) </math> は、任意の実数値 ''t'' に関する次のような[[矩形関数]]である。
**<math>\mathrm{rect}(t) = \sqcap(t) = \begin{cases}
0           & \mbox{if } |t| > \frac{1}{2} \\[3pt]
\frac{1}{2} & \mbox{if } |t| = \frac{1}{2} \\[3pt]
1           & \mbox{if } |t| < \frac{1}{2}
\end{cases} </math>
* <math>\operatorname{tri}(t) </math> は任意の実数値 ''t'' に関する次のような[[三角形関数]]である。
**<math>\operatorname{tri}(t) = \land (t) = 
\begin{cases}
1 + t; & - 1 \leq t \leq 0 \\
1 - t; & 0 < t \leq 1 \\
0 & \mbox{otherwise} 
\end{cases} </math>

{| class="wikitable"
|-
! 時間領域 <br /> <math> x[n] \, </math>
! 周波数領域 <br /><math>  X(\omega) \, </math>
! 備考
|-
| <math>\delta [n] \!</math>
| <math>1 \!</math>
| 
|-
| <math>\delta [n - M] \!</math>
| <math>e^{-i \omega M} \!</math>
| ''M'' は整数
|-
|<math>\sum_{m = -\infty}^{\infty} \delta[n - M m] \,</math>
|<math>\sum_{m = -\infty}^{\infty} e^{-i \omega M m} = \frac{1}{M}\sum_{k = -\infty}^{\infty} \delta \left( \frac{\omega}{2\pi} - \frac{k}{M} \right) \,</math>
| ''M'' は整数
|-
| <math>u[n]\!</math>
| <math>\frac{1}{1-e^{-i \omega}} \!</math>
|
|-
| <math>e^{-ian} \!</math>
| <math> 2\pi \delta (\omega + a)  \, </math>
| ''a'' は実数
|-
| <math>\cos (a n) \!</math>
| <math>\pi \left[ \delta (\omega - a) + \delta (\omega + a) \right]</math>
| ''a'' は実数
|-
| <math>\sin (a n) \!</math>
| <math>\frac{\pi}{i} \left[ \delta (\omega - a) - \delta ( \omega + a) \right]</math>
| ''a'' は実数
|-
| <math> \mathrm{rect} \left[ { ( n - M/2 ) \over M  } \right] </math>
| <math> { \sin[ \omega (M+1) / 2 ] \over \sin( \omega / 2 ) } \,  e^{ -i \omega M / 2 }</math>
| ''M'' は整数
|-
| <math>\operatorname{sinc} [(a + n)]</math>
| <math>e^{i a \omega} \!</math>
| ''a'' は実数
|-
| <math>W\cdot \operatorname{sinc}^2(W n)\,</math>
| <math>\operatorname{tri} \left( { \omega \over 2\pi W } \right)</math>
| real number ''W''<br /><math>0 < W \le 0.5</math>
|-
| <math>W\cdot \operatorname{sinc} [ W (n + a)]</math>
| <math>\operatorname{rect} \left( { \omega \over 2\pi W } \right) \cdot e^{j a \omega}</math>
| ''W'', ''a'' は実数<br /><math>0 < W \le 1</math>
|-
| <math> 
\begin{cases}
0 & n=0 \\
\frac{(-1)^n}{n} & \mbox{elsewhere}
\end{cases}
</math>
| <math>j \omega</math>
| [[微分回路]]フィルタとして機能する
|-
| <math>\frac{W}{(n + a)} \left\{ \cos [ \pi W (n+a)] - \operatorname{sinc} [ W (n+a)] \right\}</math>
| <math>j \omega \cdot \operatorname{rect} \left( { \omega \over \pi W } \right) e^{j a \omega}</math>
| ''W'',''a'' は実数<br /><math>0 < W \le 1</math>
|-
| <math>\frac{1}{\pi n^2} [(-1)^n - 1]</math>
| <math>| \omega | \!</math>
|
|-
| <math>
\begin{cases}
0; & n \mbox{ odd} \\
\frac{2}{\pi n} ; & n \mbox{ even}
\end{cases}
</math>
| <math>
\begin{cases}
j & \omega < 0 \\
0 & \omega = 0 \\
-j & \omega > 0
\end{cases}
</math>
|[[ヒルベルト変換]]
|-
| <math>\frac{C (A + B)}{2 \pi} \cdot \operatorname{sinc} \left[ \frac{A - B}{2\pi} n \right] \cdot \operatorname{sinc} \left[ \frac{A + B}{2\pi} n \right]</math>
| [[画像:Trapezoid_signal.png|250px]]
| ''A'', ''B'' は実数 <br /> ''C'' は複素数
|}

== 脚注 ==
{{脚注ヘルプ}}
=== 注釈 ===
{{notelist2}}
=== 出典 ===
{{reflist}}

== 参考文献 ==
* {{cite book | author=Alan V. Oppenheim and Ronald W. Schafer | title = Discrete-Time Signal Processing | edition = 2nd Edition | publisher = Prentice Hall Signal Processing Series | date = 1999年 | id = ISBN 0-13-754920-2}}
* {{cite book | author = William McC. Siebert | title = Circuits, Signals, and Systems | date = 1986年 | publisher = [[MIT Press]] | location = MIT Electrical Engineering and Computer Science Series.  Cambridge, MA }}
* {{cite book | author = Boaz Porat | title = A Course in Digital Signal Processing | publisher = John Wiley and Sons | id = ISBN 0-471-14961-6 | pages = pp. 27-29 and 104-105 }}
* {{Citation|和書|title=知識の森|last=越田|first=俊介|year=2008|url=https://www.ieice-hbkb.org/files/ad_base/view_pdf.html?p=/files/01/01gun_09hen_01m.pdf|publisher=電子情報通信学会|pages=1-10|chapter=1群-9編-1章 ディジタル信号処理の基礎理論}}

== 関連項目 ==
* [[フーリエ変換]]
* [[離散信号]]

{{デジタル信号処理}}

{{DEFAULTSORT:りさんしかんふりえへんかん}}
[[Category:フーリエ解析]]
[[Category:信号処理]]
[[Category:ジョゼフ・フーリエ]]
[[Category:数学に関する記事]]