高速フーリエ変換.pptVIP

ディジタル信号処理Digital Signal Processing 第１７講 ＦＦＴ Fast Fourier transform 高速フーリエ変換 高速フーリエ変換 高速フーリエ変換（ Fast Fourier Transform, FFT）とは、離散フーリエ変換 (Discrete Fourier Transform, DFT) を計算機上で高速に計算するアルゴリズム。 　　逆変換をIFFT (Inverse FFT) という。 高速フーリエ変換といえば一般的には1965年、ジェイムズ?クーリー (J. W. Cooley) とジョン?テューキー (J. W. Tukey) が発見したとされ，Cooley-Tukey型FFTアルゴリズムを呼ぶ。しかし、1805年ごろにガウスが同様のアルゴリズムを独自に発見していたといわれている。 信州大学工学部?井澤裕司先生URL:http://laputa.cs.shinshu-u.ac.jp　　/~yizawa/InfSys1/basic/chap7　　/index.htm わかりやすい記事が載っています 先頭ページは 　　　http://laputa.cs.shinshu-u.ac.jp 　　　　　/~yizawa/InfSys1/basic 　　　　　/index.htm 原理の簡単な説明(Wikipedia) 離散フーリエ係数は、exp( -2 π ｊ / N ) を使うと、次のように表せる。 5.4　高速フーリエ変換　（FFT : Fast Fourier Transform） 5.4.1 基数２のFFT ■ビット反転の関係 5.4.2　混合基数FFT 5.4.3　高速フーリエ逆変換 計算の高速化計算量を減せば早く答が得られる 式をそのまま行うと膨大な計算量になる 計算量を大幅に減少させるアルゴリズムを高速フーリエ変換（fast Fourier transform : FFT ）という 基本的な考え方： 　①　データ数 Ｎ を ２p にする 　②　計算順序を変える ③　偶数番のデータと奇数番のデータ (n=1,2,3,????,N/2-1) x0(n)=x(2n)　　 ???偶数番 　　x1(n)=x(2n+1)　???奇数番 高速フーリエ変換の計算手数 データの並べ替えは無視する 積和演算の回数は，r個の積和演算をq回行い，さらにこれをp回繰り返すことになる N＝q?r＝２pであるから，Nだけで表すと演算回数＝N?log2N 演算回数比η＝FFT／DFT＝ N?log2N／N2 Nが大きいときFFTの高価は顕著に表れる 高速フーリエ変換のメリット次の例で考えてみる Ｎ＝８ 鋸歯状波 離散時間信号 x(0)=0,x(1)=1,x(2)=2,x(3)=3, 　x(4)=4,x(5)=5,x(6)=6,x(7)=7　 離散フーリエ変換の式（３．２）式??標準式 X(k)の計算 X(k)=Σx(n)?exp(-j?2?π?k?n / N)を展開して Ｘ(k)={ｘ(0)?exp(-j?2?π?k?0 / 8) 　+ｘ(1)?exp(-j?2?π?k?１ / 8) 　+ｘ(2)?exp(-j?2?π?k?2 / 8) 　+ｘ(3)?exp(-j?2?π?k?3 / 8) 　+ｘ(4)?exp(-j?2?π?k?4 / 8) 　+ｘ(5)?exp(-j?2?π?k?5 / 8) 　+ｘ(6)?exp(-j?2?π?k?6 / 8) 　+ｘ(7)?exp(-j?2?π?k?7 / 8)} Ｃ言語では，複素計算ができないので， 　　exp(-j2πkn/N)=cos(2πkn/N)-jsin(2πkn/N) 　と実数部，虚数部別々に計算しなければならない Ｘ(k)にオイラーの公式を当てはめ Ｘ(k)={ｘ(0)?exp(-j2πk0/8)+ｘ(1)?exp(-j2πk１/8)+ｘ(2)?exp(-j2πk2/8) 　+ｘ(3)?exp(-j2πk3/8)+ｘ(4)?exp(-j2πk4/8)+ｘ(5)?exp(-j2πk5/8) 　+ｘ(6)?exp(-j2πk6/8)+ｘ(7)?exp(-j2πk7/8)} 　　= [ｘ(0)?{cos(2πk0/8)-j?sin(2πk0/8)} 　　+ｘ(1)?{cos(2πk1/8)-j?sin(2πk1/8)} 　　+ｘ(2)?{cos(2πk2/8)-j?sin(2πk2/8)} 　　+ｘ(3)?{cos(2πk3/8)-j?sin(2