1秒ごとの重心周波数の変化を表示したい

Question

音声データの1秒ごとの重心周波数の変化を表示したいです。

現在音声データの解析を行なっています。
使用環境はPythonです。

1秒間でのフーリエ変換後のデータから重心周波数が求められるのではないかと考えましたが、コードの書き方がわかりません。

ネットで検索もしたのですが、まず重心周波数というワード自体があまり表示されないです。

サンプルプログラム、参考URLでも構いません。
わかる方に教えていただきたいです。

Comments

[Orlofsky]

日本語の「重心周波数」でググって適切な情報がhitしなかたら、英訳「Center of Gravity Frequency？」で当たってみては？
直接コードは見つからなくて参考程度にはなるかも？

Spectral centroid スペクトル重心
https://en.wikipedia.org/wiki/Spectral_centroid

[hirohiro1999]

なるほど！
探してみます！！！

[meg_]

「重心周波数」とはどのような計算で求められるのでしょうか？

[hirohiro1999]

(周波数α×振幅a＋周波数β×振幅b＋・・・＋周波数γ×振幅n)÷n
だと思っています！
後ろの記号は適当ですが全ての周波数とその大きさを掛け合わせて足したものを平均した数値じゃ無いでしょうか？
違うかもしれませんが。。。

[meg_]

下記は参考になりませんか？
https://jorublog.site/python-voice-analysis/#%E5%91%A8%E6%B3%A2%E6%95%B0%E6%88%90%E5%88%86%E3%82%92%E8%A1%A8%E7%A4%BA

[hirohiro1999]

見てみます！