ゆき社長

シーゲンガーのお勉強　ゲームプログラマ、ゲーマー、色々！

FFT 始めました！ part.2

---- いろいろ ----

Part 1

ってことで　前回は、FFTWの導入までおわったかんじっ

だけどまた　　例のページ　ぱくるのだよ

とりあえず、FFTができたっぽいけど、数字じゃよくわからないわけで、グラフに表示しましょう

ってゆーのと、正しいデータなのか 大まかにチェックできるようにしましょう　ってことで

サンプリングデータを 正弦波にして 図形にプロットしようと思います！！

例のページパクリですが

gnuprot よりDownload

バージョンは　 4.6.3 にしましたわ　gp463-win32-setup.exe　でインストールしたし。

で プログラムの最後の部分に csvファイル出力、gpファイル出力ルーチンを

加えます



		// fftw plot




		FILE *fp;


        fp = fopen("fftw.csv","w");


        for (i = 0; i< SIZE; i++)


                fprintf(fp,"%d,%f,%f\n", i, out_buf[i][0] * rsize, out_buf[i][1] * rsize);


        fclose(fp);




		fp = fopen ("fftw.gp","w");


        fprintf(fp,"plot \"fftw.csv\" using 1:2 w l,\"fftw.csv\" using 1:3 w l\n");


        fclose(fp);

って　この辺までは 例のサイトパクリなんで 飛ばしますが

グラフが表示できたはずです

とりまここまでで！

次は　大きなサイズのサイン波をフーリエ変換して　時間測定したいです