옥타브 대역 분석: 수학적 및 공학적 근거

옥타브-밴드 분석은 로그 주파수 축에서 일정 비율 대역폭(constant-percentage bandwidth)을 사용하여 상세 스펙트럼을 표준화된 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의 핵심)

수평 축을 로그 스케일로 바꾸면, 이 밴드들은 동일한 폭/동일한 간격으로 보이며, 이것이 바로 "상대 대역폭 일정(constant percentage bandwidth)"의 의미입니다.

이 두 그림은 핵심 개념을 포착합니다. 옥타브 대역 분석은 Hz에서의 동일 간격이 아니라, 로그 주파수 스케일에서의 동일 간격을 사용합니다.

표준 및 용어: IEC/ANSI/ISO 시스템은 실제로 무엇을 규정하는가?

실무에서 "1/3 옥타브 분석을 한다"는 것은 단순히 대역 경계를 정하는 것 이상의 제약을 받습니다. 표준은 중심 주파수의 정의(정확값 vs. 공칭값), 옥타브 비 정의(10진 기반 vs. 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의 경우가 왜 "반 스텝 이동된(half-step shifted)" 것처럼 보일까? 직관적으로, 중심 주파수 그리드는 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)이 선택되면, 전체 대역 집합(중심과 경계)은 수학적으로 고정됩니다.

정확값 vs. 공칭값: 하나의 대역에 두 개의 "중심 주파수"가 등장하는 이유

"정확한" 중심 주파수는 수학적으로 일관된 정의와 필터 설계를 위해 사용되고, "공칭" 값은 라벨링과 보고에 사용됩니다. [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))을 규정합니다. 주요 동기는 다음과 같습니다:

• Alignment with decimal preferred-number series (ISO 266 is tied to R10). [2]
• 국제적 일관성: IEC 61260-1:2014는 10진 기반을 규정하며, 2진 기반 설계는 기준 주파수에서 멀어질수록 규격을 유지하기가 더 어렵다고 명시합니다. [1]

10진 기반에서 1/3 옥타브는 10^(1/10) ≈ 1.258925(1/10 데케이드로도 해석 가능)에 해당하며, 이로써 한 데케이드당 10개의 1/3 옥타브 밴드라는 깔끔한 매핑이 형성됩니다.

"10 one-third-octave bands = 1 decade": why this matters

With base-10 one-third-octave spacing, each step multiplies frequency by r = 10^(1/10). 따라서:

• 연속된 10개의 1/3 옥타브 대역은 주파수를 정확히 10배(한 데케이드) 증가시킵니다.
• 이는 ISO 266/R10 관행과 일치하며, 표, 플로팅 및 커뮤니케이션을 단순화합니다.
• 표준화에서는 순수한 수학적 엄밀성만큼이나 가독성과 일관성도 중시됩니다.

그림 3: Base-10 one-third-octave spacing-10 equal ratio steps per decade (×10 in frequency)

ANSI S1.11 / ANSI/ASA S1.11: 허용도 클래스와 과도 신호에 대한 주의 사항

ANSI S1.11(및 이후 IEC 61260-1과 정렬된 ANSI/ASA 채택본)은 필터 셋과 분석기에 대한 성능 요구사항을 규정하며, 허용도 클래스(판에 따라 보통 클래스 0/1/2)를 포함합니다. [3][4]

ANSI 문서에는 실무적인 주의 사항이 있습니다. 과도 신호의 경우, 규격을 만족하는 서로 다른 구현이라도 서로 다른 결과를 낼 수 있다는 점입니다. [3] 이는 과도 응답 분석에서 시간 응답(군지연, 링잉, 평균 시간 상수 등)이 중요하다는 점을 강조합니다.

클래스/마스크/유효 대역폭은 실제로 무엇을 제어하는가?

"1/3 옥타브 대역을 사용했다"는 것은 공칭 대역 경계만을 의미하지 않습니다. 표준은 다음 요소들을 제한함으로써, 서로 다른 계측기/알고리즘이 비교 가능한 결과를 내도록 하는 것을 목표로 합니다:

• 주파수 간격: 중심 주파수 시퀀스와 경계 정의(10진 기반, 정확값/공칭값, f1/f2).
• 진폭 응답 허용도(마스크): 통과 대역 근처의 허용 리플과 중심에서 멀어진 곳에서 요구되는 감쇠량.
• 광대역 소음에 대한 에너지 일관성: 구현 간에 대역 레벨이 비교 가능하도록 유효 대역폭에 대한 제약.

유효 대역폭이 중요한 이유는 실제 필터가 이상적인 벽형(brick-wall) 특성을 가지지 않기 때문입니다. 광대역 소음의 경우, 출력 에너지는 ∫|H(f)|^2 S(f)df에 의존합니다. 통과 대역 리플, 스커트(skirt), 롤오프 차이로 인해 체계적인 오프셋이 발생할 수 있습니다. 표준은 이러한 오프셋을 허용 가능한 범위 안에 두기 위해 유효 대역폭을 제한합니다. [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는 옥타브 비입니다. 그러면:

For base-10 one-third octave (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진 기반의 경우:

so

and you can invert as

그림 4: 일반적인 분수 옥타브 대역폭(10진 정의)에 대한 Q 인자

"1/3 옥타브"의 두 가지 의미: 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를 대역 내에서 적분

Octave-band level is essentially the integral of power spectral density over a frequency band. 음압 p(t)의 경우:

진동(속도/가속도)의 경우에도, 단위와 기준 값만 다를 뿐 동일한 논리가 적용됩니다.

핵심 포인트: dB는 로그 척도이므로, 모든 합산 또는 평균은 먼저 선형 전력/평균제곱 영역에서 수행해야 합니다.

두 가지 이산 구현: 필터 뱅크 RMS vs. FFT/PSD 빈닝

필터 뱅크 방법: y_b(t)=BandPass_b{x(t)}를 구한 뒤, 시간 평균(옵션)을 포함하여 mean(y_b^2)를 해당 대역의 평균제곱으로 계산합니다.

FFT/PSD 빈닝 방법: S_pp(f)를(예: 피리어도그램/웰치 방법으로) 추정한 뒤, [f1,f2] 내의 빈을 수치 적분/합산합니다.

긴 정Stationary 신호의 경우, 평균화된 결과는 두 방법이 매우 근접할 수 있습니다. 과도 신호, 스위프, 짧은 이벤트의 경우에는 두 방법의 결과가 종종 다릅니다.

사용 중인 스펙트럼이 무엇인지 명확히 하십시오: 크기, 파워, 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를 곱합니다(음의 주파수에 있던 전력이 양의 쪽으로 접혀 들어오기 때문입니다). 소프트웨어마다 이러한 관행을 처리하는 방식이 다르므로, 결과를 비교하기 전에 정의를 맞추어야 합니다.

창 함수 보정: 코히어런트 게인(톤) vs. ENBW(노이즈)는 서로 다르다

창 함수를 사용하면 스펙트럼 누설은 줄어들지만 스케일링이 변합니다:

• 톤 진폭의 경우: 코히어런트 게인(CG)으로 보정하며, 종종 CG = sum(w)/N을 사용합니다.
• 광대역 노이즈/PSD의 경우: 등가 노이즈 대역폭(ENBW)으로 보정하며, 예를 들어 ENBW = fs·sum(w²)/(sum(w))²입니다. [9]

CG는 피크 진폭을 제어하고, ENBW는 평균 노이즈 플로어 면적을 제어합니다. 옥타브 대역 레벨은 에너지 통계량이므로 ENBW에 더 민감합니다.

윈도우	코히어런트 게인(CG)	ENBW(빈 수 기준)
사각	1.000	1.000
Hann	0.500	1.500
해밍	0.540	1.363
Blackman	0.420	1.727

부분 빈 가중치: 대역 경계가 FFT 빈에 정렬되지 않을 때 무엇을 할 것인가

대역 경계가 빈 주파수에 정확히 위치하는 경우는 드뭅니다. PSD가 폭 Δf를 갖는 각 빈 내부에서 거의 상수라고 보고, 경계에 걸친 빈은 중첩 비율에 따라 가중치를 부여합니다:

이렇게 하면 N이나 대역 경계가 변하더라도 더 매끄럽고 물리적으로 일관된 대역 레벨을 얻을 수 있습니다.

그림 5: 대역 경계가 FFT 빈과 일치하지 않을 때의 부분 빈 가중치 개념도

통합 공식의 통일적 관점: 두 방법 모두 ∫|H_b(f)|² S_xx(f) df를 계산한다

필터 뱅크와 PSD 빈닝 두 방법은 모두 다음과 같이 쓸 수 있습니다:

벽형(brick-wall) 빈닝은 |H_b|²가 [f1,f2] 내부에서는 1, 외부에서는 0인 경우에 해당합니다. 표준을 준수하는 실제 필터는 롤오프와 리플을 가지며, 이 때문에 표준에서 마스크와 유효 대역폭을 제한합니다.

대역 집계: 1/3 옥타브에서 1 옥타브로 합산하고, 전체 레벨을 형성하기

이상적인 분할과 에너지 보존이 이루어진다는 전제하에서:

• 연속된 세 개의 1/3 옥타브 대역을 합치면, 하나의 전체 옥타브 대역을 근사할 수 있습니다.
• 커버되는 범위 내 모든 대역 에너지를 합하면 전체 에너지가 됩니다.

항상 에너지 도메인에서 합산하십시오.

대역 레벨 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)는 불안정하고 계산 비용이 큽니다:

• 랜덤 노이즈에 대한 추정량 분산 증가(각 대역에 포착되는 에너지가 적어짐).
• 계산량 증가와 보고 작업 부담 증가.
• 규제나 등급 산정에서 요구하는 것보다 과도한 세부 정보를 제공하는 경우가 많습니다.

One-third octave is the classic compromise: enough resolution for engineering insight, stable enough for standardized measurements, and broadly supported by instruments and software.

정신음향학: 중간 주파수 대역에서의 임계 대역폭은 1/3 옥타브에 근접

많은 정신음향학 자료에서 가청 범위 전체에 약 24개의 임계 대역이 있다고 설명하며, 중간 주파수 영역에서 임계 대역폭은 종종 1/3 옥타브 대역폭과 유사합니다. [7][8] 이 때문에 1/3 옥타브는 지각된 소리와 관련된 문제를 다룰 때 자연스러운 중간 표현이 되며, Bark/ERB 스케일보다 더 표준화되어 있다는 장점도 있습니다.

직접적인 표준/응용 수요: 많은 워크플로가 1/3 옥타브 입·출력을 요구

Once major standards define inputs/outputs in 1/3 octave, ecosystems (instruments, software, reporting templates) converge around it. 예시:

• 건축 음향 등급: ISO 717-1은 단일 수치 성능지수 계산을 위해 1/3 옥타브 밴드를 참조합니다. [5]
• 실내 음향 파라미터(예: 잔향 시간)는 일반적으로 옥타브/1/3 옥타브 밴드(ISO 3382 시리즈)로 보고됩니다. [6]

추가 10진 기반 이점: R10 표, 데케이드당 10개 대역, 가독성

• 10 bands per decade: multiplying frequency by 10 corresponds to exactly 10 one-third-octave steps (very clean for log plots).
• R10 선호수: 1.00, 1.25, 1.60, 2.00, 2.50, 3.15, 4.00, 5.00, 6.30, 8.00 (×10^n)은 널리 알려져 있고 전달하기 쉽습니다.
• Compared with base-2, decimal labeling is less awkward and cross-standard ambiguity is reduced.

옥타브 대역 분석은 일반적으로 FFT 빈닝 또는 필터 뱅크 두 가지 방법 중 하나를 사용하여 구현됩니다. 계속 읽어 보십시오 -> 옥타브 대역 분석 가이드: FFT 빈닝 vs. 필터 뱅크

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 (mentions 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 (mentions midrange close to 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