No.1
- 回答日時:
質問の意味がよく分からないのですが、例えば、Excelの「ツール」→「分析ツール」→「フーリエ解析」で、ある波形データをFFT演算した後に出てくるデータ列に対応する最小周波数と最大周波数を知りたいということでしょうか?
Yes→時間間隔をΔt [s] として取った波形 A(1)~A(N) をFFTにかけて、B(1)~B(N)というデータ列が得られたとき、B(1)に対応する周波数 f1 [Hz]は、f1 = 1/(N*Δt)、B(2)の周波数は f2 = 2/(N*Δt)、・・・、B(N)の周波数は fn = 1/(Δt) となります。
例えば、Δt = 1ms でとった1024点からなる波形をFFTにかけたとき、周波数範囲は 1/1.024 Hz、2/1.024 Hz、・・・、1000 Hz になります。波形取り込み時間(N*Δt)が 1s なら f1 = 1Hz から 1Hz きざみのスペクトルになります。0.1Hzの分解能のスペクトルを得たいのなら、10秒間の波形データを取る必要があります。一方、最大周波数は、1/Δt ですので、10kHzまでのスペクトルを得たいのなら、Δt = 0.1 ms とする必要があります。
FFTだと N が2の倍数しか許されないので、(N*Δt)を 1s にしようとすると、Δt がキリのいい数値になりません(N = 1024 だと Δt = 1/1024 [s])。私個人としては、そういうのが気に入らないので、FFTでなくDFTを使って、N = 1000、Δt = 1ms とか、キリのいい数値にして計算したりします(N が10000未満ならExcel VBAで DFT 計算させてもそんなに時間がかからないので)。趣味の問題ですが。
Excelの分析ツールのフーリエ変換の結果は複素数で出てきますので、関数 IMABS(FFT後の数値)で絶対値に変換すれば振幅スペクトルになります。振幅をdB単位で出したいなら20*log10(IMABS(FFT後の数値))としてください。関数 IMABS()は「ツール」→「アドイン」で「分析ツール」にチェックしないと使えなかったと思います。
ありがとうございます。
今FFTアナライザのようなものを作っています。
指定されたサンプリング件数と周波数レンジを使って、サンプリングしたデータをFFTし、リニアスペクトルのグラフを表示するものです。
http://www.onosokki.co.jp/HP-WK/c_support/newrep …
↑のページに書いてあることを参考しています。
質問も↑のページにあることを実現すれば良いかどうかを問う内容です。
No.2
- 回答日時:
まだ詳しくURL先を見ていないので詳しいこといえませんが
2.56というのは、このFFTのサンプリング点の数なので
自分で処理をするとなると勝手が違ってくると思います。
No.3ベストアンサー
- 回答日時:
>X軸の最大値が周波数レンジになるはずですが…。
この質問は参考サイトの周波数レンジと、自分で造ったプログラムまたはエクセルなどによる解析結果の最大値が異なることに対する質問では似ないでしょうか?
FFTでは、ナイキストによりサンプリング周波数の半分までの周波数しか解析できません。おそらく最大値というのはその値になっていると思います。
参照サイトの例を用いると、サンプリング周波数は25.6kHzですので、FFTで解析すると、その半分の周波数の12.3kHzまで算出することができるのです。
しかし、小野測器などの専用FFTアナライザーでは、ちょっと取扱を変えています(他のメーカーでも私が使用したことがあるFFTアナライザーはみんなそうなっていました)。
本来1/2まで解析できるのですが、それを変えて1/2.56までしか表示しないようにして、それを周波数レンジとしています。
2048個のデータからは1024個の解析ライン数ができるのですが、800個のラインで打ち切っているようです。
なぜそうしているのか、その理由はよく知りませんが、801以上の解析ラインに対するデータの信頼性の問題ではないかと思っています。
この取扱の違いが、一般的なFFTプログラムを用いて出されるX軸の最大値と小野測器のいう周波数レンジの違いとなって現れます。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 フーリエ変換後の負の周波数成分の扱いについて 4 2022/09/03 10:18
- 物理学 打撃力は広範囲の連続周波数スペクトルを有するとあります。 どういうことなのでしょうか。 連続周波数ス 6 2022/04/21 15:50
- その他(自然科学) 原子の線スペクトルについて 1 2022/05/07 09:10
- 物理学 風力発電での音 1 2023/04/16 08:55
- 物理学 風車音の周波数の特徴とその理由。 2 2023/03/18 17:13
- 工学 CR発振回路 C1=C2=C=0.001【uF】 R1=R2=R3=16【kΩ】 R4は出力電圧が正 2 2023/05/10 19:13
- 物理学 風車音の周波数の特徴とその理由。 3 2023/03/18 16:42
- 工学 セシウムの線スペクトルについて 1 2022/05/04 09:23
- 物理学 測定値と理論値の誤差について 交流回路の実験でRL回路、RC回路、RLC直列回路の周波数を上げた時の 1 2022/05/22 23:37
- 物理学 物理の波のグラフについてですが、 「この波のグラフを1/2周期(半周期)だけ進める」 「+x軸方向に 1 2022/08/28 18:03
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
近似直線について
-
MOSFETのgm-Vgs特性
-
半値幅の測り方
-
分解能の値が大きいと小さいも...
-
1/3オクターブバンドについて質...
-
拍手の間隔を解析で求めるには?
-
フーリエ変換、FFT演算について
-
有限要素法解析のエラーについて
-
プランク定数の実験をしたので...
-
画像の式はフーリエ級数展開を...
-
静的・動的の意味
-
有効数字…?
-
yの二乗をXで微分したら2y・y' ...
-
時間を100進法であらわしたい。
-
合成関数の微分を使う時と、使...
-
穴が開く? 空く? 明く?
-
「強度」は高い?強い?
-
植込みボルト製図
-
sin75°×sin15°の値を教えてくだ...
-
積分定数Cとは一体なんですか?
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報