def create_test_time_data(nsamples):
"""
creates test data for a single moving monopole emitting white noise
Parameters
----------
speed : float
source velocity in m/s.
Returns
-------
None.
"""
n1 = WNoiseGenerator(sample_freq=SFREQ, numsamples=nsamples, seed=SEED)
p1 = MovingPointSource(signal=n1,
mics=MGEOM,
trajectory=TRAJ,
conv_amp=CONV_AMP)
wh5 = WriteH5(source=p1, name=FNAME)
print(50 * "#")
print(f"create {FNAME} ...")
print(f"num samples: {nsamples}, pass-by time: {t_passby}s")
print(50 * "#")
wh5.save()
#===============================================================================
# define the different source signals
#===============================================================================
if sys.version_info > (3, ):
long = int
nsamples = long(sfreq * tmax)
n1 = WNoiseGenerator(sample_freq=sfreq, numsamples=nsamples)
s1 = SineGenerator(sample_freq=sfreq, numsamples=nsamples, freq=freq)
s2 = SineGenerator(sample_freq=sfreq, numsamples=nsamples, freq=freq, \
phase=pi)
#===============================================================================
# define the moving source and one fixed source
#===============================================================================
p0 = MovingPointSource(signal=s1, mpos=m, trajectory=tr1)
#t = p0 # use only moving source
p1 = PointSource(signal=n1, mpos=m, loc=(0, R, Z))
t = Mixer(source=p0, sources=[
p1,
]) # mix both signals
#t = p1 # use only fix source
# uncomment to save the signal to a wave file
#ww = WriteWAV(source = t)
#ww.channels = [0,14]
#ww.save()
#===============================================================================
# fixed focus frequency domain beamforming
#===============================================================================
# define environment and grid
#===============================================================================
ufe = UniformFlowEnvironment(ma=0.5, fdv=(0, 1, 0))
g = RectGrid(x_min=-3.0,
x_max=+3.0,
y_min=-3.0,
y_max=+3.0,
z=Z,
increment=0.2)
#===============================================================================
# define some sources
#===============================================================================
mp = MovingPointSource(signal=s1, env=ufe, mpos=m, trajectory=tr1)
ps = PointSource(signal=s1, mpos=m, loc=(0, R, Z))
pd = PointSourceDipole(signal=s1,
mpos=m,
direction=(0.5, 0, 0),
loc=(0, -2, Z))
un = UncorrelatedNoiseSource(signal=n1, mpos=mg)
mix = SourceMixer(sources=[ps, pd])
#===============================================================================
# 3rd Octave Filter
#===============================================================================
fi = FiltFiltOctave(source=un, band=freq, fraction='Third octave')
#==============================================================================