c言語でDFTのプログラムを作成しました。
以下にソースを載せます。
#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include<time.h>
#define PI 3.141592653589793
#define N 64 //データ数
DFT(double result[]){
int i,k;
double A[N],B[N],T=0; //A[N]:実数部,B[N]:虚数部
double a,b;
for(k=0;k<N;++k){
a=b=0;
for(i=0;i<N;++i){
a=a+result[i]*cos(2*PI*i*k/N);
b=b+(-1.0)*result[i]*sin(2*PI*i*k/N);
}
A[k]=a/N;
B[k]=b/N;
}
for(i=0;i<N;++i){
printf("[%f秒]:Re:%f,Im:%f\n",T,A[i],B[i]); //変換後の値を表示
T=T+(0.1/N);
}
}
main(){
int i;
double T=0;
double Sampdata;
double result[N];
for(i=0;i<N;++i){
Sampdata=5*sin(20*PI*T); //0~0.1秒間をN個にサンプリング
result[i]=Sampdata; //サンプリングデータを代入
T=T+(0.1/N);
}
clock_t start,end; //処理時間計測開始
start=clock();
DFT(result);
end=clock();
printf("%.2f秒かかりました\n",(double)(end-start)/CLOCKS_PER_SEC); //処理時間表示
}
元信号には5sin(20πt)の値を入れています。この信号は周期は0.1secです。
これでフーリエ変換を行うとデータ数N/2を中心に対称なデータが出てくるのですが、処理が終わるのが早い気がするんです。
例えば2^15個のデータで実行しても2分もかからずに処理が終わってしまいます。一応、対称性が出てるとはいえ、終わるのが早すぎる気がするのですが、おかしい所があれば教えていただけると嬉しいです。
よろしくお願いします。
No.1ベストアンサー
- 回答日時:
計算に使っているコンピュータの性能はどれくらいなのでしょう?
このプログラムで一番時間がかかると思われるのは
> for(i=0;i<N;++i){
> printf("[%f秒]:Re:%f,Im:%f\n",T,A[i],B[i]); //変換後の値を表示>
> T=T+(0.1/N);
> }
この出力の部分です。
手許のcore i7 2.66GHz
gcc -O0 (最適化無し)
DFT(result); を2^15回ループ
結果出力部分をコメントアウト
で22秒程でした。
CPUの速度、コンパイラの最適化等でもっと速くなります。
2分かからないくらいなら、別に変では無いと思います。
心配なら、この計算結果を、別の方法で計算したもの(実績のあるDFT計算ライブラリ、Excel等の別ソフト等)と比較してみては?
この回答への補足
回答ありがとうございます。
私の周りでは2^15でDFTを行った場合に30分も時間がかかったという意見があったのですがこんなにかかるものなのでしょうか?
私のと処理の仕方が違うかもしれませんが、それでも2^15のデータの処理に30分もかかるというのが自分と違い、不安になってこの質問をしました。
その周りの人にもプログラムを教えてもらえるような状況ではないのですがどうでしょうか?
しつこいですがお願いします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# プログラミングの授業の課題です 1 2023/01/17 22:15
- C言語・C++・C# 並列プログラミングのπ計算について 1 2022/07/16 22:30
- C言語・C++・C# 10個の実数に対する降順ソート結果を出力するプログラムを作りたいのですが、以下のプログラムをどう直せ 1 2022/07/09 22:16
- C言語・C++・C# LU分解法のピボット選択機能実装について(C言語・gcc-9) 1 2022/07/22 15:20
- C言語・C++・C# LU分解法のピボッティングについて(C言語/gcc-9) 3 2022/07/11 23:10
- C言語・C++・C# C 言語の Gauss Jordan 法について 2 2022/12/28 11:16
- C言語・C++・C# バイナリファイルをコピーするのにかかる時間を測りたいのですが実行するとFatel error:gli 2 2022/11/03 01:10
- C言語・C++・C# c言語の問題です 3 2023/01/10 16:15
- C言語・C++・C# C言語 プログラミング 4 2022/05/22 11:53
- C言語・C++・C# このプログラミングの問題を教えてほしいです。 キーボードからデータ数nとn個のデータを入力し、平均値 3 2022/12/19 22:51
関連するカテゴリからQ&Aを探す
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
プログラムでの数字につく”f”の...
-
float型とdouble型の変数の違い...
-
doubleの変数にintとintの割り...
-
C言語 関数プロトタイプ宣言の...
-
C言語で台形公式を使った二重積...
-
関数におけるif文とreturn文に...
-
C言語(プログラミング)関連の質...
-
C言語を実行すると-infが出てき...
-
C 開放してるのにエラー(doubl...
-
数値を指数部と仮数部に分離したい
-
マチンの公式による円周率のプ...
-
至急です! マクロ定義で #defi...
-
C# 分秒表示ついて
-
c言語で、繰り返し文の中で、0....
-
C言語でdouble型の小数点の引き...
-
学校の課題で2次方程式のプログ...
-
ニュートン法
-
C言語で表記についの質問です
-
c言語の問題
-
2分法で方程式の複数の解を自...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
プログラムでの数字につく”f”の...
-
float型とdouble型の変数の違い...
-
doubleの変数にintとintの割り...
-
C言語を実行すると-infが出てき...
-
C 開放してるのにエラー(doubl...
-
至急です! マクロ定義で #defi...
-
c言語で、繰り返し文の中で、0....
-
関数におけるif文とreturn文に...
-
C言語 関数プロトタイプ宣言の...
-
C言語初心者 構造体 課題について
-
C言語の型による処理速度の違い
-
Cで3乗根を求める方法
-
C言語で-23乗を取り扱うには
-
2分法で方程式の複数の解を自...
-
doubleは常に%lfとするべきなのか
-
c言語のコンパイルエラー canno...
-
C言語で直角三角形の斜辺を求め...
-
C言語のプログラムで#include<m...
-
int とdoubleの比較
-
C++で外積
おすすめ情報