Skip to content

Power Module

Sound power level calculations.

Functions

lw_iso3746 

lw_iso3746(
    LpAi: ndarray,
    LpAiB: ndarray,
    S: float,
    alpha: ndarray,
    surfaces: ndarray,
) -> float

Calculate sound power level according to ISO 3746:2010.

PARAMETER DESCRIPTION
LpAi

Sound pressure levels of the source \(L_{pAi}\).

TYPE: ndarray

LpAiB

Background noise sound pressure levels \(L_{pAiB}\).

TYPE: ndarray

S

Area in square meters of the measurement surface \(S\).

TYPE: float

alpha

Absorption coefficients of the room \(\alpha\).

TYPE: ndarray

surfaces

Room surfaces.

TYPE: ndarray

RETURNS DESCRIPTION
float

Sound power level \(L_{w}\).
Source code in acoustic_toolbox/power.py 
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
def lw_iso3746(
    LpAi: np.ndarray,
    LpAiB: np.ndarray,
    S: float,
    alpha: np.ndarray,
    surfaces: np.ndarray,
) -> float:
    r"""Calculate sound power level according to ISO 3746:2010.

    Args:
      LpAi: Sound pressure levels of the source $L_{pAi}$.
      LpAiB: Background noise sound pressure levels $L_{pAiB}$.
      S: Area in square meters of the measurement surface $S$.
      alpha: Absorption coefficients of the room $\alpha$.
      surfaces: Room surfaces.

    Returns:
      Sound power level $L_{w}$.
    """
    LpA = 10.0 * np.log10(np.sum(10.0 ** (0.1 * LpAi)) / LpAi.size)
    LpAB = 10.0 * np.log10(np.sum(10.0 ** (0.1 * LpAiB)) / LpAiB.size)
    deltaLpA = LpA - LpAB

    if deltaLpA > 10.0:
        k_1a = 0.0
    elif 3.0 <= deltaLpA <= 10.0:
        k_1a = -10.0 * np.log10(1.0 - 10.0 ** (-0.1 * deltaLpA))
    else:
        # This should alert to user because poor condition of the measurement.
        k_1a = 3.0

    S0 = 1.0
    Sv = np.sum(surfaces)
    alpha_mean = np.average(alpha, axis=0, weights=surfaces)
    A = alpha_mean * Sv

    k_2a = 10.0 * np.log10(1.0 + 4.0 * S / A)

    LpA_mean = LpA - k_1a - k_2a
    L_WA = LpA_mean + 10.0 * np.log10(S / S0)
    return L_WA

:::