fft_row

【目的】

fft_row - 行毎の高速フーリエ変換

【形式】

Y = fft_row(X)               Y = fft_row(X,k)
   (CoMatrix|CoArray) Y;        (CoMatrix|CoArray) Y;
   (Matrix|Array) X;            (Matrix|Array) X;
                                Integer k;

【詳細】

fft_row(X)は，2 のべき乗の長さのベクトル X について，基底が 2 の 高速フーリエ変換を計算する。fft_row(X,k) は k 点の高速フーリエ変換 をする。X の長さが k より短いとき，ゼロが後ろに付け加えられ， X の長さが k より長いとき，k 番目以降が切り捨てられる。 X が行列のとき，行毎に高速フーリエ変換を求める。

【例題】

正規白色雑音の中に隠された50Hzと120Hzの信号を見つけ出す。

t = [0:0.001:0.6];
x = Array([[sin(2*PI*50*t)][sin(2*PI*120*t)]]);
y = x + randn(2,length(t));
Y = fft_row(y,512);
Pyy = Re(Y * conj(Y))/512;
f = 1000*[0:255]/512;
mgplot(1,f,Pyy(:,1:256));

【参照】

fft(2.66), fft_col(2.67)