# set up microphone at (0,0,0)
m = MicGeom()
m.mpos_tot = array([[0, 0, 0]])
# create noise source
n1 = WNoiseGenerator(sample_freq=sfreq, numsamples=10 * sfreq, seed=1)
t = PointSource(signal=n1, mics=m, loc=(1, 0, 1))
# create power spectrum
f = PowerSpectra(time_data=t, window='Hanning', overlap='50%', block_size=4096)
# get spectrum data
spectrum_data = real(f.csm[:, 0,
0]) # get power spectrum from cross-spectral matrix
freqs = f.fftfreq() # FFT frequencies
# use barspectrum from acoular.tools to create third octave plot data
(f_borders, p, f_center) = barspectrum(spectrum_data, freqs, band, bar=True)
(f_borders_, p_, f_center_) = barspectrum(spectrum_data,
freqs,
band,
bar=False)
# create figure with barspectra
figure(figsize=(20, 6))
title("Powerspectrum")
plot(f_borders, L_p(p), label="bar=True")
plot(f_borders_, L_p(p_), label="bar=False")
xlim(f_borders[0] * 2**(-1. / 6), f_borders[-1] * 2**(1. / 6))
ylim(50, 90)
mpos = m,
loc = (1, 0, 1))
f = PowerSpectra(time_data = t,
window = 'Hanning',
overlap = '50%',
block_size = 4096)
###################################################################
### Plotting ###
from pylab import figure,plot,show,xlim,ylim,xscale,xticks,xlabel,ylabel,grid,real
from acoular import L_p
band = 3 # octave: 1 ; 1/3-octave: 3
(f_borders, p, f_center) = barspectrum(real(f.csm[:,0,0]), f.fftfreq(), band)
label_freqs = [str(int(_)) for _ in f_center]
figure(figsize=(20, 6))
plot(f_borders,L_p(p))
xlim(f_borders[0]*2**(-1./6),f_borders[-1]*2**(1./6))
ylim(40,90)
xscale('symlog')
xticks(f_center,label_freqs)
xlabel('f in Hz')
ylabel('SPL in dB')
m = MicGeom()
m.mpos_tot = array([[0, 0, 0]])
t = PointSource(signal=n1, mpos=m, loc=(1, 0, 1))
f = PowerSpectra(time_data=t,
window='Hanning',
overlap='50%',
block_size=4096)
###################################################################
### Plotting ###
from pylab import figure, plot, show, xlim, ylim, xscale, xticks, xlabel, ylabel, grid, real
from acoular import L_p
band = 3 # octave: 1 ; 1/3-octave: 3
(f_borders, p, f_center) = barspectrum(real(f.csm[:, 0, 0]), f.fftfreq(),
band)
label_freqs = [str(int(_)) for _ in f_center]
figure(figsize=(20, 6))
plot(f_borders, L_p(p))
xlim(f_borders[0] * 2**(-1. / 6), f_borders[-1] * 2**(1. / 6))
ylim(40, 90)
xscale('symlog')
xticks(f_center, label_freqs)
xlabel('f in Hz')
ylabel('SPL in dB')