オクターブ帯域分析:その数学的・工学的根拠

オクターブ帯域分析は、対数周波数軸上の定割合帯域幅を用いて、詳細なスペクトルを標準化された1/1オクターブ帯域および1/3オクターブ帯域に変換します. 本記事では、CPB(定割合帯域幅)の数学的な基礎、IEC 61260-1およびANSI S1.11がオクターブ帯域をそのように定義している理由、そして実務上どのように帯域レベルが計算されるのか(FFTビニング vs. フィルタバンクRMS)を解説します. 目的は、音響、NVH、各種規格に準拠した測定において、再現性が高く相互比較可能な結果を得ることです.

オクターブ帯域分析とは何か、そしてそれはどのような問題を解決するのか?

オクターブ帯域分析は、周波数軸を対数スケール上で分割し、バンドパス帯域に区切るスペクトル解析手法のファミリーです. 各帯域は、その上下の遮断周波数間の比が一定(定割合帯域幅:CPB)になるように設計されています. 各帯域内では、細かな線スペクトルの詳細は無視し、その帯域内の総エネルギー／RMS(または電力)に着目します.

言い換えると、「1 Hzごとに何が起きているか」ではなく、「等しい相対帯域幅にわたってエネルギーがどのように分布しているか」を見るものです.

この表現は、人間の聴覚や、多くの工学システムに自然に対応しており、それらの周波数分解能は固定Hzスケールよりも相対的(対数)スケールに近いことが多くあります.

• これは、多くの規格で要求される一般的な報告フォーマットでもあり、室内音響パラメータ、防音性能評価、環境騒音、機械騒音、風騒音／路面騒音などで、1/3オクターブ帯域がよく用いられます.

線形Hzから対数周波数へ:なぜCPBは「工学の言語」に近いのか

等幅の周波数ビン(例:10 Hzごと)でエネルギーを積算すると、スペクトル全体で一貫しない挙動になります:

• 低周波数では、10 Hzビンは広すぎて、詳細をぼかしてしまうことがあります.
• 高周波数では、10 Hzビンは狭すぎて、ランダム雑音に対する分散が大きくなり、推定値が不安定になります.

対照的に、CPBでは帯域幅は周波数とともに増加します(Δf ∝ f). 各帯域は同程度の相対変化をカバーするため、安定性と再現性が向上し、標準化試験にとって重要です.

視覚的な直感:線形軸では帯域幅が増大し、対数軸では一様になる

図1:同じ1/3オクターブ帯域を線形周波数軸上にプロットした場合-高周波数側では帯域幅が大きく見える

各水平線分は1/3オクターブ帯域[f1, f2]を表し、短い縦線が帯域中心周波数fmです. 線形軸では、高周波数側の帯域ほど広く見えます.

図2:同じ帯域を対数周波数軸上にプロットした場合-帯域が等間隔に並ぶ(これがCPBの本質)

横軸を対数スケールにすると、これらの帯域は等しい幅／等しい間隔に見えるようになります.これこそが「定割合帯域幅」の意味するところです.

これらの図は中核となる考え方を示しています.すなわち、オクターブ帯域分析はHzでの等間隔ではなく、対数周波数スケール上の等間隔ステップを用いるということです.

規格と用語:IEC／ANSI／ISOの体系は実際に何を規定しているのか?

実務上、「1/3オクターブ分析を行う」ということは、単に帯域端だけで決まるわけではありません. 規格は(明示的あるいは強く示唆する形で)、中心周波数の定義方法(厳密値と公称値)、オクターブ比の定義(ベース10かベース2か)、フィルタの許容差／クラス、さらには帯域レベルを形成する際の測定・平均化の流儀まで規定しています.

IEC 61260-1:2014の要点:ベース10の比、基準周波数、および中心周波数の式

IEC 61260-1:2014は、オクターブ帯域および分数オクターブ帯域フィルタに関する主要な仕様です. この規格はベース10の設計を採用しており、オクターブの周波数比としてG = 10^(3/10) ≈ 1.99526(2に非常に近いが、厳密には2ではない)を用います. 基準周波数はfr = 1000 Hzです. 規格では、正確な中帯域(中心)周波数を与える式を示すとともに、帯域端周波数の幾何平均が中心周波数に等しくなるよう規定しています. [1]

主要な式(規格から並べ替えたもの):[1]

分数の分母bが奇数(例:1, 3, 5, ...)の場合:

If b is even (e.g., 2, 4, 6, ...):

常に次が成り立ちます:

なぜbが偶数の場合は「半ステップずれて」見えるのか? 直感的には、中心周波数のグリッドはlog(f)上で等間隔になっています. When b is even, IEC chooses a half-step offset relative to fr so that band edges align more neatly in common reporting conventions. 実装上の堅牢な方法は、規格の式を用いて正確なfm列を生成し、その後f1 = fm / G^(1/(2b))、f2 = fm * G^(1/(2b))により帯域端を計算し、最後に通常の公称周波数で帯域にラベル付けすることです.

OpenTestでデータを表示する(IEC 61260-1オクターブ帯域分析) ->

帯域端、中心周波数、および帯域幅指定子b

規格では一般に1/bを「帯域幅指定子」として用い、1/1が1オクターブ、1/3が1/3オクターブなどと表します. [1] いったん(G, b, fr)が選ばれれば、帯域集合全体(中心と端)は数学的に一意に定まります.

厳密値と公称値:同じ帯域に二つの「中心周波数」が現れる理由

「厳密」中心周波数は、数学的に一貫した定義やフィルタ設計に用いられ、「公称」値はラベル付けや報告に用いられます. [1] ISO 266:1997は、ISO 3の推奨数列(R10)に基づき、1000 Hzを基準とした音響測定用の推奨周波数を定義しています. [2]

その結果、厳密な等比数列は、通常次のような見慣れた公称値でラベル付けされます:

20, 25, 31.5, 40, 50, 63, 80, 100, 125, 160, ..., 1k, 1.25k, 1.6k, 2k, 2.5k, 3.15k, ..., 20k.

実装上のヒント:まず厳密な中心周波数から帯域端を計算し、公称値への丸め／表示はその後に行います. これにより、規格から外れていってしまうことを防げます.

ベース10 vs ベース2:なぜ規格は厳密な2:1オクターブにこだわらないのか

「オクターブ」はしばしば2:1として捉えられますが、IEC 61260-1はG=2ではなく、ベース10(G=10^(3/10))を採用しています. 主な理由には次のようなものがあります:

• 10進の推奨数列との整合(ISO 266はR10に結び付けられています). [2]
• 国際的な一貫性:IEC 61260-1:2014はベース10を規定しており、ベース2設計は基準周波数から離れた領域で適合性を維持しにくいことが示されています. [1]

ベース10では、1/3オクターブは10^(1/10) ≈ 1.258925(1/10デケードとも解釈可能)に対応し、1デケードあたり正確に10個の三分の一オクターブ帯域というきれいな対応が得られます.

「10個の三分の一オクターブ帯域＝1デケード」:これが重要な理由

ベース10の三分の一オクターブ間隔では、各ステップで周波数はr = 10^(1/10)倍されます. したがって:

• 連続する10個の三分の一オクターブ帯域で、周波数は正確に10倍(1デケード)になります.
• これはISO 266／R10の慣習と一致しており、表の作成やプロット、コミュニケーションを簡素化します.
• 標準化では、純粋な数学的厳密さと同じくらい、可読性と一貫性が重視されます.

図3:ベース10の三分の一オクターブ間隔-1デケード(周波数×10)あたり等比ステップが10個

ANSI S1.11／ANSI/ASA S1.11:許容クラスと過渡信号に関する注意

ANSI S1.11(およびIEC 61260-1と整合した後続のANSI/ASA版)は、フィルタセットやアナライザの性能要件を規定しており、許容クラス(版によってクラス0/1/2など)も含まれます. [3][4]

ANSI文書には実務上の注意として、過渡信号では、規格に適合した異なる実装同士でも異なる結果になりうることが明記されています. [3] This highlights that time response (group delay, ringing, averaging time constants) matters for transient analysis.

クラス／マスク／有効帯域幅は実際には何を制御しているのか?

「1/3オクターブ帯域を使った」と言うことは、公称の帯域端を使ったというだけの意味ではありません. 規格は、次の点を制約することで、異なる機器／アルゴリズム間で結果が比較可能になることを目指しています:

• 周波数の配置:中心周波数の列と帯域端の定義(ベース10、厳密値／公称値、f1／f2).
• 振幅応答の許容差(マスク):通過帯域近傍で許容されるリップルと、中心周波数から離れた領域で必要な減衰量.
• 広帯域雑音に対するエネルギーの一貫性:実装間で帯域レベルが比較可能になるよう、有効帯域幅に制約を設けます.

有効帯域幅が重要なのは、実際のフィルタが理想的なブリックウォールではないからです. 広帯域雑音については、出力エネルギーは∫|H(f)|^2 S(f)dfに依存します. 通過帯域リップルやスカート部、ロールオフの違いによって、系統的なオフセットが生じる可能性があります. 規格は有効帯域幅を制約することで、そのようなオフセットが許容範囲内に収まるようにしています. [1][3][4]

過渡信号に関する注意は矛盾ではありません.マスクは主に定常状態の周波数領域での挙動を制約する一方で、過渡応答は位相／群遅延、リンギング、時間平均に強く依存するためです. [3]

数学:帯域定義、帯域幅、Q、および帯域インデックス付け

CPBと対数軸上の等間隔配置

CPBは、対数周波数における等幅配置と同値です. If u = log(f), then every band spans a fixed Δu. 多くのスペクトル(例えば1/f型)は、対数周波数で見るとより滑らかで統計的に安定して見えます.

幾何平均定義から導かれる帯域端式(一般的な1/b形)

IECでは、中心周波数を帯域端の幾何平均fm = sqrt(f1 f2)として定義します. [1] 1/bオクターブ帯域では、帯域端の比は通常f2/f1 = G^(1/b)となり、ここでGはオクターブ比です. このとき:

ベース10の三分の一オクターブ(b=3)の場合:G=10^(3/10). 隣接する中心周波数の比はr = G^(1/3) = 10^(1/10) ≈ 1.258925であり、帯域端の倍率はk = 10^(1/20) ≈ 1.122018です.

Qファクタと分解能:オクターブ分析は定Q分析である

Q = fm / (f2 − f1)と定義します. CPB帯域ではΔf = f2 − f1はfmとともにスケールするため、Qは周波数には依存せず、bとGのみに依存します.

簡易リファレンス(ベース10、fr=1000 Hz):

分数オクターブ	帯域比 f2/f1	相対帯域幅 Δf/fm	Q = fm/Δf
1/1	1.995262	0.704592	1.419
1/2	1.412538	0.347107	2.881
1/3	1.258925	0.230768	4.333
1/6	1.122018	0.115193	8.681
1/12	1.059254	0.057573	17.369

解釈:1/3オクターブではQ≈4.33であり、各帯域の幅は中心周波数に対しておよそ23%です. より細かい帯域(1/6, 1/12など)は分解能が高くなりますが、ランダム雑音に対する分散が大きくなり、一般により長い平均化時間が必要になります.

帯域番号付け(整数インデックス)と式による列挙

実装では、整数の帯域インデックスxを用いることがよくあります. IECでは、このxが中心周波数の式fm = fr * G^(x/b)の中に直接現れます. [1] これにより、対象とする周波数範囲をカバーする全帯域を安定して列挙でき、隙間のない、規格に一貫した帯域端が保証されます.

ベース10の場合:

したがって

and you can invert as

図4:代表的な分数オクターブ帯域(ベース10の定義)に対するQファクタ

「三分の一オクターブ」の二つの意味:ベース2 vs ベース10-混同しないこと

文献によっては2ベースを用い、隣接する中心周波数が2^(1/3)になるとするものもあります. IEC 61260-1および現代の音響分野の多くの実務ではベース10が用いられ、隣接する中心周波数は10^(1/10)となります. 簡単なチェックとして、公称中心周波数が1.0k → 1.25k → 1.6k → 2.0k(R10スタイル)のような並びであれば、ベース10である可能性が高いと言えます.

帯域レベルの数学的定義:PSD積分からdB報告まで

連続周波数の見方:帯域内でPSDを積分する

オクターブ帯域レベルは、本質的にはある周波数帯域にわたるパワースペクトル密度の積分です. 音圧p(t)の場合:

振動(速度／加速度)の場合も、単位と基準量が異なるだけで同様の考え方が適用されます.

重要な点:dBは対数量であるため、和や平均は必ず線形の電力／平均二乗量の領域で行ったうえで、最後にdBへ変換しなければなりません.

二つの離散実装:フィルタバンクRMS vs FFT／PSDビニング

フィルタバンク法:y_b(t)=BandPass_b{x(t)}とし、その平均二乗mean(y_b^2)を帯域平均二乗値として計算します(必要に応じて時間平均を実施).

FFT／PSDビニング法:S_pp(f)を(ピリオドグラム／Welch法などで)推定し、[f1,f2]内のビンを数値積分／加算します.

長く定常的な信号で平均をとれば、これら二つの手法の結果は非常に近くなり得ます. 一方、過渡、スイープ、短時間事象では、両者はしばしば異なる結果を生じます.

自分が扱っているスペクトルが何かを明確にする:振幅、電力、PSD(およびdB/Hz)

• 振幅スペクトル|X(f)|:振幅単位(例:Pa)で表され、トーン／高調波に有用です.
• パワースペクトル|X(f)|²:平均二乗単位(Pa²)です.
• パワースペクトル密度(PSD):Hzあたりの平均二乗量(Pa²/Hz)で表され、雑音に最もよく用いられます.

オクターブ帯域レベルは帯域平均二乗／電力を表すため、どの表現から出発したとしても、最終的にはPa²(あるいは同等の量)で積分／加算する必要があります.

周波数分解能と片側スペクトル:Δf, 0..fs/2, および「×2」ルール

FFTのビン間隔はΔf = fs/Nです. 典型的な離散近似は次のようになります:

片側スペクトル(0..fs/2)を用いる場合、エネルギー保存のため、DCおよびナイキストを除くすべてのビンを通常2倍します(負の周波数側の電力を正の側に折り畳むため). ソフトウェアによってこれらの取り扱いは異なるため、結果を比較する前に定義をそろえておきましょう.

窓補正:コヒーレントゲイン(トーン)とENBW(雑音)は別物

窓関数はスペクトルリークを低減しますが、スケーリングを変化させます:

• トーン振幅に対しては、コヒーレントゲイン(CG)で補正し、しばしばCG = sum(w)/Nが用いられます.
• 広帯域雑音／PSDに対しては、等価雑音帯域幅(ENBW)で補正し、例えばENBW = fs·sum(w²)/(sum(w))²とします. [9]

CGはピーク振幅を、ENBWは平均的なノイズフロアの面積を制御します. オクターブ帯域レベルはエネルギー統計量であり、ENBWの影響をより強く受けます.

ウィンドウ	コヒーレントゲイン(CG)	ENBW(ビン数)
矩形	1.000	1.000
ハン	0.500	1.500
ハミング	0.540	1.363
ブラックマン	0.420	1.727

部分ビン重み付け:帯域端がFFTビンにそろわない場合の扱い

帯域端がビンの周波数にちょうど一致することはほとんどありません. PSDを、幅Δfの各ビン内でほぼ一定とみなし、境界ビンについては帯域との重なり率で重み付けします:

これにより、Nや帯域端が変化しても、より滑らかで物理的に一貫した帯域レベルが得られます.

図5:帯域端がFFTビンに整合しない場合の部分ビン重み付けの模式図

統一的な式:両手法は∫|H_b(f)|² S_xx(f) dfを計算している

フィルタバンク法とPSDビニング法はいずれも、次の形で表現できます:

ブリックウォール型ビニングは、|H_b|²が[f1,f2]の内側で1、外側で0である場合に対応します. 規格に準拠した真のフィルタにはロールオフとリップルが存在するため、規格はマスクと有効帯域幅に制約を課しているのです.

帯域集約:1/3オクターブから1オクターブを合成し、全体レベルを形成する

理想的な分割とエネルギー収支のもとでは:

• 隣接する三つの三分の一オクターブ帯域を組み合わせることで、一つのフルオクターブ帯域を近似できます.
• カバーしている全帯域のエネルギーを合計すれば、総エネルギーが得られます.

必ずエネルギー領域で合成してください.

L_iをdB表記の帯域レベルとすると、エネルギーはE_i = 10^(L_i/10)となります. このとき:

IEC 61260-1は、分数オクターブの結果を合成して、より広い帯域レベルを形成できることを明示しています. [1]

有効帯域幅:なぜ規格がそれを規定するのか

実際のフィルタは理想的な矩形特性ではありません. 白色雑音(PSDが一定S0)の場合、出力平均二乗値は次の式で表されます:

ピンクノイズ(PSD ~ 1/f)のような非白色スペクトルに対しては、標準的な工学的雑音スペクトル全体で結果の比較性を保つため、重み付けを伴う正規化有効帯域幅が規格で定義される場合があります. [1]

実務上の含意として、FFTによる「ハードビニング」は、B_eff = (f2 − f1)をもつブリックウォールフィルタを暗黙に仮定していることになります. 規格に準拠したオクターブフィルタにはスカートが存在するため、B_effはわずかに(かつクラスによって)異なり得ます. 結果をそろえたい場合は、標準の|H(f)|²を周波数領域で近似するか、手法の違いを明示的に文書化してください.

なぜ1/3オクターブが好まれるのか(数学＋聴覚＋工学的トレードオフ)

情報密度が「ちょうどよい」:1オクターブより細かく、極細分より安定

一つのオクターブ帯域だけでは粗すぎてスペクトル形状を隠してしまう一方で、非常に細かい分数(例:1/12, 1/24)は不安定でコストが高くなり得ます:

• ランダム雑音に対する推定分散が大きくなる(各帯域が捉えるエネルギーが少ない).
• 計算量が増え、報告の負担も大きくなります.
• 多くの場合、規制や評価スキームが必要とする以上の詳細になってしまいます.

三分の一オクターブは古典的な折衷案であり、工学的な洞察を得るのに十分な分解能を持ちつつ、標準化された測定に耐える安定性があり、多くの測定器やソフトウェアで広くサポートされています.

聴覚心理学:中域の臨界帯域は1/3オクターブに近い

多くの聴覚心理学の文献では、可聴域全体におよそ24個の臨界帯域があると記述されており、中周波数帯では臨界帯域幅がおおむね1/3オクターブ帯域幅に近いとされています. [7][8] これにより、1/3オクターブは知覚音に関わる問題に対する自然な中間表現となり、Bark／ERBスケールほど特殊ではない一方で、より標準化された表現として利用できます.

規格／アプリケーションからの直接的な牽引:多くのワークフローが1/3オクターブ入出力を必須としている

主要な規格が三分の一オクターブで入力／出力を定義すると、測定器、ソフトウェア、レポート様式(テンプレート)といったエコシステム全体がそれに収束していきます. 例:

• 建築音響評価:ISO 717-1は、一値指標の算出に三分の一オクターブ帯域を参照しています. [5]
• 室内音響パラメータ(例:残響時間)は、ISO 3382シリーズにおいて、オクターブ／三分の一オクターブ帯域で報告されるのが一般的です. [6]

追加のベース10の利点:R10表、10帯域／デケード、可読性

• 10帯域／デケード:周波数を10倍することが、正確に10回の三分の一オクターブステップに対応します(対数プロットにとって非常に扱いやすい).
• R10推奨数:1.00, 1.25, 1.60, 2.00, 2.50, 3.15, 4.00, 5.00, 6.30, 8.00 (×10^n) は広く認知されており、伝達しやすい値です.
• ベース2と比べて、10進ラベリングの方が扱いにくさが少なく、規格間のあいまいさも減ります.

オクターブ帯域分析は、通常FFTビニングまたはフィルタバンクのいずれかの方法で実装されます. 続きを読む -> オクターブ帯域分析ガイド:FFTビニングとフィルタバンクの比較

OpenTestは両方の手法を統合しています.今すぐダウンロードして使い始めましょう -> または、下のフォームに入力してライブデモの予約を行ってください.

www.opentest.comで、さらに多くの機能やアプリケーション事例をご覧ください.

参考

[1] IEC 61260-1:2014 PDFサンプル(iTeh):https://cdn.standards.iteh.ai/samples/13383/3c4ae3e762b540cc8111744cb8f0ae8e/IEC-61260-1-2014.pdf

[2] ISO 266:1997, Acoustics - Preferred frequencies(ISO):https://www.iso.org/obp/ui/

[3] ANSI S1.11-2004 プレビューPDF(ASA/ANSI):https://webstore.ansi.org/preview-pages/ASA/preview_ANSI%2BS1.11-2004.pdf

[4] ANSI/ASA S1.11-2014/Part 1 / IEC 61260-1:2014 プレビュー:https://webstore.ansi.org/preview-pages/ASA/preview_ANSI%2BASA%2BS1.11-2014%2BPart%2B1%2BIEC%2B61260-1-2014%2B%28R2019%29.pdf

[5] ISO 717-1:2020 abstract(one-third-octave usageに言及):https://www.iso.org/standard/77435.html

[6] ISO 3382-2:2008 アブストラクト(室内音響パラメータ):https://www.iso.org/standard/36201.html

[7] Ansys Help: Bark scale and critical bands(midrangeがthird octaveに近いことに言及):https://ansyshelp.ansys.com/public/Views/Secured/corp/v252/en/Sound_SAS_UG/Sound/UG_SAS/bark_scale_and_critical_bands_179506.html

[8] Simon Fraser University Sonic Studio Handbook: Critical Band and Critical Bandwidth:https://www.sfu.ca/sonic-studio-webdav/cmns/Handbook5/handbook/Critical_Band.html

[9] MathWorks: ENBW定義の例:https://www.mathworks.com/help/signal/ref/enbw.html