def zap(msec=10, samplingfreq=44100, amplitude=0.95, shape='sine'): """create a bipolar pulse"""### f = 1000.0 / float(msec) if shape in ['sine', 'sin', 'square']: carrier = Basic.wavegen(freq_hz=f, samplingfreq_hz=samplingfreq, duration_msec=msec, container=wav(), waveform=numpy.sin) stim = Basic.ampmod(carrier, freq_hz=f) else: raise ValueError, 'unrecognized shape "%s"' % shape if shape == 'square': y = stim.y y[numpy.where(y < 0)] = -1 y[numpy.where(y > 0)] = +1 stim.autoscale(amplitude) return stim
def amstim(msec=1000,
modfreq=20,
carrierfreq=200,
samplingfreq=1000,
amplitude=1.0,
carriershape=numpy.sin,
modshape=numpy.sin):
"""create an amplitude-modulated periodic wave""" ###
carrier = Basic.wavegen(freq_hz=carrierfreq,
samplingfreq_hz=samplingfreq,
duration_msec=msec,
container=wav(),
waveform=carriershape)
stim = Basic.ampmod(carrier, freq_hz=modfreq, waveform=modshape)
stim.autoscale(amplitude)
return stim
def zap(msec=10, samplingfreq=44100, amplitude=0.95, shape='sine'):
"""create a bipolar pulse""" ###
f = 1000.0 / float(msec)
if shape in ['sine', 'sin', 'square']:
carrier = Basic.wavegen(freq_hz=f,
samplingfreq_hz=samplingfreq,
duration_msec=msec,
container=wav(),
waveform=numpy.sin)
stim = Basic.ampmod(carrier, freq_hz=f)
else:
raise ValueError, 'unrecognized shape "%s"' % shape
if shape == 'square':
y = stim.y
y[numpy.where(y < 0)] = -1
y[numpy.where(y > 0)] = +1
stim.autoscale(amplitude)
return stim
def amstim(msec=1000, modfreq=20, carrierfreq=200, samplingfreq=1000, amplitude=1.0, carriershape=numpy.sin, modshape=numpy.sin): """create an amplitude-modulated periodic wave"""### carrier = Basic.wavegen(freq_hz=carrierfreq, samplingfreq_hz=samplingfreq, duration_msec=msec, container=wav(), waveform=carriershape) stim = Basic.ampmod(carrier, freq_hz=modfreq, waveform=modshape) stim.autoscale(amplitude) return stim
