よろしくお願いします。
1024日分(1日1回測定)の体重のデータがあったとします。
これをEXCELでFFTを行い,結果をIMABS関数で複素数の絶対値のみを求めました。
グラフは対象形になるので512までの結果をグラフにしました。
ここで質問なのですが,横軸(周波数軸)の解釈をどのようにしたらよいのかが分からなくなってしまいました。
単純に1日1回なので横軸は(日)を考え,150日付近にピークがありましたから,「150日周期ぐらいの体重変動がある」と考えたのですが,この解釈は誤りでしょうか。
厳密な解釈をすれば,様々あると思うのですが,「およそ●●周期の変動」という概要が知りたいと思っています。ご指導をお願いします。
No.1ベストアンサー
- 回答日時:
横軸は周波数軸なので周波数を表します。
一番最後(512番目)がサンプリング周波数の半分(ナイキスト周波数)になるので、1日1回測定の場合は1/2[1/日]、要するに2日周期の変動ですね。
150番目あたりだと、150/1024[1/日]なのでおおよそ7日周期の変動でしょう。日常生活からしてもありそうな周期ですね。
No.2
- 回答日時:
詳しい中身は分かりませんが、「窓関数」はどのように設定されましたか?
これによっては、変な結果が出ます。
http://laputa.cs.shinshu-u.ac.jp/~yizawa/InfSys1 …
http://www.logical-arts.jp/?p=124
また「グラフは対称形になるので」「512までの結果」の意味がよくわかりません。
おそらく「サンプリング定理」から、「1024個のデータを、周波数 512Hz までの周期関数」としてフーリエ解析した、という意味でしょうか。
「グラフ」の横軸が「周波数(1/日)」なのか「周期(日)」なのか、ご質問分からはよく分かりません。
データに「日」という単位を付けて入力したとは思えませんので、「横軸が周期(日)」で「150日付近にピークがある」ということはなさそうです。
横軸が「周波数 512Hz までの周波数」の場合は、「150Hz」ということでしょう。データ全体を大きな「1つの周期関数」とみなしているので、150Hzは
1024(日)× (150/512) = 300(日)
という意味ではないかと思います。
かなりの誤差を含んで「年周期」ということではないでしょうか。
No.3
- 回答日時:
No.2です。
「150Hz」という呼び方は適切ではないかもしれません。データ全体を大きな「1つの周期関数」とみなして、最大周波数512Hzで決まる「最小周期」
1024(日)× (1/512) = 2(日)
に対して、その何倍かを示す数値といった方が適切かと思います。
No.4
- 回答日時:
No.2&3です。
何度もすみません。No.2&3は、一つの「可能性」ということであって、あて推量です。
常識的に考えれば、No.1さんのように、「150」という数値を
「150Hz」→ 1024(日)× (1/150) ≒ 6.83(日)
で約7日周期と考えるのが、一番妥当かと思います。
実施した処理内容と、その表示グラフの内容を見て判断してください。
ありがとうございました。私の不十分な情報を基に判断してくださり感謝です。
実際のデータを見るとおよそ7日周期で増減を繰り返しておるのでNo.1の方の考えが妥当だと思われます。
いずれにしてもありがとうございました。
ちなみにハニングウインドを掛けないとおかしな計算になるという御指摘ですが,Excelで行う場合はどのようにしたらよいのでしょうか。(今回は窓関数を掛けずに処理しました。)
No.5
- 回答日時:
No.4です。
>ちなみにハニングウインドを掛けないとおかしな計算になるという御指摘ですが,Excelで行う場合はどのようにしたらよいのでしょうか。
エクセルでFFTを行ったことがないのでよくわかりませんが、窓関数に関する設定項目が何もないようなので、始めからハニング窓などが標準で設定されているのかもしれません。
いろいろ調べてみたら、ご質問の「グラフは対象形になる」というのは、「エイリアシング」のことなのですね。
ということは、明らかに横軸は「周波数」ですね。
No.2&3は「妄想」ですので忘れてください。
ありがとうございました。日々の体重変化のようなデータでもFFTによって周期性が見いだせることが分かれば十分です。
あっという間に解決して感謝です。
No.6
- 回答日時:
・窓関数
フーリエ変換は一般にデータ列が周期データとみなして計算をします。
その関係で、最初のデータと最後のデータにギャップがあるとき、窓関数を使わないと最後のデータと最初のデータのつなぎ目で急激な変化があったように計算され、大きな高周波成分が生成されるのです。
それを避けるために窓関数を使うのですが、それはそれでデータを歪めるので、使った方が良いとも言い難いですね。
・グラフは対象形になる
これは、グラフが対称形になるですね。
一般に実データのFFTパワースペクトルは、サンプリング周波数の半分を軸として左右対称になります。複素数でみると実部は同じで虚部が双対になるわけですが。
対称形の間違いです。御指摘ありがとうございました。
左右対称は実際にグラフを書いて試してみてたら対称形のグラフになるので納得していましたが,互いにa+ib と a-ib のように共役の複素数の関係になっているのですね。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 風力発電での音 1 2023/04/16 08:55
- 工学 CR発振回路 C1=C2=C=0.001【uF】 R1=R2=R3=16【kΩ】 R4は出力電圧が正 2 2023/05/10 19:13
- 筋トレ・加圧トレーニング 減量期に体脂肪率の変化がないことが気になっています 2 2022/08/18 21:24
- 物理学 風車が回転する理由 5 2023/05/08 08:03
- Visual Basic(VBA) Excel のユーザー定義関数でソルバーが動作しない 1 2022/09/05 19:51
- 数学 離散フーリエ逆変換が周波数分割数をNにできる理由について 4 2022/09/18 12:56
- その他(コンピューター・テクノロジー) PIC16F1シリーズマイコンのNCO機能について 1 2023/04/18 08:41
- 統計学 学業成績に関する重回帰分析の見方について 4 2022/06/06 17:19
- 物理学 風車音の周波数の特徴とその理由。 3 2023/03/18 16:42
- C言語・C++・C# C言語初心者 ポインタについて、お助けください、、 2 2023/03/15 23:50
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
サンプル数の異なる2群間にお...
-
応答で層別した場合の検定について
-
EXCELにてローパスフィルタを作...
-
心理機能診断をしたのですが、...
-
下の対数表示のグラフから低域...
-
エクセルのグラフから半値幅を...
-
脳波(EEG)の周波数は0.5~60Hzで...
-
検量線の決定係数について
-
極値をもつ時と持たない時、単...
-
ポアソン回帰でのカウントデー...
-
検定統計量の値がマイナス
-
統計学 カイ二乗検定とt検定の...
-
数3の複素数です 12分の25πが12...
-
t検定の使いどころ
-
片対数グラフで…
-
両対数グラフでの直線式の求め方
-
高校 数学 aを実数の定数とする...
-
最小二乗法を反比例の式を元に...
-
t検定の記述の仕方を教えてくだ...
-
パイロットサンプルって何ですか?
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
EXCELにてローパスフィルタを作...
-
サンプル数の異なる2群間にお...
-
心理機能診断をしたのですが、...
-
検量線の決定係数について
-
下の対数表示のグラフから低域...
-
エクセルのグラフから半値幅を...
-
x^2+y^2はどのような分布をする?
-
検定統計量の値がマイナス
-
統計について
-
【統計】有意に「高い」?「低...
-
統計学の基本的なことについて...
-
最小二乗法を反比例の式を元に...
-
心理学の統計について
-
ポアソン回帰でのカウントデー...
-
[Excel] リストからの無作為抽出
-
対数グラフの対数とは・・・?
-
極値をもつ時と持たない時、単...
-
対数目盛の読み方を教えてください
-
溶解度の問題 理科
-
正規分布について
おすすめ情報
早速大変よく分かる解説をくださり感謝いたします。
1024回の測定は全観測時間 1(日/回)×1024(回)=1024日分 したがって512番目の目盛りは 512/1024 =0.5(回/日) =2(日/回) 要するに2日周期の変動
ということは,150番目の目盛りは 150/1024 =0.146(回/日)=6.8(日/回)
他のサイトに書かれているナイキスト周波数の意味がこれでよく分かりました。ちなみに実際のデータも土日に体重増加 平日に減少に転じる変化を繰り返しております。