FFT

入力の高速フーリエ変換 (FFT)

ライブラリ:
DSP System Toolbox / Transforms

説明

FFT ブロックは、N 次元の入力配列 u の最初の次元全体において高速フーリエ変換 (FFT) を計算します。ブロックは 2 つの可能な FFT 実装のどちらかを使用します。FFTW ライブラリを基にした実装か、基数 2 のアルゴリズムのコレクションを基にした実装を選択できます。ブロックで実装を選択できるように、[自動] を選択できます。FFT 実装の詳細については、アルゴリズムを参照してください。

ユーザー指定の FFT 長 (P と等しくない) では、ゼロパディングまたは切り捨て、あるいは長さを法とするデータラッピングが FFT 演算の前に発生します。P ≤ M の FFT の場合、次のようになります。

y = fft(u,M) % P ≤ M

ラッピング:

y(:,L) = fft(datawrap(u(:,L),M)) % P > M; L = 1,...,N

切り捨て:

y (:,L) = fft(u,M) % P > M; L = 1,...,N

ヒント

入力の長さ P が FFT 長 M より大きい場合、FFT 出力での振幅が増加することがあります。振幅が増加するのは、FFT ブロックが M を法とするデータラッピングを使用してすべての使用可能な入力サンプルを保持するためです。

このような振幅の増加を避けるため、入力サンプル P の長さを FFT 長 M まで切り捨てることができます。これを行うには、Pad ブロックをモデル内の FFT ブロックの前に配置します。

例

時間領域データから周波数領域への変換

この例では、FFT ブロックを使用して時間領域データを周波数領域に変換する方法について説明します。

スクリプトを開く

端子

入力

すべて展開する

Port_1 — 入力信号
ベクトル | 行列 | N 次元配列

FFT を計算する入力信号。ブロックは N 次元の入力信号の最初の次元に沿って FFT を計算します。

ブロックが FFT を計算する方法の詳細については、説明およびアルゴリズムを参照してください。

出力

すべて展開する

Port_1 — 入力の FFT
ベクトル | 行列 | N 次元配列

FFT。N 次元の入力配列の最初の次元全体で計算されます。ブロックの出力に整数または小数点データ型が含まれている場合、常に符号付きです。

L 番目の出力チャネル y(k,L) の k 番目のエントリは、L 番目の入力チャネルにおける M 点離散フーリエ変換 (DFT) の k 番目の点と等価です。

$y (k, L) = \sum_{p = 1}^{P} u (p, L) e^{- j 2 π (p - 1) (k - 1) / M} k = 1, \dots, M$

ブロックが FFT を計算する方法の詳細については、説明およびアルゴリズムを参照してください。

データ型: single | double | int8 | int16 | int32 | fixed point
複素数のサポート: あり

パラメーター

すべて展開する

メイン

FFT の実装 — FFT の実装
`自動` (既定値) | `Radix-2` | `FFTW`

このパラメーターを [FFTW] に設定して、任意の長さの入力信号をサポートします。ブロックは FFTW 実装で生成されたコードを MATLAB^® を実行可能なホストコンピューターのみに制限します。

ビット反転した処理、固定小数点データまたは浮動小数点データ、または Simulink^® Coder™ を使った移植可能な C コード生成では、このパラメーターを [Radix-2] に設定します。M 行 N 列の入力行列の次元 M は、2 のべき乗でなければなりません。他の入力サイズで操作するには、Pad ブロックを使ってこれらの次元を 2 のべき乗にパディングまたは切り捨てをするか、可能な場合は FFTW 実装を選択します。[Radix-2] モードで使用されるアルゴリズムの詳細については、基数 2 の実装を参照してください。

このパラメーターを [自動] に設定して、ブロックが FFT 実装を選択できるようにします。2 のべき乗でない変換の長さの浮動小数点の場合、FFTW アルゴリズムは自動的に選択されます。それ以外の場合は、Radix-2 アルゴリズムが自動的に選択されます。2 のべき乗でない変換の長さの場合、ブロックは生成したコードを MATLAB のホストコンピューターに制限します。

ビット反転順での出力 — ビット反転順での出力
`off` (既定値) | `on`

入力要素の順序に対する出力チャネル要素の順序を指定します。このチェックボックスをオンにすると、出力チャネル要素が入力順序に対してビット反転した順序で表示されます。このチェックボックスをオフにすると、出力チャネル要素が入力順序に対して線形順に表示されます。

メモ

FFT ブロックは出力をビット反転した順序で計算します。FFT ブロックの出力を線形順にするには、追加のビット反転演算が必要になります。多くの状況では [ビット反転順での出力] チェックボックスをオンにすると FFT ブロックの速度を上げることができます。

出力の順序の詳細については、Linear and Bit-Reversed Output Orderを参照してください。

依存関係

このパラメーターを有効にするには、[FFT 実装] を [自動] または [Radix-2] に設定します。

出力を FFT 長で除算 — 出力を FFT 長で除算
`off` (既定値) | `on`

このパラメーターを選択すると、ブロックは FFT 長によって FFT の出力を分割します。このオプションは、FFT の出力をその入力と同じ振幅範囲内に留める必要がある場合に役立ちます。これは特に、固定小数点データ型での作業に役立ちます。

FFT 長を入力次元から継承 — FFT 長を入力次元から継承
`on` (既定値) | `off`

入力次元から FFT 長を継承する場合に選択します。このチェックボックスをオンにする場合、入力の長さは 2 のべき乗でなければなりません。

依存関係

このチェックボックスをオフにする場合、長さを指定する [FFT 長] パラメーターが使用できるようになります。

FFT 長 — FFT 長
`64` (既定値) | 整数

FFT 長を 2 以上の整数で指定します。

[FFT 実装] パラメーターを Radix-2 に設定する場合または [ビット反転順での出力] チェックボックスをオンにする場合、この値は 2 のべき乗でなければなりません。

依存関係

このパラメーターを有効にするには、[FFT 長を入力次元から継承] チェックボックスをオフにします。

FFT 長が入力長よりも短い場合は入力データをラップ — 入力のラップまたは切り捨て
`on` (既定値) | `off`

FFT 長に応じて入力をラップするか、切り捨てるかを選択します。このパラメーターをオンにすると、FFT 長が入力の長さよりも短い場合、FFT 演算の前に長さを法とするデータラッピングが発生します。このチェックボックスをオフにすると、FFT 演算の前に入力データが FFT 長に切り捨てられます。

依存関係

このパラメーターを有効にするには、[FFT 長を入力次元から継承] チェックボックスをオフにします。

データ型

丸めモード — 丸め方法
`負方向` (既定値) | `正方向` | `最も近い偶数方向` | `最も近い正の整数方向` | `最も近い整数方向` | `最も簡潔` | `ゼロ方向`

固定小数点演算の丸めモードを選択します。

制限

正弦表の値はこのパラメーターに従わず、常に [最も近い正の整数方向] に丸められます。

[丸めモード] パラメーターは、以下のすべての条件に一致する場合、数値結果には影響しません。

[乗算出力] のデータ型が次である。[継承: 内部ルールによる継承]。
[アキュムレータ] のデータ型が次である。[継承: 内部ルールによる継承]。

これらのデータ型設定を使用すると、ブロックは完全精度モードで動作します。

整数オーバーフローで飽和 — 整数オーバーフローで飽和
`off` (既定値) | `on`

このパラメーターをオンにすると、ブロックは固定小数点演算の結果を飽和させます。このパラメーターをオフにすると、ブロックは固定小数点演算の結果をラップします。saturate および wrap の詳細については、固定小数点演算のオーバーフローモードを参照してください。

制限

[整数オーバーフローで飽和] のパラメーターは、以下のすべての条件に一致する場合、数値結果には影響しません。

[乗算出力] のデータ型が次である。[継承: 内部ルールによる継承]。
[アキュムレータ] のデータ型が次である。[継承: 内部ルールによる継承]。

これらのデータ型設定を使用すると、ブロックは完全精度モードで動作します。

正弦表 — 正弦表の値のデータ型
`継承: 入力と同じ語長` (既定値) | `fixdt(1,16)`

正弦表の値の語長を指定する方法を選択します。正弦表の値の小数部の長さは常に語長から 1 を引いた値に等しくなります。このパラメーターは以下のように設定できます。

データ型継承ルール。例: [継承: 入力と同じ語長]
有効なデータ型として評価する式。例: fixdt(1,16)

[データ型アシスタントを表示] ボタンをクリックして、[正弦表] パラメーターの設定を行うための [データ型アシスタント] を表示します。

詳細については、データ型アシスタントを利用したデータ型の指定 (Simulink)を参照してください。

制限

正弦表の値は [丸めモード] パラメーターおよび [整数オーバーフローで飽和] パラメーターには従わず、常に飽和して [最も近い正の整数方向] に丸められます。

乗算出力 — 乗算出力のデータ型
`継承: 内部ルールによる継承` (既定値) | `継承: 入力と同じ` | `fixdt(1,16,0)`

乗算出力のデータ型を指定します。このブロックの乗算出力データ型の使い方を示す図は、固定小数点および乗算のデータ型を参照してください。このパラメーターは以下のように設定できます。

データ型継承ルール。例: [継承: 内部ルールによる継承]。このルールの詳細については、内部ルールによる継承を参照してください。
有効なデータ型として評価する式。例: fixdt(1,16,0)

[データ型アシスタントを表示] ボタンをクリックして、[乗算出力] パラメーターの設定を行うための [データ型アシスタント] を表示します。

詳細については、データ型アシスタントを利用したデータ型の指定 (Simulink)を参照してください。

アキュムレータ — アキュムレータのデータ型
`継承: 内部ルールによる継承` (既定値) | `継承: 入力と同じ` | `継承: 乗算出力と同じ` | `fixdt(1,16,0)`

アキュムレータのデータ型を指定します。このブロックのアキュムレータデータ型の使い方を示す図は、固定小数点を参照してください。このパラメーターは以下のように設定できます。

データ型継承ルール。例: [継承: 内部ルールによる継承]。このルールの詳細については、内部ルールによる継承を参照してください。
有効なデータ型として評価する式。例: fixdt(1,16,0)

[データ型アシスタントを表示] ボタンをクリックして、[アキュムレータ] パラメーターの設定を行うための [データ型アシスタント] を表示します。

詳細については、データ型アシスタントを利用したデータ型の指定 (Simulink)を参照してください。

出力 — 出力データ型
`継承: 内部ルールによる継承` (既定値) | `継承: 入力と同じ` | `fixdt(1,16,0)`

出力データ型を指定します。このブロックの出力データ型の使い方を示す図は、固定小数点を参照してください。このパラメーターは以下のように設定できます。

データ型継承ルール。例: [継承: 内部ルールによる継承]。
[継承: 内部ルールによる継承] を選択すると、ブロックは出力の語長と小数部の長さを自動的に計算します。理想的な出力の語長と小数部の長さを計算するためにブロックが使用する方程式は、[出力を FFT 長で除算] チェックボックスの設定によって異なります。
- [出力を FFT 長で除算] チェックボックスをオンにした場合、理想的な出力の語長と小数部の長さは入力の語長と小数部の長さと同じになります。
- [出力を FFT 長で除算] チェックボックスをオフにした場合、ブロックは次の方程式に従って理想的な出力の語長と小数部の長さを計算します。
  $W L_{i d e a l o u t p u t} = W L_{i n p u t} + floor (\log_{2} (F F T l e n g t h - 1)) + 1$
  $F L_{i d e a l o u t p u t} = F L_{i n p u t}$
これらの理想的な結果を使用して、内部ルールによって、使用しているハードウェアに適した語長と小数部の長さが選択されます。詳細については、内部ルールによる継承を参照してください。
有効なデータ型として評価する式。例: fixdt(1,16,0)

[データ型アシスタントを表示] ボタンをクリックして、[出力] パラメーターの設定を行うための [データ型アシスタント] を表示します。

詳細については、信号のデータ型の制御 (Simulink)を参照してください。

出力の [最小値] — ブロックが出力しなければならない最小値
`[]` (既定値) | スカラー

ブロックが出力する最小値を指定します。既定値は [] (指定なし) です。Simulink ソフトウェアは、以下を行う際にこの値を使用します。

シミュレーション範囲のチェック (信号範囲の指定 (Simulink)を参照)
固定小数点データ型の自動スケーリング

出力の [最大値] — ブロックが出力しなければならない最大値
`[]` (既定値) | スカラー

ブロックが出力する最大値を指定します。既定値は [] (指定なし) です。Simulink ソフトウェアは、以下を行う際にこの値を使用します。

シミュレーション範囲のチェック (信号範囲の指定 (Simulink)を参照)
固定小数点データ型の自動スケーリング

固定小数点ツールによる変更に対してデータ型の設定をロックする — 固定小数点ツールによるデータ型のオーバーライドの回避
`off` (既定値) | `on`

ブロックダイアログボックスで指定するデータ型が固定小数点ツールによってオーバーライドされないようにするには、このパラメーターを選択します。

ブロックの特性

データ型	`double` \| `fixed point` \| `integer` \| `single`
直達	`いいえ`
多次元信号	`はい`
可変サイズの信号	`はい^a`
ゼロクロッシング検出	`いいえ`
^a 可変サイズの信号は、[Inherit FFT length from input dimensions] チェックボックスが選択されている場合にのみサポートされます。

アルゴリズム

すべて展開する

FFTW 実装

FFTW 実装は、シミュレーションとコード生成の両方で、2 のべき乗および 2 のべき乗でない変換の長さのサポートを含む、最適化された FFT 計算を提供します。FFTW 実装を使用して生成されたコードは MATLAB を実行可能なコンピューターでのみ実行できます。入力データ型は浮動小数点でなければなりません。

基数 2 の実装

基数 2 の実装は、ビット反転した処理や固定小数点データまたは浮動小数点データをサポートし、ブロックで Simulink Coder を使用した移植可能な C コードを生成できるようにします。M 行 N 列の入力行列の次元 M は、2 のべき乗でなければなりません。他の入力サイズで操作するには、Pad ブロックを使ってこれらの次元を 2 のべき乗にパディングまたは切り捨てをするか、可能な場合は FFTW 実装を選択します。

基数 2 を選択すると、ブロックは次のアルゴリズムの 1 つ以上を実装します。

バタフライ演算
double 信号のアルゴリズム
半分長のアルゴリズム
基数 2 の時間間引き (DIT) アルゴリズム
基数 2 の周波数間引き (DIF) アルゴリズム

実数信号と複素数信号に対する基数 2 のアルゴリズム

実数/複素数の入力	出力順序	FFT 計算に使用するアルゴリズム
複素数	線形	ビット反転した演算と基数 2 の DIT
複素数	ビット反転	基数 2 の DIF
実数	線形	半分長および double 信号のアルゴリズムを併用した、ビット反転した演算と基数 2 の DIT
実数	ビット反転	半分長および double 信号のアルゴリズムを併用した基数 2 の DIF

DFT を計算する前に実数値のシーケンスから複素数値のシーケンスを形成することにより、実数入力信号での FFT アルゴリズムの効率性が向上する可能性があります。2N+1 の実数入力チャネルがある場合、FFT ブロックは double 信号のアルゴリズムを最初の 2N 入力チャネルに、半分長のアルゴリズムを最後の奇数チャネルに適用して、これらの複素数値のシーケンスを形成します。

固定小数点データ型の実数入力信号では、すべての入力チャネルが同一であっても、最後の奇数チャネルの出力で異なる数値結果が表示されることがあります。この数値の相違は、double 信号のアルゴリズムと半分長のアルゴリズムの相違に起因します。

この数値の相違は、次の 2 つの方法で除去することができます。

固定小数点入力信号に完全精度の演算を使用する
入力データ型を浮動小数点に変更する

double 信号のアルゴリズムの詳細については、[2]、475 ページの「Efficient Computation of the DFT of Two Real Sequences」を参照してください。半分長のアルゴリズムの詳細については、[2]、476 ページの「Efficient Computation of the DFT of a 2N-Point Real Sequence」を参照してください。

三角関数値のテーブルに対する基数 2 の最適化

特定の状況では、ブロックの基数 2 のアルゴリズムは回転因子のすべての可能な三角関数値を計算します。

$e^{j \frac{2 π k}{K}}$

ここで、K は M か N のいずれか大きい方の値であり、また $k = 0, \dots, K - 1$ です。ブロックはこれらの値をテーブルに格納し、シミュレーション中に取得します。固定小数点と浮動小数点のテーブルエントリ数は、次の表にまとめられています。

N 点 FFT のテーブルエントリ数
浮動小数点	3N/4
固定小数点	N

参照

[1] Orfanidis, S. J. Introduction to Signal Processing. Upper Saddle River, NJ: Prentice Hall, 1996, p. 497.

[2] Proakis, John G. and Dimitris G. Manolakis. Digital Signal Processing, 3rd ed. Upper Saddle River, NJ: Prentice Hall, 1996.

[3] FFTW (https://www.fftw.org)

[4] Frigo, M. and S. G. Johnson, “FFTW: An Adaptive Software Architecture for the FFT,”Proceedings of the International Conference on Acoustics, Speech, and Signal Processing, Vol. 3, 1998, pp. 1381-1384.

拡張機能

C/C++ コード生成
Simulink® Coder™ を使用して C および C++ コードを生成します。

使用上の注意および制限:

次の条件に該当する場合、このブロックから生成される実行可能ファイルは、MATLAB に付属のプリビルドされたダイナミックライブラリファイル (.dll ファイル) に依存します。
- [FFT の実装] が [FFTW] に設定されている。
- [FFT 長を入力次元から継承] がオフで、[FFT 長] が 2 のべき乗でない値に設定されている。
このブロックから生成されたコードとすべての関連ファイルを zip 圧縮ファイルにパッケージ化するには、関数 packNGo を使用します。この zip ファイルを使用して、MATLAB がインストールされていない他の開発環境にプロジェクトを移動して解凍し、リビルドすることができます。詳細については、How To Run a Generated Executable Outside MATLABを参照してください。
FFT 長が 2 のべき乗の場合は、このブロックからスタンドアロンの C コードおよび C++ コードを生成できます。

固定小数点の変換
Fixed-Point Designer™ を使用して固定小数点システムの設計とシミュレーションを行います。

次の図は、FFT ブロックで固定小数点信号に使用されるデータ型を示します。パラメーターで説明されているように、FFT ダイアログボックスの図に表示される [正弦表]、[アキュムレータ]、[乗算出力] および [出力] データ型を設定できます。

FFT ブロックへの入力は、まず出力データ型にキャストされてから出力バッファーに格納されます。各バタフライ段階の後、アキュムレータデータ型で信号が処理され、バタフライの最終出力が出力データ型にキャストし直されます。ブロックは、時間間引き FFT での各バタフライ段階の前と、周波数間引き FFT での各バタフライ段階の後で、回転因数で乗算されます。

乗算器への入力は両方とも複素数であるため、乗算器の出力がアキュムレータデータ型で表示されます。実行される複素数の乗算の詳細については、乗算のデータ型を参照してください。

メモ

ブロックの入力が固定小数点の場合、すべての内部データ型は符号付き固定小数点になります。

バージョン履歴

R2006a より前に導入

参考

FFT

説明

例

時間領域データから周波数領域への変換

端子

入力

Port_1 — 入力信号 ベクトル | 行列 | N 次元配列

出力

Port_1 — 入力の FFT ベクトル | 行列 | N 次元配列

パラメーター

メイン

FFT の実装 — FFT の実装 自動 (既定値) | Radix-2 | FFTW

ビット反転順での出力 — ビット反転順での出力 off (既定値) | on

依存関係

出力を FFT 長で除算 — 出力を FFT 長で除算 off (既定値) | on

FFT 長を入力次元から継承 — FFT 長を入力次元から継承 on (既定値) | off

依存関係

FFT 長 — FFT 長 64 (既定値) | 整数

依存関係

FFT 長が入力長よりも短い場合は入力データをラップ — 入力のラップまたは切り捨て on (既定値) | off

依存関係

データ型

丸めモード — 丸め方法 負方向 (既定値) | 正方向 | 最も近い偶数方向 | 最も近い正の整数方向 | 最も近い整数方向 | 最も簡潔 | ゼロ方向

制限

整数オーバーフローで飽和 — 整数オーバーフローで飽和 off (既定値) | on

制限

正弦表 — 正弦表の値のデータ型 継承: 入力と同じ語長 (既定値) | fixdt(1,16)

制限

乗算出力 — 乗算出力のデータ型 継承: 内部ルールによる継承 (既定値) | 継承: 入力と同じ | fixdt(1,16,0)

アキュムレータ — アキュムレータのデータ型 継承: 内部ルールによる継承 (既定値) | 継承: 入力と同じ | 継承: 乗算出力と同じ | fixdt(1,16,0)

出力 — 出力データ型 継承: 内部ルールによる継承 (既定値) | 継承: 入力と同じ | fixdt(1,16,0)

出力の [最小値] — ブロックが出力しなければならない最小値 [] (既定値) | スカラー

出力の [最大値] — ブロックが出力しなければならない最大値 [] (既定値) | スカラー

固定小数点ツールによる変更に対してデータ型の設定をロックする — 固定小数点ツールによるデータ型のオーバーライドの回避 off (既定値) | on

ブロックの特性

アルゴリズム

FFTW 実装

基数 2 の実装

参照

拡張機能

C/C++ コード生成 Simulink® Coder™ を使用して C および C++ コードを生成します。

固定小数点の変換 Fixed-Point Designer™ を使用して固定小数点システムの設計とシミュレーションを行います。

バージョン履歴

参考

オブジェクト

関数

ブロック

トピック

Port_1 — 入力信号
ベクトル | 行列 | N 次元配列

Port_1 — 入力の FFT
ベクトル | 行列 | N 次元配列

FFT の実装 — FFT の実装
`自動` (既定値) | `Radix-2` | `FFTW`

ビット反転順での出力 — ビット反転順での出力
`off` (既定値) | `on`

出力を FFT 長で除算 — 出力を FFT 長で除算
`off` (既定値) | `on`

FFT 長を入力次元から継承 — FFT 長を入力次元から継承
`on` (既定値) | `off`

FFT 長 — FFT 長
`64` (既定値) | 整数

FFT 長が入力長よりも短い場合は入力データをラップ — 入力のラップまたは切り捨て
`on` (既定値) | `off`

丸めモード — 丸め方法
`負方向` (既定値) | `正方向` | `最も近い偶数方向` | `最も近い正の整数方向` | `最も近い整数方向` | `最も簡潔` | `ゼロ方向`

整数オーバーフローで飽和 — 整数オーバーフローで飽和
`off` (既定値) | `on`

正弦表 — 正弦表の値のデータ型
`継承: 入力と同じ語長` (既定値) | `fixdt(1,16)`

乗算出力 — 乗算出力のデータ型
`継承: 内部ルールによる継承` (既定値) | `継承: 入力と同じ` | `fixdt(1,16,0)`

アキュムレータ — アキュムレータのデータ型
`継承: 内部ルールによる継承` (既定値) | `継承: 入力と同じ` | `継承: 乗算出力と同じ` | `fixdt(1,16,0)`

出力 — 出力データ型
`継承: 内部ルールによる継承` (既定値) | `継承: 入力と同じ` | `fixdt(1,16,0)`

出力の [最小値] — ブロックが出力しなければならない最小値
`[]` (既定値) | スカラー

出力の [最大値] — ブロックが出力しなければならない最大値
`[]` (既定値) | スカラー

固定小数点ツールによる変更に対してデータ型の設定をロックする — 固定小数点ツールによるデータ型のオーバーライドの回避
`off` (既定値) | `on`

C/C++ コード生成
Simulink® Coder™ を使用して C および C++ コードを生成します。

固定小数点の変換
Fixed-Point Designer™ を使用して固定小数点システムの設計とシミュレーションを行います。