高速フーリエ変換からパワースペクトルを求める (ID:52785)

名前

ホームページ(ブログ、Twitterなど)のURL (省略可)

本文

>> 振幅スペクトルの２乗がパワースペクトルになるはずなんですよ・・

通じてなかったかな？
x2[i]*x2[i] + y2[i]*y2[i]
実部と虚部をそれぞれ２乗して加える。
複素数の２乗でなく、複素数のノルムの２乗です。
複素数の２乗(の実部)だと x*x - y*y になってしまうが、これは関係ない。

log(|fft|)   ... 振幅スペクトル
log(|fft|^2) ... パワースペクトル

結局 log() だから値は２倍になるだけだけどね。
↓ここをクリック。
http://www.onosokki.co.jp/HP-WK/c_support/faq/cf5200/CF5200_time_spect.htm

>> 得られた値をエクセルで書きたいです。
コンソールアプリを作成して、
データの出力をタブで区切ればいいということですね。
で出力を data.txt としたら、
これを１度ノートパッドで開き、
全選択＝＞エクセルに貼り付けです。
あなたのプログラムでは

printf("%4d | %6.3f %6.3f | %6.3f %6.3f | %6.3f %6.3f\n",
        i, x1[i], y1[i], x2[i], y2[i],log10(x2[i]*x2[i]+y2[i]y2[i]));

を
 printf("%4d\t%6.3f\t%6.3f\t%6.3f\t%6.3f\t%6.3f\n",
        i, x1[i], y1[i], x2[i], y2[i],log10(x2[i]*x2[i]+y2[i]*y2[i]));

と直せばよいですね。
ファイルに出力するなら

FILE * fout;
if( fout = fopen( "data.txt", "w" ))
{
  for( i = 0; i < N; i ++ )
  {
     fprintf( fout, "%4d\t%6.3f\t%6.3f\t%6.3f\t%6.3f\t%6.3f\n",
        i, x1[i], y1[i], x2[i], y2[i],log10(x2[i]*x2[i]+y2[i]*y2[i]));
  }
  fclose( fout );
}

とすればよい。

>> ソースを教えてもらうこと
はできないこともないと思いますが、
あなた自身、どこまでできているのかな。
...って、じゃあ、こうしよう。

一番最初のソースを作ったところまではできているとして、

まずは、[なーめ 2003/12/06(土) 04:42:49] のところで示した

make_bitrev()
の関数を

static void make_sintbl(n, sintbl)

の上に貼り付けてください。
次に

/* #define PI 3.14159265358979323846 */

を

#define PI 3.14159265358979323846

に直してください。
次にクラシックなスタイルのＣ言語を
Visual C++ で使える ANSIスタイルにします。

static void make_sintbl(n, sintbl)
int n;
float *sintbl;
{

を

static void make_sintbl( int n, float * sintbl)
{

に直してください。

int fft(n, x, y)
int n;
float *x;
float *y;
{

を

int fft( int n, float * x, float * y )
{

に直してください。

int main(argc, argv)
int argc;
char *argv[];
{

を

int main( int argc, char ** argv )
{

に直してください。

for (i = 0; i < N; i++)
        printf("%4d | %6.3f %6.3f | %6.3f %6.3f | %6.3f %6.3f\n",
        i, x1[i], y1[i], x2[i], y2[i],log10(x2[i]*x2[i]+y2[i]y2[i]));

を

に直してください。
次もクラシックＣでは型チェックが甘かったものが
ANSIで厳しくなったものです。
修正：

sintbl = malloc((n + n4) * sizeof(double));
 bitrev = malloc(n * sizeof(int));

を

sintbl = (double *)malloc((n + n4) * sizeof(double));
 bitrev = (int *)malloc(n * sizeof(int));

に直してください。
main() の中の

if (fft(N, x2, y2)) return EXIT_FAILURE;
を
	fclose( fp );
	if (fft(N, x2, y2)) return EXIT_FAILURE;

に直してください。

これでコンパイルすると、

'double' から 'float' に変換しました。

というメッセージがいっぱいでるので、
全ての float を double に直してください
ソースコードで float を検索し、もれなく行ってください。

これでとりあえずコンパイルできます。
コンパイルは Visual C++ の標準的な設定では F7 キーを押します。

まだ、直すべきところはありますが、
とりあえずここまでやってください。

ここまでで、エラーが出たら、

--------------------構成: m04 - Win32 Debug--------------------
コンパイル中...
m04.cpp
A:\Projects\test\math\m04\m04.cpp(134) : error C2146: 構文エラー : ')' が、識別子 'y2' の前に必要です。
I:\Projects\test\math\m04\m04.cpp(134) : error C2059: 構文エラー : ')'
cl.exe の実行エラー

のように貼り付けてください。
では、また「後」で。(^^;;

←解決時は質問者本人がここをチェックしてください。

戻る