nsamples = freq * duration datafile = 'cry_n0000001.wav' micgeofile = path.join(path.split(bpath)[0], 'xml', 'array_64_8mic.xml') h5savefile = 'cry_n0000001.wav' m = MicGeom(from_file=micgeofile) n = WNoiseGenerator(sample_freq=sfreq, numsamples=nsamples, seed=1) #1pascal p = PointSource(signal=n, mpos=m, loc=(-0.1, -0.1, 0.3)) p1 = Mixer(source=p) wh5 = WriteH5(source=p, name=h5savefile) wh5.save() #definition the different source signal r = 52.5 nsamples = long(sfreq * 0.3) 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 a circular array of 8 microphones m = MicGeom('array_64_8mic.xml') m.mpos_tot = array([(r*sin(2*pi*i+pi/8), r*cos(2*pi*i+pi/8), 0) \ for i in linspace(0.0, 1.0, 8, False)]).T t = MaskedTimeSamples(name=datafile) f = PowerSpectra(time_data=t, window='Hanning', overlap='50%', block_size=128, \ ind_low=1,ind_high=30) g = RectGrid(x_min=0 - .2, x_max=0.3, y_min=-0.13, y_max=0.3, z=0,
m.mpos_tot = array([(r*sin(2*pi*i+pi/4), r*cos(2*pi*i+pi/4), 0) \ for i in linspace(0.0, 1.0, 28, False)]).T mg_file = path.join(path.split(acoular.__file__)[0], 'xml', 'array_64.xml') mg = MicGeom(from_file=mg_file) #=============================================================================== # define the different 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) p1 = PNoiseGenerator(seed=1, numsamples=nsamples) #=============================================================================== # 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) #===============================================================================