• ベストアンサー

1/3オクターブバンド分析の和の求め方

1~200Hzまで0.01Hz間隔で振動レベルがわかっている場合、80Hzを中心とする1/3オクターブバンドレベルの大きさはどのようにして求めるのですか? まず、80Hzのバンドを中心としたフィルターを通して、その後の各周波数の振動レベルの和を求めていくと考えているのですが、うまくいきません。 教えていただけないでしょうか? 各周波数ごとの和をとった場合、周波数刻み幅が大きいほど値が大きくなってしまうと思うのですが、周波数間隔で割るなどの処理が必要なのでしょうか?

質問者が選んだベストアンサー

  • ベストアンサー
  • ufon
  • ベストアンサー率64% (37/57)
回答No.1

中心周波数80Hzの周波数範囲は71Hzから90 Hzですから、中心周波数80Hzの振動レベルは71Hzから90Hzの19個の振動レベル L71、L72、・・・ L90 dB のdB和として次式で求められます。 10log[10^(L71/10)+10^(L72/10)+・・・+10^(L90/10)] dB  また71Hzの振動レベル L71は70.01~71.00Hzの100個の振動レベルのdB和として上式で求められます。 72、73・・・90Hzのレベルについても同じです。 計算は以上のとおりですが、1~200Hzの範囲を中心周波数80Hzの1/3オクターブバンドフィルターを通した結果得られるのが、上記計算結果です。                     以上

saterain20
質問者

お礼

ありがとうございました。

関連するQ&A

  • オクターブバンドの周波数帯の考え方

    例えば、 ●中心周波数500Hzで1オクターブバンドだと 周波数の領域は 500*2^(-1/2)~500~500*2^(1/2) ●中心周波数500Hzで1/3オクターブバンドだと 周波数の領域は 500*2^(-1/6)~500~500*2^(1/6) ●中心周波数1500Hzで1/2オクターブバンドだと 周波数の領域は 1500*2^(-1/4)~1500~1500*2^(1/4) というふうになる気がするのですが、 合っているでしょうか? どなたか答えあわせをお願いします。

  • 1/3オクターブバンド分析

    200Hzでサンプリングした振動波形(加速度データ)のCSVデータがあります。この振動波形データを1/3オクターブバンド分析したいのですが、参考となる書籍が見つからず困っています。これまで調べた書籍、ネットの情報から、まずFFTをかけて、周波数(Hz)と加速度(m/s^2)のデータに変換。このデータを1/3オクターブバンド毎に2乗和して平方根する?・・・・そのあと逆フーリエ変換?・・・ぜんぜんわかっていない状況です。 なにか参考になる良い書籍があればご教授お願いします。

  • 1/3オクターブ分析

    ある振動波形のFFT結果として、周波数(Hz)と加速度(m/s^2)があります。このデータをExcelを用いて1/3オクターブ分析の値を計算上(数学的な算術結果)求めたいのですが、計算式(定義式?)をご教示いただけませんでしょうか?何卒宜しくお願いいたします。

  • オクターブバンド分析とは

    床の上を人が歩くことに伴う振動測定(鉛直振動)を行い、オフライン処理(ソフトウェア上)で測定データから、1/3オクターブバンド分析による加速度の最大値(0-P)が求めたい結果です。 そもそもオクターブバンド分析というものは、全体波形(振動波形)のどのタイミングで処理を行っているものなのでしょうか?振動波形の中で最大値を示す箇所(もしくは、オーバオールが最大の位置)で処理を行いたいのですが、関連書物を見てもいまいちピンとくるものがありませんでした。 私自身の勉強不足ではありますが、何卒ご教示宜しくお願いいたします。

  • 1/3オクターブバンド分析について

    とある音のラウドネスを求めるために1/3オクターブバンド分析を行うプログラムを組みたいのですが、詳しい式や過程が載っているサイトや書籍がなかなか見つかりません。なので基本的なところから説明を頂きたいと思っています。 また、ここでの質問の中に周波数のほかに加速度を用いるといった記述があったのですがそれは音響分野でも同様なのでしょうか。 かなり初歩的な質問なのかもしれませんがこれらにご回答いただけると嬉しいです。 よろしくお願いします。

  • バンドパスIIRフィルタの低域におけるエラー

    バンドパスIIRフィルタの低域におけるエラー 申し訳ないのですが教えて下さい。  波形データのオクターブバンド分析結果を出す為に、ScilabやOctaveなどMatlabクローンのソフトを使ってバンドパスIIRフィルタを使用して各バンドでフィルタを通過させるのですが、ナイキスト周波数に対して通過域が低域にありすぎるとエラーを起こすのか、結果が出てきません。サンプリング周波数の値を変えて遮断周波数との差を縮めればそのようなエラーは起きません。例えば、480Hzサンプリングでも、1Hzと2Hzの中心周波数のバンドパスフィルタは値が帰っていなかったり、波形がおかしかったりします。  フィルタ次数は4で行っていて、1などに下げればそのナイキスト周波数と遮断帯域の差ではエラーを起こさなくなりますが、もっと差が広がった場所では起きますし、JISで定めるオクターブフィルターの遮断性能を考えると次数を下げる訳にも行きません。  本家Matlabはもっていないので判りませんが、会社のLabviewですと問題はありません。解決策はデサンプリングしてエラーを起こす帯域ではサンプリング周波数を落とすしか無いのかな?と思うのですが、実際はどうなのでしょうか?お教え頂けると幸いです。 宜しくお願い致します。

  • 1オクターブはなぜ12半音なのですか

    音楽に関して色々疑問がでてきたので質問するに至ります。 ドレミファソラシドという音階が昔からあり、周波数の測定ができるようになってから 黒鍵が追加された、ということが調べてわかったのですが どうして12半音になったのかが理解できません。 ドレミファソラシド(2 2 1 2 2 2 1)のミとファ、シとドの間隔に合わせて黒鍵を追加した、ということでしょうか 調べているうちに半音につき2の12乗根(≒1.05946309436)倍,周波数が高くなっている、だとか オクターブ上がるということは周波数が倍になっているだとかがでてきてしまって こんがらがってしまいました。 A110Hzから(1.05946309436)をかけていったところ 220.0000000017562Hz、440.000000007025Hz と端数がでてきてしまいさらにわけがわからなくなりました。 こんな文章で伝わるでしょうか・・・ できれば色々こんがらがってしまっているので簡単に説明していただけると助かります。 もし回答ですら理解できないようであれば1から楽典の勉強を始めようと思っています。

  • デジタル信号処理でのたたみ込み和

    ここに質問を書くのは初めてです。 間違ったこと、してはいけないことがあった場合はご指摘ください。 質問です。 与えられた信号を解析して、そのパワースペクトルを図示する問題がでたのですが、 データ長が異なる場合のたたみ込み演算を行う方法が分かりません。 16000Hzで16000点の信号と標本化周波数の16000Hzでフィルタ次数101点の低域通過フィルタが与えられています。 この二つのたたみ込み和を周波数軸上でおこなって、そのパワースペクトルを図示したいのですがどうしたら良いでしょうか? 因みに、C言語で書かなければなりません。 ただ、計算方法が分かれば、プログラミングはできると思うので、「データ長が異なる2つのデータのたたみ込み和を周波数軸上でおこなう方法」を優先で教えて下さい。 質問文で分かるかもしれませんが、いまいち、理解できていないので、易しめに教えていただけると助かります。

  • オクターブの換算について

    オクターブの換算について 1/12oct.から1/3oct.への変換でわからなくて困っています。 1/3oct.;1000Hzの値を1/12oct.の値から求める方法を教えてください。 下記に例として1/12oct.で800~1250Hzまでのレベル(a,b,c,...)を記入していますが、 ここで、1000Hzのところの1/3oct.のレベルを求めたいです。 800Hzおよび、1250Hzについては別途求めることとします。 1/12 Level 1/3 ----------------- 800 a 800 850 b 900 c 950 d 1000 e 1000 1060 f 1120 g 1180 h 1250 i 1250 ----------------- よろしくお願いします。

  • 任意の周波数からオクターブ幅がすぐに分かるようにしたいのですが

    1000Hzの1オクターブ上が2000Hzというのは分かるのですが、ではエクセルのA1セルに1000、A2セルに2000と入れたらA3セルに1、と答えが出る数式はどういうふうになるのでしょうか? 任意の周波数間(例えば5500Hz~20000Hz)のオクターブを即、分かるようにしたいのですが、どういう数式をA3セルに入れるといいのでしょうか。 音楽理論を習った方なら初歩の質問かもしれませんが、ひとつよろしく回答をお願い致します。