def test_third2oct_2darray_axis0():
levels = np.array([[1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0],
[10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0],
[100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0]])
generated = third2oct(levels, axis=0)
real = np.array([100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0])
assert_array_almost_equal(generated, real)
def test_third2oct_2darray():
levels = np.array([[100, 95, 80, 55, 65, 85, 75, 70, 90, 95, 105, 110],
[100, 95, 80, 55, 65, 85, 75, 70, 90, 95, 105, 110]])
generated = third2oct(levels, axis=1)
real = np.array([[101.22618116, 85.04751156, 90.17710468, 111.29641738],
[101.22618116, 85.04751156, 90.17710468, 111.29641738]])
assert_array_almost_equal(generated, real)
def test_third2oct_2darray():
levels = np.array([[100, 95, 80, 55, 65, 85, 75, 70, 90, 95, 105, 110],
[100, 95, 80, 55, 65, 85, 75, 70, 90, 95, 105, 110]])
generated = third2oct(levels, axis=1)
real = np.array([[101.22618116, 85.04751156, 90.17710468, 111.29641738],
[101.22618116, 85.04751156, 90.17710468, 111.29641738]])
assert_array_almost_equal(generated, real)
def test_third2oct_2darray_axis1():
levels = np.array([[1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0],
[10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0],
[100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0]])
generated = third2oct(levels, axis=1)
real = np.array([[5.77121255, 5.77121255, 5.77121255], [14.77121255, 14.77121255, 14.77121255],
[104.77121255, 104.77121255, 104.77121255]])
assert_array_almost_equal(generated, real)
def test_third2oct_2darray_axis0():
levels = np.array([[ 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0],
[ 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0],
[100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0]
])
generated = third2oct(levels, axis=0)
real = np.array([100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0])
assert_array_almost_equal(generated, real)
def test_third2oct_2darray_axis1():
levels = np.array(
[[1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0],
[10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0],
[100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0]])
generated = third2oct(levels, axis=1)
real = np.array([[5.77121255, 5.77121255, 5.77121255],
[14.77121255, 14.77121255, 14.77121255],
[104.77121255, 104.77121255, 104.77121255]])
assert_array_almost_equal(generated, real)
def test_third2oct_3darray_axis0():
# Array of ones with shape (3,4,5)
levels = np.array([[[1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.]],
[[1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.]],
[[1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.]]])
generated = third2oct(levels, axis=0)
real = np.array([[5.77121255, 5.77121255, 5.77121255, 5.77121255, 5.77121255],
[5.77121255, 5.77121255, 5.77121255, 5.77121255, 5.77121255],
[5.77121255, 5.77121255, 5.77121255, 5.77121255, 5.77121255],
[5.77121255, 5.77121255, 5.77121255, 5.77121255, 5.77121255]])
assert_array_almost_equal(generated, real)
def test_third2oct_3darray_axis0():
# Array of ones with shape (3,4,5)
levels = np.array([[[1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.],
[1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.]],
[[1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.],
[1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.]],
[[1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.],
[1., 1., 1., 1., 1.], [1., 1., 1., 1., 1.]]])
generated = third2oct(levels, axis=0)
real = np.array(
[[5.77121255, 5.77121255, 5.77121255, 5.77121255, 5.77121255],
[5.77121255, 5.77121255, 5.77121255, 5.77121255, 5.77121255],
[5.77121255, 5.77121255, 5.77121255, 5.77121255, 5.77121255],
[5.77121255, 5.77121255, 5.77121255, 5.77121255, 5.77121255]])
assert_array_almost_equal(generated, real)
def test_third2oct():
levels = np.array([10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0])
generated = third2oct(levels)
real = np.array([14.77121255, 14.77121255, 14.77121255])
assert_array_almost_equal(generated, real)
def test_third2oct():
levels = np.array([10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0])
generated = third2oct(levels)
real = np.array([14.77121255, 14.77121255, 14.77121255])
assert_array_almost_equal(generated, real)
print "reference SPL:", acd.REFERENCE_PRESSURE print "equivalent SPL:",acd.equivalent_sound_pressure_level(sigy) print "peak SPL:", acd.peak_sound_pressure(sigy) """ IEC 61672 """ import acoustics.standards as acs print "third octave center frequency vector:", acs.iec_61672_1_2013.NOMINAL_OCTAVE_CENTER_FREQUENCIES print "reference frequency:", acs.iec_61672_1_2013.REFERENCE_FREQUENCY print "slow time constant:", acs.iec_61672_1_2013.SLOW """ third octave band """ import acoustics.bands as acb print "some thrid octave frequency:", acb.third(100, 500) print "third octave to octave:",acb.third2oct([10,10,10,10,10,10]) """ third octave analysis """ import acoustics.signal as ass a,AX = ass.third_octaves(sigy, fs) a,AZ1 = ass.third_octaves(sigz1, fs) a,AZ2 = ass.third_octaves(sigz2, fs) """ show graphics using pylab """ import pylab as plt plt.hold(True) x = acoustics.bands.THIRD_OCTAVE_CENTER_FREQUENCIES
