fft

iddata オブジェクトの高速フーリエ変換 (FFT)

構文

datf = fft(data)

datf = fft(data,N)

datf = fft(data,N,'complex')

説明

datf = fft(data) は高速フーリエ変換 (FFT) アルゴリズムを使用して、時間領域 iddata オブジェクトデータを周波数領域 iddata オブジェクト datf に変換します。data は実数値信号を含み、定数のサンプル時間 Ts をもちます。datf は、0 からナイキスト周波数までの均等に分布した周波数値をもちます。ナイキスト周波数は、信号長が偶数の場合は含まれますが、信号長が奇数の場合は含まれません。信号強度とノイズレベルを保持するために、FFT は各変換を信号長の平方根で除算して正規化されます。

例

datf = fft(data,N) は、変換の長さ N を指定します。既定では、変換の長さは信号長によって決まります。N を含めると、FFT の変換において、data の信号の方が N より短い場合はゼロでパディングされ、それ以外の場合は信号が切り捨てられます。したがって、実数信号の場合の周波数の数は、N が偶数の場合は (N/2)+1、N が奇数の場合は (N+1)/2 になります。data に複数の実験が含まれる場合は、N を対応する長さの行ベクトルにすることができます。

例

datf = fft(data,N,'complex') は、負の周波数を含めるように指定します。実数信号の場合、既定では datf に非負の周波数のみが含まれますが、複素数値信号の場合は常に負の周波数が含まれます。実数の場合に負の周波数を強制するには、'complex' を追加します。

例

すべて折りたたむ

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

ライブスクリプトを開く

300 のサンプルが含まれる時間領域データ z1 を読み込んでプロットします。

load iddata1 z1;
plot(z1)

データを周波数領域に変換してプロットします。

datf = fft(z1);
plot(datf)

変換の長さを指定したデータの変換

ライブスクリプトを開く

データを読み込みます。fft を使用し、変換の長さ N を 100 にして時間領域データ z1 を周波数領域に変換します。

load iddata1 z1;
datf = fft(z1, 100);
plot(datf)

負の周波数をもつ実数データの変換

ライブスクリプトを開く

データを読み込みます。負の周波数を含める zf2 を作成するときに 'complex' を使用して、2 周波数領域オブジェクト zf1 および zf2 を作成します。

結果を比較します。

load iddata1 z1
zf1 = fft(z1,300);
zf2 = fft(z1,300,'complex');
h = plot(zf2,zf1);
legend('zf2','zf1')
opt = getoptions(h);
opt.FreqScale = 'linear';
opt.PhaseMatching = 'on';
setoptions(h,opt)
xlim([-32 32])