• ベストアンサー

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

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

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

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

 1/3オクターブバンドの加速度はバンド内周波数の加速度の自乗の和の平方根になります。  例えば中心周波数10Hz(9~11Hz)の1/3オクターブバンドの加速度m/s^2は、バンド内の周波数9、10 、11Hzの加速度の自乗の和の平方根になります。1/3オクターブ中心周波数ごとに順次この計算をすると、 分析結果がえられます。ただこれだけです。

gufo02
質問者

補足

ufonさん 回答ありがとうございます。しかし、基本的なことがわかっていないので、段階を追った分析方法をきちんと知りたいのです。 波形データを 1/3オクターブバンドの幅でバンドパスフィルタをかける。かけたデータを ufonさんの回答のように処理すればよいのでしょうか? 詳しく書かれた参考書などはないでしょうか? できれはVBA等でプログラムを作って処理したいと思っています。

その他の回答 (2)

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

回答への補足に対する回答です。 1/3オクターブバンドごとのバンドレベルを求めるには、波形データを1/3オクターブバンドの幅でバンドパスフィルタをかけ、得られた1/3オクターブバンドごとの加速度の実効値aを求め、これから振動加速度レベル20log[a/10^(-5)]dBを求めると、これが1/3オクターブバンドレベルになります。 オクターブバンドレベルを求めるには、三つの1/3オクターブバンドレベルを10log[10^(L1/10)+10^(L2/10)+10^(L3/10)]として求めることができます。

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

補足質問に対す回答です。 波形データを1/3オクターブバンドの幅でバンドパスフィルタをかけると、得られるデータは1/3オクターブバンドごとの加速度になります。これが1/3オクターブ分析結果です。かけたデータを処理する必要はありません。得られた1/3オクターブバンドごとの加速度を加速度レベルに換算したりする場合は処理が必要ですが。 基本的なことなので、特に詳しく書かれた参考書などは見当たりません。 

gufo02
質問者

補足

ufon さん引き続き回答ありがとうございます。 こちらの説明不足でしたが、もとめたいのは、1/3オクターブバンドごとのオクターブバンドレベルで、オクターブバンドスペクトルを作りたいというのが目的です。 すいません、ほんとに基本的なことからわかっていないもので・・・

関連するQ&A

  • 1/3オクターブ分析

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

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

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

  • 1/3オクターブバンドについて質問です。

    1/3オクターブバンドについて質問です。 加速度信号をFFTアナライザを用いて、パワースペクトルを計測しました。 そのパワースペクトルのデータを用いて、計算上で1/3オクターブバンドに変換することは可能でしょうか? いろいろ調べたのですが、残念ながらよく分かりませんでした。 どなたか、解る方がいましたら、ご回答お願い致します。

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

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

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

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

  • 断続切削の振動

    旋盤で断続切削中の振動波形をみようと、市販の加速度計とロガーで計測を行いました。 加速度計は、加工に近い刃物台に設置。サンプリングは100μ(10000Hz)で行い、フーリエ変換した結果より特定周波数を見出そうとしております。 振動波形を見る限り、減衰波形などもとれているので計測状態は問題ないとみております。 質問したいのは、フーリエ変換して出てくる周波数に断続切削で出ると想定していた周波数領域が現れないことです。 詳しくは、変換後には振動減衰時の波形周波数(1.5kHzあたり)が大きく出ており、断続の影響で現れる波形(100Hz)が全く出ておりません。 これは、この加工状態で問題となるのは1.5kHzであり、断続の100Hzは問題視する必要なしと捉えていいものでしょうか? 加工の音的には、断続の音が聞こえてくるのですが...

  • 非定常波形のフィルタリング

    構造物の振動計測で加速度波形を測定します。アナログ波形を200Hzほどでデジタルサンプリングします。計測対象周波数帯域は5Hz以上なのですが、測定データには2Hz以下の低周波成分が混じります。高価な計測機器であればハイパスフィルタで3Hz未満をカットしてしまえばいいのですが、あいにく手持ちの機器には10Hz以上を通すハイパスの設定しかありません。 そこでAD変換をしてサンプリングした後に後処理で3Hz以下をカットしたいと思います。手始めに全体をFFTして低周波成分の係数を除いて逆変換してみました。ところが低周波ノイズは非定常で、すなわち時刻暦波形の中で低周波ノイズが入っている部分と入っていない部分があります。それを全体一括で扱うと言うのは無理があります。ハードウェアでアナログフィルタをかけるように波形の中のその場その場でフィルタをかけて低周波成分が混じっている部分のみカットするということが出来ればうれしいです。 1.こういう場合、どのような手法が取られるのでしょうか? 2.あるいは10万円程度以下のソフトでこういった処理の出来るソフトはありませんか?   振動に限らず、音声や通信など他分野のソフトでも結構です。 よろしくお願いします。

  • A/D変換器のサンプリング周期誤差について

    A/D変換器のサンプリング周期誤差について 簡易FFTアナライザーの作成を考えております。 そこで、A/D変換に考えている装置が1msでのサンプリング周期で ±10usのサンプリング周期誤差があります。 この装置からの入力データで簡易FFTアナライザーの作成は可能でしょうか。 (たとえは、300Hzぐらいの波形を測定する時に、サンプリング ノイズでスペクトル  が見えなくなり、使用できない等) 測定周波数はDC~400Hz程度で1024点程度のデータ点数を考えています。

  • FFTの結果からフーリエ係数は求められるものでしょうか?

    FFTの結果からフーリエ係数は求められるものでしょうか? 実測値の波形データを関数で表したいのです。 フーリエ展開でなんとかならないかと考えたのですが、 計算の速いFFTを使うと、逆変換で波形データに戻せるものの、 関数として表す事が出来ない事に、 FFT後に気付いてしまいました。 数学は専門外で、入門書と格闘しながらやっているのですが、 行き詰ってしまいました。

  • フーリエ変換後の解析について

    フーリエ変換後の解析について ある振動の波形を測定していたとします。 No.1 最初は特に何も振動を加えず、ありのままを測定したとします。 No.2 そして次に周期的に弱い振動を与えたとします。 1と2の結果得られた振動波形を見比べた結果、ほとんど違いは見られなかった場合 FFT解析を行う意味はあるのでしょうか? 波形にほとんど違いがなかったということは本当に微妙な振動が加わっただけということですよね? そんな微妙な違いしかない結果をFFTしたところでほとんどスペクトルにも差がでないと思うのですが・・・。 それともたとえ微妙であっても周期的に加えてる以上どこかしらにピークがでてくるものなのでしょうか?